Genomic regions underlying susceptibility to bovine tuberculosis 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.

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

Exploring the genetic factors behind the discrepancy in resistance to bovine tuberculosis between African zebu cattle and European taurine cattle

Caused by the pathogenic agent Mycobacterium bovis, bovine tuberculosis (bTB) is a major concern in cattle breeding due to both its zoonotic potential and economic impact. Greater resistance to this disease has been reported in certain African zebu breeds compared to European taurine breeds. However the genetic basis for the lower susceptibility to bTB infection observed in zebu cattle remains poorly explored. This study was conducted on whole genome sequencing data of three bTB infection-resistant African zebu breeds and two bTB infection-susceptible taurine breeds to decipher the genetic background. A set of four selection signature statistics based on linkage disequilibrium, site frequency spectrum, and population differentiation were used on SNPs whereas between population variance based VST and t-test were used on CNVs. As a complement, genes from previous literature reported as candidate genes for bTB resistance were also inspected to identify genetic variations. Interestingly, the resulting nine candidate genes had deleterious missense variants (SHC3, IFNGR1, TLR2, TLR6, IL1A, LRRK2, EP300 and IRAK4) or a CNV difference (CD48) segregating between the groups. The genes found in the study play a role in immune pathways activated during Mycobacterium infection, contributing to the proliferation of immune cells and the granuloma formation, ultimately modulating the outcome of the infectious event. In particular, a deleterious variant in the LRRK2 gene, whose deficiency has been linked to improved prognosis upon tuberculosis infection, was found in the bTB infection-resistant zebu breeds. Therefore, these genes constitute credible candidates in explaining the discrepancy in Mycobacterium bovis infection susceptibility among different breed.

Selective breeding can contribute to bovine tuberculosis control and eradication

Bovine tuberculosis (bTB) persists in many countries having a significant impact on public health and livestock industry finances. The incidence and prevalence of new cases in parts of the UK and elsewhere over the past decades warrant intensified efforts towards achieving Officially Tuberculosis Free (OTF) status in the respective regions. Genetic selection aiming to identify and remove inherently susceptible animals from breeding has been proposed as an additional measure in ongoing programmes towards controlling the disease. The presence of genetic variation among individual animals in their capacity to respond to Mycobacterium bovis exposure has been documented and heritability estimates of 0.06-0.18 have been reported. Despite their moderate magnitude, these estimates suggest that host resistance to bTB is amenable to improvement with selective breeding. Although relatively slow, genetic progress can be constant, cumulative and permanent, thereby complementing ongoing disease control measures. Importantly, mostly no antagonistic genetic correlations have been found between bTB resistance and other animal traits suggesting that carefully incorporating the former in breeding decisions should not adversely affect bovine productivity. Simulation studies have demonstrated the potential impact of genetic selection on reducing the probability of a breakdown to occur or the duration and severity of a breakdown that has already been declared. Furthermore, research on the bovine genome has identified multiple genomic markers and genes associated with bTB resistance. Nevertheless, the combined outcomes of these studies suggest that host resistance to bTB is a complex, polygenic trait, with no single gene alone explaining the inherent differences between resistant and susceptible animals. Such results support the development of accurate genomic breeding values that duly capture the collective effect of multiple genes to underpin selective breeding programmes. In addition to improving host resistance to bTB, scientists and practitioners have considered the possibility of reducing host infectivity. Ongoing studies have suggested the presence of genetic variation for infectivity and confirmed that bTB eradication would be accelerated if selective breeding considered both host resistance and infectivity traits. In conclusion, research activity on bTB genetics has generated knowledge and insights to support selective breeding as an additional measure towards controlling and eradicating the disease.

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.

Mycobacterial Infection of Precision-Cut Lung Slices Reveals Type 1 Interferon Pathway Is Locally Induced by Mycobacterium bovis but Not M. tuberculosis in a Cattle Breed

Tuberculosis exacts a terrible toll on human and animal health. While Mycobacterium tuberculosis (Mtb) is restricted to humans, Mycobacterium bovis (Mb) is present in a large range of mammalian hosts. In cattle, bovine TB (bTB) is a noticeable disease responsible for important economic losses in developed countries and underestimated zoonosis in the developing world. Early interactions that take place between mycobacteria and the lung tissue early after aerosol infection govern the outcome of the disease. In cattle, these early steps remain poorly characterized. The precision-cut lung slice (PCLS) model preserves the structure and cell diversity of the lung. We developed this model in cattle in order to study the early lung response to mycobacterial infection. In situ imaging of PCLS infected with fluorescent Mb revealed bacilli in the alveolar compartment, in adjacent or inside alveolar macrophages, and in close contact with pneumocytes. We analyzed the global transcriptional lung inflammation signature following infection of PCLS with Mb and Mtb in two French beef breeds: Blonde d'Aquitaine and Charolaise. Whereas, lungs from the Blonde d'Aquitaine produced high levels of mediators of neutrophil and monocyte recruitment in response to infection, such signatures were not observed in the Charolaise in our study. In the Blonde d'Aquitaine lung, whereas the inflammatory response was highly induced by two Mb strains, AF2122 isolated from cattle in the UK and Mb3601 circulating in France, the response against two Mtb strains, H37Rv, the reference laboratory strain, and BTB1558, isolated from zebu in Ethiopia, was very low. Strikingly, the type I interferon pathway was only induced by Mb but not Mtb strains, indicating that this pathway may be involved in mycobacterial virulence and host tropism. Hence, the PCLS model in cattle is a valuable tool to deepen our understanding of early interactions between lung host cells and mycobacteria. It revealed striking differences between cattle breeds and mycobacterial strains. This model could help in deciphering biomarkers of resistance vs. susceptibility to bTB in cattle as such information is still critically needed for bovine genetic selection programs and would greatly help the global effort to eradicate bTB.

Genetic Analyses and Genome-Wide Association Studies on Pathogen Resistance of Bos taurus and Bos indicus Cattle Breeds in Cameroon

Autochthonous taurine and later introduced zebu cattle from Cameroon differ considerably in their resistance to endemic pathogens with little to no reports of the underlying genetic make-up. Breed history and habitat variations are reported to contribute significantly to this diversity worldwide, presumably in Cameroon as well, where locations diverge in climate, pasture, and prevalence of infectious agents. In order to investigate the genetic background, the genotypes of 685 individuals of different Cameroonian breeds were analysed by using the BovineSNP50v3 BeadChip. The variance components including heritability were estimated and genome-wide association studies (GWAS) were performed. Phenotypes were obtained by parasitological screening and categorised in Tick-borne pathogens (TBP), gastrointestinal nematodes (GIN), and onchocercosis (ONC). Estimated heritabilities were low for GIN and TBP (0.079 (se = 0.084) and 0.109 (se = 0.103) respectively) and moderate for ONC (0.216 (se = 0.094)). Further than revealing the quantitative nature of the traits, GWAS identified putative trait-associated genomic regions on five chromosomes, including the chromosomes 11 and 18 for GIN, 20 and 24 for TBP, and 12 for ONC. The results imply that breeding for resistant animals in the cattle population from Northern Cameroon might be possible for the studied pathogens; however, further research in this field using larger datasets will be required to improve the resistance towards pathogen infections, propose candidate genes or to infer biological pathways, as well as the genetic structures of African multi-breed populations.

Development of polymorphic markers in the immune gene complex loci of cattle

The addition of cattle health and immunity traits to genomic selection indices holds promise to increase individual animal longevity and productivity, and decrease economic losses from disease. However, highly variable genomic loci that contain multiple immune-related genes were poorly assembled in the first iterations of the cattle reference genome assembly and underrepresented during the development of most commercial genotyping platforms. As a consequence, there is a paucity of genetic markers within these loci that may track haplotypes related to disease susceptibility. By using hierarchical assembly of bacterial artificial chromosome inserts spanning 3 of these immune-related gene regions, we were able to assemble multiple full-length haplotypes of the major histocompatibility complex, the leukocyte receptor complex, and the natural killer cell complex. Using these new assemblies and the recently released ARS-UCD1.2 reference, we aligned whole-genome shotgun reads from 125 sequenced Holstein bulls to discover candidate variants for genetic marker development. We selected 124 SNPs, using heuristic and statistical models to develop a custom genotyping panel. In a proof-of-principle study, we used this custom panel to genotype 1,797 Holstein cows exposed to bovine tuberculosis (bTB) that were the subject of a previous GWAS study using the Illumina BovineHD array. Although we did not identify any significant association of bTB phenotypes with these new genetic markers, 2 markers exhibited substantial effects on bTB phenotypic prediction. The models and parameters trained in this study serve as a guide for future marker discovery surveys particularly in previously unassembled regions of the cattle genome.

Detection of genomic regions underlying milk production traits in Valle del Belice dairy sheep using regional heritability mapping

The aim of this study was to identify genomic regions underlying milk production traits in the Valle del Belice dairy sheep using regional heritability mapping (RHM). Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY), collected over a period of 6 years (2006–2012) on 481 Valle del Belice ewes, were used for the analysis. Animals were genotyped with the Illumina 50k SNP chip. Variance components, heritabilities and repeatabilities within and across lactations were estimated, fitting parity, litter size, season of lambing and fortnights in milk, as fixed; and additive genetic, permanent environment within and across lactations, flock by test‐day interaction and residual as random effects. For the RHM analysis, the model included the same fixed and random effects as before, plus an additional regional genomic additive effect (specific for the region being tested) as random. While the whole genomic additive effect was estimated using the genomic relationship matrix (GRM) constructed from all SNPs, the regional genomic additive effect was estimated from a GRM matrix constructed from the SNPs within each region. Heritability estimates ranged between 0.06 and 0.15, with repeatabilities being between 0.14 and 0.24 across lactations and between 0.23 and 0.39 within lactation for all milk production traits. A substantial effect of flock‐test‐day on milk production traits was also estimated. Significant genomic regions at either genome‐wide (p < .05) or suggestive (i.e., one false positive per genome scan) level were identified on chromosome (OAR) 2, 3 and 20 for F% and on OAR3 for P%, with the regions on OAR3 in common between the two traits. Our results confirmed the role of LALBA and AQP genes, on OAR3, as candidate genes for milk production traits in sheep.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

High-throughput genotype-based population structure analysis of selected buffalo breeds

India is considered as the home tract of some of the best buffalo breeds. However, the genetic structure of the Indian river buffalo is poorly understood. Hence, there is a need to characterize the populations and understand the genetic structure of various buffalo breeds for selection and to design breeding strategies. In this study, we have analyzed genetic variability and population structure of seven buffalo breeds from their respective geographical regions using Axiom Buffalo Genotyping Array. Diversity, as measured by expected heterozygosity, ranged from 0.364 in Surti to 0.384 in Murrah breed, and pair-wise F_ST values revealed the lowest genetic distance between Murrah and Nili-Ravi (0.0022), while the highest between Surti and Pandharpuri (0.030). Principal component analysis and structure analysis unveiled the differentiation of Surti, Pandharpuri, and Jaffarabadi in first two principal components and at K = 4, respectively, while remaining breeds were grouped together as a separate single cluster and admixed. Murrah and Mehsana showed early linkage disequilibrium (LD) decay, while Surti breed showed late decay. In LD blocks to quantitative trait locis (QTLs) concordance analysis, 4.65% of concordance was observed with 873 LD blocks overlapped with 2,330 QTLs. Overall, total 4,090 markers were identified from all LD blocks for six types of traits. Results of this study indicated that these single-nucleotide polymorphism (SNP) markers could differentiate phenotypically distinct breeds like Surti, Pandharpuri, and Jaffarabadi but not others. So, there is a need to develop SNP chip based on SNP markers identified by sequence information of local breeds.

Genetic Diversity of Cameroon Cattle and a Putative Genomic Map for Resistance to Bovine Tuberculosis

Bovine Tuberculosis (bTB) caused by Mycobacterium bovis is a livestock disease of global economic and public health importance. There are currently no effective vaccines available for livestock and so control relies on animal level surveillance and pasteurization of dairy products. A new alternative control approach is to exploit the genetic variability of the host; recent studies have demonstrated that breeding European taurine cattle, such as Holsteins for increased resistance to bTB is feasible. The utility of such an approach is still unknown for African cattle populations. This study aims to assess genetic variation in bTB resistance and the underlying genomic architecture of cattle in Cameroon. We conducted a cross-sectional study of 2,346 slaughter cattle in Cameroon. Retropharyngeal lymph node samples were collected and cultured on Lowenstein Jensen media and the BACTEC MGIT 960 system, and M. bovis was identified using the Hain® Genotype kits. A total of 153 cattle were positive for M. bovis and were archived along with a random selection of negative samples. In this study, we genotyped archived samples from 212 cattle. Their genomic diversity was characterized using PCA, hierarchical clustering and admixture analysis. We assessed genetic variation in bTB resistance using heritability analysis and compared quantitative trait loci. Previous research on this study population have shown that Fulani cattle are more susceptible to bTB than mixed breeds. However, here we show that these apparent phenotypic differences in breeds are not reflected by clear genomic differences. At the genetic level, both the Fulani and mixed cattle show similar patterns of admixture with evidence of both taurine and indicine ancestry. There was little European taurine introgression within the studied population. Hierarchical clustering showed clusters of cattle that differed in their susceptibility to bTB. Our findings allude to bTB resistance being polygenic in nature. This study highlights the potential for genetic control of bTB in Africa and the need for further research into the genetics of bTB resistance within African cattle populations.

Analysis of Genetic Variation in the Bovine SLC11A1 Gene, Its Influence on the Expression of NRAMP1 and Potential Association With Resistance to Bovine Tuberculosis

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is a chronic zoonotic disease where host genetics is thought to contribute to susceptibility or resistance. One of the genes implicated is the SLC11A1 gene, that encodes for the natural resistance-associated macrophage protein 1 (NRAMP1). The aim of this study was to identify SLC11A1 polymorphisms and to investigate any resulting functional differences in NRAMP1 expression that might be correlated with resistance/susceptibility to M. bovis infection. Sequencing of the SLC11A1 gene in cDNA isolated from Brown Swiss, Holstein Friesian, and Sahiwal cattle identified five single nucleotide polymorphisms (SNPs) in the coding region, but only one of these (SNP4, c.1066C>G, rs109453173) was present in all three cattle breeds and therefore warranted further investigation. Additionally, variations of 10, 11, and 12 GT repeats were identified in a microsatellite (MS1) in the SLC11A1 3′UTR. Measurement of NRAMP1 expression in bovine macrophages by ELISA showed no differences between cells generated from the different breeds. Furthermore, variations in the length of the MS1 microsatellite did not impact on NRAMP1 protein expression as analyzed by luciferase reporter assay. However, further analysis of the ELISA data identified that the presence of the alternative G allele at SNP4 was associated with increased expression of NRAMP1 in bovine macrophages. Since NRAMP1 has been shown to influence the survival of intracellular pathogens such as M. bovis through the sequestering of iron, it is possible that cattle expressing the alternative G allele might have an increased resistance to bTB through increased NRAMP1 expression in their macrophages.

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.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Simple Summary

Bovine tuberculosis is an infectious disease of cattle caused by Mycobacterium bovis characterized by the formation of tubercles in any organ or tissue. Bovine tuberculosis represents a significant veterinary and public health problem in many parts of the world. It is zoonotic, transmitted to humans through consumption of infected milk and other cattle products. Although many factors influence infection and progression of the disease, there must be an important host genetic component that explains why some animals get sick and others remain healty. We present evidence of genetic variants associated with resistance to tuberculosis in Mexican Holstein dairy cattle using a case-control approach with a selective DNA pooling. Here, we identified novel quantitative trait loci regions harboring genes involved in Mycobacterium spp. immune response. This is a first screening about resistance to tuberculosis infection on Mexican Holstein cattle based on a dense single nucleotide polymorphism chip. The identified genes belong to both, the already known, and the undisclosed quantitative trait loci regions.

Abstract

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

Impact of Genetic Selection for Increased Cattle Resistance to Bovine Tuberculosis on Disease Transmission Dynamics

Bovine tuberculosis (bTB) poses a challenge to animal health and welfare worldwide. Presence of genetic variation in host resistance to Mycobacterium bovis infection makes the trait amenable to improvement with genetic selection. Genetic evaluations for resistance to infection in dairy cattle are currently available in the United Kingdom (UK), enabling genetic selection of more resistant animals. However, the extent to which genetic selection could contribute to bTB eradication is unknown. The objective of this study was to quantify the impact of genetic selection for bTB resistance on cattle-to-cattle disease transmission dynamics and prevalence by developing a stochastic genetic epidemiological model. The model was used to implement genetic selection in a simulated cattle population. The model considered various levels of selection intensity over 20 generations assuming genetic heterogeneity in host resistance to infection. Our model attempted to represent the dairy cattle population structure and current bTB control strategies in the UK, and was informed by genetic and epidemiological parameters inferred from data collected from UK bTB infected dairy herds. The risk of a bTB breakdown was modeled as the percentage of herds where initially infected cows (index cases) generated secondary cases by infecting herd-mates. The model predicted that this risk would be reduced by half after 4, 6, 9, and 15 generations for selection intensities corresponding to genetic selection of the 10, 25, 50, and 70% most resistant sires, respectively. In herds undergoing bTB breakdowns, genetic selection reduced the severity of breakdowns over generations by reducing both the percentage of secondary cases and the duration over which new secondary cases were detected. Selection of the 10, 25, 50, and 70% most resistant sires reduced the percentage of secondary cases to <1% in 4, 5, 7, and 11 generations, respectively. Similarly, the proportion of long breakdowns (breakdowns in which secondary cases were detected for more than 365 days) was reduced by half in 2, 2, 3, and 4 generations, respectively. Collectively, results suggest that genetic selection could be a viable tool that can complement existing management and surveillance methods to control and ultimately eradicate bTB.

Epidemiology of Mycobacterium bovis and Mycobacterium tuberculosis in animals: Transmission dynamics and control challenges of zoonotic TB in Ethiopia

Mycobacterium tuberculosis complex is the cause of tuberculosis (TB) in humans and other animals. Specifically, Mycobacterium bovis (M. bovis) and Mycobacterium tuberculosis (M. tuberculosis) are highly pathogenic mycobacteria that may infect different animal species and are the sources of TB in humans. The objective of this paper was to review the epidemiology of M. bovis and M. tuberculosis in animals. The review also highlighted the transmission dynamics of M. bovis and M. tuberculosis in humans and animals and control challenges of zoonotic TB in Ethiopia. The literature review focused on scientific peer-reviewed articles from studies exclusively conducted in Ethiopia that were published from 1998 to 2017. Husbandry system, breed and herd size have significant role in the epidemiology of bovine tuberculosis (BTB) in Ethiopia. The information presented reveals that different strains of M. bovis are widely distributed in domestic animals predominantly in the Ethiopian cattle and the main strain was found to be SB1176. In addition, the isolation of M. tuberculosis from domestic animals in different settings signifies the circulation of the agent between humans and animals in Ethiopia. The life styles of the Ethiopian communities, close contact with domestic animals and/or the habit of consuming raw animal products, are suggested as the main factors for transmission of M. bovis and M. tuberculosis between human and animal which may have impact on the TB control program in human. In Ethiopia, a human TB control program has been widely implemented, however, the role of animal in the transmission of the causative agent has been neglected which could be one of the challenges for an effective control program. This warrants the need for incorporating animal TB control programs using "One Health" approach for effective TB control for both human and animal.