Comparative functional genomics and the bovine macrophage response to strains of the Mycobacterium genus

Mycobacterial infections are major causes of morbidity and mortality in cattle and are also potential zoonotic agents with implications for human health. Despite the implementation of comprehensive animal surveillance programs, many mycobacterial diseases have remained recalcitrant to eradication in several industrialized countries. Two major mycobacterial pathogens of cattle are Mycobacterium bovis and Mycobacterium avium subspecies paratuberculosis (MAP), the causative agents of bovine tuberculosis (BTB) and Johne's disease (JD), respectively. BTB is a chronic, granulomatous disease of the respiratory tract that is spread via aerosol transmission, while JD is a chronic granulomatous disease of the intestines that is transmitted via the fecal-oral route. Although these diseases exhibit differential tissue tropism and distinct complex etiologies, both M. bovis and MAP infect, reside, and replicate in host macrophages - the key host innate immune cell that encounters mycobacterial pathogens after initial exposure and mediates the subsequent immune response. The persistence of M. bovis and MAP in macrophages relies on a diverse series of immunomodulatory mechanisms, including the inhibition of phagosome maturation and apoptosis, generation of cytokine-induced necrosis enabling dissemination of infection through the host, local pathology, and ultimately shedding of the pathogen. Here, we review the bovine macrophage response to infection with M. bovis and MAP. In particular, we describe how recent advances in functional genomics are shedding light on the host macrophage-pathogen interactions that underlie different mycobacterial diseases. To illustrate this, we present new analyses of previously published bovine macrophage transcriptomics data following in vitro infection with virulent M. bovis, the attenuated vaccine strain M. bovis BCG, and MAP, and discuss our findings with respect to the differing etiologies of BTB and JD

GOexpress: an R/Bioconductor package for the identification and visualisation of robust gene ontology signatures through supervised learning of gene expression data

Background: Identification of gene expression profiles that differentiate experimental groups is critical for discovery and analysis of key molecular pathways and also for selection of robust diagnostic or prognostic biomarkers. While integration of differential expression statistics has been used to refine gene set enrichment analyses, such approaches are typically limited to single gene lists resulting from simple two-group comparisons or time-series analyses. In contrast, functional class scoring and machine learning approaches provide powerful alternative methods to leverage molecular measurements for pathway analyses, and to compare continuous and multi-level categorical factors. Results: We introduce GOexpress, a software package for scoring and summarising the capacity of gene ontology features to simultaneously classify samples from multiple experimental groups. GOexpress integrates normalised gene expression data (e.g., from microarray and RNA-seq experiments) and phenotypic information of individual samples with gene ontology annotations to derive a ranking of genes and gene ontology terms using a supervised learning approach. The default random forest algorithm allows interactions between all experimental factors, and competitive scoring of expressed genes to evaluate their relative importance in classifying predefined groups of samples. Conclusions: GOexpress enables rapid identification and visualisation of ontology-related gene panels that robustly classify groups of samples and supports both categorical (e.g., infection status, treatment) and continuous (e.g., time-series, drug concentrations) experimental factors. The use of standard Bioconductor extension packages and publicly available gene ontology annotations facilitates straightforward integration of GOexpress within existing computational biology pipelines.Department of Agriculture, Food and the MarineEuropean Commission - Seventh Framework Programme (FP7)Science Foundation IrelandUniversity College Dubli

RNA-seq transcriptional profiling of peripheral blood leukocytes from cattle infected with Mycobacterium bovis

Bovine tuberculosis, caused by infection with 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 gene expression microarrays and RNA-sequencing (RNA-seq), 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 analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix(®) GeneChip(®) Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity(®) Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top-ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA

Meta-Analysis of Mitochondrial DNA Reveals Several Population Bottlenecks during Worldwide Migrations of Cattle

Several studies have investigated the differentiation of mitochondrial DNA in Eurasian, African and American cattle as well as archaeological bovine material. A global survey of these studies shows that haplogroup distributions are more stable in time than in space. All major migrations of cattle have shifted the haplogroup distributions considerably with a reduction of the number of haplogroups and/or an expansion of haplotypes that are rare or absent in the ancestral populations. The most extreme case is the almost exclusive colonization of Africa by the T1 haplogroup, which is rare in Southwest Asian cattle. In contrast, ancient samples invariably show continuity with present-day cattle from the same location. These findings indicate strong maternal founder effects followed by limited maternal gene flow when new territories are colonized. However, effects of adaptation to new environments may also play a rol

Interrogation of modern and ancient genomes reveals the complex domestic history of cattle

The analysis of mitochondrial and nuclear DNA sequence polymorphisms from modern cattle populations has had a profound impact on our understanding of the events surrounding the domestication of cattle. From these studies, it has been possible to distinguish between pre- and post-domestic genetic differentiation, supporting previous assertions from archaeological studies and, in some cases, revealing novel aspects of the demographic history of cattle. Analyses of genetic material retrieved from the remains of extinct ancestral wild cattle have also added valuable layers of information pertaining to cattle domestic origins; however, information from these investigations have, in general, been limited to small, variable portions of the mitochondrial genome owing to technical challenges associated with the retrieval and amplification of ancient DNA. In recent years, however, new high-throughput, massively parallel genomics technology platforms, such as single-nucleotide polymorphism (SNP) genotyping arrays and next-generation sequencing (NGS), have provided a new impetus to the studies of genetic variation in extant and ancient cattle. Arrays of SNP have facilitated high-resolution genetic surveys of global cattle populations and detection of ancient and recent genomic selective sweeps. Next-generation sequencing analyses of modern and ancient cattle hold great promise for identifying and cataloging of pre- and post-domestication patterns of genomic variation and correlating this with natural and artificial selection processes

Transcriptional adaptations following exercise in Thoroughbred horse skeletal muscle highlights molecular mechanisms that lead to muscle hypertrophy

<p>Abstract</p> <p>Background</p> <p>Selection for exercise-adapted phenotypes in the Thoroughbred racehorse has provided a valuable model system to understand molecular responses to exercise in skeletal muscle. Exercise stimulates immediate early molecular responses as well as delayed responses during recovery, resulting in a return to homeostasis and enabling long term adaptation. Global mRNA expression during the immediate-response period has not previously been reported in skeletal muscle following exercise in any species. Also, global gene expression changes in equine skeletal muscle following exercise have not been reported. Therefore, to identify novel genes and key regulatory pathways responsible for exercise adaptation we have used equine-specific cDNA microarrays to examine global mRNA expression in skeletal muscle from a cohort of Thoroughbred horses (<it>n = </it>8) at three time points (before exercise, immediately post-exercise, and four hours post-exercise) following a single bout of treadmill exercise.</p> <p>Results</p> <p>Skeletal muscle biopsies were taken from the <it>gluteus medius </it>before (T₀), immediately after (T₁) and four hours after (T₂) exercise. Statistically significant differences in mRNA abundance between time points (T₀<it>vs </it>T₁and T₀<it>vs </it>T₂) were determined using the empirical Bayes moderated <it>t</it>-test in the Bioconductor package Linear Models for Microarray Data (LIMMA) and the expression of a select panel of genes was validated using real time quantitative reverse transcription PCR (qRT-PCR). While only two genes had increased expression at T₁(<it>P </it>< 0.05), by T₂932 genes had increased (<it>P </it>< 0.05) and 562 genes had decreased expression (<it>P </it>< 0.05). Functional analysis of genes differentially expressed during the recovery phase (T₂) revealed an over-representation of genes localized to the actin cytoskeleton and with functions in the MAPK signalling, focal adhesion, insulin signalling, mTOR signaling, p53 signaling and Type II diabetes mellitus pathways. At T₁, using a less stringent statistical approach, we observed an over-representation of genes involved in the stress response, metabolism and intracellular signaling. These findings suggest that protein synthesis, mechanosensation and muscle remodeling contribute to skeletal muscle adaptation towards improved integrity and hypertrophy.</p> <p>Conclusions</p> <p>This is the first study to characterize global mRNA expression profiles in equine skeletal muscle using an equine-specific microarray platform. Here we reveal novel genes and mechanisms that are temporally expressed following exercise providing new knowledge about the early and late molecular responses to exercise in the equine skeletal muscle transcriptome.</p

Transcriptional profiling of cattle infected with Trypanosoma congolense highlights gene expression signatures underlying trypanotolerance and trypanosusceptibility

<p>Abstract</p> <p>Background</p> <p>African animal trypanosomiasis (AAT) caused by tsetse fly-transmitted protozoa of the genus <it>Trypanosoma </it>is a major constraint on livestock and agricultural production in Africa and is among the top ten global cattle diseases impacting on the poor. Here we show that a functional genomics approach can be used to identify temporal changes in host peripheral blood mononuclear cell (PBMC) gene expression due to disease progression. We also show that major gene expression differences exist between cattle from trypanotolerant and trypanosusceptible breeds. Using bovine long oligonucleotide microarrays and real time quantitative reverse transcription PCR (qRT-PCR) validation we analysed PBMC gene expression in naïve trypanotolerant and trypanosusceptible cattle experimentally challenged with <it>Trypanosoma congolense </it>across a 34-day infection time course.</p> <p>Results</p> <p>Trypanotolerant N'Dama cattle displayed a rapid and distinct transcriptional response to infection, with a ten-fold higher number of genes differentially expressed at day 14 post-infection compared to trypanosusceptible Boran cattle. These analyses identified coordinated temporal gene expression changes for both breeds in response to trypanosome infection. In addition, a panel of genes were identified that showed pronounced differences in gene expression between the two breeds, which may underlie the phenomena of trypanotolerance and trypanosusceptibility. Gene ontology (GO) analysis demonstrate that the products of these genes may contribute to increased mitochondrial mRNA translational efficiency, a more pronounced B cell response, an elevated activation status and a heightened response to stress in trypanotolerant cattle.</p> <p>Conclusion</p> <p>This study has revealed an extensive and diverse range of cellular processes that are altered temporally in response to trypanosome infection in African cattle. Results indicate that the trypanotolerant N'Dama cattle respond more rapidly and with a greater magnitude to infection compared to the trypanosusceptible Boran cattle. Specifically, a subset of the genes analyzed by real time qRT-PCR, which display significant breed differences, could collectively contribute to the trypanotolerance trait in N'Dama.</p

Tracing Cattle Breeds with Principal Components Analysis Ancestry Informative SNPs

The recent release of the Bovine HapMap dataset represents the most detailed survey of bovine genetic diversity to date, providing an important resource for the design and development of livestock production. We studied this dataset, comprising more than 30,000 Single Nucleotide Polymorphisms (SNPs) for 19 breeds (13 taurine, three zebu, and three hybrid breeds), seeking to identify small panels of genetic markers that can be used to trace the breed of unknown cattle samples. Taking advantage of the power of Principal Components Analysis and algorithms that we have recently described for the selection of Ancestry Informative Markers from genomewide datasets, we present a decision-tree which can be used to accurately infer the origin of individual cattle. In doing so, we present a thorough examination of population genetic structure in modern bovine breeds. Performing extensive cross-validation experiments, we demonstrate that 250-500 carefully selected SNPs suffice in order to achieve close to 100% prediction accuracy of individual ancestry, when this particular set of 19 breeds is considered. Our methods, coupled with the dense genotypic data that is becoming increasingly available, have the potential to become a valuable tool and have considerable impact in worldwide livestock production. They can be used to inform the design of studies of the genetic basis of economically important traits in cattle, as well as breeding programs and efforts to conserve biodiversity. Furthermore, the SNPs that we have identified can provide a reliable solution for the traceability of breed-specific branded products