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.

Structural variation and introgression from wild populations in East Asian cattle genomes confer adaptation to local environment.

BACKGROUND: Structural variations (SVs) in individual genomes are major determinants of complex traits, including adaptability to environmental variables. The Mongolian and Hainan cattle breeds in East Asia are of taurine and indicine origins that have evolved to adapt to cold and hot environments, respectively. However, few studies have investigated SVs in East Asian cattle genomes and their roles in environmental adaptation, and little is known about adaptively introgressed SVs in East Asian cattle. RESULTS: In this study, we examine the roles of SVs in the climate adaptation of these two cattle lineages by generating highly contiguous chromosome-scale genome assemblies. Comparison of the two assemblies along with 18 Mongolian and Hainan cattle genomes obtained by long-read sequencing data provides a catalog of 123,898 nonredundant SVs. Several SVs detected from long reads are in exons of genes associated with epidermal differentiation, skin barrier, and bovine tuberculosis resistance. Functional investigations show that a 108-bp exonic insertion in SPN may affect the uptake of Mycobacterium tuberculosis by macrophages, which might contribute to the low susceptibility of Hainan cattle to bovine tuberculosis. Genotyping of 373 whole genomes from 39 breeds identifies 2610 SVs that are differentiated along a "north-south" gradient in China and overlap with 862 related genes that are enriched in pathways related to environmental adaptation. We identify 1457 Chinese indicine-stratified SVs that possibly originate from banteng and are frequent in Chinese indicine cattle. CONCLUSIONS: Our findings highlight the unique contribution of SVs in East Asian cattle to environmental adaptation and disease resistance.

Deletion in KARLN intron 5 and predictive relationship with bovine tuberculosis and brucellosis infection phenotype.

Kalirin (gene: KALRN) is a Rho-GEF kinase linked to neurodegenerative diseases in humans. Unexpectedly, various polymorphisms in KALRN gene were previously associated with resistance to bacterial infections in ruminants. In this study, we evaluated the effect of the rs384223075 (RS-075) deletion in KALRN intron 5 on the occurrence of Mycobacterium bovis and Brucella abortus infections in cattle. We performed two separate case-control association analyses: one for bovine tuberculosis (bTB) using 308 Holstein and Jersey cows from three herds with prevalence between 5 and 15% for this infection; and another for brucellosis using 140 Holstein and beef crossbred cows from two herds with high prevalence for brucellosis (> 30%). In the bTB analysis, the RS-075 deletion frequency was higher among cases than controls (p = 0.0001), and the absence of the RS-075 deletion allele was associated with negative PPD-skin test results (p = 0.0009) at genotype level. On the contrary, RS-075 was not associated with Brucella spp. serological status (p = 0.72) but, unexpectedly, the deletion allele was more frequent among controls than cases in the beef crossbred herd (0.31 vs. 0.14, p = 0.02). In concordance with this observation, in vitro assays showed that the RS-075 deletion could be linked to an enhanced cellular response to bacterial antigens and unspecific stimulation in mononuclear cells derived from beef crossbred cows, specifically the reactive nitrogen species production (p = 0.008) and proliferation capacity (p = 0.018). This study is consistent with other reports that support an important role of the KALRN gene and its polymorphisms in the host response to intracellular pathogens.

High-resolution transcriptomics of bovine purified protein derivative-stimulated peripheral blood from cattle infected with Mycobacterium bovis across an experimental time course.

OBJECTIVES: Improved bovine tuberculosis (bTB) diagnostics with higher sensitivity and specificity are urgently required. A better understanding of the peripheral blood transcriptional response of Mycobacterium bovis-infected animals after bovine purified protein derivative (PPD-b) stimulation of whole blood-an important component of current bTB diagnostics-will provide new information for development of better diagnostics. METHODS: RNA sequencing (RNA-seq) was used to study the peripheral blood transcriptome after stimulation with PPD-b across four time points (-1 wk pre-infection, and +1 wk, +2 wk, and +10 wk post-infection) from a 14-week M. bovis infection time course experiment with ten age-matched Holstein-Friesian cattle. RESULTS: In vitro PPD-b stimulation of peripheral blood from M. bovis-infected and non-infected cattle elicited a strong transcriptional response. Comparison of PPD-b stimulated, and unstimulated samples revealed higher expression of genes encoding cytokine receptors, transcription factors, and interferon-inducible proteins. Lower expression was seen for genes encoding proteins involved in antimicrobial activity, C-type lectin receptors, inhibition of signal transduction, and genes encoding metal ion transporters. CONCLUSIONS: A transcriptional signature associated with the peripheral blood response to PPD-b stimulation consisting of 170 genes was identified exclusively in the post-infection time points. Therefore, this represents a panel of potential biomarkers of M. bovis infection.

Defining the Genes Required for Survival of Mycobacterium bovis in the Bovine Host Offers Novel Insights into the Genetic Basis of Survival of Pathogenic Mycobacteria.

Tuberculosis has severe impacts on both humans and animals. Understanding the genetic basis of survival of both Mycobacterium tuberculosis, the human-adapted species, and Mycobacterium bovis, the animal-adapted species, is crucial to deciphering the biology of both pathogens. There are several studies that identify the genes required for survival of M. tuberculosis in vivo using mouse models; however, there are currently no studies probing the genetic basis of survival of M. bovis in vivo. In this study, we utilize transposon insertion sequencing in M. bovis AF2122/97 to determine the genes required for survival in cattle. We identify genes encoding established mycobacterial virulence functions such as the ESX-1 secretion system, phthiocerol dimycocerosate (PDIM) synthesis, mycobactin synthesis, and cholesterol catabolism that are required in vivo. We show that, as in M. tuberculosis H37Rv, phoPR is required by M. bovis AF2122/97 in vivo despite the known defect in signaling through this system. Comparison to studies performed in species that are able to use carbohydrates as an energy source, such as M. bovis BCG and M. tuberculosis, suggests that there are differences in the requirement for genes involved in cholesterol import (mce4 operon) and oxidation (hsd). We report a good correlation with existing mycobacterial virulence functions but also find several novel virulence factors, including genes involved in protein mannosylation, aspartate metabolism, and glycerol-phosphate metabolism. These findings further extend our knowledge of the genetic basis of survival in vivo in bacteria that cause tuberculosis and provide insight for the development of novel diagnostics and therapeutics. IMPORTANCE This is the first report of the genetic requirements of an animal-adapted member of the Mycobacterium tuberculosis complex (MTBC) in a natural host. M. bovis has devastating impacts on cattle, and bovine tuberculosis is a considerable economic, animal welfare, and public health concern. The data highlight the importance of mycobacterial cholesterol catabolism and identify several new virulence factors. Additionally, the work informs the development of novel differential diagnostics and therapeutics for TB in both human and animal populations.

Association of NOS2A gene polymorphisms with susceptibility to bovine tuberculosis in Chinese Holstein cattle.

Bovine tuberculosis (bTB) is a global zoonotic disease that has detrimental economic impacts worldwide. The NOS2A gene plays a key role in immunological control of many infectious diseases. However, research on the association between NOS2A polymorphisms and bTB infection in Holstein cattle reared on the Yunnan-Guizhou plateau of China is scarce. This study investigated a possible linkage between NOS2A polymorphisms and risk of developing bTB in Chinese Holstein cattle. The NOS2A gene was genotyped in 144 bTB-infected Holstein cows and 139 healthy controls were genotyped through nucleotide sequencing. Ten single-nucleotide polymorphisms (SNPs) were detected, six of which were associated with susceptibility/resistance patterns of bTB. Furthermore, the C/T genotypes of 671 and 2793, and T/T genotype of E22 (+15) were significantly associated with susceptibility risk; the G/A genotype of 2857, T/T genotype of E9 (+65), and C/C genotype of E9 (+114) probably increased resistance to bTB. In addition, the haplotypes of NOS2A-2 and NOS2A-9 were risk factors for bTB susceptibility, while the NOS2A-5 and NOS2A-8 haplotypes were contributing protective variants against tuberculosis. There is a significant association between variation in SNPs of NOS2A and tuberculosis susceptibility/resistance pattern. These findings suggest that substitution of genetic selection would be helpful for eradicating bTB. However, further investigation is required to study the underlying mechanism through which NOS2A polymorphisms affect bTB infection.

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.

Exploring the Potential of Interferon Gamma Gene as Major Immune Responder for Bovine Tuberculosis in River Buffalo.

Bovine tuberculosis (bTB) is a widespread zoonotic infection targeting the livestock sector, especially in developing countries, and posing a risk to humans and animal populations. Its recent prevalence in river buffaloes has been estimated as higher as 33.7%. In emergent countries like Pakistan, there is likeliness of human-livestock interfaces extensively and lacking of effective preventive measures that illustrate the risk of spreading the infection at a remarkable rate. The river buffalo (Bubalus bubalis) is an upkeep host of Mycobacterium bovis and is responsible for disease transmission among buffaloes and other livestock species. In this study, potential molecular biomarkers in the Interferon-gamma gene (IFNg) were identified after genomic screening of river buffaloes. Unique genomic loci in river buffalo proved the novelty of the genomic structure of this phenomenal animal but also highlighted its significance in natural immunity against the Mycobacterium. A total of eight single nucleotide polymorphisms were identified in the coding region of IFNg. The SNPs in the exonic region were all transitions, i.e., the conversion of purines to purines. These SNPs were analyzed for Hardy Weinberg Equilibrium, chi2 test, gene diversity, and protein structural conformation. Pathway analysis in tuberculosis revealed that IFNg inhibits the antigen-presenting cells (APC) through JAK and STAT pathways. Network analysis of IFNg proteins in both species showed strong associations among the immunity-related proteins (interleukins, tissue necrosis factors) and receptors of interferons. The identified polymorphic sites might be novel-potentiated markers for the selection of animals with superior immune response against bTB and can be exploited as promising genomic sites for breeding the resistant animal herds to combat Mycobacterium infection in a long run.

HMEJ-based safe-harbor genome editing enables efficient generation of cattle with increased resistance to tuberculosis.

The CRISPR/Cas9 system has been used in a wide range of applications in the production of gene-edited animals and plants. Most efforts to insert genes have relied on homology-directed repair (HDR)-mediated integration, but this strategy remains inefficient for the production of gene-edited livestock, especially monotocous species such as cattle. Although efforts have been made to improve HDR efficiency, other strategies have also been proposed to circumvent these challenges. Here we demonstrate that a homology-mediated end-joining (HMEJ)-based method can be used to create gene-edited cattle that displays precise integration of a functional gene at the ROSA26 locus. We found that the HMEJ-based method increased the knock-in efficiency of reporter genes by eightfold relative to the traditional HDR-based method in bovine fetal fibroblasts. Moreover, we identified the bovine homology of the mouse Rosa26 locus that is an accepted genomic safe harbor and produced three live-born gene-edited cattle with higher rates of pregnancy and birth, compared with previous work. These gene-edited cattle exhibited predictable expression of the functional gene natural resistance-associated macrophage protein-1 (NRAMP1), a metal ion transporter that should and, in our experiments does, increase resistance to bovine tuberculosis, one of the most detrimental zoonotic diseases. This research contributes to the establishment of a safe and efficient genome editing system and provides insights for gene-edited animal breeding.

Long-term molecular surveillance provides clues on a cattle origin for Mycobacterium bovis in Portugal.

Animal tuberculosis (TB), caused by Mycobacterium bovis, is maintained in Portugal in a multi-host system, with cattle, red deer and wild boar, playing a central role. However, the ecological processes driving transmission are not understood. The main aim of this study was thus to contribute to the reconstruction of the spatiotemporal history of animal TB and to refine knowledge on M. bovis population structure in order to inform novel intervention strategies. A collection of 948 M. bovis isolates obtained during long-term surveillance (2002-2016, 15 years) of cattle (n = 384), red deer (n = 303) and wild boar (n = 261), from the main TB hotspot areas, was characterized by spoligotyping and 8 to 12-loci MIRU-VNTR. Spoligotyping identified 64 profiles and MIRU-VNTR distinguished 2 to 36 subtypes within each spoligotype, enabling differentiation of mixed or clonal populations. Common genotypic profiles within and among livestock and wildlife in the same spatiotemporal context highlighted epidemiological links across hosts and regions, as for example the SB0119-M205 genotype shared by cattle in Beja district or SB0121-M34 shared by the three hosts in Castelo Branco and Beja districts. These genomic data, together with metadata, were integrated in a Bayesian inference framework, identifying five ancestral M. bovis populations. The phylogeographic segregation of M. bovis in specific areas of Portugal where the disease persists locally is postulated. Concurrently, robust statistics indicates an association of the most probable ancient population with cattle and Beja, providing a clue on the origin of animal TB epidemics. This relationship was further confirmed through a multinomial probability model that assessed the influence of host species on spatiotemporal clustering. Two significant clusters were identified, one that persisted between 2004 and 2010, in Beja district, with Barrancos county at the centre, overlapping the central TB core area of the Iberian Peninsula, and highlighting a significant higher risk associated to cattle. The second cluster was predominant in the 2012-2016 period, holding the county Rosmaninhal at the centre, in Castelo Branco district, for which wild boar contributed the most in relative risk. These results provide novel quantitative insights beyond empirical perceptions, that may inform adaptive TB control choices in different regions.

A molecular strategy to optimize bovine tuberculosis post-mortem diagnosis and the exposure to Mycobacterium tuberculosis variant bovis.

Bovine tuberculosis (bTB) is a zoonotic disease caused by Mycobacterium tuberculosis var. bovis, for which the definitive diagnosis is accomplished by bacterial isolation, which has biosafety issues and requires long time. Thus, diagnostic methods with potential to be faster and more efficient can represent an advance in bTB epidemiological knowledge and decrease exposure to M. tuberculosis var. bovis. This study aimed to validate a molecular test for bTB post-mortem diagnosis, as a strategy to reduce waste in bovine production. A total of 185 tissues from animals of infected herds or with suspected lesions at abattoir were evaluated through bacterial isolation, PCR and histopathology. PCR and histopathology showed sensitivities of 45.1% and 71.2%, respectively, and specificities of 83.3% and 83.0%, respectively, when compared to bacterial isolation. The combination of both tests resulted in enhanced specificity and positive predictive values.Therefore, PCR in conjunction with histopathology may be used as screening, in which concordant results can be considered conclusive, and discordant results may be submitted to bacterial isolation.

Poor stimulation of bovine dendritic cells by Mycobacterium bovis culture supernatant and surface extract is associated with decreased activation of ERK and NF-κB and higher expression of SOCS1 and 3.

The cell envelope of pathogenic mycobacteria interfaces with the host. As such, the interaction of bacterial products localized at or released from the cell surface with the host's immune system can determine the fate of the bacterium in its host. In this study, the effects of three different types of Mycobacterium bovis cell envelope fractions-purified protein derivative, total cell wall lipids and culture supernatant and surface extract-on bovine dendritic cells were assessed. We found that the culture supernatant and surface extract fraction induced little to no production of the pro-inflammatory cytokines TNF-α and IL-12 in bovine dendritic cells. Moreover, this muted response was associated with poor activation of ERK and NF-κB, both of which are critical for the pro-inflammatory response. Furthermore, culture supernatant and surface extract treatment increased the expression of suppressor of cytokine signaling 1 and 3, both of which are negative regulators of pro-inflammatory signaling, in bovine dendritic cells. These observations taken together suggest the M. bovis culture supernatant and surface extract fraction contain immunomodulatory molecules that may aid in M. bovis pathogenesis.

BoLA-DRB3 exon2 polymorphisms among tuberculous cattle: Nucleotide and functional variability and their association with bovine tuberculosis pathology.

Bovine tuberculosis (bTB) is caused by Mycobacterium bovis and disseminated worldwide. In Argentina, the highest prevalence occurs in dairy areas. BoLA DRB3.2 is related to the adaptive immunity in mycobacterial infections. Genetic polymorphisms of this marker have been associated with resistance or susceptibility to bovine diseases. We evaluated the association between BoLA DRB3.2 polymorphisms and bTB pathology scores in dairy and beef cattle breeds of Argentina. Most bovines exhibited visible lesions compatible with tuberculosis and, furthermore, 150 (85.7%) were also positive by bacteriology. A pathology index showed a variable degree of disease, from 3 to 76 (median pathology score = 9 (IQR: 7-15)). Thirty-five BoLA DRB3.2 alleles were identified with an associated frequency from 16% to 0.3%, distributed 73% (n = 128) in heterozygosis and 27% (n = 47) in homozygosis, with 12 BoLA DRB3.2 alleles (*0101, *1101, *1501, *0201, *2707 *1001, *1002, *1201, *14011, *0501 *0902 and *0701) representing the 74.7% of the population variability. A functional analysis grouped them in 4 out of 5 clusters (A-D), suggesting a functional overlapping. Among the 90 identified genotypes, *1101/*1101, *1101/*1501 and *0101/*0101 were the most frequent (10%, 8.9% and 8.9%, respectively). No association was detected between the pathology scores and a specific DRB3.2 allele (p > .05). Animals infected with M. bovis spoligotype SB0153 showed a significantly higher pathology score than those affected by the spoligotype SB0145 (p = .018). Furthermore, the Aberdeen Angus breed exhibited highest pathological scores (p < .0001), which were associated with disseminated lesion, thus suggesting that the host component could be important to the disease progression.

A natural gene drive system influences bovine tuberculosis susceptibility in African buffalo: Possible implications for disease management.

Bovine tuberculosis (BTB) is endemic to the African buffalo (Syncerus caffer) of Hluhluwe-iMfolozi Park (HiP) and Kruger National Park, South Africa. In HiP, the disease has been actively managed since 1999 through a test-and-cull procedure targeting BTB-positive buffalo. Prior studies in Kruger showed associations between microsatellite alleles, BTB and body condition. A sex chromosomal meiotic drive, a form of natural gene drive, was hypothesized to be ultimately responsible. These associations indicate high-frequency occurrence of two types of male-deleterious alleles (or multiple-allele haplotypes). One type negatively affects body condition and BTB resistance in both sexes. The other type has sexually antagonistic effects: negative in males but positive in females. Here, we investigate whether a similar gene drive system is present in HiP buffalo, using 17 autosomal microsatellites and microsatellite-derived Y-chromosomal haplotypes from 401 individuals, culled in 2002-2004. We show that the association between autosomal microsatellite alleles and BTB susceptibility detected in Kruger, is also present in HiP. Further, Y-haplotype frequency dynamics indicated that a sex chromosomal meiotic drive also occurred in HiP. BTB was associated with negative selection of male-deleterious alleles in HiP, unlike positive selection in Kruger. Birth sex ratios were female-biased. We attribute negative selection and female-biased sex ratios in HiP to the absence of a Y-chromosomal sex-ratio distorter. This distorter has been hypothesized to contribute to positive selection of male-deleterious alleles and male-biased birth sex ratios in Kruger. As previously shown in Kruger, microsatellite alleles were only associated with male-deleterious effects in individuals born after wet pre-birth years; a phenomenon attributed to epigenetic modification. We identified two additional allele types: male-specific deleterious and beneficial alleles, with no discernible effect on females. Finally, we discuss how our findings may be used for breeding disease-free buffalo and implementing BTB test-and-cull programs.

Atypical granuloma formation in Mycobacterium bovis-infected calves.

Bovine tuberculosis is a chronic inflammatory disease that causes granuloma formation. Characterization of granulomatous lesions of Mycobacterium bovis (M. bovis) experimentally infected cattle has helped to better understand the pathogenesis of this disease. However, few studies have described granulomas found in M. bovis naturally infected cattle. The aim of this work was to examine granulomas from Holstein-Friesian cattle naturally infected with M. bovis from a dairy basin located in the central region of Mexico. Tissue samples from thirty-two cattle with lesions suggestive of tuberculosis were collected post-mortem. Fifteen of the 32 sampled animals (46.8%) were 4 months of age or younger (calves), whereas the rest (53.2%, 17/32) were over one year old (adults). Macroscopic lesions suggestive of tuberculosis were found in the mediastinal lymph node chain of all animals (32/32). From the 1,143 granulomatous lesions that were microscopically analyzed, 34.6% (396/1143) were collected from adult animals and subsequently classified according to the nomenclature suggested by Wangoo et al., 2005. Surprisingly, lesions from calf tissues showed an atypical pattern which could not be fitted into the established developmental stages of this classification. Granulomatous lesions found in calves covered most of the affected organ, histologically showed large necrotic areas with central calcification, absence of a connective tissue capsule, and few giant cells. Also, there was a higher percentage of lesions with acid-fast bacilli (AFB) when compared to studied granulomas in adults. Growth of Mycobacterium spp was detected in 11 bacteriological tissue cultures. Genotypic identification of M. bovis was performed by DNA extraction from bacterial isolates, formalin-fixed and paraffin-embedded (FFPE) tissues and samples without bacterial isolation. M. bovis was detected by PCR in 84.3% (27/32) of the studied cases; whereas other AFB were observed in tissues of the remaining sampled animals (5/32). Our results describe atypical granuloma formation in calves 4 months of age or younger, naturally infected with M. bovis. These findings contribute to better understanding the physiopathology of M. bovis infection in cattle.

The bovine alveolar macrophage DNA methylome is resilient to infection with Mycobacterium bovis.

DNA methylation is pivotal in orchestrating gene expression patterns in various mammalian biological processes. Perturbation of the bovine alveolar macrophage (bAM) transcriptome, due to Mycobacterium bovis (M. bovis) infection, has been well documented; however, the impact of this intracellular pathogen on the bAM epigenome has not been determined. Here, whole genome bisulfite sequencing (WGBS) was used to assess the effect of M. bovis infection on the bAM DNA methylome. The methylomes of bAM infected with M. bovis were compared to those of non-infected bAM 24 hours post-infection (hpi). No differences in DNA methylation (CpG or non-CpG) were observed. Analysis of DNA methylation at proximal promoter regions uncovered >250 genes harbouring intermediately methylated (IM) promoters (average methylation of 33-66%). Gene ontology analysis, focusing on genes with low, intermediate or highly methylated promoters, revealed that genes with IM promoters were enriched for immune-related GO categories; this enrichment was not observed for genes in the high or low methylation groups. Targeted analysis of genes in the IM category confirmed the WGBS observation. This study is the first in cattle examining genome-wide DNA methylation at single nucleotide resolution in an important bovine cellular host-pathogen interaction model, providing evidence for IM promoter methylation in bAM.

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.

Comparative Study of the Molecular Basis of Pathogenicity of M. bovis Strains in a Mouse Model.

It is widely accepted that different strains of Mycobacterium tuberculosis have variable degrees of pathogenicity and induce different immune responses in infected hosts. Similarly, different strains of Mycobacterium bovis have been identified but there is a lack of information regarding the degree of pathogenicity of these strains and their ability to provoke host immune responses. Therefore, in the current study, we used a mouse model to evaluate various factors involved in the severity of disease progression and the induction of immune responses by two strains of M. bovis isolated from cattle. Mice were infected with both strains of M. bovis at different colony-forming unit (CFU) via inhalation. Gross and histological findings revealed more severe lesions in the lung and spleen of mice infected with M. bovis N strain than those infected with M. bovis C68004 strain. In addition, high levels of interferon-γ (IFN-γ), interleukin-17 (IL-17), and IL-22 production were observed in the serum samples of mice infected with M. bovis N strain. Comparative genomic analysis showed the existence of 750 single nucleotide polymorphisms and 145 small insertions/deletions between the two strains. After matching with the Virulence Factors Database, mutations were found in 29 genes, which relate to 17 virulence factors. Moreover, we found an increased number of virulent factors in M. bovis N strain as compared to M. bovis C68004 strain. Taken together, our data reveal that variation in the level of pathogenicity is due to the mutation in the virulence factors of M. bovis N strain. Therefore, a better understanding of the mechanisms of mutation in the virulence factors will ultimately contribute to the development of new strategies for the control of M. bovis infection.

Identifying genotype specific elevated-risk areas and associated herd risk factors for bovine tuberculosis spread in British cattle.

Bovine tuberculosis (bTB) is a chronic zoonosis with major health and economic impact on the cattle industry. Despite extensive control measures in cattle and culling trials in wildlife, the reasons behind the expansion of areas with high incidence of bTB breakdowns in Great Britain remain unexplained. By balancing the importance of cattle movements and local transmission on the observed pattern of cattle outbreaks, we identify areas at elevated risk of infection from specific Mycobacterium bovis genotypes. We show that elevated-risk areas (ERAs) were historically more extensive than previously understood, and that cattle movements alone are insufficient for ERA spread, suggesting the involvement of other factors. For all genotypes, we find that, while the absolute risk of infection is higher in ERAs compared to areas with intermittent risk, the statistically significant risk factors are remarkably similar in both, suggesting that these risk factors can be used to identify incipient ERAs before this is indicated by elevated incidence alone. Our findings identify research priorities for understanding bTB dynamics, improving surveillance and guiding management to prevent further ERA expansion.