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.

Individualized Proteogenomics Reveals the Mutational Landscape of Melanoma Patients in Response to Immunotherapy.

Immune checkpoint inhibitors are used to restore or augment antitumor immune responses and show great promise in the treatment of melanoma and other types of cancers. However, only a small percentage of patients are fully responsive to immune checkpoint inhibition, mostly due to tumor heterogeneity and primary resistance to therapy. Both of these features are largely driven by the accumulation of patient-specific mutations, pointing to the need for personalized approaches in diagnostics and immunotherapy. Proteogenomics integrates patient-specific genomic and proteomic data to study cancer development, tumor heterogeneity and resistance mechanisms. Using this approach, we characterized the mutational landscape of four clinical melanoma patients. This enabled the quantification of hundreds of sample-specific amino acid variants, among them many that were previously not reported in melanoma. Changes in abundance at the protein and phosphorylation site levels revealed patient-specific over-represented pathways, notably linked to melanoma development (MAPK1 activation) or immunotherapy (NLRP1 inflammasome). Personalized data integration resulted in the prediction of protein drug targets, such as the drugs vandetanib and bosutinib, which were experimentally validated and led to a reduction in the viability of tumor cells. Our study emphasizes the potential of proteogenomic approaches to study personalized mutational landscapes, signaling networks and therapy options.

Proteogenomics Reveals Perturbed Signaling Networks in Malignant Melanoma Cells Resistant to BRAF Inhibition.

Analysis of nucleotide variants is a cornerstone of cancer medicine. Although only 2% of the genomic sequence is protein coding, mutations occurring in these regions have the potential to influence protein structure or modification status and may have severe impact on disease aetiology. Proteogenomics enables the analysis of sample-specific nonsynonymous nucleotide variants with regard to their effect at the proteome and phosphoproteome levels. Here, we developed a proof-of-concept proteogenomics workflow and applied it to the malignant melanoma cell line A375. Initially, we studied the resistance to serine/threonine-protein kinase B-raf (BRAF) inhibitor (BRAFi) vemurafenib in A375 cells. This allowed identification of several oncogenic nonsynonymous nucleotide variants, including a gain-of-function variant on aurora kinase A (AURKA) at F31I. We also detected significant changes in abundance among (phospho)proteins, which led to reactivation of the MAPK signaling pathway in BRAFi-resistant A375 cells. Upon reconstruction of the multiomic integrated signaling networks, we predicted drug therapies with the potential to disrupt BRAFi resistance mechanism in A375 cells. Notably, we showed that AURKA inhibition is effective and specific against BRAFi-resistant A375 cells. Subsequently, we investigated amino acid variants that interfere with protein posttranslational modification (PTM) status and potentially influence A375 cell signaling irrespective of BRAFi resistance. Mass spectrometry (MS) measurements confirmed variant-driven PTM changes in 12 proteins. Among them was the runt-related transcription factor 1 (RUNX1) displaying a variant on a known phosphorylation site S(Ph)276L. We confirmed the loss of phosphorylation site by MS and demonstrated the impact of this variant on RUNX1 interactome.

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.

Quantitative Proteomics Links the Intermediate Filament Nestin to Resistance to Targeted BRAF Inhibition in Melanoma Cells.

Targeted inhibition of mutated kinases using selective MAP kinase inhibitors in malignant melanoma often results in temporary improvement of clinical symptoms followed by rapid development of resistance. To gain insights in molecular processes that govern resistance, we performed SILAC-based quantitative proteomics profiling of vemurafenib-resistant and -sensitive melanoma cells. Among downregulated proteins in vemurafenib-resistant cell lines we detected multiple proteins involved in cytoskeletal organization and signaling, including the intermediate filament nestin, which was one of the most downregulated proteins. Previous studies showed that nestin is expressed in various types of solid tumors and its abundance correlates with malignant phenotype of transformed cells. However, the role of nestin in cancer cells regarding acquired resistance is still poorly understood. We performed CRISPR/Cas9 knockout of the nestin gene (NES) in vemurafenib-sensitive cells and showed that loss of nestin leads to increased cellular proliferation and colony formation upon treatment with BRAFV600E and MEK inhibitors. Moreover, nestin depletion led to increased invasiveness and metalloproteinase activity like the phenotype of melanoma cells with acquired resistance to the BRAF inhibitor. Finally, phosphoproteome analysis revealed that nestin depletion influenced signaling through integrin and PI3K/AKT/mTOR pathways and led to increased focal adhesion kinase abundance and phosphorylation. Taken together, our results reveal that nestin is associated with acquired vemurafenib resistance in melanoma cells.

In-depth analysis of Bacillus subtilis proteome identifies new ORFs and traces the evolutionary history of modified proteins.

Bacillus subtilis is a sporulating Gram-positive bacterium widely used in basic research and biotechnology. Despite being one of the best-characterized bacterial model organism, recent proteomics studies identified only about 50% of its theoretical protein count. Here we combined several hundred MS measurements to obtain a comprehensive map of the proteome, phosphoproteome and acetylome of B. subtilis grown at 37 °C in minimal medium. We covered 75% of the theoretical proteome (3,159 proteins), detected 1,085 phosphorylation and 4,893 lysine acetylation sites and performed a systematic bioinformatic characterization of the obtained data. A subset of analyzed MS files allowed us to reconstruct a network of Hanks-type protein kinases, Ser/Thr/Tyr phosphatases and their substrates. We applied genomic phylostratigraphy to gauge the evolutionary age of B. subtilis protein classes and revealed that protein modifications were present on the oldest bacterial proteins. Finally, we performed a proteogenomic analysis by mapping all MS spectra onto a six-frame translation of B. subtilis genome and found evidence for 19 novel ORFs. We provide the most extensive overview of the proteome and post-translational modifications for B. subtilis to date, with insights into functional annotation and evolutionary aspects of the B. subtilis genome.

RNA Sequencing (RNA-Seq) Reveals Extremely Low Levels of Reticulocyte-Derived Globin Gene Transcripts in Peripheral Blood From Horses (Equus caballus) and Cattle (Bos taurus).

RNA-seq has emerged as an important technology for measuring gene expression in peripheral blood samples collected from humans and other vertebrate species. In particular, transcriptomics analyses of whole blood can be used to study immunobiology and develop novel biomarkers of infectious disease. However, an obstacle to these methods in many mammalian species is the presence of reticulocyte-derived globin mRNAs in large quantities, which can complicate RNA-seq library sequencing and impede detection of other mRNA transcripts. A range of supplementary procedures for targeted depletion of globin transcripts have, therefore, been developed to alleviate this problem. Here, we use comparative analyses of RNA-seq data sets generated from human, porcine, equine, and bovine peripheral blood to systematically assess the impact of globin mRNA on routine transcriptome profiling of whole blood in cattle and horses. The results of these analyses demonstrate that total RNA isolated from equine and bovine peripheral blood contains very low levels of globin mRNA transcripts, thereby negating the need for globin depletion and greatly simplifying blood-based transcriptomic studies in these two domestic species.

Proteome Response of a Metabolically Flexible Anoxygenic Phototroph to Fe(II) Oxidation.

The oxidation of Fe(II) by anoxygenic photosynthetic bacteria was likely a key contributor to Earth's biosphere prior to the evolution of oxygenic photosynthesis and is still found in a diverse range of modern environments. All known phototrophic Fe(II) oxidizers can utilize a wide range of substrates, thus making them very metabolically flexible. However, the underlying adaptations required to oxidize Fe(II), a potential stressor, are not completely understood. We used a combination of quantitative proteomics and cryogenic transmission electron microscopy (cryo-TEM) to compare cells of Rhodopseudomonas palustris TIE-1 grown photoautotrophically with Fe(II) or H2 and photoheterotrophically with acetate. We observed unique proteome profiles for each condition, with differences primarily driven by carbon source. However, these differences were not related to carbon fixation but to growth and light harvesting processes, such as pigment synthesis. Cryo-TEM showed stunted development of photosynthetic membranes in photoautotrophic cultures. Growth on Fe(II) was characterized by a response typical of iron homeostasis, which included an increased abundance of proteins required for metal efflux (particularly copper) and decreased abundance of iron import proteins, including siderophore receptors, with no evidence of further stressors, such as oxidative damage. This study suggests that the main challenge facing anoxygenic phototrophic Fe(II) oxidizers comes from growth limitations imposed by autotrophy, and, once this challenge is overcome, iron stress can be mitigated using iron management mechanisms common to diverse bacteria (e.g., by control of iron influx and efflux).IMPORTANCE The cycling of iron between redox states leads to the precipitation and dissolution of minerals, which can in turn impact other major biogeochemical cycles, such as those of carbon, nitrogen, phosphorus and sulfur. Anoxygenic phototrophs are one of the few drivers of Fe(II) oxidation in anoxic environments and are thought to contribute significantly to iron cycling in both modern and ancient environments. These organisms thrive at high Fe(II) concentrations, yet the adaptations required to tolerate the stresses associated with this are unclear. Despite the general consensus that high Fe(II) concentrations pose numerous stresses on these organisms, our study of the large-scale proteome response of a model anoxygenic phototroph to Fe(II) oxidation demonstrates that common iron homeostasis strategies are adequate to manage this. The bulk of the proteome response is not driven by adaptations to Fe(II) stress but to adaptations required to utilize an inorganic carbon source. Such a global overview of the adaptation of these organisms to Fe(II) oxidation provides valuable insights into the physiology of these biogeochemically important organisms and suggests that Fe(II) oxidation may not pose as many challenges to anoxygenic phototrophs as previously thought.

Parallel reaction monitoring on a Q Exactive mass spectrometer increases reproducibility of phosphopeptide detection in bacterial phosphoproteomics measurements.

Increasing number of studies report the relevance of protein Ser/Thr/Tyr phosphorylation in bacterial physiology, yet the analysis of this type of modification in bacteria still presents a considerable challenge. Unlike in eukaryotes, where tens of thousands of phosphorylation events likely occupy more than two thirds of the proteome, the abundance of protein phosphorylation is much lower in bacteria. Even the state-of-the-art phosphopeptide enrichment protocols fail to remove the high background of abundant unmodified peptides, leading to low signal intensity and undersampling of phosphopeptide precursor ions in consecutive data-dependent MS runs. Consequently, large-scale bacterial phosphoproteomic datasets often suffer from poor reproducibility and a high number of missing values. Here we explore the application of parallel reaction monitoring (PRM) on a Q Exactive mass spectrometer in bacterial phosphoproteome analysis, focusing especially on run-to-run sampling reproducibility. In multiple measurements of identical phosphopeptide-enriched samples, we show that PRM outperforms data-dependent acquisition (DDA) in terms of detection frequency, reaching almost complete sampling efficiency, compared to 20% in DDA. We observe a similar trend over multiple heterogeneous phosphopeptide-enriched samples and conclude that PRM shows a great promise in bacterial phosphoproteomics analyses where reproducible detection and quantification of a relatively small set of phosphopeptides is desired. SIGNIFICANCE: Bacterial phosphorylated peptides occur in low abundance compared to their unmodified counterparts, and are therefore rarely reproducibly detected in shotgun (DDA) proteomics measurements. Here we show that parallel reaction monitoring complements DDA analyses and makes detection of known, targeted phosphopeptides more reproducible. This will be of significance in replicated MS measurements that have a goal to reproducibly detect and quantify phosphopeptides of interest.

Comparative 'omics analyses differentiate Mycobacterium tuberculosis and Mycobacterium bovis and reveal distinct macrophage responses to infection with the human and bovine tubercle bacilli.

Members of the Mycobacterium tuberculosis complex (MTBC) are the causative agents of tuberculosis in a range of mammals, including humans. A key feature of MTBC pathogens is their high degree of genetic identity yet distinct host tropism. Notably, while Mycobacterium bovis is highly virulent and pathogenic for cattle, the human pathogen M. tuberculosis is attenuated in cattle. Previous research also suggests that host preference amongst MTBC members has a basis in host innate immune responses. To explore MTBC host tropism, we present in-depth profiling of the MTBC reference strains M. bovis AF2122/97 and M. tuberculosis H37Rv at both the global transcriptional and the translational level via RNA-sequencing and SWATH MS. Furthermore, a bovine alveolar macrophage infection time course model was used to investigate the shared and divergent host transcriptomic response to infection with M. tuberculosis H37Rv or M. bovis AF2122/97. Significant differential expression of virulence-associated pathways between the two bacilli was revealed, including the ESX-1 secretion system. A divergent transcriptional response was observed between M. tuberculosis H37Rv and M. bovis AF2122/97 infection of bovine alveolar macrophages, in particular cytosolic DNA-sensing pathways at 48 h post-infection, and highlights a distinct engagement of M. bovis with the bovine innate immune system. The work presented here therefore provides a basis for the identification of host innate immune mechanisms subverted by virulent host-adapted mycobacteria to promote their survival during the early stages of infection.

Circulating microRNAs as Potential Biomarkers of Infectious Disease.

microRNAs (miRNAs) are a class of small non-coding endogenous RNA molecules that regulate a wide range of biological processes by post-transcriptionally regulating gene expression. Thousands of these molecules have been discovered to date, and multiple miRNAs have been shown to coordinately fine-tune cellular processes key to organismal development, homeostasis, neurobiology, immunobiology, and control of infection. The fundamental regulatory role of miRNAs in a variety of biological processes suggests that differential expression of these transcripts may be exploited as a novel source of molecular biomarkers for many different disease pathologies or abnormalities. This has been emphasized by the recent discovery of remarkably stable miRNAs in mammalian biofluids, which may originate from intracellular processes elsewhere in the body. The potential of circulating miRNAs as biomarkers of disease has mainly been demonstrated for various types of cancer. More recently, however, attention has focused on the use of circulating miRNAs as diagnostic/prognostic biomarkers of infectious disease; for example, human tuberculosis caused by infection with Mycobacterium tuberculosis, sepsis caused by multiple infectious agents, and viral hepatitis. Here, we review these developments and discuss prospects and challenges for translating circulating miRNA into novel diagnostics for infectious disease.

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.

RNA sequencing provides exquisite insight into the manipulation of the alveolar macrophage by tubercle bacilli.

Mycobacterium bovis, the agent of bovine tuberculosis, causes an estimated $3 billion annual losses to global agriculture due, in part, to the limitations of current diagnostics. Development of next-generation diagnostics requires a greater understanding of the interaction between the pathogen and the bovine host. Therefore, to explore the early response of the alveolar macrophage to infection, we report the first application of RNA-sequencing to define, in exquisite detail, the transcriptomes of M. bovis-infected and non-infected alveolar macrophages from ten calves at 2, 6, 24 and 48 hours post-infection. Differentially expressed sense genes were detected at these time points that revealed enrichment of innate immune signalling functions, and transcriptional suppression of host defence mechanisms (e.g., lysosome maturation). We also detected differentially expressed natural antisense transcripts, which may play a role in subverting innate immune mechanisms following infection. Furthermore, we report differential expression of novel bovine genes, some of which have immune-related functions based on orthology with human proteins. This is the first in-depth transcriptomics investigation of the alveolar macrophage response to the early stages of M. bovis infection and reveals complex patterns of gene expression and regulation that underlie the immunomodulatory mechanisms used by M. bovis to evade host defence mechanisms.

Analysis of the Bovine Monocyte-Derived Macrophage Response to Mycobacterium avium Subspecies Paratuberculosis Infection Using RNA-seq.

Johne's disease, caused by infection with Mycobacterium avium subsp. paratuberculosis, (MAP), is a chronic intestinal disease of ruminants with serious economic consequences for cattle production in the United States and elsewhere. During infection, MAP bacilli are phagocytosed and subvert host macrophage processes, resulting in subclinical infections that can lead to immunopathology and dissemination of disease. Analysis of the host macrophage transcriptome during infection can therefore shed light on the molecular mechanisms and host-pathogen interplay associated with Johne's disease. Here, we describe results of an in vitro study of the bovine monocyte-derived macrophage (MDM) transcriptome response during MAP infection using RNA-seq. MDM were obtained from seven age- and sex-matched Holstein-Friesian cattle and were infected with MAP across a 6-h infection time course with non-infected controls. We observed 245 and 574 differentially expressed (DE) genes in MAP-infected versus non-infected control samples (adjusted P value ≤0.05) at 2 and 6 h post-infection, respectively. Functional analyses of these DE genes, including biological pathway enrichment, highlighted potential functional roles for genes that have not been previously described in the host response to infection with MAP bacilli. In addition, differential expression of pro- and anti-inflammatory cytokine genes, such as those associated with the IL-10 signaling pathway, and other immune-related genes that encode proteins involved in the bovine macrophage response to MAP infection emphasize the balance between protective host immunity and bacilli survival and proliferation. Systematic comparisons of RNA-seq gene expression results with Affymetrix(®) microarray data generated from the same experimental samples also demonstrated that RNA-seq represents a superior technology for studying host transcriptional responses to intracellular infection.

MicroRNA profiling of the bovine alveolar macrophage response to Mycobacterium bovis infection suggests pathogen survival is enhanced by microRNA regulation of endocytosis and lysosome trafficking.

Mycobacterium bovis, the causative agent of bovine tuberculosis, a major problem for global agriculture, spreads via an airborne route and is taken up by alveolar macrophages (AM) in the lung. Here, we describe the first next-generation sequencing (RNA-seq) approach to temporally profile miRNA expression in primary bovine AMs post-infection with M. bovis. One, six, and forty miRNAs were identified as significantly differentially expressed at 2, 24 and 48 h post-infection, respectively. The differential expression of three miRNAs (bta-miR-142-5p, bta-miR-146a, and bta-miR-423-3p) was confirmed by RT-qPCR. Pathway analysis of the predicted mRNA targets of differentially expressed miRNAs suggests that these miRNAs preferentially target several pathways that are functionally relevant for mycobacterial pathogenesis, including endocytosis and lysosome trafficking, IL-1 signalling and the TGF-ß pathway. Over-expression studies using a bovine macrophage cell-line (Bomac) reveal the targeting of two key genes in the innate immune response to M. bovis, IL-1 receptor-associated kinase 1 (IRAK1) and TGF-ß receptor 2 (TGFBR2), by miR-146. Taken together, our study suggests that miRNAs play a key role in tuning the complex interplay between M. bovis survival strategies and the host immune response.

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.

Key Hub and Bottleneck Genes Differentiate the Macrophage Response to Virulent and Attenuated Mycobacterium bovis.

Mycobacterium bovis is an intracellular pathogen that causes tuberculosis in cattle. Following infection, the pathogen resides and persists inside host macrophages by subverting host immune responses via a diverse range of mechanisms. Here, a high-density bovine microarray platform was used to examine the bovine monocyte-derived macrophage transcriptome response to M. bovis infection relative to infection with the attenuated vaccine strain, M. bovis Bacille Calmette-Guérin. Differentially expressed genes were identified (adjusted P-value ≤0.01) and interaction networks generated across an infection time course of 2, 6, and 24 h. The largest number of biological interactions was observed in the 24-h network, which exhibited scale-free network properties. The 24-h network featured a small number of key hub and bottleneck gene nodes, including IKBKE, MYC, NFKB1, and EGR1 that differentiated the macrophage response to virulent and attenuated M. bovis strains, possibly via the modulation of host cell death mechanisms. These hub and bottleneck genes represent possible targets for immuno-modulation of host macrophages by virulent mycobacterial species that enable their survival within a hostile environment.

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.

Innate cytokine profiling of bovine alveolar macrophages reveals commonalities and divergence in the response to Mycobacterium bovis and Mycobacterium tuberculosis infection.

Despite sharing >99.9% genome sequence similarity at the nucleotide level, Mycobacterium tuberculosis and Mycobacterium bovis-the causative agents of human and bovine tuberculosis, respectively-exhibit distinct host preferences. M. bovis can cause disease in both cattle and humans yet rarely transmits between immuno-competent human hosts, while M. tuberculosis is a highly successful pathogen of humans that does not sustain in animal populations. Based on the key role played by alveolar macrophages during mycobacterial infection, we hypothesised that the immunological and pathological differences observed in cattle infected with virulent M. bovis and M. tuberculosis may have a basis in innate immune mechanisms; these differences, in turn, would be reflected at the macrophage mRNA and protein level. To investigate this, we have analysed the transcriptional profile of innate immune genes in bovine alveolar macrophages following 24 and 48 h infection with the genome-sequenced strains, M. bovis AF2122/97 and M. tuberculosis H37Rv. A bespoke multiplex ELISA was also used to quantify corresponding cytokine secretion in supernatants from the same infected alveolar macrophages. All cytokines showed similar significant patterns of expression (i.e., up- or down-regulation) at both the mRNA and protein levels in infected macrophages relative to parallel non-infected controls at the two time points (P ≤ 0.05). However, significant upregulation and downregulation of several innate immune genes-including TLR2, FOS, PIK3IP1, CCL4, IL1B, IL6 and TNF-and the CCL-4 protein was observed in the M. bovis-infected macrophages relative to the M. tuberculosis-infected macrophages 48 h post-infection (P ≤ 0.05). These results support the hypothesis that the divergent virulence of M. bovis and M. tuberculosis in cattle has a basis in innate immune mechanisms, which may contribute to host preference within the M. tuberculosis complex of strains.

Profiling microRNA expression in bovine alveolar macrophages using RNA-seq.

MicroRNAs (miRNAs) are important regulators of gene expression and are known to play a key role in regulating both adaptive and innate immunity. Bovine alveolar macrophages (BAMs) help maintain lung homeostasis and constitute the front line of host defense against several infectious respiratory diseases, such as bovine tuberculosis. Little is known, however, about the role miRNAs play in these cells. In this study, we used a high-throughput sequencing approach, RNA-seq, to determine the expression levels of known and novel miRNAs in unchallenged BAMs isolated from lung lavages of eight different healthy Holstein-Friesian male calves. Approximately 80 million sequence reads were generated from eight BAM miRNA Illumina sequencing libraries, and 80 miRNAs were identified as being expressed in BAMs at a threshold of at least 100 reads per million (RPM). The expression levels of miRNAs varied over a large dynamic range, with a few miRNAs expressed at very high levels (up to 800,000RPM), and the majority lowly expressed. Notably, many of the most highly expressed miRNAs in BAMs have known roles in regulating immunity in other species (e.g. bta-let-7i, bta-miR-21, bta-miR-27, bta-miR-99b, bta-miR-146, bta-miR-147, bta-miR-155 and bta-miR-223). The most highly expressed miRNA in BAMs was miR-21, which has been shown to regulate the expression of antimicrobial peptides in Mycobacterium leprae-infected human monocytes. Furthermore, the predicted target genes of BAM-expressed miRNAs were found to be statistically enriched for roles in innate immunity. In addition to profiling the expression of known miRNAs, the RNA-seq data was also analysed to identify potentially novel bovine miRNAs. One putatively novel bovine miRNA was identified. To the best of our knowledge, this is the first RNA-seq study to profile miRNA expression in BAMs and provides an important reference dataset for investigating the regulatory roles miRNAs play in this important immune cell type.