Short communication: Uncovering quantitative trait loci associated with resistance to Mycobacterium avium ssp. paratuberculosis infection in Holstein cattle using a high-density single nucleotide polymorphism panel

Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle

Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.

The genetic architecture of complete blood counts in lactating Holstein dairy cows

Complete blood counts (CBCs) measure the abundance of individual immune cells, red blood cells, and related measures such as platelets in circulating blood. These measures can indicate the health status of an animal; thus, baseline circulating levels in a healthy animal may be related to the productive life, resilience, and production efficiency of cattle. The objective of this study is to determine the heritability of CBC traits and identify genomic regions that are associated with CBC measurements in lactating Holstein dairy cattle. The heritability of CBCs was estimated using a Bayes C0 model. The study population consisted of 388 cows with genotypes at roughly 75,000 markers and 16 different CBC phenotypes taken at one to three time points (n = 33, 131, and 224 for 1, 2, and 3 time points, respectively). Heritabilities ranged from 0.00 ± 0.00 (red cell distribution width) to 0.68 ± 0.06 (lymphocytes). A total of 96 different 1-Mb windows were identified that explained more than 1% of the genetic variance for at least one CBC trait, with 10 windows explaining more than 1% of the genetic variance for two or more traits. Multiple genes in the identified regions have functions related to immune response, cell differentiation, anemia, and disease. Positional candidate genes include RAD52 motif-containing protein 1 (RDM1), which is correlated with the degree of immune infiltration of immune cells, and C-X-C motif chemokine ligand 12 (CXCL12), which is critically involved in neutrophil bone marrow storage and release regulation and enhances neutrophil migration. Since animal health directly impacts feed intake, understanding the genetics of CBCs may be useful in identifying more disease-resilient and feed-efficient dairy cattle. Identification of genes responsible for variation in CBCs will also help identify the variability in how dairy cattle defend against illness and injury.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation

Background

Bovine paratuberculosis, or Johne’s disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions.

Results

Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection.

Conclusions

This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

Genome-wide association studies for the identification of cattle susceptible and resilient to paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) causes Johne's disease or paratuberculosis (PTB), with important animal health and economic implications. There are no therapeutic strategies to control this disease, and vaccination with inactivated vaccines is limited in many countries because it can interfere with the intradermal test used for bovine tuberculosis detection. Thus, infected animals either get culled after a positive ELISA or fecal PCR result or die due to clinical disease. In this study, we review recent studies aimed to discover genetic markers which could help to identify and select cattle less susceptible and more resilient to PTB. In recent years, the genotyping and subsequent imputation to whole-genome sequence (WGS) has allowed the identification of single-nucleotide polymorphisms (SNPs), quantitative trait loci (QTL), and candidate genes in the Bos taurus genome associated with susceptibility to MAP infection. In most of these genome-wide association studies (GWAS), phenotypes were based on ante-mortem test results including serum ELISA, milk ELISA, and detection of MAP by fecal PCR and bacteriological culture. Cattle infected with MAP display lesions with distinct severity but the associations between host genetics and PTB-associated pathology had not been explored until very recently. On the contrary, the understanding of the mechanisms and genetic loci influencing pathogen resistance, and disease tolerance in asymptomatic individuals is currently very limited. The identification of long-time asymptomatic cattle that is able to resist the infection and/or tolerate the disease without having their health and milk production compromised is important for disease control and breeding purposes.

A framework for non-preserved consensus gene module detection in Johne's disease

Johne's disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) is a major concern in dairy industry. Since, the pathogenesis of the disease is not clearly known, it is necessary to develop an approach to discover molecular mechanisms behind this disease with high confidence. Biological studies often suffer from issues with reproducibility. Lack of a method to find stable modules in co-expression networks from different datasets related to Johne's disease motivated us to present a computational pipeline to identify non-preserved consensus modules. Two RNA-Seq datasets related to MAP infection were analyzed, and consensus modules were detected and were subjected to the preservation analysis. The non-preserved consensus modules in both datasets were determined as they are modules whose connectivity and density are affected by the disease. Long non-coding RNAs (lncRNAs) and TF genes in the non-preserved consensus modules were identified to construct integrated networks of lncRNA-mRNA-TF. These networks were confirmed by protein-protein interactions (PPIs) networks. Also, the overlapped hub genes between two datasets were considered hub genes of the consensus modules. Out of 66 consensus modules, 21 modules were non-preserved consensus modules, which were common in both datasets and 619 hub genes were members of these modules. Moreover, 34 lncRNA and 152 TF genes were identified in 12 and 19 non-preserved consensus modules, respectively. The predicted PPIs in 17 non-preserved consensus modules were significant, and 283 hub genes were commonly identified in both co-expression and PPIs networks. Functional enrichment analysis revealed that eight out of 21 modules were significantly enriched for biological processes associated with Johne's disease including “inflammatory response,” “interleukin-1-mediated signaling pathway”, “type I interferon signaling pathway,” “cytokine-mediated signaling pathway,” “regulation of interferon-beta production,” and “response to interferon-gamma.” Moreover, some genes (hub mRNA, TF, and lncRNA) were introduced as potential candidates for Johne's disease pathogenesis such as TLR2, NFKB1, IRF1, ATF3, TREM1, CDH26, HMGB1, STAT1, ISG15, CASP3. This study expanded our knowledge of molecular mechanisms involved in Johne's disease, and the presented pipeline enabled us to achieve more valid results.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

Genomic loci and genetic parameters for uterine diseases in first-parity Holstein cows and associations with milk production and fertility

Based on the clinical stage (e.g., vaginal discharge) and bacterial species, several forms of uterine diseases (UD) exist and can be classified as different traits [i.e., different stages of endometritis (EM) and metritis (MET)], which may differ in their genetic background and causal physiological mechanisms. Consequently, the present study aimed to study (1) the effect of UD on 305-d lactation and fertility, (2) the estimation of heritabilities for UD traits using pedigree- and SNP-based relationships, and (3) genome-wide associations to detect significant SNP markers and to infer candidate genes for UD traits. The data set contained herd manager and veterinarian recorded UD traits of 14,810 first-lactating genotyped Holstein cows from 63 large-scale contract herds. Binary defined UD traits (healthy or diseased) according to the clinical stage were endometritis catarrhalis (EM I), endometritis mucopurulenta (EM II), endometritis purulenta (EM III), pyometra (EM IV), endometritis (EM_SOD; superordinate diagnosis = no specific clinical stage defined), and MET. The binary defined trait UDall included all EM and MET diagnoses. The prevalence of UDall was 26.7%. The effect of UD on 305-d lactation and fertility was estimated via linear and generalized linear mixed models. We applied linear single-trait animal models and threshold models to estimate pedigree- and SNP-based heritabilities for UD traits, and bivariate linear models for genetic correlation estimations between UDall with 305-d lactation and fertility traits. A diagnosis for UDall had significant unfavorable effects on the female fertility traits calving interval, interval from calving to first service, days open, and nonreturn rate after 90 d, but was unrelated to 305-d lactation records for production traits milk yield, protein yield, and fat yield. Heritabilities for UDall and EM stages were close to zero, displaying maximal values of 0.05 for pedigree and 0.07 for SNP-based relationship matrices. For MET, pedigree- and SNP-based heritabilities were <0.001 and 0.07, respectively. Genetic correlations ranged from 0.20 to 0.31 between UDall with 305-d milk, protein, and fat yield, and from 0.17 to 0.40 with fertility traits. The GWAS revealed 5 SNP on bovine chromosomes (BTA) 1, 8, 10, 23 for UDall, 5 SNP on BTA 26 for EM I, 1 SNP on BTA 19 for EM II, 4 SNP on BTA 2, 18, 20, 25 for EM III, and 4 SNP on BTA 4, 16, 20 for EM IV above the significance threshold. For EM_SOD, we identified 15 significantly associated SNP on 4 chromosomes, and 4 significant SNP on BTA 3, 20, 22, 28 for MET. Marker associations for UD traits were annotated to 24 potential candidate genes using the ENSEMBL database. Six of these genes were previously reported to be involved in uterine defense mechanisms or in endometritis. Further detected genes contribute to immune response mechanisms during bacterial infections. Different SNP significantly influenced different UD stages, explaining the inter-individual variations in clinical severity of uterine infections.

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

Johne's disease (JD), also known as paratuberculosis, is a severe production-limiting disease with significant economic and welfare implications for the global cattle industry. Caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP), JD manifests as chronic enteritis in infected cattle. In addition to the economic losses and animal welfare issues associated with JD, MAP has attracted public health concerns with potential association with Crohn's disease, a human inflammatory bowel disease. The lack of effective treatment options, such as a vaccine, has hampered JD control resulting in its increasing global prevalence. The disease was first reported in 1895, but in recognition of its growing economic impact, extensive recent research facilitated by a revolution in technological approaches has led to significantly enhanced understanding of the immunological, genetic, and pathogen factors influencing disease pathogenesis. This knowledge has been derived from a variety of diverse models to elucidate host-pathogen interactions including in vivo and in vitro experimental infection models, studies measuring immune parameters in naturally-infected animals, and by studies conducted at the population level to enable the estimation of genetic parameters, and the identification of genetic markers and quantitative trait loci (QTL) putatively associated with susceptibility or resistance to JD. The main objectives of this review are to summarize these recent developments from an immunogenetics perspective and attempt to extract the principal and common findings emerging from this wealth of recent information. Based on these analyses, and in light of emerging technologies such as gene-editing, we conclude by discussing potential future avenues for effectively mitigating JD in cattle.

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

Population differentiated copy number variation of Bos taurus, Bos indicus and their African hybrids

Background

CNV comprises a large proportion in cattle genome and is associated with various traits. However, there were few population-scale comparison studies on cattle CNV.

Results

Here, autosome-wide CNVs were called by read depth of NGS alignment result and copy number variation regions (CNVRs) defined from 102 Eurasian taurine (EAT) of 14 breeds, 28 Asian indicine (ASI) of 6 breeds, 22 African taurine (AFT) of 2 breeds, and 184 African humped cattle (AFH) of 17 breeds. The copy number of every CNVRs were compared between populations and CNVRs with population differentiated copy numbers were sorted out using the pairwise statistics V_ST and Kruskal-Wallis test. Three hundred sixty-two of CNVRs were significantly differentiated in both statistics and 313 genes were located on the population differentiated CNVRs.

Conclusion

For some of these genes, the averages of copy numbers were also different between populations and these may be candidate genes under selection. These include olfactory receptors, pathogen-resistance, parasite-resistance, heat tolerance and productivity related genes. Furthermore, breed- and individual-level comparison was performed using the presence or copy number of the autosomal CNVRs. Our findings were based on identification of CNVs from short Illumina reads of 336 individuals and 39 breeds, which to our knowledge is the largest dataset for this type of analysis and revealed important CNVs that may play a role in cattle adaption to various environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07808-7.

Genome-Wide Scanning for Signatures of Selection Revealed the Putative Genomic Regions and Candidate Genes Controlling Milk Composition and Coat Color Traits in Sahiwal Cattle

Background

In the evolutionary time scale, selection shapes the genetic variation and alters the architecture of genome in the organisms. Selection leaves detectable signatures at the genomic coordinates that provide clues about the protein-coding regions. Sahiwal is a valuable indicine cattle adapted to tropical environments with desirable milk attributes. Insights into the genomic regions under putative selection may reveal the molecular mechanisms affecting the quantitative and other important traits. To understand this, the present investigation was undertaken to explore signatures of selection in the genome of Sahiwal cattle using a medium-density genotyping INDUS chip.

Result

De-correlated composite of multiple selection signals (DCMS), which combines five different univariate statistics, was computed in the dataset to detect the signatures of selection in the Sahiwal genome. Gene annotations, Quantitative Trait Loci (QTL) enrichment, and functional analyses were carried out for the identification of significant genomic regions. A total of 117 genes were identified, which affect a number of important economic traits. The QTL enrichment analysis highlighted 14 significant [False Discovery Rate (FDR)-corrected p-value ≤ 0.05] regions on chromosomes BTA 1, 3, 6, 11, 20, and 21. The top three enriched QTLs were found on BTA 6, 20, and 23, which are associated with exterior, health, milk production, and reproduction traits. The present study on selection signatures revealed some key genes related with coat color (PDGFRA, KIT, and KDR), facial pigmentation (LEF), milk fat percent (MAP3K1, HADH, CYP2U1, and SGMS2), sperm membrane integrity (OSTC), lactation persistency (MRPS30, NNT, CCL28, HMGCS1, NIM1K, ZNF131, and CCDC152), milk yield (GHR and ZNF469), reproduction (NKX2-1 and DENND1A), and bovine tuberculosis susceptibility (RNF144B and PAPSS1). Further analysis of candidate gene prioritization identified four hub genes, viz., KIT, KDR, MAP3K1, and LEF, which play a role in coat color, facial pigmentation, and milk fat percentage in cattle. Gene enrichment analysis revealed significant Gene ontology (GO) terms related to breed-specific coat color and milk fat percent.

Conclusion

The key candidate genes and putative genomic regions associated with economic traits were identified in Sahiwal using single nucleotide polymorphism data and the DCMS method. It revealed selection for milk production, coat color, and adaptability to tropical climate. The knowledge about signatures of selection and candidate genes affecting phenotypes have provided a background information that can be further utilized to understand the underlying mechanism involved in these traits in Sahiwal cattle.

Functional Variants Surrounding Endothelin 2 Are Associated With Mycobacterium avium Subspecies paratuberculosis Infection

Bovine paratuberculosis, caused by Mycobacterium avium subspecies paratuberculosis (MAP), continues to impact the dairy industry through increased morbidity, mortality, and lost production. Although genome-wide association analyses (GWAAs) have identified loci associated with susceptibility to MAP, limited progress has been made in identifying mutations that cause disease susceptibility. A 235-kb region on Bos taurus chromosome 3 (BTA3), containing a 70-kb haplotype block surrounding endothelin 2 (EDN2), has previously been associated with the risk of MAP infection. EDN2 is highly expressed in the gut and is involved in intracellular calcium signaling and a wide array of biological processes. The objective of this study was to identify putative causal mutations for disease susceptibility in the region surrounding EDN2 in Holstein and Jersey cattle. Using sequence data from 10 Holstein and 10 Jersey cattle, common variants within the 70-kb region containing EDN2 were identified. A custom SNP genotyping array fine-mapped the region using 221 Holstein and 51 Jersey cattle and identified 17 putative causal variants (P < 0.01) located in the 5′ region of EDN2 and a SNP in the 3′ UTR (P = 0.00009) associated with MAP infection. MicroRNA interference assays, mRNA stability assays, and electrophoretic mobility shift assays were performed to determine if allelic changes at each SNP resulted in differences in EDN2 stability or expression. Two SNPs [rs109651404 (G/A) and rs110287192 (G/T)] located within the promoter region of EDN2 displayed differential binding affinity for transcription factors in binding sequences harboring the alternate SNP alleles. The luciferase reporter assay revealed that the transcriptional activity of the EDN2 promoter was increased (P < 0.05) with the A allele for rs109651404 and the G allele for rs110287192. These results suggest that the variants rs109651404 and rs110287192 are mutations that alter transcription and thus may alter susceptibility to MAP infection in Holstein and Jersey cattle.

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.

Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds

Background

Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions.

Results

Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C).

Conclusions

We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Ancestral Haplotype Mapping for GWAS and Detection of Signatures of Selection in Admixed Dairy Cattle of Kenya

Understanding the genetic structure of adaptation and productivity in challenging environments is necessary for designing breeding programs that suit such conditions. Crossbred dairy cattle in East Africa resulting from over 60 years of crossing exotic dairy breeds with indigenous cattle plus inter se matings form a highly variable admixed population. This population has been subject to natural selection in response to environmental stresses, such as harsh climate, low-quality feeds, poor management, and strong disease challenge. Here, we combine two complementary sets of analyses, genome-wide association (GWA) and signatures of selection (SoS), to identify genomic regions that contribute to variation in milk yield and/or contribute to adaptation in admixed dairy cattle of Kenya. Our GWA separates SNP effects due to ancestral origin of alleles from effects due to within-population linkage disequilibrium. The results indicate that many genomic regions contributed to the high milk production potential of modern dairy breeds with no region having an exceptional effect. For SoS, we used two haplotype-based tests to compare haplotype length variation within admixed and between admixed and East African Shorthorn Zebu cattle populations. The integrated haplotype score (iHS) analysis identified 16 candidate regions for positive selection in the admixed cattle while the between population Rsb test detected 24 divergently selected regions in the admixed cattle compared to East African Shorthorn Zebu. We compare the results from GWA and SoS in an attempt to validate the most significant SoS results. Only four candidate regions for SoS intersect with GWA regions using a low stringency test. The identified SoS candidate regions harbored genes in several enriched annotation clusters and overlapped with previously found QTLs and associations for different traits in cattle. If validated, the GWA and SoS results indicate potential for SNP-based genomic selection for genetic improvement of smallholder crossbred cattle.

Association of genetic polymorphisms related to Johne's disease with estimated breeding values of Holstein sires for milk ELISA test scores

BACKGROUND

Johne's disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls' estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle.

RESULTS

Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls' EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively.

CONCLUSIONS

The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs.

Identification of the ABCC4, IER3, and CBFA2T2 candidate genes for resistance to paratuberculosis from sequence-based GWAS in Holstein and Normande dairy cattle

Background

Bovine paratuberculosis is a contagious disease, caused by Mycobacterium avium subsp. paratuberculosis (MAP), with adverse effects on animal welfare and serious economic consequences. Published results on host genetic resistance to MAP are inconsistent, mainly because of difficulties in characterizing the infection status of cows. The objectives of this study were to identify quantitative trait loci (QTL) for resistance to MAP in Holstein and Normande cows with an accurately defined status for MAP.

Results

From MAP-infected herds, cows without clinical signs of disease were subjected to at least four repeated serum ELISA and fecal PCR tests over time to determine both infected and non-infected statuses. Clinical cases were confirmed using PCR. Only cows that had concordant results for all tests were included in further analyses. Positive and control cows were matched within herd according to their birth date to ensure a same level of exposure to MAP. Cows with accurate phenotypes, i.e. unaffected (control) or affected (clinical or non-clinical cases), were genotyped with the Illumina BovineSNP50 BeadChip. Genotypes were imputed to whole-genome sequences using the 1000 Bull Genomes reference population (run6). A genome-wide association study (GWAS) of MAP status of 1644 Holstein and 649 Normande cows, using either two (controls versus cases) or three classes of phenotype (controls, non-clinical and clinical cases), revealed three regions, on Bos taurus (BTA) chromosomes 12, 13, and 23, presenting significant effects in Holstein cows, while only one of those was identified in Normande cows (BTA23). The most significant effect was found on BTA13, in a short 8.5-kb region. Conditional analyses revealed that only one causal variant may be responsible for the effects observed on each chromosome with the ABCC4 (BTA12), CBFA2T2 (BTA13), and IER3 (BTA23) genes as good functional candidates.

Conclusions

A sequence-based GWAS on cows for which resistance to MAP was accurately defined, was able to identify candidate variants located in genes that were functionally related to resistance to MAP; these explained up to 28% of the genetic variance of the trait. These results are very encouraging for efforts towards implementation of a breeding strategy aimed at improving resistance to paratuberculosis in Holstein cows.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Characterization of the bovine salivary gland transcriptome associated with Mycobacterium avium subsp. paratuberculosis experimental challenge

BACKGROUND

Mycobacterium avium subsp. paratuberculosis (MAP), the etiologic agent of Johne's disease is spread between cattle via the fecal-oral route, yet the functional changes in the salivary gland associated with infection remain uncharacterized. In this study, we hypothesized that experimental challenge with MAP would induce stable changes in gene expression patterns in the salivary gland that may shed light on the mucosal immune response as well as the regional variation in immune capacity of this extensive gland. Holstein-Friesian cattle were euthanized 33 months' post oral challenge with MAP strain CIT003 and both the parotid and mandibular salivary glands were collected from healthy control (n = 5) and MAP exposed cattle (n = 5) for histopathological and transcriptomic analysis.

RESULTS

A total of 205, 21, 61, and 135 genes were significantly differentially expressed between control and MAP exposed cattle in dorsal mandibular (M1), ventral mandibular (M2), dorsal parotid (P1) and ventral parotid salivary glands (P2), respectively. Expression profiles varied between the structurally divergent parotid and mandibular gland sections which was also reflected in the enriched biological pathways identified. Changes in gene expression associated with MAP exposure were detected with significantly elevated expression of BoLA DR-ALPHA, BOLA-DRB3 and complement factors in MAP exposed cattle. In contrast, reduced expression of genes such as polymeric immunoglobin receptor (PIGR), TNFSF13, and the antimicrobial genes lactoferrin (LF) and lactoperoxidase (LPO) was detected in MAP exposed animals.

CONCLUSIONS

This first analysis of the transcriptomic profile of salivary glands in cattle adds an important layer to our understanding of salivary gland immune function. Transcriptomic changes associated with MAP exposure have been identified including reduced LF and LPO. These critical antimicrobial and immunoregulatory proteins are known to be secreted into saliva and their downregulation may contribute to disease susceptibility. Future work will focus on the validation of their expression levels in saliva from additional cattle of known infection status as a potential strategy to augment disease diagnosis.