Functional analysis of the Mycobacterium bovis AF2122/97 PhoPR system.

The PhoPR system is a master regulator in Mycobacterium tuberculosis. A key difference between M. tuberculosis and Mycobacterium bovis is a G71I substitution in the M. bovis PhoR orthologue. Functional studies of the M. bovis PhoPR system have generated conflicting findings, with some research suggesting that the M. bovis PhoPR is defective while others indicate it is functional. We sought to revisit the functionality of the M. bovis PhoPR system. To address this, we constructed a phoPR mutant in the reference strain M. bovis AF2122/97. We employed a combination of growth assays and transcriptomics analyses to assess the phenotype of the mutant vs wild type and complemented strains. We found that the M. bovis AF2122/97 ΔphoPR mutant showed a growth defect on solid and liquid media compared to the wild type and complemented strains. The transcriptome of the M. bovis AF2122/97 ΔphoPR mutant was also altered as compared to wild type, including differential expression of genes involved in lipid metabolism and secretion. Our work provides further insight into the activity of PhoPR in M. bovis and underlines the importance of the PhoPR system as a master regulator of gene expression in the Mycobacterium tuberculosis complex.

Recent selection and introgression facilitated high-altitude adaptation in cattle.

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle.

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Unveiling the Role of Endoplasmic Reticulum Stress Pathways in Canine Demodicosis.

Canine demodicosis is a prevalent skin disease caused by overpopulation of a commensal species of Demodex mite, yet its precise cause remains unknown. Research suggests that T-cell exhaustion, increased immunosuppressive cytokines, induction of regulatory T cells and increased expression of immune checkpoint inhibitors may contribute to its pathogenesis. This study aimed to gain a deeper understanding of the molecular changes occurring in canine demodicosis using mass spectrometry and pathway enrichment analysis. The results indicate that endoplasmic reticulum stress promotes canine demodicosis through regulation of three linked signalling pathways: eIF2, mTOR, and eIF4 and p70S6K. These pathways are involved in the modulation of Toll-like receptors, most notably TLR2, and have been shown to play a role in the pathogenesis of skin diseases in both dogs and humans. Moreover, these pathways are also implicated in the promotion of immunosuppressive M2 phenotype macrophages. Immunohistochemical analysis, utilising common markers of dendritic cells and macrophages, verified the presence of M2 macrophages in canine demodicosis. The proteomic analysis also identified immunological disease, organismal injury and abnormalities and inflammatory response as the most significant underlying diseases and disorders associated with canine demodicosis. This study demonstrates that Demodex mites, through ER stress, unfolded protein response and M2 macrophages contribute to an immunosuppressive microenvironment, thereby assisting in their proliferation.

Transcriptomics analysis of the bovine endometrium during the perioestrus period.

During the oestrous cycle, the bovine endometrium undergoes morphological and functional changes, which are regulated by alterations in the levels of oestrogen and progesterone and consequent changes in gene expression. To clarify these changes before and after oestrus, RNA-seq was used to profile the transcriptome of oestrus-synchronized beef heifers. Endometrial samples were collected from 29 animals, which were slaughtered in six groups beginning 12 h after the withdrawal of intravaginal progesterone releasing devices until seven days post-oestrus onset (luteal phase). The groups represented proestrus, early oestrus, metoestrus and early dioestrus (luteal phase). Changes in gene expression were estimated relative to gene expression at oestrus. Ingenuity Pathway Analysis (IPA) was used to identify canonical pathways and functional processes of biological importance. A total of 5,845 differentially expressed genes (DEGs) were identified. The lowest number of DEGs was observed at the 12 h post-oestrus time point, whereas the greatest number was observed at Day 7 post-oestrus onset (luteal phase). A total of 2,748 DEGs at this time point did not overlap with any other time points. Prior to oestrus, Neurological disease and Organismal injury and abnormalities appeared among the top IPA diseases and functions categories, with upregulation of genes involved in neurogenesis. Lipid metabolism was upregulated before oestrus and downregulated at 48h post-oestrus, at which point an upregulation of immune-related pathways was observed. In contrast, in the luteal phase the Lipid metabolism and Small molecule biochemistry pathways were upregulated.

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens.

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

Preferential differential gene expression within the WC1.1+ γδ T cell compartment in cattle naturally infected with Mycobacterium bovis.

Bovine tuberculosis (bTB), caused by infection with Mycobacterium bovis, continues to cause significant issues for the global agriculture industry as well as for human health. An incomplete understanding of the host immune response contributes to the challenges of control and eradication of this zoonotic disease. In this study, high-throughput bulk RNA sequencing (RNA-seq) was used to characterise differential gene expression in γδ T cells - a subgroup of T cells that bridge innate and adaptive immunity and have known anti-mycobacterial response mechanisms. γδ T cell subsets are classified based on expression of a pathogen-recognition receptor known as Workshop Cluster 1 (WC1) and we hypothesised that bTB disease may alter the phenotype and function of specific γδ T cell subsets. Peripheral blood was collected from naturally M. bovis-infected (positive for single intradermal comparative tuberculin test (SICTT) and IFN-γ ELISA) and age- and sex-matched, non-infected control Holstein-Friesian cattle. γδ T subsets were isolated using fluorescence activated cell sorting (n = 10-12 per group) and high-quality RNA extracted from each purified lymphocyte subset (WC1.1+, WC1.2+, WC1- and γδ-) was used to generate transcriptomes using bulk RNA-seq (n = 6 per group, representing a total of 48 RNA-seq libraries). Relatively low numbers of differentially expressed genes (DEGs) were observed between most cell subsets; however, 189 genes were significantly differentially expressed in the M. bovis-infected compared to the control groups for the WC1.1+ γδ T cell compartment (absolute log2 FC ≥ 1.5 and FDR P adj. ≤ 0.1). The majority of these DEGs (168) were significantly increased in expression in cells from the bTB+ cattle and included genes encoding transcription factors (TBX21 and EOMES), chemokine receptors (CCR5 and CCR7), granzymes (GZMA, GZMM, and GZMH) and multiple killer cell immunoglobulin-like receptor (KIR) proteins indicating cytotoxic functions. Biological pathway overrepresentation analysis revealed enrichment of genes with multiple immune functions including cell activation, proliferation, chemotaxis, and cytotoxicity of lymphocytes. In conclusion, γδ T cells have important inflammatory and regulatory functions in cattle, and we provide evidence for preferential differential activation of the WC1.1+ specific subset in cattle naturally infected with M. bovis.

Transcriptome profiling of differentiating adipose-derived stem cells across species reveals new genes regulating adipogenesis.

Adipose-derived stem cells (ADSCs) that are enriched in adipose tissue with multilineage differentiation potential have become an important tool in therapeutic research and tissue engineering. Certain breeds of sheep exhibit a unique fat tail trait such that tail tissue accounts for approximately 10 % of body weight and can provide an excellent source of ADSCs. Here, we describe isolation of primary ADSCs from ovine embryonic fat tail tissues that displayed high self-renewal capacity, multilineage differentiation and excellent adipogenic ability. Through transcriptome analysis covering ADSCs differentiating into adipocytes, 37 transcription factors were involved in early transcriptional events that initiate a regulatory cascade of adipogenesis; the entire adipogenic activity consists of a reduction in proliferation ability and upregulation of genes related to lipid generation and energy metabolism, as well as several genes associated with myogenesis. Furthermore, Comparative transcriptome analysis across species (sheep, human, and mouse) revealed enhanced basal metabolic ability in differentiating ovine ADSCs, which may relate to the excellent adipogenic capability of these cells. We also identified a small evolutionarily conserved gene set, consisting of 21 and 22 genes exhibiting increased and decreased expression, respectively. Almost half (20) of these genes have not previously been reported to regulate adipogenesis in mammals. In this study, we identified important regulators that trigger ovine adipocyte differentiation, main biological pathways involved in adipogenesis as well as the evolutionarily conserved genes governing adipogenic process across species. Our study provides a novel excellent biomaterial and novel genes regulating adipogenesis for cellular transplantation therapy and investigations of fat metabolism.

Integrative genomics analysis highlights functionally relevant genes for equine behaviour.

Behavioural plasticity enables horses entering an exercise training programme to adapt with reduced stress. We characterised SNPs associated with behaviour in yearling Thoroughbred horses using genomics analyses for two phenotypes: (1) handler-assessed coping with early training events [coping] (n = 96); and (2) variation in salivary cortisol concentration at the first backing event [cortisol] (n = 34). Using RNA-seq derived gene expression data for amygdala and hippocampus tissues from n = 2 Thoroughbred stallions, we refined the SNPs to those with functional relevance to behaviour by cross-referencing to the 500 most highly expressed genes in each tissue. The SNPs of high significance (q < 0.01) were in proximity to genes (coping - GABARAP, NDM, OAZ1, RPS15A, SPARCL1, VAMP2; cortisol - CEBPA, COA3, DUSP1, HNRNPH1, RACK1) with biological functions in social behaviour, autism spectrum disorder, suicide, stress-induced anxiety and depression, Alzheimer's disease, neurodevelopmental disorders, neuroinflammatory disease, fear-induced behaviours and alcohol and cocaine addiction. The strongest association (q = 0.0002) was with NDN, a gene previously associated with temperament in cattle. This approach highlights functionally relevant genes in the behavioural adaptation of Thoroughbred horses that will contribute to the development of genetic markers to improve racehorse welfare.

Inbreeding depression and durability in the North American Thoroughbred horse.

The proportion of the genome containing runs of homozygosity (ROH) affects production traits in livestock populations. In European and Australasian Thoroughbreds inbreeding, quantified using ROH (FROH ), is associated with the probability of ever racing. Here, we measured FROH using 333 K SNP genotypes from 768 Thoroughbred horses born in North America to evaluate the effect of inbreeding on racing traits in that region. Among North American horses, FROH was not associated (p = 0.518) with the probability of ever racing but was significantly associated with the number of race starts (p = 0.002). Among raced horses, those with a 10% higher FROH than the mean inbreeding coefficient were predicted to have 3.5 fewer race starts compared to horses with a mean inbreeding coefficient. Considering the trend of increasing inbreeding and a decline in the average number of race starts per runner in North America, mitigating inbreeding in the population could positively influence racing durability.

Selection signatures for local and regional adaptation in Chinese Mongolian horse breeds reveal candidate genes for hoof health.

BACKGROUND: Thousands of years of natural and artificial selection since the domestication of the horse has shaped the distinctive genomes of Chinese Mongolian horse populations. Consequently, genomic signatures of selection can provide insights into the human-mediated selection history of specific traits and evolutionary adaptation to diverse environments. Here, we used genome-wide SNPs from five distinct Chinese Mongolian horse populations to identify genomic regions under selection for the population-specific traits, gait, black coat colour, and hoof quality. Other global breeds were used to identify regional-specific signatures of selection. RESULTS: We first identified the most significant selection peak for the Wushen horse in the region on ECA23 harbouring DMRT3, the major gene for gait. We detected selection signatures encompassing several genes in the Baicha Iron Hoof horse that represent good biological candidates for hoof health, including the CSPG4, PEAK1, EXPH5, WWP2 and HAS3 genes. In addition, an analysis of regional subgroups (Asian compared to European) identified a single locus on ECA3 containing the ZFPM1 gene that is a marker of selection for the major domestication event leading to the DOM2 horse clade. CONCLUSIONS: Genomic variation at these loci in the Baicha Iron Hoof may be leveraged in other horse populations to identify animals with superior hoof health or those at risk of hoof-related pathologies. The overlap between the selection signature in Asian horses with the DOM2 selection peak raises questions about the nature of horse domestication events, which may have involved a prehistoric clade other than DOM2 that has not yet been identified.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds.

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

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.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations.

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

Inbreeding depression and the probability of racing in the Thoroughbred horse.

Small effective population sizes and active inbreeding can lead to inbreeding depression due to deleterious recessive mutations exposed in the homozygous state. The Thoroughbred racehorse has low levels of population genetic diversity, but the effects of genomic inbreeding in the population are unknown. Here, we quantified inbreeding based on runs of homozygosity (ROH) using 297 K SNP genotypes from 6128 horses born in Europe and Australia, of which 13.2% were unraced. We show that a 10% increase in inbreeding (FROH) is associated with a 7% lower probability of ever racing. Moreover, a ROH-based genome-wide association study identified a haplotype on ECA14 which, in its homozygous state, is linked to a 32.1% lower predicted probability of ever racing, independent of FROH. The haplotype overlaps a candidate gene, EFNA5, that is highly expressed in cartilage tissue, which when damaged is one of the most common causes of catastrophic musculoskeletal injury in racehorses. Genomics-informed breeding aiming to reduce inbreeding depression and avoid damaging haplotype carrier matings will improve population health and racehorse welfare.

Transcriptome Profiling of Developing Ovine Fat Tail Tissue Reveals an Important Role for MTFP1 in Regulation of Adipogenesis.

Fat-tail sheep exhibit a unique trait whereby substantial adipose tissue accumulates in the tail, a phenotype that is advantageous in many agroecological environments. In this study, we conducted histological assays, transcriptome analysis and functional assays to examine morphogenesis, characterize gene expression, and elucidate mechanisms that regulate fat tail development. We obtained the microstructure of tail before and after fat deposition, and demonstrated that measurable fat deposition occurred by the 80-day embryo (E80) stage, earlier than other tissues. Transcriptome profiling revealed 1,058 differentially expressed genes (DEGs) with six markedly different expression trends. GSEA enrichment and other downstream analyses showed important roles for genes and pathways involving in metabolism and that mitochondrial components were specifically overexpressed in the fat tail tissue of the 70-day embryo (E70). One hundred and eighty-three genes were further identified by leading edge gene analysis, among which, 17 genes have been reported in previous studies, including EEF1D, MTFP1, PPP1CA, PDGFD. Notably, the MTFP1 gene was highly correlated with the expression of other genes and with the highest enrichment score and gene expression change. Knockdown of MTFP1 in isolated adipose derived stem cells (ADSCs) inhibited cell proliferation and migration ability, besides, promoted the process of adipogenesis in vitro.

Timing of Transcriptomic Peripheral Blood Mononuclear Cell Responses of Sheep to Fasciola hepatica Infection Differs From Those of Cattle, Reflecting Different Disease Phenotypes.

Infection with the zoonotic trematode Fasciola hepatica, common in many regions with a temperate climate, leads to delayed growth and loss of productivity in cattle, while infection in sheep can have more severe effects, potentially leading to death. Previous transcriptomic analyses revealed upregulation of TGFB1, cell death and Toll-like receptor signalling, T-cell activation, and inhibition of nitric oxide production in macrophages in response to infection. However, the differences between ovine and bovine responses have not yet been explored. The objective of this study was to further investigate the transcriptomic response of ovine peripheral blood mononuclear cells (PBMC) to F. hepatica infection, and to elucidate the differences between ovine and bovine PBMC responses. Sixteen male Merino sheep were randomly assigned to infected or control groups (n = 8 per group) and orally infected with 120 F. hepatica metacercariae. Transcriptomic data was generated from PBMC at 0, 2 and 16 weeks post-infection (wpi), and analysed for differentially expressed (DE) genes between infected and control animals at each time point (analysis 1), and for each group relative to time 0 (analysis 2). Analysis 2 was then compared to a similar study performed previously on bovine PBMC. A total of 453 DE genes were found at 2 wpi, and 2 DE genes at 16 wpi (FDR < 0.1, analysis 1). Significantly overrepresented biological pathways at 2 wpi included role of PKR in interferon induction and anti-viral response, death receptor signalling and RIG-I-like receptor signalling, which suggested that an activation of innate response to intracellular nucleic acids and inhibition of cellular apoptosis were taking place. Comparison of analysis 2 with the previous bovine transcriptomic study revealed that anti-inflammatory response pathways which were significantly overrepresented in the acute phase in cattle, including IL-10 signalling, Th2 pathway, and Th1 and Th2 activation were upregulated only in the chronic phase in sheep. We propose that the earlier activation of anti-inflammatory responses in cattle, as compared with sheep, may be related to the general absence of acute clinical signs in cattle. These findings offer scope for "smart vaccination" strategies for this important livestock parasite.

Transcriptomic Analysis of Ovine Hepatic Lymph Node Following Fasciola hepatica Infection - Inhibition of NK Cell and IgE-Mediated Signaling.

Fasciola hepatica is a trematode parasite responsible for major economic losses in livestock production, and is also a food-borne zoonotic agent in developing rural regions. For years, the immunoregulatory mechanisms employed by the parasite have hampered efforts to develop a successful vaccine candidate. Given that a comprehensive understanding of the immune response to infection is needed, we investigated the gene expression changes in ovine hepatic lymph nodes after experimental infection with F. hepatica. Lymph nodes from uninfected and infected animals were processed for RNA sequencing (RNA-seq) at 16 weeks post-infection. Comparison of groups revealed 5,132 differentially-expressed genes (DEGs). An inhibition of pro-inflammatory pathways, which has previously been described during fasciolosis, was evident in our data. However, other signals previously identified in ruminant peripheral blood mononuclear cells (PBMC) or liver tissue, such as activation of TGF-ß or apoptosis-related pathways were not detected. We found inhibition of some key immunological pathways, including natural killer (NK) cell activity and IgE-mediated signaling. These may point to additional some as yet unrecognized mechanisms employed by the parasite to evade the host immune response. Understanding these, and leveraging information from this and other omics studies, will be important for the development of future vaccine prototypes against this parasite.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course.

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at -1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the -1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-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.