Allele substitution and dominance effects of CD166/ALCAM gene polymorphisms for endoparasite resistance and test-day traits in a small cattle population using logistic regression analyses

Selective Sweeps in Cattle Genomes in Response to the Influence of Urbanization and Environmental Contamination

A genomic study was conducted to identify the effects of urbanization and environmental contaminants with heavy metals on selection footprints in dairy cattle populations reared in the megacity of Bengaluru, South India. Dairy cattle reared along the rural-urban interface of Bengaluru with/without access to roughage from public lakeshores were selected. The genotyped animals were subjected to the cross-population-extended haplotype homozygosity (XP-EHH) methodology to infer selection sweeps caused by urbanization (rural, mixed, and urban) and environmental contamination with cadmium and lead. We postulated that social-ecological challenges contribute to mechanisms of natural selection. A number of selection sweeps were identified when comparing the genomes of cattle located in rural, mixed, or urban regions. The largest effects were identified on BTA21, displaying pronounced peaks for selection sweeps for all three urbanization levels (urban_vs_rural, urban_vs_mixed and rural_vs_mixed). Selection sweeps are located in chromosomal segments in close proximity to the genes lrand rab interactor 3 (RIN3), solute carrier family 24 member 4 (SLC24A4), tetraspanin 3 (TSPAN3), and proline-serine-threonine phosphatase interacting protein 1 (PSTPIP1). Functional enrichment analyses of the selection sweeps for all three comparisons revealed a number of gene ontology (GO) and KEGG terms, which were associated with reproduction, metabolism, and cell signaling-related functional mechanisms. Likewise, a number of the chromosomal segments under selection were observed when creating cattle groups according to cadmium and lead contaminations. Stronger and more intense positive selection sweeps were observed for the cadmium contaminated group, i.e., signals of selection on BTA 16 and BTA19 in close proximity to genes regulating the somatotropic axis (growth factor receptor bound protein 2 (GRB2) and cell ion exchange (chloride voltage-gated channel 6 (CLCN6)). A few novel, so far uncharacterized genes, mostly with effects on immune physiology, were identified. The lead contaminated group revealed sweeps which were annotated with genes involved in carcass traits (TNNC2, SLC12A5, and GABRA4), milk yield (HTR1D, SLCO3A1, TEK, and OPCML), reproduction (GABRA4), hypoxia/stress response (OPRD1 and KDR), cell adhesion (PCDHGC3), inflammatory response (ADORA2A), and immune defense mechanism (ALCAM). Thus, the findings from this study provide a deeper insight into the genomic regions under selection under the effects of urbanization and environmental contamination.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN)

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.

Genome-Wide Association Study Using Whole-Genome Sequence Data for Fertility, Health Indicator, and Endoparasite Infection Traits in German Black Pied Cattle

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.

Proofs for genotype by environment interactions considering pedigree and genomic data from organic and conventional cow reference populations

The aim of the present study was to prove genotype by environment interactions (G × E) for production, longevity, and health traits considering conventional and organic German Holstein dairy cattle subpopulations. The full data set included 141,778 Holstein cows from 57 conventional herds and 7,915 cows from 9 organic herds. The analyzed traits were first-lactation milk yield and fat percentage (FP), the length of productive life (LPL) and the health traits mastitis, ovarian cycle disorders, and digital dermatitis in first lactation. A subset of phenotyped cows was genotyped and used for the implementation of separate cow reference populations. After SNP quality controls, the cow reference sets considered 40,830 SNP from 19,700 conventional cows and the same 40,830 SNP from 1,282 organic cows. The proof of possible G × E was made via multiple-trait model applications, considering same traits from the conventional and organic population as different traits. In this regard, pedigree (A), genomic (G) and combined relationship (H) matrices were constructed. For the production traits, heritabilities were very similar in both organic and conventional populations (i.e., close to 0.70 for FP and close to 0.40 for milk yield). For low heritability health traits and LPL, stronger heritability fluctuations were observed, especially for digital dermatitis with 0.05 ± 0.01 (organic, A matrix) to 0.33 ± 0.04 (conventional, G matrix). Quite large genetic correlations between same traits from the 2 environments were estimated for production traits, especially for high heritability FP. For LPL, the genetic correlation was 0.67 (A matrix) and 0.66 (H matrix). The genetic correlation between LPL organic with LPL conventional was 0.94 when considering the G matrix, but only 213 genotyped cows were included. For health traits, genetic correlations were throughout lower than 0.80, indicating possible G × E. Genetic correlations from the different matrices A, G, and H for health and production traits followed the same pattern, but the estimates from G for health traits were associated with quite large standard errors. In genome-wide association studies, significantly associated SNP for production traits overlapped in the conventional and organic population. In contrast, for low heritability LPL and health traits, significantly associated SNP and annotated potential candidate genes differed in both populations. In this regard, significantly associated SNP for mastitis from conventional cows were located on Bos taurus autosomes 6 and 19, but on Bos taurus autosomes 1, 10, and 22 in the organic population. For the remaining health traits and LPL, different potential candidate genes were annotated, but the different genes reflect similar physiological pathways. We found evidence of G × E for low heritability functional traits, suggesting different breeding approaches in organic and conventional populations. Nevertheless, for a verification of results and implementation of alternative breeding strategies, it is imperative to increase the organic cow reference population.