Genetic variation of toll-like receptor genes and infection by Mycobacterium avium ssp. paratuberculosis in Holstein-Friesian cattle

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Association of Variants in Innate Immune Genes TLR4 and TLR5 with Reproductive and Milk Production Traits in Czech Simmental Cattle

Bovine genes TLR4 and TLR5, which encode antibacterial toll-like receptors, were screened for polymorphisms in Czech Red Pied (Czech Simmental) cattle to identify variants associated with reproduction, udder health, and milk production traits. Variants were discovered by hybrid resequencing of 164 bulls using HiSeq X-Ten and PacBio technologies and then individually genotyped. Nominal p-values < 0.05 for associations were detected in 18 combinations between 14 polymorphisms and 15 traits using one-way analysis of variance (ANOVA). The TLR4 variants g.610C>T (rs43578094) and g.10310T>G (rs8193072) in reference AC000135.1 were strictly associated with the index of early reproductive disorders and maternal calving ease, respectively, at false discovery rate (FDR) < 0.05. A highly permissive false discovery rate cutoff of 0.6 separated seventeen combinations in both genes comprising eight positives. In the case of the TLR4 variant g.9422T>C (rs8193060), indications were obtained for the association with as many as four reproductive traits: incidence of cystic ovaries, early reproductive disorders, calving ease, and production longevity. The permissive FDR interpretation for the TLR5 data indicated associations with cyst incidence and early reproduction disorders with maternal calving ease. Moreover, three TLR5 polymorphisms correlated with milk production traits. The discrepancy of the observed associations with the predicted impacts of the SNPs on protein function points to the role of haplotypes. Nevertheless, this question should be resolved on a larger scale. The observed associations are endorsed by independent evidence from the published functional roles in other species and by the published QTL mapping data.

Genetic polymorphism in European and African sheep breeds reared in Hungary based on 48 SNPs associated with resistance to gastrointestinal parasite infection using KASP-PCR technique

This pilot study used an alternative and economically efficient technique, the Kompetitive Allele-Specific Polymerase Chain Reaction (KASP-PCR) to examine 48 SNPs from 11 parasite-resistance genes found on 8 chromosomes in 110 animals from five sheep breeds reared in Hungary; Hungarian Tsigai, White Dorper, Dorper, Ile de France, and Hungarian Merino. Allele and genotype frequencies, fixation index, observed heterozygosity, expected heterozygosity, F statistic, and their relationship with the Hardy-Weinberg equilibrium (WHE) and the polymorphic information content (PIC) were determined, followed by principal component analysis (PCA). As much as 32 SNPs out of the 48 initially studied were successfully genotyped. A total of 9 SNPs, 4 SNPs in TLR5, 1 SNP in TLR8, and 4 SNPs in TLR2 genes, were polymorphic. The variable genotype and allele frequency of the TLRs gene indicated genetic variability among the studied sheep breeds, with the Hungarian Merino exhibiting the most polymorphisms, while Dorper was the population with the most SNPs departing from the HWE. According to the PIC value, the rs430457884-TLR2, rs55631273-TLR2, and rs416833129-TLR5 were found to be informative in detecting polymorphisms among individuals within the populations, whereas the rs429546187-TLR5 and rs424975389-TLR5 were found to have a significant influence in clustering the population studied. This study reported a moderate level of genetic variability and that a low to moderate within-breed diversity was maintained in the studied populations.

Role of Toll-Like Receptor 4 in Mycobacterium avium subsp. paratuberculosis Infection of Bovine Mammary Epithelial (MAC-T) Cells In Vitro

Toll-like receptor 4 (TLR4) encodes an innate immune cell pattern-recognition receptor implicated in the recognition of Mycobacterium avium subsp. paratuberculosis (MAP), the causative agent of Johne's disease in ruminants. Polymorphisms in TLR4 have been associated with susceptibility to MAP infection. In this study, a previously developed TLR4 knockout (TLR4KO) bovine mammary epithelial (MAC-T) cell line and wild-type MAC-T cells (WT) were infected with live MAP for 72 h to identify potential immunoregulatory miRNAs, inflammatory genes, and cytokines/chemokines impacted by MAP infection in the presence/absence of TLR4. Cytokines/chemokines production in culture supernatants was measured by multiplexing immunoassay. Total RNA was extracted from the remaining MAC-T cells, and quantitative PCR was performed to determine the expression of inflammatory genes and selected bovine miRNAs. Results showed that the levels of tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), CXCL8, CXCL10, CCL4, and CCL3 were significantly induced in WT MAC-T cells during MAP infection. However, TLR4KO MAC-T cells had greater secretion of CCL3, IL-6, vascular endothelial growth factor (VEGF-α), and TNF-α and decreased secretion of CXCL10 and CCL2. Moreover, the expression of inflammatory genes was induced in TLR4KO cells. The expression of miRNAs (miR133b, miR-92a, and miR-184) was increased in WT MAC-T cells post-MAP infection; however, there was no significant induction of these miRNAs in TLR4KO cells, which suggests they are involved in regulating the innate immune response to MAP infection. Target gene function analysis further suggests that miR-92a may be involved in TLR and interleukin signaling and miR-133b and miR-184 may be involved in other signaling pathways. These findings support the involvement of TLR4 in the regulation of innate immune response to MAP. IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent for paratuberculosis or Johne's disease (JD) in ruminants, a disease clinically very similar to Crohn's disease in humans. Polymorphisms in the bovine Toll-like receptor genes (TLR1, TLR2, and TLR4) have been shown to affect MAP recognition and host innate immune response and have been associated with increased susceptibility of cattle to paratuberculosis. Our results demonstrated that knocking out the TLR4 gene in bovine MAC-T cells enhanced inflammation in response to MAP. These findings show divergent roles for TLR4 in Escherichia coli lipopolysaccharide and mycobacterial infections, and this may have important consequences for the treatment of these inflammatory diseases and for genetic selection to improve disease resistance. It advances our understanding of the role of TLR4 in the context of MAP infection.

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.

CRISPR-Cas9-mediated knockout of TLR4 modulates Mycobacterium avium ssp. paratuberculosis cell lysate-induced inflammation in bovine mammary epithelial cells

Toll-like receptor 4 (TLR4) is a pattern-recognition receptor involved in the recognition of microbial pathogens and host alarmins. Ligation to TLR4 initiates a signaling cascade that leads to inflammation. Polymorphisms in bovine TLR4 have been associated with Mycobacterium avium ssp. paratuberculosis (MAP) susceptibility and resistance, the cause of Johne's disease, and milk somatic cell score, a biomarker of mastitis. Although the contribution of TLR4 to recognition of bacterial lipopolysaccharide (LPS) has been well characterized, its role in MAP recognition is less certain. Clustered regularly interspaced short palindromic repeats-Cas9 mediated gene editing was performed to generate TLR4 knockout (KO) mammary epithelial cells to determine if TLR4 expression is involved in the initiation of the host inflammatory response to MAP cell lysate (5 and 10 µg/mL) and Escherichia coli LPS (5 µg/mL). The absence of TLR4 in KO cells resulted in enhanced expression of key inflammatory genes (TNFA and IL6), anti-inflammatory genes (IL10 and SOCS3), and supernatant cytokine and chemokine levels (TNF-α, IL-6, IL-10, CCL3) in response to the MAP cell lysate (10 µg/mL). However, in response to LPS, the KO cells showed reduced expression of key inflammatory genes (TNFA, IL1A, IL1B, and IL6) and supernatant cytokine levels (TNF-α, IL-6, CCL2, IL-8) as compared with unedited cells. Overall, these results confirm that TLR4 is essential for eliciting inflammation in response to LPS; however, exacerbated gene and protein expression in TLR4 KO cells in response to MAP cell lysate suggests a different mechanism of infection and host response for MAP, at least in terms of how it interacts with TLR4. These novel findings show potential divergent roles for TLR4 in mycobacterial infections, and this may have important consequences for the therapeutic control of inflammation in cattle.

Survey of Candidate Single-Nucleotide Polymorphisms in SLC11A1, TLR4, NOD2, PGLYRP1, and IFNγ in Ankole Longhorn Cattle in Central Region of Uganda to Determine Their Role in Mycobacterium avium Subspecies paratuberculosis Infection Outcome

Mycobacterium avium ssp. paratuberculosis (MAP) is the cause of Johne's disease (JD) in a wide range of domestic and wild ruminants. Single-nucleotide polymorphisms (SNPs) in several genes including solute-like carrier 11A1 (SLC11A1), interferon gamma (IFNγ), Toll-like receptor 4 (TLR4), nucleotide-binding oligomerization domain 2 gene (NOD2), and bovine peptidoglycan recognition protein 1 (PGLYRP1) have been implicated in influencing the infection outcome of MAP in cattle. We have carried out a survey in a population of Ankole cattle from three districts in the central region of Uganda including Isingiro, Lyantonde, and Rakai to determine the role played by several SNPs on the above genes in the infection outcome of local cattle in Uganda. Nine hundred fifty-five heads of cattle obtained from 93 herds were tested using ELISA. Thirty-five ELISA-positive cattle and 35 negative herd mates from a total of 955 cattle tested for MAP were genotyped using iPLEX MassARRAY genotyping systems to detect the presence of a total of 13 SNPS in five different genes (SLC11A1, IFNγ, TLR4, NOD2, and PGLYRP1). The cow-level prevalence of MAP infection in Ankole Longhorn cattle in the three districts was 3.98% (35/955), while the herd-level prevalence was 27.9% and within-herd prevalence was 12 ± 1.5% (95% CI = 9.1-14.8%). The genotypes and allele frequencies of the MAP-positive cattle were compared with those of their ELISA-negative herd mates to determine the significance of the polymorphisms. The results showed that SNPs rs109915208, rs110514940, and rs110905610 on SLC11A1, c.480G>A and c.625C>A on PGLYRP1, and c.2021C>T on TLR4 were monomorphic in both seropositive and seronegative cattle and therefore had no influence on the infection outcome. The remaining SNPs studied in the five genes [SLC11A1: rs109614179; TLR4: rs29017188 (c.226G>C), c.2021C>T; NOD2: rs110536091, rs111009394; PGLYRP1: c.102G>C, c.480G>A, c.625C>A; IFNγ: rs110853455] were polymorphic, but their allele and genotype frequencies did not show any significant difference between the seropositive and seronegative cattle. No significant difference was observed for any haplotype at the gene level.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case–control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Correlations between single nucleotide polymorphisms in bovine CD209, SLC11A1, SP110 and TLR2 genes and estimated breeding values for several traits in Spanish Holstein cattle

Little is known about the correlations between the genetic susceptibility/resistance to Mycobacterium avium subsp. paratuberculosis (Map) infection and the estimated breeding values for type, production and functional traits. Previously, we identified 70 combinations of five single nucleotide polymorphisms (SNPs) in four bovine innate immune genes (SLC11A1, SP110, TLR2, CD209) that are associated with the genetic risk of paratuberculosis (PTB) or Johne's disease progression, which can be graded as low (LOWIN), latent (LATIN), or patent (PATIN) risk. Other possible combinations of these 5 SNPs were grouped in the average group (AVERIN). In the current study, differences in estimated breeding values (EBVs) for several traits were analyzed using linear models in a large cohort of Holstein cows (N = 15656) genotyped across Spain in 2016 or 2017. After the assignment of each genotyped cow to a risk group, cows within the PATIN risk group (N = 1448) had a superior combined genetic index (2797.57), type genetic index (524.62), milk yield (653.92 kg), protein yield (21.77 kg), fat yield (24.82 kg) and economic merit index (125 Euros) compared with the other three risk groups. Statistically significant differences in the longevity scores between the cows that were included in the PATIN risk group (108.85) and the LOWIN (107.82) and AVERIN (107.92) groups were also observed. The associations between the genetic risk groups and PTB diagnostic results were validated in a population of 99 cows from a Spanish farm with a high prevalence of PTB. Significant differences in ELISA readings between the PATIN (65.49 %) and the AVERIN (15.97 %), LATIN (2.11 %), and LOWIN (3.27 %) groups were observed. In addition, significant differences in Map DNA copies/gram of feces were observed between the PATIN and the other three risk groups. These results together with the substantial economic impact of PTB in dairy cattle support the selection of the animals with less susceptibility to PTB in the Spanish breeding program.

Association of TLR gene variants in a Czech Red Pied cattle population with reproductive traits

The bovine genes TLR1, TLR2 and TLR6, which encode Toll-like receptors, key components of the innate immune system, were screened for polymorphisms in Czech Red Pied (Czech Simmental) cattle, and the different variants present in the population were tested for association with reproductive and fitness traits. Diversity was investigated in a group of 164 bulls using hybrid resequencing of pooled amplicons with PacBio technology and of pooled genomic DNA using HiSeq X-Ten technology. The validated single nucleotide polymorphisms (SNPs) were genotyped in individual animals using the primer extension technique. The association of genotypic classes of 16 polymorphisms with six phenotypic traits were estimated with one-way analysis of variance (ANOVA) and with restricted maximum likelihood (REML) algorithm. The evaluated traits included the incidence of cystic ovaries, index of early reproductive disorders, paternal and maternal indicators of calving ease, production longevity and calf vitality index. The estimated breeding values were used for combined trait quantification. Early traits, namely, cystic ovaries and early reproductive disorders, were not associated with any of the tested polymorphisms according to the general ANOVA test. By contrast, five variants of all three genes were associated with calving ease, both paternal and maternal. The production longevity correlated with two variants of TLR1 and the calf vitality index correlated with the 1044 T > C (rs68268249) polymorphism in TLR2. The false discovery rate (FDR) according to Benjamini-Hochberg was favourable for the calving ease trait (0.221) and maternal calving ease (0.214), which allows to consider the observed associations real, regardless of the error arising from the multiple comparisons. These results were supported by REML only partially, probably in view of the additivity assumption. Two mechanisms of action on calving are conceivable, either via infection resistance or via the involvement of TLR2 in signalling in the myometrium. The known formation of heterodimers by the TLR1, -2 and -6 products might be responsible for the shared pattern of action in these genes. The association of the calf vitality index with TLR2 variation might reflect the increased role of infections in calves compared to adult animals.

Potential of TLR-gene diversity in Czech indigenous cattle for resistance breeding as revealed by hybrid sequencing

A production herd of Czech Simmental cattle (Czech Red Pied, CRP), the conserved subpopulation of this breed, and the ancient local breed Czech Red cattle (CR) were screened for diversity in the antibacterial toll-like receptors (TLRs), which are members of the innate immune system. Polymerase chain reaction (PCR) amplicons of TLR1, TLR2, TLR4, TLR5, and TLR6 from pooled DNA samples were sequenced with PacBio technology, with 3- 5 ×  coverage per gene per animal. To increase the reliability of variant detection, the gDNA pools were sequenced in parallel with the Illumina X-ten platform at low coverage ( 60 × per gene). The diversity in conserved CRP and CR was similar to the diversity in conserved and modern CRP, representing 76.4 % and 70.9 % of its variants, respectively. Sixty-eight (54.4 %) polymorphisms in the five TLR genes were shared by the two breeds, whereas 38 (30.4 %) were specific to the production herd of CRP; 4 (3.2 %) were specific to the broad CRP population; 7 (5.6 %) were present in both conserved populations; 5 (4.0 %) were present solely for the conserved CRP; and 3 (2.4 %) were restricted to CR. Consequently, gene pool erosion related to intensive breeding did not occur in Czech Simmental cattle. Similarly, no considerable consequences were found from known bottlenecks in the history of Czech Red cattle. On the other hand, the distinctness of the conserved populations and their potential for resistance breeding were only moderate. This relationship might be transferable to other non-abundant historical cattle breeds that are conserved as genetic resources. The estimates of polymorphism impact using Variant Effect Predictor and SIFT software tools allowed for the identification of candidate single-nucleotide polymorphisms (SNPs) for association studies related to infection resistance and targeted breeding. Knowledge of TLR-gene diversity present in Czech Simmental populations may aid in the potential transfer of variant characteristics from other breeds.

Genetic association of polymorphisms in bovine TLR2 and TLR4 genes with Mycobacterium avium subspecies paratuberculosis infection in Indian cattle population

Toll like receptors (TLRs) are pattern recognition molecules involved in cellular recognition of Mycobacterium avium subspecies paratuberculosis (MAP), the infectious agent causing Paratuberculosis (PTB), a notified disease of domestic and wild ruminants. The present study was undertaken to investigate the presence of single nucleotide polymorphisms (SNPs) in TLR2 and TLR4 gene and to evaluate association of these SNPs with occurrence of PTB in Indian cattle. A total of 213 cattle, were subjected to multiple diagnostic tests viz. Johnin PPD, ELISA test (Indigenous and Parachek kit method), fecal microscopy and fecal culture for detection of MAP infection. Based on screening results 51 animals each were assigned to case and control population. Two SNPs viz. rs55617172, rs41830058 in TLR2 gene and two SNPs viz. rs8193046, rs8193060 in TLR4 gene and were genotyped by PCR-RFLP method. All SNPs were found to be polymorphic except rs41830058 in the case-control population. Both SNPs in TLR4 gene but none in TLR2 genes were significantly associated with the occurrence of PTB in our population. The genotypes in SNP rs8193046 and SNP rs8193060 were significantly (P < 0.01) different in case-control population. These findings suggest that SNPs rs8193046 and rs8193060 are likely a potential marker against MAP infection and a selection programme eliminating AG genotype for rs8193046 and CT genotype for rs8193060 might be beneficial in conferring resistance to MAP infection in Indian cattle population.

Genome-wide association study of Mycobacterium avium subspecies Paratuberculosis infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50 K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50 K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P < 5 × 10^− 5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1 Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5 Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5385-3) contains supplementary material, which is available to authorized users.

Toll-like receptor 4 gene polymorphisms influence milk production traits in Chinese Holstein cows

The research reported in this Research Communication aimed to describe the influence of Toll-like receptor 4 gene polymorphisms on milk production traits in Chinese Holstein cows. Toll-like receptor 4 (TLR4) is an important member of the toll-like receptor gene family that is widely found in various organisms. Since TLR4 can identify molecular patterns from various pathogenic microorganisms and induce natural and acquired immunity, it plays an important role in disease resistance in dairy cows. Two single nucleotide polymorphisms (SNPs) of TLR4 (c.-226 G > C and c.2021 C > T) that were previously found to be associated with health traits were genotyped using Sequenom MassARRAY (Sequenom Inc., San Diego, CA) for Chinese Holstein cows (n = 866). The associations between SNPs or their haplotypes and milk production traits and somatic cell count were analyzed by the generalized linear model procedure of Statistics Analysis System software (SAS). The c.-226 G > C and c.2021 C > T showed low linkage disequilibrium (r2 = 0·192). There was no association between these two SNPs and SCC, but significant effects were found for SNP c.-226 G > C on test-day milk yield, fat content, protein content, and total solid and milk urea nitrogen (P T and the SNP haplotypes on test-day milk yield, fat content, protein content, lactose content and total solids (P C was located within several potential transcription factor binding sites, including transcription factor AP-2. The polymorphisms c.-226 G > C and c.2021 C > T had significant effects on the milk production for Chinese Holstein, and these SNP could be used for molecular marker-assisted selection of milk production.

Association between combinations of genetic polymorphisms and epidemiopathogenic forms of bovine paratuberculosis

Control of major mycobacterial diseases affecting livestock is a challenging issue that requires different approaches. The use of genetic markers for improving resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle has been explored as a promising population strategy We performed paratuberculosis epidemiopathogenic phenotypic and genotypic characterization involving 24 SNPs in six candidate genes (NOD2, CD209, SLC11A1, SP110, TLR2 and TLR4) on 502 slaughtered Friesian cows. In the current study, we investigate whether recently proposed paratuberculosis (PTB) epidemiopathogenic (EP) forms (apparently free-AF, latent-LAT and patent-PAT) could be associated with some combination of these 24 SNPs. Best EP form grouping was obtained using a combination of 5 SNPs in four genes (CD209: rs210748127; SLC11A1: rs110090506; SP110: rs136859213 and rs110480812; and TLR2: rs41830058). These groups were defined according to the level of infection progression risk to patent epidemiopathogenic forms and showed the following distributions: LOWIN (low) with 39 (8%) cases (94.9% AF/5.1% LAT/0% PAT); LATIN (low) with 17 (3%) cases (5.9% AF/94.1% LAT/0% PAT); AVERIN (average) with 413 (82%) cases (52.1% AF/38.5% LAT/9.4% PAT) and PATIN (patent) with 33 (7%) cases (36.4% AF/24.2% LAT/39.4% PAT). Age of slaughter was significantly higher for LATIN (88.3 months) compared to AVERIN (65.3 months; p = 0.0007) and PATIN (59.1 months; p = 0.0004), and for LOWIN (73.9 months) compared to PATIN (p = 0.0233), and nearly significant compared to AVERIN (p = 0.0572) These results suggest that some selected genetic polymorphisms have a potential use as markers of PTB EP forms and thus add a new tool for the control of this widespread infection.

Diversity of the TLR4 Immunity Receptor in Czech Native Cattle Breeds Revealed Using the Pacific Biosciences Sequencing Platform

The allelic variants of immunity genes in historical breeds likely reflect local infection pressure and therefore represent a reservoir for breeding. Screening to determine the diversity of the Toll-like receptor gene TLR4 was conducted in two conserved cattle breeds: Czech Red and Czech Red Pied. High-throughput sequencing of pooled PCR amplicons using the PacBio platform revealed polymorphisms, which were subsequently confirmed via genotyping techniques. Eight SNPs found in coding and adjacent regions were grouped into 18 haplotypes, representing a significant portion of the known diversity in the global breed panel and presumably exceeding diversity in production populations. Notably, the ancient Czech Red breed appeared to possess greater haplotype diversity than the Czech Red Pied breed, a Simmental variant, although the haplotype frequencies might have been distorted by significant crossbreeding and bottlenecks in the history of Czech Red cattle. The differences in haplotype frequencies validated the phenotypic distinctness of the local breeds. Due to the availability of Czech Red Pied production herds, the effect of intensive breeding on TLR diversity can be evaluated in this model. The advantages of the Pacific Biosciences technology for the resequencing of long PCR fragments with subsequent direct phasing were independently validated.

Molecular markers for resistance against infectious diseases of economic importance

Huge livestock population of India is under threat by a large number of endemic infectious (bacterial, viral, and parasitic) diseases. These diseases are associated with high rates of morbidity and mortality, particularly in exotic and crossbred cattle. Beside morbidity and mortality, economic losses by these diseases occur through reduced fertility, production losses, etc. Some of the major infectious diseases which have great economic impact on Indian dairy industries are tuberculosis (TB), Johne's disease (JD), mastitis, tick and tick-borne diseases (TTBDs), foot and mouth disease, etc. The development of effective strategies for the assessment and control of infectious diseases requires a better understanding of pathogen biology, host immune response, and diseases pathogenesis as well as the identification of the associated biomarkers. Indigenous cattle (Bos indicus) are reported to be comparatively less affected than exotic and crossbred cattle. However, genetic basis of resistance in indigenous cattle is not well documented. The association studies of few of the genes associated with various diseases, namely, solute carrier family 11 member 1, Toll-like receptors 1, with TB; Caspase associated recruitment domain 15, SP110 with JD; CACNA2D1, CD14 with mastitis and interferon gamma, BoLA--DRB3.2 alleles with TTBDs, etc., are presented. Breeding for genetic resistance is one of the promising ways to control the infectious diseases. High host resistance is the most important method for controlling such diseases, but till today no breed is total immune. Therefore, work may be undertaken under the hypothesis that the different susceptibility to these diseases are exhibited by indigenous and crossbred cattle is due to breed-specific differences in the dealing of infected cells with other immune cells, which ultimately influence the immune response responded against infections. Achieving maximum resistance to these diseases is the ultimate goal, is technically possible to achieve, and is permanent. Progress could be enhanced through introgression of resistance genes to breeds with low resistance. The quest for knowledge of the genetic basis for infectious diseases in indigenous livestock is strongly warranted.

Association of <i>TLR4</i> and <i>CARD15/NOD2</i> polymorphisms with SCC in Holstein–Friesian cattle

Abstract. Mastitis is one of the most important dairy cattle diseases which results in economic losses in dairy production. Mastitis cases can be classified as subclinical or clinical. All forms of mastitis lead to changes in milk composition and induce an increase in somatic cell count (SCC). SCC is a very important and basic indicator of udder health. An increase in SCC is usually caused by the immune response to the invasion of pathogens contributing to mastitis. The aim of this study was to investigate associations between the polymorphisms of selected genes (TLR4 and CARD15/NOD2) whose products are involved in the identification of pathogen-associated molecular patterns (PAMPs) during the innate immune response to infection, and immunity to mastitis expressed as SCC. The genes under study were also examined for epistatic effects as well as effects of interactions with parity and stages of lactation. In all the studied classes, allele G of TLR4 had a favourable additive effect with negative values, contributing to a lower lnSCC. Allele A of CARD15/NOD2 had a desirable additive effect which varied with time and the changing internal environment during lactation. With regard to the dominance effect, allele A of CARD15/NOD2 was found to be significantly associated with a higher SCC in milk in the first lactation and in the third stage of each single lactation. Moreover, statistically significant epistatic effects were found, in particular additive–additive and dominance–additive interactions were favourably associated with SCC which was lower than expected in the case of no epistasis.

Genetic association analysis of paratuberculosis forms in holstein-friesian cattle

A genetic susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infections in ruminants has been longtime suspected to exist. Recently, natural infections in cattle have been reclassified into latent and patent forms based on histopathological findings and their associations with immunological and microbiological variables. This study aims to explore whether these newly defined phenotypes are associated with twenty-four single-nucleotide polymorphisms (SNPs) in six bovine candidate genes: nucleotide-binding oligomerization domain 2 (NOD2), solute carrier family 11 member A1 (SLC11A1), nuclear body protein SP110 (SP110), toll-like receptors (TLRs) 2 and 4, and CD209 (also known as DC-SIGN, dendritic cell-specific ICAM3-grabbing nonintegrin). SNPs were genotyped for 772 Holstein-Friesian animals (52.6% apparently free; 38.1% latent; 9.3% patent) by TaqMan OpenArray technology. Genotypic-phenotypic associations were assessed by logistic regression analysis adjusted for age at slaughter, under five models (codominant, dominant, recessive, overdominant, and log-additive), and corrected for multiple testing. The rs208222804 C allele (CD209 gene) was found to be associated with latent paratuberculosis (log-additive model: P < 0.0034 after permutation procedure; OR = 0.64, 95% CI = 0.48–0.86). No significant association was detected between any SNP and the patent phenotype. Consequently, CD209 gene may play a key role in the pathogenesis of bovine paratuberculosis.

Impact of host genetics on susceptibility and resistance to Mycobacterium avium subspecies Paratuberculosis infection in domestic ruminants

Johne's disease or Paratuberculosis has emerged as major infectious disease of animals in general and domestic livestock in particular on global basis. There have been major initiatives in developed countries for the control of this incurable malady of animals and human beings alike (inflammatory bowel disease or Crohn's disease). Disease has not received similar attention due to inherent complexities of disease, diagnosis and control, in resource poor counties around the world. However, the rich genetic diverstiy of the otherwise low productive animal population offers opportunity for the control of Johne's disease and improve per animal productivity. Present review aims to gather and compile information available on genetics or resistance to Johne's disease and its future exploitation by resource poor countries rich in animal diversity. This review will also help to create awareness and share knowledge and experience on prevalence and opportunities for control of Johne's disease in the livestock population to boost per animal productivity among developing and poor countries of the world. Breeding of animals for disease resistance provides good, safe, effective and cheaper way of controlling Johne's disease in animals, with especial reference to domestic livestock of developing and poor countries. Study will help to establish better understanding of the correlation between host cell factors and resistance to MAP infection which may have ultimately help in the control of Johne's disease in future.

Polymorphism of cytokine and innate immunity genes associated with bovine brucellosis in cattle

Genetic susceptibility to brucellosis is multifactorial, and it is known that impairment of the immune system could contribute to risk for getting brucellosis. The aim of the study was to find association of bovine brucellosis with 20 SNPs pertaining to bovine cytokine (IFNG, IFNGR1, IFNGR2, TNFA) and innate immunity (SLC11A1, TLR1, TLR4, and TLR9) genes using PCR-RFLP genotyping technique and it was observed that SLC11A1 (+1066 C/G), TLR1 (+1446 C/A), TLR1 (+1380 G/A), TLR4 (+10 C/T) and TLR4 (+399 C/T) loci were significantly (P≤0.05) associated with bovine brucellosis. The odds ratios (OR) of CG and CC genotypes versus GG genotype were 0.31 (0.12-0.82; 95% CI) and 0.18 (0.03-1.06; 95% CI) at SLC11A1 (+1066 C/G) locus in cases of brucellosis affected cattle. For TLR1 (+1380 G/A) locus, the OR for AG and AA genotypes versus GG genotypes were 0.15 (0.05-0.44; 95% CI) and 0.26 (0.04-1.47; 95% CI) which indicated that proportion of GG homozygote was significantly higher in brucellosis affected animals as compared to control. At TLR1 (+1446 C/A) locus the OR of AC genotype versus CC genotype was 0.24 (0.08-0.68; 95% CI) which revealed that relative proportion CC genotypes was significantly higher in case population. The TLR4 (+10 C/T) locus had three genotypes (TT, CT and CC) where OR of CT and CC genotypes versus TT genotype were near to zero. The OR of CT genotypes versus CC genotypes was 8.25 (0.94-71.92; 95% CI) at TLR4 (+399 C/T) locus and indicated that CT genotype had higher odds of bovine brucellosis than control animals.

Functional polymorphisms in Toll-like receptor genes for innate immunity in farm animals

The exploitation of the genetic factors affecting the health status of farm animals represents an alternative approach to controlling the diseases caused by microbial pathogens. The determination of innate immunity based on the genotype of the germplasm cells is a constraint for specificity but becomes an advantage in breeding schemes. The structural deviations among Toll-like receptors (TLRs), as the most frequently studied innate immunity components, have been documented at all levels, i.e., interspecific, inter- and intravarietal, in the main farm species. The current computational methods facilitate the prediction of the functional consequences of the observed mutations. Subsequently, these predictions can be verified through immunological responsiveness and population-wide association studies. The frequency and haplotype grouping of individual polymorphisms are used to track the origin and selection coefficient as independent indicators of functional changes. The Toll-like receptor variants associated with mastitis and mycobacterial infection have been identified in cattle, consequently, the targeting of these proteins in breeding could contribute to disease control. The range of infections affected by TLR polymorphisms suggests that the improvement of innate resistance is feasible in more species. Thus, the traditional breeds and wild populations should be regarded as the resources of genetic variability accessible for these purposes.

Toll like receptor 9 (TLR9) polymorphism G520R in sheep is associated with seropositivity for Small Ruminant Lentivirus

Infectious diseases of sheep are of major economic importance causing direct and indirect losses. Among the major sheep infectious agents are Small Ruminant Lentivirus, Chlamydophila abortus and Mycobacterium avium subsp. paratuberculosis infections, mainly due to their worldwide distribution and economic impact that they cause. Based on the differential susceptibility to infectious diseases between and within breeds and on the recent findings regarding the putative involvement of TLR9 in disease susceptibility, the aim of this study was to evaluate the levels of nucleotide variation of TLR9 and its mediator MyD88 in three sheep flocks originated from different breeds and assess their possible association with seropositivity/seronegativity for different infectious agents. The analysis indicated that the change of G to R at codon 520 of TLR9 polypeptide shows a significant association with Small Ruminant Lentivirus seropositivity. This amino-acid substitution, which can result in polarity change, might influence structure and function of LRR17, interfering with ligand binding and thus could be used in studies investigating susceptibility/resistance to Small Ruminant Lentivirus infections in sheep.

Altered Toll-like receptor 9 signaling in Mycobacterium avium subsp. paratuberculosis-infected bovine monocytes reveals potential therapeutic targets

Mycobacterium avium subsp. paratuberculosis is the causative agent of Johne's disease in cattle. The complex, multifaceted interaction of M. avium subsp. paratuberculosis with its host includes dampening the ability of infected cells to respond to stimuli that promote M. avium subsp. paratuberculosis clearance. By disrupting host defenses, M. avium subsp. paratuberculosis creates an intracellular environment that favors the establishment and maintenance of infection. Toll-like receptors (TLRs) are important sensors that initiate innate immune responses to microbial challenge and are also immunotherapeutic targets. For example, TLR9 contributes to host defense against M. avium subsp. paratuberculosis, and its agonists (CpG oligodeoxynucleotides [ODNs]) are under investigation for treatment of Johne's disease and other infections. Here we demonstrate that M. avium subsp. paratuberculosis infection changes the responsiveness of bovine monocytes to TLR9 stimulation. M. avium subsp. paratuberculosis inhibits classical TLR9-mediated responses despite a 10-fold increase in TLR9 expression and maintained uptake of CpG ODNs. Other TLR9-mediated responses, such as oxidative burst, which occur through noncanonical signaling, remain functional. Kinome analysis verifies that classic TLR9 signaling is blocked by M. avium subsp. paratuberculosis infection and that signaling instead proceeds through a Pyk2-mediated mechanism. Pyk2-mediated signaling does not hinder infection, as CpG ODNs fail to promote M. avium subsp. paratuberculosis clearance. Indeed, Pyk2 signaling appears to be an important aspect of M. avium subsp. paratuberculosis infection, as Pyk2 inhibitors significantly reduce the number of intracellular M. avium subsp. paratuberculosis bacteria. The actions of M. avium subsp. paratuberculosis on TLR9 signaling may represent a strategy to generate a host environment which is better suited for infection, revealing potential new targets for therapeutic intervention.

Systems biology analysis of gene expression during in vivo Mycobacterium avium paratuberculosis enteric colonization reveals role for immune tolerance

Survival and persistence of Mycobacterium avium subsp. paratuberculosis (MAP) in the intestinal mucosa is associated with host immune tolerance. However, the initial events during MAP interaction with its host that lead to pathogen survival, granulomatous inflammation, and clinical disease progression are poorly defined. We hypothesize that immune tolerance is initiated upon initial contact of MAP with the intestinal Peyer's patch. To test our hypothesis, ligated ileal loops in neonatal calves were infected with MAP. Intestinal tissue RNAs were collected (0.5, 1, 2, 4, 8 and 12 hrs post-infection), processed, and hybridized to bovine gene expression microarrays. By comparing the gene transcription responses of calves infected with the MAP, informative complex patterns of expression were clearly visible. To interpret these complex data, changes in the gene expression were further analyzed by dynamic Bayesian analysis, and genes were grouped into the specific pathways and gene ontology categories to create a holistic model. This model revealed three different phases of responses: i) early (30 min and 1 hr post-infection), ii) intermediate (2, 4 and 8 hrs post-infection), and iii) late (12 hrs post-infection). We describe here the data that include expression profiles for perturbed pathways, as well as, mechanistic genes (genes predicted to have regulatory influence) that are associated with immune tolerance. In the Early Phase of MAP infection, multiple pathways were initiated in response to MAP invasion via receptor mediated endocytosis and changes in intestinal permeability. During the Intermediate Phase, perturbed pathways involved the inflammatory responses, cytokine-cytokine receptor interaction, and cell-cell signaling. During the Late Phase of infection, gene responses associated with immune tolerance were initiated at the level of T-cell signaling. Our study provides evidence that MAP infection resulted in differentially regulated genes, perturbed pathways and specifically modified mechanistic genes contributing to the colonization of Peyer's patch.

Evolution of the bovine TLR gene family and member associations with Mycobacterium avium subspecies paratuberculosis infection

Members of the Toll-like receptor (TLR) gene family occupy key roles in the mammalian innate immune system by functioning as sentries for the detection of invading pathogens, thereafter provoking host innate immune responses. We utilized a custom next-generation sequencing approach and allele-specific genotyping assays to detect and validate 280 biallelic variants across all 10 bovine TLR genes, including 71 nonsynonymous single nucleotide polymorphisms (SNPs) and one putative nonsense SNP. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and specialized beef and dairy breeds could not be differentiated despite an average polymorphism density of 1 marker/158 bp. Collectively, 160 tagSNPs and two tag insertion-deletion mutations (indels) were sufficient to predict 100% of the variation at 280 variable sites for both Bos subspecies and their hybrids, whereas 118 tagSNPs and 1 tagIndel predictively captured 100% of the variation at 235 variable sites for B. t. taurus. Polyphen and SIFT analyses of amino acid (AA) replacements encoded by bovine TLR SNPs indicated that up to 32% of the AA substitutions were expected to impact protein function. Classical and newly developed tests of diversity provide strong support for balancing selection operating on TLR3 and TLR8, and purifying selection acting on TLR10. An investigation of the persistence and continuity of linkage disequilibrium (r²≥0.50) between adjacent variable sites also supported the presence of selection acting on TLR3 and TLR8. A case-control study employing validated variants from bovine TLR genes recognizing bacterial ligands revealed six SNPs potentially eliciting small effects on susceptibility to Mycobacterium avium spp paratuberculosis infection in dairy cattle. The results of this study will broadly impact domestic cattle research by providing the necessary foundation to explore several avenues of bovine translational genomics, and the potential for marker-assisted vaccination.

Altered patterns of toll-like receptor gene expression in cull cows infected with Mycobacterium avium subsp. paratuberculosis

Johne's disease caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic enteric disease of cattle. The mechanism how MAP can co-exist in the gastro-intestinal tract despite a massive infiltration of immune cells is not known. Toll-like receptors (TLRs) are known to play an important role in both innate and acquired immune responses but it is unclear what role different TLRs play in response to MAP. In this study, 38 cull cows from herds infected with MAP were classified into four groups, based on MAP culture from gut tissues and histopathological lesion scores. The expression of TLR1, 2 and 4 mRNA from MAP antigen-stimulated mesenteric lymph node (MLN) cultures and peripheral blood mononuclear cells (PBMCs) and in the MLN and ileum tissues of these animals was determined. MAP antigen-specific expression of TLR1 in MLN and PBMC was significantly lower in the MAP-infected groups than the non-infected control group, suggesting that in MAP-infected animals there is impairment in the up-regulation of TLR1 in response to MAP antigen. TLR4 expression in MLN tissues was significantly higher in the severely infected group than the control group suggesting up-regulation of endogenous TLR4 expression at a site of MAP infection in animals severely affected with Johne's disease. A preliminary screening of TLR1, 2 and 4 in the cull cows revealed the presence of polymorphisms in TLR1 and TLR2. In summary, one mechanism how MAP may subvert the immune system is that there is an apparent lack of recognition of MAP antigens as foreign by TLR1 in MAP-infected cows.