Search

Mastitis in cattle is often caused by microorganism infections in the udder. The three most common pathogens are esculin-positive streptococci (SC+), coagulase-negative staphylococci (CNS), and Escherichia coli (E. coli). In a previous study, 10 SNPs were associated with somatic cell score and mastitis in diverse Holstein populations. We tested these SNPs for their effects on individual pathogen presence. Milk and pathogen samples of 3076 Holstein cows were collected from four farms. Samples were excluded if multiple pathogens were present at the same time. Records of the same pathogen within 14 days of each other were counted as one infection. This resulted in 1129 pathogen-positive samples. Cases and controls were in ratios of 20:80 for SC+, 8:92 for CNS, and 11:89 for E. coli. The lasso, backward, and forward methods were used to narrow down SNPs associated with pathogen presence. The suitability of the SNPs to separate the samples into cases or controls for each pathogen was indicated using ROC curves. The Cochran-Armitage (CAT) and the Jonckheere-Terpstra (JTT) tests evaluated the influence of the SNPs on pathogen presence. Finally, a generalised linear mixed model (GLMM) including fixed environmental effects and a random sire effect was fitted to the binary trait of pathogen presence to test for association. In total, six out of the 10 investigated SNPs showed associations with pathogen presence based on the forward method: Two SNPs each for SC+ (rs41588957, rs41257403) and CNS (rs109934030, rs109441194), and three for E. coli (rs109934030, rs41634110, rs41636878). The CAT and GTT tests linked four SNPs (rs41588957, rs41634110, rs109441194, rs41636878) to pathogen presence, two of which were confirmed with the GLMM (rs41634110, rs109441194), with effects on CNS and E. coli. The SNPs linked to CNS and those linked to E. coli explained 13.2% and 13.8% of the variance, compared to 19% and 18.4%, respectively, of the full model with all 10 SNPs. Half of the SNP genotypes previously linked to lower SCS also decreased the probability for pathogen presence and might therefore be targets not just for lower SCS but for a better pathogen resistance. Trial Registration: Not applicable, no new data were collected for this study., (© 2024 The Author(s). Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

Since their domestication, horses have accompanied mankind, and humans have constantly shaped horses according to their needs through stallion-centered breeding. Consequently, the male-specific portion of the Y chromosome (MSY) is extremely uniform in modern horse breeds. The majority of stallions worldwide carry MSY haplotypes (HT) attributed to an only ~1,500-y-old, so-called, "Crown" haplogroup. The predominance of the Crown in modern horse breeds is thought to represent a footprint of the vast impact of stallions of "Oriental origin" in the past millennium. Here, we report the results of a fine-scaled MSY haplotyping of large datasets of patrilines comprising 1,517 males of 189 modern horse breeds, covering a broad phenotypic and geographic spectrum. We can disentangle the multilayered influence of Oriental stallions over the last few hundred years, exposing the intense linebreeding and the wide-ranging impact of Arabian, English Thoroughbred, and Coldblood sires. Iberian and New World horse breeds contain a wide range of diversified Crown lineages. Their broad HT spectrum illustrates the spread of horses of Oriental origin via the Iberian Peninsula after the Middle Ages, which is commonly referred to as the "Spanish influence." Our survey also revealed a second major historical dissemination of horses from Western Asia, attributed to the expansion of the Ottoman Empire. Our analysis shows that MSY analysis can uncover the complex history of horse breeds and can be used to establish the paternal ancestry of modern horse breeds., Competing Interests: Competing interests statement:The authors declare no competing interest.

The beneficial effects of feeding probiotic Bacillus subtilis DSM 32315 (BS) and Bacillus velezensis CECT 5940 (BV) to chickens in vivo are well-documented, with potential immune modulation as a key mechanism. In this study, we investigated the direct interactions of chicken peripheral blood mononuclear cells (PBMCs) with BS or BV in vitro through whole transcriptome profiling and cytokine array analysis. Transcriptome profiling revealed 20 significantly differentially expressed genes (DEGs) in response to both Bacillus treatments, with twelve DEGs identified in BS-treated PBMCs and eight in BV-treated PBMCs. Pathway analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated significant regulation of immune-related pathways by both BS and BV. Notably, BS treatment upregulated genes associated with immune cell surface markers (CD4, CD25, CD28), anti-inflammatory cytokine interleukin-10 (IL-10), and C-C motif chemokine ligand 5 (CCL5), while downregulating the gene encoding pro-inflammatory IL-16. BV treatment similarly affected genes associated with immune cell surface markers, IL-16, and CCL5, with no impact on the gene encoding IL-10. Both treatments induced higher expression of the gene encoding the avian β-defensin 1 (AvBD1). The results of this in vitro study indicate an immunomodulatory effect of BS and BV in chicken PBMCs by regulating genes involved in anti-inflammatory, bacteriostatic, protective, and pro-inflammatory responses. Consequently, BS and BV may serve to augment the immune system's capacity to defend against infection by modulating immune responses and cytokine expression. Thus, the administration of these probiotics holds promise for reducing reliance on antimicrobials in farming practices., Competing Interests: Author GL was employed by the company Evonik Operations GmbH. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Larsberg, Sprechert, Hesse, Falker-Gieske, Loh, Brockmann and Kreuzer-Redmer.)

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and metabolic syndrome. The sublines BFMI861-S1 and BFMI861-S2 differ in weight despite high genetic similarity and a shared obesity-related locus. This study focused on identifying additional body weight quantitative trait loci (QTLs) by analyzing weekly weight measurements in a male population of the advanced intercross line BFMI861-S1 x BFMI861-S2. QTL analysis, utilizing 200 selectively genotyped mice (GigaMUGA) and 197 males genotyped for top SNPs, revealed a genome-wide significant QTL on Chr 15 (68.46 to 81.40 Mb) for body weight between weeks 9 to 20. Notably, this QTL disappeared (weeks 21 to 23) and reappeared (weeks 24 and 25) coinciding with a diet change. Additionally, a significant body weight QTL on Chr 16 (3.89 to 22.79 Mb) was identified from weeks 6 to 25. Candidate genes, including Gpt, Cbx6, Apol6, Apol8, Sun2 (Chr 15) and Trap1, Rrn3, Mapk1 (Chr 16), were prioritized. This study unveiled two additional body weight QTLs, one of which is novel and responsive to diet changes. These findings illuminate genomic regions influencing weight in BFMI and emphasize the utility of time series data in uncovering novel genetic factors., (© 2024. The Author(s).)

When resequencing animal genomes, some short reads cannot be mapped to the reference genome and are usually discarded. In this study, unmapped reads from 302 German Black Pied cattle were analyzed to identify potential pathogenic DNA. These unmapped reads were assembled and blasted against NCBI's database to identify bacterial and viral sequences. The results provided evidence for the presence of pathogens. We found sequences of Bovine parvovirus 3 and Mycoplasma species. These findings emphasize the information content of unmapped reads for gaining insight into bacterial and viral infections, which is important for veterinarians and epidemiologists., (© 2023. L’Institut National de Recherche en Agriculture, Alimentation et Environnement (INRAE).)

Background: The gut microbiota is modulated by a combination of diet, host genetics, and sex effects. The magnitude of these effects and interactions among them is important to understanding inter-individual variability in gut microbiota. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed along with several QTLs for metabolic traits. In the current study, we searched for genetic variants underlying differences in the gut microbiota in response to American and ketogenic diets, which are high in fat and vary in carbohydrate composition, between C57BL/6 J (B6) and FVB/NJ (FVB) mouse strains., Results: Genetic mapping of microbial features revealed 18 loci under the QTL model (i.e., marginal effects that are not specific to diet or sex), 12 loci under the QTL by diet model, and 1 locus under the QTL by sex model. Multiple metabolic and microbial features map to the distal part of Chr 1 and Chr 16 along with eigenvectors extracted from principal coordinate analysis of measures of β-diversity. Bilophila, Ruminiclostridium 9, and Rikenella (Chr 1) were identified as sex- and diet-independent QTL candidate keystone organisms, and Parabacteroides (Chr 16) was identified as a diet-specific, candidate keystone organism in confirmatory factor analyses of traits mapping to these regions. For many microbial features, irrespective of which QTL model was used, diet or the interaction between diet and a genotype were the strongest predictors of the abundance of each microbial trait. Sex, while important to the analyses, was not as strong of a predictor for microbial abundances., Conclusions: These results demonstrate that sex, diet, and genetic background have different magnitudes of effects on inter-individual differences in gut microbiota. Therefore, Precision Nutrition through the integration of genetic variation, microbiota, and sex affecting microbiota variation will be important to predict response to diets varying in carbohydrate composition. Video Abstract., (© 2023. BioMed Central Ltd., part of Springer Nature.)

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity., (© 2023. The Author(s).)

Genetic evaluations of local cattle breeds are hampered due to small reference groups or biased due to the utilization of SNP effects estimated in other large populations. Against this background, there is a lack of studies addressing the possible advantage of whole-genome sequences (WGS) or consideration of specific variants from WGS data in genomic predictions for local breeds with small population size. Consequently, the aim of this study was to compare genetic parameters and accuracies of genomic estimated breeding values (GEBV) for 305-d production traits, fat-to protein ratio (FPR), and somatic cell score (SCS) at the first test date after calving and confirmation traits of the endangered German Black Pied cattle (DSN) breed using 4 different marker panels: (1) the commercial 50K Illumina BovineSNP50 BeadChip, (2) a customized 200K chip designed for DSN (DSN200K) which considers the most important variants for DSN from WGS, (3) randomly generated 200K chips based on WGS data, and (4) a WGS panel. The same number of animals was considered for all marker panel analyses (i.e., 1,811 genotyped or sequenced cows for conformation traits, 2,383 cows for lactation production traits, and 2,420 cows for FPR and SCS). Mixed models for the estimation of genetic parameters directly included the respective genomic relationship matrix from the different marker panels plus the trait-specific fixed effects. For the calculation of GEBV accuracies, we applied repeated random subsampling validation. In the process of separate cross-validations per trait, we created a validation set including 20% of cows with masked phenotypes, and a training set comprising 80% of the cows. The cows were selected randomly in a procedure with 10 replicates considering replacements in the different scenarios. The accuracy was defined as the correlation between the direct GEBV and the phenotypes with subtracted corresponding fixed effects for the cows in the validation set. For FPR and SCS, as well as for lactation production traits, heritabilities were largest based on WGS data, but the increase compared with the 50K or DSN200K applications was quite small in the range from 0.01 to 0.03. Also, for most of the conformation traits, heritabilities were largest based on WGS and DSN200K data, but the increase was in the range of the corresponding standard error. Accordingly, GEBV accuracies for most of the studied traits were highest based on WGS data or when utilizing the DSN200K chip, but the accuracy differences across the marker panels were quite small and nonsignificant. In conclusion, WGS data and the DSN200K chip only contributed to minor improvements in genomic predictions, still justifying the use of the commercial 50K chip. Nevertheless, WGS and the 200KDSN chip harbor breed-specific variants, which are valuable for studying causal genetic mechanisms in the endangered DSN population., (The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

German Black Pied (DSN) is considered an ancestral population of the Holstein breed. The goal of the current study was to fine-map genomic loci for milk production traits and to provide sequence variants for selection. We studied genome-wide associations for milk-production traits in 2160 DSN cows. Using 11.7 million variants from whole-genome sequencing of 304 representative DSN cattle, we identified 1980 associated variants (-log 10 ( p ) ≥ 7.1) in 13 genomic loci on 9 chromosomes. The highest significance was found for the MGST1 region affecting milk fat content (-log 10 ( p ) = 11.93, MAF = 0.23, substitution effect of the minor allele (ß MA ) = -0.151%). Different from Holstein, DGAT1 was fixed (0.97) for the alanine protein variant for high milk and protein yield. A key gene affecting protein content was CSN1S1 (-log 10 ( p ) = 8.47, MAF = 049, ß MA = -0.055%) and the GNG2 region (-log 10 ( p ) = 10.48, MAF = 0.34, ß MA = 0.054%). Additionally, we suggest the importance of FGF12 for protein and fat yield, HTR3C for milk yield, TLE4 for milk and protein yield, and TNKS for milk and fat yield. Selection for favored alleles can improve milk yield and composition. With respect to maintaining the dual-purpose type of DSN, unfavored linkage to genes affecting muscularity has to be investigated carefully, before the milk-associated variants can be applied for selection in the small population.

Banning antibiotic growth promotors and other antimicrobials in poultry production due to the increasing antimicrobial resistance leads to increased feeding of potential alternatives such as probiotics. However, the modes of action of those feed additives are not entirely understood. They could act even with a direct effect on the immune system. A previously established animal-related in vitro system using primary cultured peripheral blood mononuclear cells (PBMCs) was applied to investigate the effects of immune-modulating feed additives. Here, the immunomodulation of different preparations of two probiotic Bacillus strains, B. subtilis DSM 32315 (BS), and B. amyloliquefaciens CECT 5940 (BA) was evaluated. The count of T-helper cells and activated T-helper cells increased after treatment in a ratio of 1:3 (PBMCs: Bacillus ) with vital BS (CD4+: p < 0.05; CD4+CD25+: p < 0.01). Furthermore, vital BS enhanced the proliferation and activation of cytotoxic T cells (CD8+: p < 0.05; CD8+CD25+: p < 0.05). Cell-free probiotic culture supernatants of BS increased the count of activated T-helper cells (CD4+CD25+: p < 0.1). UV-inactivated BS increased the proportion of cytotoxic T cells significantly (CD8+: p < 0.01). Our results point towards a possible involvement of secreted factors of BS in T-helper cell activation and proliferation, whereas it stimulates cytotoxic T cells presumably through surface contact. We could not observe any effect on B cells after treatment with different preparations of BS. After treatment with vital BA in a ratio of 1:3 (PBMCs: Bacillus ), the count of T-helper cells and activated T-helper cells increased (CD4+: p < 0.01; CD4+CD25+: p < 0.05). Cell-free probiotic culture supernatants of BA as well as UV-inactivated BA had no effect on T cell proliferation and activation. Furthermore, we found no effect of BA preparations on B cells. Overall, we demonstrate that the two different Bacillus strains enhanced T cell activation and proliferation, which points towards an immune-modulating effect of both strains on chicken immune cells in vitro. Therefore, we suggest that administering these probiotics can improve the cellular adaptive immune defense in chickens, thereby enabling the prevention and reduction of antimicrobials in chicken farming., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F ST ), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F Hom ) and genomic inbreeding (F RoH ) through runs of homozygosity (RoHs). Region-wide F ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F Hom and F RoH values in DSN were among the lowest. Regions with high F ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Neumann, Korkuć, Arends, Wolf, May, König and Brockmann.)

Body size is one of the main selection indices in chicken breeding. Although often investigated, knowledge of the underlying genetic mechanisms is incomplete. The aim of the current study was to identify genomic regions associated with body size differences between Asian Game and Asian Bantam type chickens. In this study, 94 and 107 chickens from 4 Asian Game and 5 Asian Bantam type breeds, respectively, were genotyped using the chicken 580K single nucleotide polymorphism (SNP) array. A genome-wide association study (GWAS) and principal component analyses (PCA) were performed to identify genomic regions associated with body size related-traits such as wing length, shank length, shank thickness, keel length, and body weight. Hierarchical clustering of genotype data showed a clear genetic difference between the investigated Asian Game and Asian Bantam chicken types. GWAS identified 16 genomic regions associated with wing length (2, FDR ≤ 0.018), shank thickness (6, FDR ≤ 0.008), keel length (5, FDR ≤ 0.023), and body weight (3, FDR ≤ 0.041). PCA showed that the first principal component (PC1) separated the 2 chicken types and significantly correlated with the measured body size related-traits (P ≤ 2.24e-40). SNPs contributing significantly to PC1 were subjected to a more detailed investigation. This analysis identified 11 regions potentially associated with differences in body size related-traits. A region on chromosome 4 (GGA4) (17.3-21.3 Mb) was detected in both analyses GWAS and PCA. This region harbors 60 genes. Among them are myotubularin 1 (MTM1) and secreted frizzled-related protein 2 (SFPR2) which can be considered as potential candidate genes for body size related-traits. Our results clearly show that the investigated Asian Game type chicken breeds are genetically different from the Asian Bantam breeds. A region on GGA4 between 17.3 and 21.3 Mb was identified which contributes to the phenotypic difference, though further validation of candidate genes is necessary., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Intrauterine exosomes have been identified to be involved in the embryo development and implantation. The aim of this study was to explore the role of miRNAs in intrauterine exosomes in bovine pregnancy. Intrauterine exosomes were collected from uterine flushing fluids of three donor and three recipient Xianan cows 7 days after fertilization. Intrauterine exosomes miRNAs were extracted and the exosomal miRNAs expression levels were analyzed. Sixty miRNAs differed significantly in their amounts between donors and recipients ( p-value < 0.05, |log2(FoldChange)| > 1). Twenty-two miRNAs were upregulated and 38 downregulated in the group of donor cows. The bta-miR-184 was the most significant ( P Benjamini-Hochberg < 0.001). A total of 9,775 target genes were predicted using the 60 miRNAs. GO and KEGG analysis showed that the target genes were enriched in several biological processes or pathways associated with embryo implantation and endometrial development, such as cell adhesion, cell junction, focal adhesion, and Rap1 signaling pathway. Our findings suggest that, in cattle early pregnancy stage, these differently expressed miRNAs in intrauterine exosomes involved in embryo implantation and endometrial development, which may exert a significant effect and influence the uterine microenvironment for embryo implantation. These results could provide reference for screening and exploring the intrauterine exosomal miRNA affecting embryo implantation., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Zhai, Shi, Chu, Chen, Feng, Zhang, Qi, Arends, Brockmann, Wang and Lyu.)

The milk casein genes in goats, are highly polymorphic genes with numerous synonymous and non-synonymous mutations. So far, 20 protein variants have been reported in goats for alpha-S1-casein, eight for beta-casein, 14 for alpha-S2-casein, and 24 for kappa-casein. This review provides a comprehensive overview on identified milk casein protein variants in goat and non-coding DNA sequence variants with some affecting the expression of the casein genes. The high frequency of some casein protein variants in different goat breeds and geographical regions might reflect specific breeding goals with respect to milk processing characteristics, properties for human nutrition and health, or adaptation to the environment. Because protein names, alongside the discovery of protein variants, go through a historical process, we linked old protein names with new ones that reveal more genetic variability. The haplotypes across the cluster of the four genetically linked casein genes are recommended as a valuable genetic tool for discrimination between breeds, managing genetic diversity within and between goat populations, and breeding strategies. The enormous variation in the casein proteins and genes is crucial for producing milk and dairy products with different properties for human health and nutrition, and for genetic improvement depending on local breeding goals., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Rahmatalla, Arends and Brockmann.)

This study examined (1) if fatty acids in bovine hair are influenced by dietary energy levels and (2) if the relationship between energy availability and fatty acids in hair persists across breeds and farms. Sixty-two and 59 Fleckvieh (Simmental), and 55 German Holstein cows from three farms, respectively, were fed two levels of energy concentration of roughage (6.1 and 6.5 MJ net energy for lactation/kg dry matter) and two levels of concentrate supply (150 and 250 g/kg energy-corrected milk). The average body weight was 727 kg (Simmental) and 668 kg (Holstein). The average lactation number was 3.1. Hair samples were taken in lactation weeks 4 and 8. In Simmental cows, a lower energy deficit due to a relatively higher energy intake from high energy concentration of the roughage was associated with higher C18:2 n -6 and C18:3 n -3 contents in hair at week 8. In cows from all three farms, higher energy intake between lactation weeks 2 and 6 correlated with higher content of C18:2 n -6 in hair samples taken in lactation weeks 4 and 8. No correlation was found for C12:0. These results provide the first evidence that increased energy intake increases the contents of C18:2 n -6 in hair.

Dietary fiber (DF) is receiving increasing attention, and its importance in pig nutrition is now acknowledged. Although DF for pigs was frowned upon for a long time because of reductions in energy intake and digestibility of other nutrients, it has become clear that feeding DF to pigs can affect their well-being and health. This review aims to summarize the state of knowledge of studies on DF in pigs, with an emphasis on the underlying mode of action, by considering research using DF in sows as well as suckling and weaned piglets, and fattening pigs. These studies indicate that DF can benefit the digestive tracts and the health of pigs, if certain conditions or restrictions are considered, such as concentration in the feed and fermentability. Besides the chemical composition and the impact on energy and nutrient digestibility, it is also necessary to evaluate the possible physical and physiologic effects on intestinal function and intestinal microbiota, to better understand the relation of DF to animal health and welfare. Future research should be designed to provide a better mechanistic understanding of the physiologic effects of DF in pigs.

The Bardet-Biedl Syndrome 7 ( Bbs7 ) gene was identified as the most likely candidate gene causing juvenile obesity in the Berlin Fat Mouse Inbred (BFMI) line. Bbs7 expression is significantly lower in the brain, adipose tissue, and liver of BFMI mice compared to lean C57BL/6NCrl (B6N) mice. A DNA sequence comparison between BFMI and B6N revealed 16 sequence variants in the Bbs7 promoter region. Here, we tested if these mutations contribute to the observed differential expression of Bbs7 . In a cell-based dual-luciferase assay, we compared the effects of the BFMI and the B6N haplotypes of different regions of the Bbs7 promotor on the reporter gene expression. A single-nucleotide polymorphism (SNP) was identified causing a significant reduction in the reporter gene expression. This SNP (rs29947545) is located in the 5' UTR of Bbs7 at Chr3:36.613.350. The SNP is not unique to BFMI mice but also occurs in several other mouse strains, where the BFMI allele is not associated with lower Bbs7 transcript amounts. Thus, we suggest a compensatory mutation in the other mouse strains that keeps Bbs7 expression at the normal level. This compensatory mechanism is missing in BFMI mice and the cell lines tested.

The Berlin Fat Mouse Inbred (BFMI) line is a model for juvenile obesity. Previous studies on crosses between BFMI and C57Bl/6N (B6N) have identified a recessive defect causing juvenile obesity on chromosome 3 (jObes1). Bbs7 was identified as the most likely candidate gene for the observed effect. Comparative sequence analysis showed a 1578 bp deletion in intron 8 of Bbs7 in BFMI mice. A CTCF-element is located inside this deletion. To investigate the functional effect of this deletion, it was introduced into B6N mice using CRISPR/Cas9. Two mice containing the target deletion were obtained (B6N Bbs7 emI8∆1 and Bbs7 emI8∆2 ) and were subsequently mated to BFMI and B6N to generate two families suitable for complementation. Inherited alleles were determined and body composition was measured by quantitative magnetic resonance. Evidence for a partial complementation (13.1-15.1%) of the jObes1 allele by the CRISPR/Cas9 modified B6N Bbs7 emI8∆1 and Bbs7 emI8∆2 alleles was found. Mice carrying the complementation alleles had a 23-27% higher fat-to-lean ratio compared to animals which have a B6N allele (P (Bbs7emI8∆1)  = 4.25 × 10 -7 ; P (Bbs7emI8∆2)  = 3.17 × 10 -5 ). Consistent with previous findings, the recessive effect of the BFMI allele was also seen for the B6N Bbs7 emI8∆1 and Bbs7 emI8∆2 alleles. However, the effect size of the B6N Bbs7 emI8∆1 and Bbs7 emI8∆2 alleles was smaller than the BFMI allele, and thus showed only a partial complementation. Findings suggest additional variants near Bbs7 in addition to or interacting with the deletion in intron 8., (© 2021. The Author(s).)

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with liver weight, liver triglycerides, and body weight using the obese BFMI sub-line BFMI861-S1. BFMI861-S1 mice are insulin resistant and store ectopic fat in the liver. In generation 10, 58 males and 65 females of the advanced intercross line (AIL) BFMI861-S1xB6N were phenotyped under a standard diet over 20 weeks. QTL analysis was performed after genotyping with the MiniMUGA Genotyping Array. Whole-genome sequencing and gene expression data of the parental lines was used for the prioritization of positional candidate genes. Three QTLs associated with liver weight, body weight, and subcutaneous adipose tissue (scAT) weight were identified. A highly significant QTL on chromosome (Chr) 1 (157-168 Mb) showed an association with liver weight. A QTL for body weight at 20 weeks was found on Chr 3 (34.1-40 Mb) overlapping with a QTL for scAT weight. In a multiple QTL mapping approach, an additional QTL affecting body weight at 16 weeks was identified on Chr 6 (9.5-26.1 Mb). Considering sequence variants and expression differences, Sec16b and Astn1 were prioritized as top positional candidate genes for the liver weight QTL on Chr 1; Met and Ica1 for the body weight QTL on Chr 6. Interestingly, all top candidate genes have previously been linked with metabolic traits. This study shows once more the power of an advanced intercross line for fine mapping. QTL mapping combined with a detailed prioritization approach allowed us to identify additional and plausible candidate genes linked to metabolic traits in the BFMI861-S1xB6N AIL. By reidentifying known candidate genes in a different crossing population the causal link with specific traits is underlined and additional evidence is given for further investigations., (© 2022. The Author(s).)

The improvement of milk production of indigenous Sudanese cattle such as Bos indicus Butana and its cross with Holstein is a major goal of the Sudanese government to ensure sufficient healthy nutrition in the country. In this study, we investigated the K232A polymorphism of diacylglycerol acyltransferase (DGAT1), a well-known modulator of milk production in other breeds. We determined allele frequencies and the allele effects on milk production. Therefore, 93 purebred Butana and 203 Butana × Holstein crossbred cattle were genotyped using competitive allele-specific PCR assays. Association analysis was performed using a linear mixed model in R. In purebred Butana cattle, the lysine DGAT1 protein variant K232, which is found to be associated with higher fat and protein contents, as well as higher fat yield was highly frequent at 0.929, while its frequency in Butana × Holstein crossbred cattle was 0.394. Significant effects were found on milk yield (P = 7.6 × 10 -20 ), fat yield (P = 2.2 × 10 -17 ), protein yield (P = 2.0 × 10 -19 ) and lactose yield (P = 4.0 × 10 -18 ) in crossbred cattle. As expected, the protein variant K232 was disadvantageous since it was decreasing milk, protein, and lactose yields by 1.741 kg, 0.063 kg and 0.084 kg, respectively. No significant effects were found for milk fat, protein, and lactose contents. The high frequency of the lysine DGAT1 protein variant K232 in Butana cattle could contribute to their high milk fat content in combination with low milk yield. In Butana × Holstein crossbred cattle, the DGAT1 marker can be used for effective selection and thus genetic improvement of milk production., (© 2022. The Author(s).)

Background: The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with impaired glucose metabolism using the obese lines BFMI861-S1 and BFMI861-S2, which are genetically closely related, but differ in several traits. BFMI861-S1 is insulin resistant and stores ectopic fat in the liver, whereas BFMI861-S2 is insulin sensitive., Methods: In generation 10, 397 males of an advanced intercross line (AIL) BFMI861-S1 × BFMI861-S2 were challenged with a high-fat, high-carbohydrate diet and phenotyped over 25 weeks. QTL-analysis was performed after selective genotyping of 200 mice using the GigaMUGA Genotyping Array. Additional 197 males were genotyped for 7 top SNPs in QTL regions. For the prioritization of positional candidate genes whole genome sequencing and gene expression data of the parental lines were used., Results: Overlapping QTL for gonadal adipose tissue weight and blood glucose concentration were detected on chromosome (Chr) 3 (95.8-100.1 Mb), and for gonadal adipose tissue weight, liver weight, and blood glucose concentration on Chr 17 (9.5-26.1 Mb). Causal modeling suggested for Chr 3-QTL direct effects on adipose tissue weight, but indirect effects on blood glucose concentration. Direct effects on adipose tissue weight, liver weight, and blood glucose concentration were suggested for Chr 17-QTL. Prioritized positional candidate genes for the identified QTL were Notch2 and Fmo5 (Chr 3) and Plg and Acat2 (Chr 17). Two additional QTL were detected for gonadal adipose tissue weight on Chr 15 (67.9-74.6 Mb) and for body weight on Chr 16 (3.9-21.4 Mb)., Conclusions: QTL mapping together with a detailed prioritization approach allowed us to identify candidate genes associated with traits of the metabolic syndrome. In addition, we provided evidence for direct and indirect genetic effects on blood glucose concentration in the insulin-resistant mouse line BFMI861-S1., (© 2021. The Author(s).)

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

The Bos indicus zebu cattle Butana is the most commonly used indigenous dairy cattle breed in Sudan. In the last years, high-yielding Holstein dairy cattle were introgressed into Butana cattle to improve their milk yield and simultaneously keep their good adaption to extreme environmental conditions. With the focus on the improvement of milk production, other problems arose such as an increased susceptibility to mastitis. Thus, genetic selection for mastitis resistance should be considered to maintain healthy and productive cows. In this study, we tested 10 single nucleotide polymorphisms (SNPs) which had been associated with somatic cell score (SCS) in Holstein cattle for association with SCS in 37 purebred Butana and 203 Butana × Holstein crossbred cattle from Sudan. Animals were genotyped by competitive allele-specific PCR assays and association analysis was performed using a linear mixed model. All 10 SNPs were segregating in the crossbred Butana × Holstein populations, but only 8 SNPs in Sudanese purebred Butana cattle. The SNP on chromosome 13 was suggestively associated with SCS in the Butana × Holstein crossbred population (rs109441194, 13:79,365,467, P BF  = 0.054) and the SNP on chromosome 19 was significantly associated with SCS in both populations (rs41257403, 19:50,027,458, Butana: P BF  = 0.003, Butana × Holstein: P BF  = 6.2 × 10 -16 ). The minor allele of both SNPs showed an increase in SCS. Therefore, selection against the disadvantageous minor allele could be used for genetic improvement of mastitis resistance in the studied populations. However, investigations in a bigger population and across the whole genome are needed to identify additional genomic loci., (© 2022. The Author(s).)

While direct additive and dominance effects on complex traits have been mapped repeatedly, additional genetic factors contributing to the heterogeneity of complex traits have been scarcely investigated. To assess genetic background effects, we investigated transmission ratio distortions (TRDs) of alleles from parent to offspring using an advanced intercross line (AIL) of an initial cross between the mouse inbred strains C57BL/6NCrl (B6N) and BFMI860-12 [Berlin Fat Mouse Inbred (BFMI)]. A total of 341 males of generation 28 and their respective 61 parents and 66 grandparents were genotyped using Mega Mouse Universal Genotyping Arrays. TRDs were investigated using allele transmission asymmetry tests, and pathway overrepresentation analysis was performed. Sequencing data were used to test for overrepresentation of nonsynonymous SNPs (nsSNPs) in TRD regions. Genetic incompatibilities were tested using the Bateson-Dobzhansky-Muller two-locus model. A total of 62 TRD regions were detected, many in close proximity to the telocentric centromere. TRD regions contained 44.5% more nsSNPs than randomly selected regions (182 vs 125.9 ± 17.0, P < 1 × 10-4). Testing for genetic incompatibilities between TRD regions identified 29 genome-wide significant incompatibilities between TRD regions [P(BF) < 0.05]. Pathway overrepresentation analysis of genes in TRD regions showed that DNA methylation, epigenetic regulation of RNA, and meiotic/meiosis regulation pathways were affected independent of the parental origin of the TRD. Paternal BFMI TRD regions showed overrepresentation in the small interfering RNA biogenesis and in the metabolism of lipids and lipoproteins. Maternal B6N TRD regions harbored genes involved in meiotic recombination, cell death, and apoptosis pathways. The analysis of genes in TRD regions suggests the potential distortion of protein-protein interactions influencing obesity and diabetic retinopathy as a result of disadvantageous combinations of allelic variants in Aass, Pgx6, and Nme8. Using an AIL significantly improves the resolution at which we can investigate TRD. Our analysis implicates distortion of protein-protein interactions as well as meiotic drive as the underlying mechanisms leading to the observed TRD in our AIL. Furthermore, genes with large amounts of nsSNPs located in TRD regions are more likely to be involved in pathways that are related to the phenotypic differences between the parental strains. Genes in these TRD regions provide new targets for investigating genetic adaptation, protein-protein interactions, and determinants of complex traits such as obesity., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Knowledge about the modes of action of immunomodulating compounds such as pathogens, drugs, or feed additives, e.g., probiotics, gained through controlled but animal-related in vitro systems using primary cultured peripheral blood mononuclear cells (PBMCs) will allow the development of targeted nutrition strategies. Moreover, it could contribute to the prevention of infectious diseases and the usage of antimicrobials, and further promote the health of the animals. However, to our knowledge, a protocol for the isolation of PBMCs with reduced thrombocyte count from chicken blood and subsequent cell culture over several days to assess the effects of immunomodulating compounds is not available. Therefore, we established an optimized protocol for blood sampling and immune cell isolation, culture, and phenotyping for chicken PBMCs. For blood sampling commercial Na-citrate tubes revealed the highest count of vital cells compared to commercial Li-heparin ( p < 0.01) and K3EDTA ( p < 0.05) tubes. Using combined dextran and ficoll density gradient separation, the thrombocyte count was significantly reduced ( p < 0.01) compared to slow-speed centrifugation with subsequent ficoll. For cell culture, the supplementation of RPMI-1640 medium with 10% chicken serum resulted in the lowest relative cell count of thrombocytes compared to fetal calf serum (FCS) ( p < 0.05). To validate the ability of the cell culture system to respond to stimuli, concanavalin A (conA) was used as a positive control. The optimized protocol allows the isolation and cultivation of vital PBMCs with reduced thrombocyte count from chicken blood for subsequent investigation of the modes of action of immunomodulating compounds.

Background: German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement., Results: Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%., Conclusion: The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future., (© 2021. The Author(s).)

This genome-wide association study (GWAS) aimed to identify sequence variants (SVs) and candidate genes associated with fertility and health in endangered German Black Pied cattle (DSN) based on whole-genome sequence (WGS) data. We used 304 sequenced DSN cattle for the imputation of 1797 genotyped DSN to WGS. The final dataset included 11,413,456 SVs of 1886 cows. Cow traits were calving-to-first service interval (CTFS), non-return after 56 days (NR56), somatic cell score (SCS), fat-to-protein ratio (FPR), and three pre-corrected endoparasite infection traits. We identified 40 SVs above the genome-wide significance and suggestive threshold associated with CTFS and NR56, and three important potential candidate genes ( ARHGAP21 , MARCH11 , and ZNF462 ). For SCS, most associations were observed on BTA 25. The GWAS revealed 61 SVs, a cluster of 10 candidate genes on BTA 13, and 7 pathways for FPR, including key mediators involved in milk fat synthesis. The strongest associations for gastrointestinal nematode and Dictyocaulus viviparus infections were detected on BTA 8 and 24, respectively. For Fasciola hepatica infections, the strongest associated SVs were located on BTA 4 and 7. We detected 200 genes for endoparasite infection traits, related to 16 pathways involved in host immune response during infection.

Background/objectives: There is a growing appreciation for individual responses to diet. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed. In this study, we searched for genetic variants underlying differences in the responses to American and ketogenic diets between C57BL/6J (B6) and FVB/NJ (FVB) mouse strains., Results: Genetic mapping of fat and lean mass gain revealed QTLs on Chromosome (Chr) 1 at 191.6 Mb (Fmgq1) (P < 0.001, CI = 180.2-194.4 Mb), Chr5 at 73.7 Mb (Fmgq2, Lmgq1) (P < 0.001, CI = 66.1-76.6 Mb), and Chr7 at 40.5 Mb (Fmgq3) (P < 0.01, CI = 36.6-44.5 Mb). Analysis of serum HDL cholesterol concentration identified a significant (P < 0.001, CI = 160.6-176.1 Mb) QTL on Chr1 at 168.6 Mb (Hdlq1). Causal network inference suggests that HDL cholesterol and fat mass gain are both linked to Fmgq1., Conclusions: Strong sex effects were identified at both Fmgq2 and Lmgq1, which are also diet-dependent. Interestingly, Fmgq2 and Fmgq3 affect fat gain directly, while Fmgq1 influences fat gain directly and via an intermediate change in serum cholesterol. These results demonstrate how precision nutrition will be advanced through the integration of genetic variation and sex in physiological responses to diets varied in carbohydrate composition.

German Black Pied cattle (DSN) is an endangered population of about 2,550 dual-purpose cattle in Germany. Having a milk yield of about 2,500 kg less than the predominant dairy breed Holstein, the preservation of DSN is supported by the German government and the EU. The identification of the genomic loci affecting milk production in DSN can provide a basis for selection decisions for genetic improvement of DSN in order to increase market chances through the improvement of milk yield. A genome-wide association analysis of 30 milk traits was conducted in different lactation periods and numbers. Association using multiple linear regression models in R was performed on 1,490 DSN cattle genotyped with BovineSNP50 SNP-chip. 41 significant and 20 suggestive SNPs affecting milk production traits in DSN were identified, as well as 15 additional SNPs for protein content which are less reliable due to high inflation. The most significant effects on milk yield in DSN were detected on chromosomes 1, 6, and 20. The region on chromosome 6 was located nearby the casein gene cluster and the corresponding haplotype overlapped the CSN3 gene (casein kappa). Associations for fat and protein yield and content were also detected. High correlation between traits of the same lactation period or number led to some SNPs being significant for multiple investigated traits. Half of all identified SNPs have been reported in other studies, previously. 15 SNPs were associated with the same traits in other breeds. The other associated SNPs have been reported previously for traits such as exterior, health, meat and carcass, production, and reproduction traits. No association could be detected between DGAT1 and other known milk genes with milk production traits despite the close relationship between DSN and Holstein. The results of this study confirmed that many SNPs identified in other breeds as associated with milk traits also affect milk traits in dual-purpose DSN cattle and can be used for further genetic analysis to identify genes and causal variants that affect milk production in DSN cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Korkuć, Arends, May, König and Brockmann.)

Genetic variations in the four casein genes CSN1S1 , CSN2 , CSN1S2 , and CSN3 have obtained substantial attention since they affect the milk protein yield, milk composition, cheese processing properties, and digestibility as well as tolerance in human nutrition. Furthermore, milk protein variants are used for breed characterization, biodiversity, and phylogenetic studies. The current study aimed at the identification of casein protein variants in five domestic goat breeds from Sudan (Nubian, Desert, Nilotic, Taggar, and Saanen) and three wild goat species [ Capra aegagrus aegagrus (Bezoar ibex), Capra nubiana (Nubian ibex), and Capra ibex (Alpine ibex)]. High-density capture sequencing of 33 goats identified in total 22 non-synonymous and 13 synonymous single nucleotide polymorphisms (SNPs), of which nine non-synonymous and seven synonymous SNPs are new. In the CSN1S1 gene, the new non-synonymous SNP ss7213522403 segregated in Alpine ibex. In the CSN2 gene, the new non-synonymous SNPs ss7213522526, ss7213522558, and ss7213522487 were found exclusively in Nubian and Alpine ibex. In the CSN1S2 gene, the new non-synonymous SNPs ss7213522477, ss7213522549, and ss7213522575 were found in Nubian ibex only. In the CSN3 gene, the non-synonymous SNPs ss7213522604 and ss7213522610 were found in Alpine ibex. The identified DNA sequence variants led to the detection of nine new casein protein variants. New variants were detected for alpha S1 casein in Saanen goats ( CSN1S1 ∗ C1), Bezoar ibex ( CSN1S1 ∗ J), and Alpine ibex ( CSN1S1 ∗ K), for beta and kappa caseins in Alpine ibex ( CSN2 ∗ F and CSN3 ∗ X), and for alpha S2 casein in all domesticated and wild goats ( CSN1S2 ∗ H), in Nubian and Desert goats ( CSN1S2 ∗ I), or in Nubian ibex only ( CSN1S2 ∗ J and CSN1S2 ∗ K). The results show that most novel SNPs and protein variants occur in the critically endangered Nubian ibex. This highlights the importance of the preservation of this endangered breed. Furthermore, we suggest validating and further characterizing the new casein protein variants., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Rahmatalla, Arends, Said Ahmed, Hassan, Krebs, Reissmann and Brockmann.)

The laboratory mouse is the most widely used animal model for biomedical research, due in part to its well-annotated genome, wealth of genetic resources, and the ability to precisely manipulate its genome. Despite the importance of genetics for mouse research, genetic quality control (QC) is not standardized, in part due to the lack of cost-effective, informative, and robust platforms. Genotyping arrays are standard tools for mouse research and remain an attractive alternative even in the era of high-throughput whole-genome sequencing. Here, we describe the content and performance of a new iteration of the Mouse Universal Genotyping Array (MUGA), MiniMUGA, an array-based genetic QC platform with over 11,000 probes. In addition to robust discrimination between most classical and wild-derived laboratory strains, MiniMUGA was designed to contain features not available in other platforms: (1) chromosomal sex determination, (2) discrimination between substrains from multiple commercial vendors, (3) diagnostic SNPs for popular laboratory strains, (4) detection of constructs used in genetically engineered mice, and (5) an easy-to-interpret report summarizing these results. In-depth annotation of all probes should facilitate custom analyses by individual researchers. To determine the performance of MiniMUGA, we genotyped 6899 samples from a wide variety of genetic backgrounds. The performance of MiniMUGA compares favorably with three previous iterations of the MUGA family of arrays, both in discrimination capabilities and robustness. We have generated publicly available consensus genotypes for 241 inbred strains including classical, wild-derived, and recombinant inbred lines. Here, we also report the detection of a substantial number of X O and XXY individuals across a variety of sample types, new markers that expand the utility of reduced complexity crosses to genetic backgrounds other than C57BL/6, and the robust detection of 17 genetic constructs. We provide preliminary evidence that the array can be used to identify both partial sex chromosome duplication and mosaicism, and that diagnostic SNPs can be used to determine how long inbred mice have been bred independently from the relevant main stock. We conclude that MiniMUGA is a valuable platform for genetic QC, and an important new tool to increase the rigor and reproducibility of mouse research., (Copyright © 2020 by the Genetics Society of America.)

The dual-purpose German Black Pied Cattle (DSN) has become an endangered breed of approximately 2,550 registered cows in Germany. The breed is genetically related to Holstein-Friesian cattle because the old DSN breed contributed to the selection of the modern Holstein dairy cow. In dairy farms, breeders aim to improve animal health and well-being by reducing the number of mastitis cases, which would also reduce milk losses and treatment costs. On the genomic level, no markers associated with clinical mastitis have been reported in DSN. Therefore, we performed a genome-wide association study on 1,062 DSN cows using a univariate linear mixed model that included a relatedness matrix to correct for population stratification. Although the statistical power was limited by the small population size, 3 markers were significantly associated, and 2 additional markers showed a suggestive association with clinical mastitis. Those markers accounted for 1 to 3% of the variance of clinical mastitis in the examined DSN population. One marker was found in the intragenic region of NEURL1 on BTA26, and the other 4 markers in intergenic regions on BTA3, BTA6, and BTA9. Further analyses identified 23 positional candidate genes. Among them is BMPR1B, which has been previously associated with clinical mastitis in other dairy cattle breeds. The markers presented here can be used for selection for mastitis-resistant animals in the endangered DSN population, and can broadly contribute to a better understanding of mastitis determinants in dairy cattle breeds., (Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Lignocellulose is a constituent of plant cell walls and might be used as a fiber source in poultry nutrition. The current study investigated the impact of increasing dietary levels of lignocellulose on performance, nutrient digestibility, excreta DM, intestinal microbiota, and bacterial metabolites in slow growing broilers. At an age of 10 wk, 60 male broilers of an intercross line (New Hampshire × White Leghorn) were allocated to cages and fed isoenergetic and isonitrogenous diets containing 0.8% (LC1), 5% (LC2), or 10% (LC3) lignocellulose. After 23 D of feeding, broilers were killed and digesta samples of ileum and excreta analyzed for nutrient digestibility and DM. Cecal contents were analyzed for microbial composition and metabolites. Broiler performance was not affected by feeding dietary lignocellulose. LC3 fed broilers showed reduced ileal digestibility of protein compared to chickens fed LC1 (P = 0.003). Moreover, increasing levels of dietary lignocellulose reduced apparent digestibility of organic matter and gross energy (P < 0.001). Feeding of lignocellulose had no impact on the excreta DM of broilers. Increasing levels of dietary lignocellulose lowered cecal counts of Escherichia/Hafnia/Shigella (P = 0.029) and reduced the total concentration of short-chain fatty acids (P < 0.001), lactate (P < 0.05), and ammonia (P = 0.009). The molar ratio of cecal acetic acid was higher in LC3 fed broilers (P < 0.001), while the proportions of cecal propionic acid and n-butyric acid were higher in LC1 and LC2 fed chickens (P < 0.001). Correlation analyses indicated that dietary lignocellulose was negatively related to the total concentration of cecal bacterial metabolites (P < 0.001). In conclusion, the feeding of lignocellulose did not affect growth performance, but impaired nutrient digestibility of slow growing broilers. While minor changes in cecal microbial composition were detected, cecal bacterial metabolite concentrations were significantly reduced with increasing levels of dietary lignocellulose. These findings suggest that lignocellulose is not extensively degraded by bacteria residing in the large intestine of broilers., (Copyright © 2020. Published by Elsevier Inc.)

The aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual-purpose black and white (DSN) population, as well as from generated sub-populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East-DSN and West-DSN sub-populations. In addition, two further sub-populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross-population extended haplotype homozygosity methodology (XP-EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright's fixation index (F ST ). Chromosomal segments for the top 0.1 percentile of negative or positive XP-EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP-EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP-EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29-38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average F ST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis., (© 2020 Stichting International Foundation for Animal Genetics.)

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

Casein proteins were repeatedly examined for protein polymorphisms and frequencies in diverse cattle breeds. The occurrence of casein variants in Holstein Friesian, the leading dairy breed worldwide, is well known. The frequencies of different casein variants in Holstein are likely affected by selection for high milk yield. Compared to Holstein, only little is known about casein variants and their frequencies in German Black Pied cattle ("Deutsches Schwarzbuntes Niederungsrind," DSN). The DSN population was a main genetic contributor to the current high-yielding Holstein population. The goal of this study was to investigate casein (protein) variants and casein haplotypes in DSN based on the DNA sequence level and to compare these with data from Holstein and other breeds. In the investigated DSN population, we found no variation in the alpha-casein genes CSN1S1 and CSN1S2 and detected only the CSN1S1*B and CSN1S2*A protein variants. For CSN2 and CSN3 genes, non-synonymous single nucleotide polymorphisms leading to three different β and κ protein variants were found, respectively. For β-casein protein variants A 1 , A 2 , and I were detected, with CSN2*A 1 (82.7%) showing the highest frequency. For κ-casein protein variants A , B , and E were detected in DSN, with the highest frequency of CSN3*A (83.3%). Accordingly, the casein protein haplotype CSN1S1*B-CSN2*A 1 - CSN1S2*A - CSN3*A (order of genes on BTA6) is the most frequent haplotype in DSN cattle., (Copyright © 2019 Meier, Korkuć, Arends and Brockmann.)

Background: Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections., Results: A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (p Bonf  = 4.47 × 10 - 7 ), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10 - 6 on BTA 1 to 2.18 × 10 - 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections., Conclusions: We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records.

The imputation from lower density SNP chip genotypes to whole-genome sequence level is an established approach to generate high density genotypes for many individuals. Imputation accuracy is dependent on many factors and for small cattle populations such as the endangered German Black Pied cattle (DSN), determining the optimal imputation strategy is especially challenging since only a low number of high density genotypes is available. In this paper, the accuracy of imputation was explored with regard to (1) phasing of the target population and the reference panel for imputation, (2) comparison of a 1-step imputation approach, where 50 k genotypes are directly imputed to sequence level, to a 2-step imputation approach that used an intermediate step imputing first to 700 k and subsequently to sequence level, (3) the software tools Beagle and Minimac, and (4) the size and composition of the reference panel for imputation. Analyses were performed for 30 DSN and 30 Holstein Frisian cattle available from the 1000 Bull Genomes Project. Imputation accuracy was assessed using a leave-one-out cross validation procedure. We observed that phasing of the target populations and the reference panels affects the imputation accuracy significantly. Minimac reached higher accuracy when imputing using small reference panels, while Beagle performed better with larger reference panels. In contrast to previous research, we found that when a low number of animals is available at the intermediate imputation step, the 1-step imputation approach yielded higher imputation accuracy compared to a 2-step imputation. Overall, the size of the reference panel for imputation is the most important factor leading to higher imputation accuracy, although using a larger reference panel consisting of a related but different breed (Holstein Frisian) significantly reduced imputation accuracy. Our findings provide specific recommendations for populations with a limited number of high density genotyped or sequenced animals available such as DSN. The overall recommendation when imputing a small population are to (1) use a large reference panel of the same breed, (2) use a large reference panel consisting of diverse breeds, or (3) when a large reference panel is not available, we recommend using a smaller same breed reference panel without including a different related breed.

We hypothesised that a relationship would exist between hair fatty acids, especially C12:0, C14:0 and C16:0, and parameters of energy metabolism such as energy intake, energy mobilisation, and energy requirement for maintenance and milk performance. For this study, 11 primiparous German Holstein cows were available from which hair samples at weeks 6 and 8 of lactation were analysed. The average body weight of these animals was 558 ± 27 kg at calving and milk yield at 100-days in milk was 3,537 ± 529 kg. Feed intake and milk yield were measured daily. Body weight and back fat thickness were measured at calving and in weeks 2, 4, and 8 of lactation. Energy balance and energy utilisation were calculated until week 6 of lactation. Spearman's correlation coefficients were found to be significantly positive for the relationship between the percentage of C12:0 and C14:0 fatty acids in the hair in lactation week 8 and energy intake in weeks 5 and 6 (0.62 &lt; r &lt; 0.65, P &lt; 0.05). If the animals are grouped according to their energy utilisation between weeks 1 and 6 into two groups higher (n = 6) or lower (n = 5) than the median, animals of the high energy utilising group had a higher energy intake. These animals had also higher percentages of the C12:0 fatty acid in their hair fat (week 6: 4.9% vs. 3.1%, P &lt; 0.05; week 8: 4.3% vs. 2.9%, P = 0.05). Our hypothesis is supported, and this study justifies further investigation of the content of medium-chain fatty acids in hair samples as biomarkers for the metabolic status of a cow during early lactation.

Zn serves as a powerful feed additive to reduce post-weaning diarrhoea in pigs. However, the mechanisms responsible for Zn-associated effects on the adaptive immune responses following feeding of a very high dosage of Zn remain elusive. In this study, we examined the T-cell response in gut-associated lymphatic tissues of seventy-two weaned piglets. Piglets received diets with 57 mg Zn/kg (low Zn concentration, LZn), 164 mg Zn/kg (medium Zn concentration, MZn) or 2425 mg Zn/kg (high Zn concentration, HZn) mg Zn/kg feed for 1, 2 or 4 weeks. We observed that feeding the HZn diet for 1 week increased the level of activated T-helper cells (CD4+ and CD8α dim) compared with feeding MZn and LZn (P<0·05). In addition, we observed higher transcript amounts of interferon γ and T-box 21 (TBET) in the HZn group compared with the MZn and LZn groups (P<0·05). A gene set enrichment analysis revealed an over-representation of genes associated with 'cytokine signalling in immune system'. Remarkably, feeding of a very high Zn dosage led to a switch in the immune response after 2 weeks. We detected higher relative cell counts of CD4+CD25high regulatory T-helper cells (P<0·05) and a higher expression of forkhead box P3 (FOXP3) transcripts (P<0·05). After 4 weeks of feeding a high-dosage Zn diet, the relative CD4+ T-cell count (P<0·05) and the relative CD8β + T-cell count (P<0·1) were reduced compared with the MZn group. We hypothesise that after 1 week the cellular T-helper 1 response is switched on and after 2 weeks it is switched off, leading to decreased numbers of T-cells.

Long-term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome-wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185-scaffold1365-620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482-scaffold89-467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1-Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes., (© 2018 Stichting International Foundation for Animal Genetics.)

In our previous research, we identified a QTL with an interval of 3.4 Mb for growth on chicken chromosome (GGA) 4 in an advanced intercross population of an initial cross between the New Hampshire inbred line (NHI) and the White Leghorn inbred line (WL77). In the current study, an association analysis was performed in a population of purebred white layers (WLA) with White Leghorn origin. Genotypic data of 130 SNPs within the previously identified 3.4-Mb region were obtained using a 60K SNP chip. In total, 24 significant SNPs (LOD ≥ 4.44) on GGA4 were detected for daily weigh gain from 8 to 14 weeks and two SNPs (LOD ≥ 4.80) for body weight at 14 weeks. The QTL interval was reduced by 1.9 Mb to an interval of 1.5 Mb (74.6-76.1 Mb) that harbors 15 genes. Furthermore, to identify additional loci for chicken growth, a genome-wide association study (GWAS) was carried out in a WLA population. The GWAS identified an additional QTL on GGA6 for body weight at six weeks (19.8-21.2 Mb). Our findings showed that by using a WLA population we were able to further reduce the QTL confidence interval previously detected using a NHI × WL77 advanced intercross population., (© 2018 Stichting International Foundation for Animal Genetics.)

Acyl-CoA synthetases (ACSs) are responsible for acyl-CoA synthesis from nonpolar hydrophilic fatty acids and play a vital role in many metabolic processes. As a category of ACS isozymes, members of ACS family (ACSF1-3) participate in lipid metabolism; however, their expression patterns, regulatory mechanisms and effects on egg-laying performance in chicken are poorly understood. Our in vivo and in vitro studies showed that ACSF1-3 genes were extensively expressed, and their expression levels changed dynamically in the liver among different development stages. Moreover, ACSF1 expression was upregulated and ACSF2 expression was downregulated by estrogen, but ACSF3 showed no response to estrogen treatment. The regulatory effect of estrogen on ACSF1 expression was mediated via ERα. The ACSF2 was highly expressed in the liver in peak-laying hens compared with pre-laying and late-laying hens, and also highly expressed in the liver continued egg-laying hens compared with inactive egg-laying hens. It is suggested that hepatic ACSF2 expression level might relate to egg-laying performance in chicken. In conclusion, the expression of ACSF1 was upregulated by estrogen via ERα, and the expression of ACSF2 was downregulated by estrogen and might be related to egg-laying performance in chicken.