Comprehensive insights into transcriptional adaptation of intracellular mycobacteria by microbe-enriched dual RNA sequencing

BackgroundThe human pathogen Mycobacterium tuberculosis has the capacity to escape eradication by professional phagocytes. During infection, M. tuberculosis resists the harsh environment of phagosomes and actively manipulates macrophages and dendritic cells to ensure prolonged intracellular survival. In contrast to other intracellular pathogens, it has remained difficult to capture the transcriptome of mycobacteria during infection due to an unfavorable host-to-pathogen ratio.ResultsWe infected the human macrophage-like cell line THP-1 with the attenuated M. tuberculosis surrogate M. bovis Bacillus Calmette–Guérin (M. bovis BCG). Mycobacterial RNA was up to 1000-fold underrepresented in total RNA preparations of infected host cells. We employed microbial enrichment combined with specific ribosomal RNA depletion to simultaneously analyze the transcriptional responses of host and pathogen during infection by dual RNA sequencing. Our results confirm that mycobacterial pathways for cholesterol degradation and iron acquisition are upregulated during infection. In addition, genes involved in the methylcitrate cycle, aspartate metabolism and recycling of mycolic acids were induced. In response to M. bovis BCG infection, host cells upregulated de novo cholesterol biosynthesis presumably to compensate for the loss of this metabolite by bacterial catabolism.ConclusionsDual RNA sequencing allows simultaneous capture of the global transcriptome of host and pathogen, during infection. However, mycobacteria remained problematic due to their relatively low number per host cell resulting in an unfavorable bacterium-to-host RNA ratio. Here, we use a strategy that combines enrichment for bacterial transcripts and dual RNA sequencing to provide the most comprehensive transcriptome of intracellular mycobacteria to date. The knowledge acquired into the pathogen and host pathways regulated during infection may contribute to a solid basis for the deployment of novel intervention strategies to tackle infection.

[1]  H. Mollenkopf,et al.  Isolation of RNA from mycobacteria grown under in vitro and in vivo conditions. , 2000, FEMS microbiology letters.

[2]  Rob Knight,et al.  Direct sequencing of the human microbiome readily reveals community differences , 2010, Genome Biology.

[3]  S. R. Kushner,et al.  Bacterial/archaeal/organellar polyadenylation , 2011, Wiley interdisciplinary reviews. RNA.

[4]  Marlene Belfort,et al.  Multiple small RNAs identified in Mycobacterium bovis BCG are also expressed in Mycobacterium tuberculosis and Mycobacterium smegmatis , 2010, Nucleic acids research.

[5]  Dandan Huang,et al.  BSRD: a repository for bacterial small regulatory RNA , 2012, Nucleic Acids Res..

[6]  Brunero Liseo,et al.  Profiling of Mycobacterium tuberculosis gene expression during human macrophage infection: upregulation of the alternative sigma factor G, a group of transcriptional regulators, and proteins with unknown function. , 2006, Research in microbiology.

[7]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[8]  R. Brosch,et al.  Disruption of the ESX‐5 system of Mycobacterium tuberculosis causes loss of PPE protein secretion, reduction of cell wall integrity and strong attenuation , 2012, Molecular microbiology.

[9]  M. Voskuil Mycobacterium tuberculosis Cholesterol Catabolism Requires a New Class of Acyl Coenzyme A Dehydrogenase , 2013, Journal of bacteriology.

[10]  P. De Baetselier,et al.  The Mycobacterium tuberculosis cell wall component mycolic acid elicits pathogen‐associated host innate immune responses , 2005, European journal of immunology.

[11]  Julian Parkhill,et al.  The complete genome sequence of Mycobacterium bovis , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Katherine H. Huang,et al.  Efficient and robust RNA-seq process for cultured bacteria and complex community transcriptomes , 2012, Genome Biology.

[13]  Sharon L. Kendall,et al.  Cholesterol utilization in mycobacteria is controlled by two TetR-type transcriptional regulators: kstR and kstR2 , 2010, Microbiology.

[14]  S. Kaufmann,et al.  Mycobacterium tuberculosis: success through dormancy. , 2012, FEMS microbiology reviews.

[15]  Christopher M. Sassetti,et al.  igr Genes and Mycobacterium tuberculosis Cholesterol Metabolism , 2009, Journal of bacteriology.

[16]  I. Smith,et al.  Global Analysis of the Mycobacterium tuberculosis Zur (FurB) Regulon , 2006, Journal of bacteriology.

[17]  T. Myers,et al.  The Transcriptional Responses of Mycobacterium tuberculosis to Inhibitors of Metabolism , 2004, Journal of Biological Chemistry.

[18]  Didier Raoult,et al.  The rpoB gene as a tool for clinical microbiologists. , 2009, Trends in microbiology.

[19]  P. Verma,et al.  Mycobacterium tuberculosis-driven targeted recalibration of macrophage lipid homeostasis promotes the foamy phenotype. , 2012, Cell host & microbe.

[20]  William R Bishai,et al.  Altered expression of isoniazid-regulated genes in drug-treated dormant Mycobacterium tuberculosis. , 2008, The Journal of antimicrobial chemotherapy.

[21]  J. Anzola,et al.  Erratum for Rodriguez et al., Global Adaptation to a Lipid Environment Triggers the Dormancy-Related Phenotype of Mycobacterium tuberculosis , 2014, mBio.

[22]  S. Guan,et al.  Mycobacterium tuberculosis CYP125A1, a steroid C27 monooxygenase that detoxifies intracellularly generated cholest‐4‐en‐3‐one , 2010, Molecular microbiology.

[23]  Matthew R. Laird,et al.  Protein Protein Interaction Network Evaluation for Identifying Potential Drug Targets , 2009 .

[24]  J. Dziadek,et al.  ChoD and HsdD can be dispensable for cholesterol degradation in mycobacteria , 2013, The Journal of Steroid Biochemistry and Molecular Biology.

[25]  Matthew D. McMahon,et al.  Analyses of MbtB, MbtE, and MbtF Suggest Revisions to the Mycobactin Biosynthesis Pathway in Mycobacterium tuberculosis , 2012, Journal of bacteriology.

[26]  B. Abomoelak,et al.  A Novel In Vitro Multiple-Stress Dormancy Model for Mycobacterium tuberculosis Generates a Lipid-Loaded, Drug-Tolerant, Dormant Pathogen , 2009, PloS one.

[27]  Y. Poquet,et al.  A central role for aspartate in Mycobacterium tuberculosis physiology and virulence , 2013, Front. Cell. Infect. Microbiol..

[28]  B. Matthews,et al.  Peptide methionine sulfoxide reductase: structure, mechanism of action, and biological function. , 2002, Archives of biochemistry and biophysics.

[29]  G. Dougan,et al.  Sequence-Based Analysis Uncovers an Abundance of Non-Coding RNA in the Total Transcriptome of Mycobacterium tuberculosis , 2011, PLoS pathogens.

[30]  J. Anzola,et al.  Global Adaptation to a Lipid Environment Triggers the Dormancy-Related Phenotype of Mycobacterium tuberculosis , 2014, mBio.

[31]  D. Young,et al.  Non-coding RNA and its potential role in Mycobacterium tuberculosis pathogenesis , 2012, RNA biology.

[32]  A. Westermann,et al.  Dual RNA-seq of pathogen and host , 2012, Nature Reviews Microbiology.

[33]  B. Barrell,et al.  Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence , 1998, Nature.

[34]  N. Sampson,et al.  Shrinking the FadE Proteome of Mycobacterium tuberculosis: Insights into Cholesterol Metabolism through Identification of an α2β2 Heterotetrameric Acyl Coenzyme A Dehydrogenase Family , 2013, Journal of bacteriology.

[35]  C. Consolandi,et al.  An efficient rRNA removal method for RNA sequencing in GC-rich bacteria , 2013, Microbial Informatics and Experimentation.

[36]  Paul J. Hertzog,et al.  Type I Interferon Receptors: Biochemistry and Biological Functions* , 2007, Journal of Biological Chemistry.

[37]  Julian Parkhill,et al.  Genome plasticity of BCG and impact on vaccine efficacy , 2007, Proceedings of the National Academy of Sciences.

[38]  A. Munro,et al.  Structural and Biochemical Characterization of Mycobacterium tuberculosis CYP142 , 2010, The Journal of Biological Chemistry.

[39]  N. Sampson,et al.  Pathogen roid rage: Cholesterol utilization by Mycobacterium tuberculosis , 2014, Critical reviews in biochemistry and molecular biology.

[40]  D. Kirschner,et al.  Intracellular Bacillary Burden Reflects a Burst Size for Mycobacterium tuberculosis In Vivo , 2013, PLoS pathogens.

[41]  Gérald Larrouy-Maumus,et al.  Mycobacterium tuberculosis nitrogen assimilation and host colonization require aspartate , 2013, Nature chemical biology.

[42]  F. Sutterwala,et al.  Cutting Edge: Mycobacterium tuberculosis but Not Nonvirulent Mycobacteria Inhibits IFN-β and AIM2 Inflammasome–Dependent IL-1β Production via Its ESX-1 Secretion System , 2013, The Journal of Immunology.

[43]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[44]  Christopher M. Sassetti,et al.  Mycobacterial persistence requires the utilization of host cholesterol , 2008, Proceedings of the National Academy of Sciences.

[45]  C. Garzelli,et al.  Oligo(dT)-primed synthesis of cDNA by reverse transcriptase in mycobacteria. , 1998, Biochemical and biophysical research communications.

[46]  V. M. D. Martins dos Santos,et al.  Systems-level modeling of mycobacterial metabolism for the identification of new (multi-)drug targets. , 2014, Seminars in Immunology.

[47]  I. Smith,et al.  The Mycobacterium tuberculosis High-Affinity Iron Importer, IrtA, Contains an FAD-Binding Domain , 2009, Journal of bacteriology.

[48]  L. Dijkhuizen,et al.  Structural Features in the KshA Terminal Oxygenase Protein That Determine Substrate Preference of 3-Ketosteroid 9α-Hydroxylase Enzymes , 2011, Journal of bacteriology.

[49]  Yang Liu,et al.  Transcriptional Adaptation of Mycobacterium tuberculosis within Macrophages , 2003, The Journal of experimental medicine.

[50]  Y. Chai,et al.  A Gene Expression Study of the Activities of Aromatic Ring-Cleavage Dioxygenases in Mycobacterium gilvum PYR-GCK to Changes in Salinity and pH during Pyrene Degradation , 2013, PloS one.

[51]  Jesse C. J. van Dam,et al.  Integration of heterogeneous molecular networks to unravel gene-regulation in Mycobacterium tuberculosis , 2014, BMC Systems Biology.

[52]  M. Niederweis,et al.  Discovery of a Siderophore Export System Essential for Virulence of Mycobacterium tuberculosis , 2013, PLoS pathogens.

[53]  G. D. de Souza,et al.  Bacterial proteins with cleaved or uncleaved signal peptides of the general secretory pathway. , 2011, Journal of proteomics.

[54]  P. Butcher,et al.  Gene expression during host—pathogen interactions: Approaches to bacterial mRNA extraction and labelling for microarray analysis , 2002 .

[55]  Gérald Larrouy-Maumus,et al.  Mycobacterium tuberculosis Exploits Asparagine to Assimilate Nitrogen and Resist Acid Stress during Infection , 2014, PLoS pathogens.

[56]  Sailu Yellaboina,et al.  Comparative analysis of iron regulated genes in mycobacteria , 2006, FEBS letters.

[57]  M. Voskuil,et al.  Unique Roles of DosT and DosS in DosR Regulon Induction and Mycobacterium tuberculosis Dormancy , 2009, Infection and Immunity.

[58]  N. Brot,et al.  Mycobacterium tuberculosis expresses methionine sulphoxide reductases A and B that protect from killing by nitrite and hypochlorite , 2009, Molecular microbiology.

[59]  Carolyn R Bertozzi,et al.  Cholesterol catabolism by Mycobacterium tuberculosis requires transcriptional and metabolic adaptations. , 2012, Chemistry & biology.

[60]  N. Sarkar,et al.  Polyadenylation of mRNA in prokaryotes. , 1997, Annual review of biochemistry.

[61]  D. Eisenberg,et al.  Mycobacterium tuberculosis gene expression profiling within the context of protein networks. , 2006, Microbes and infection.

[62]  Jun Liu,et al.  Systematic Genetic Nomenclature for Type VII Secretion Systems , 2009, PLoS pathogens.

[63]  L. Eltis,et al.  FadD3 is an acyl‐CoA synthetase that initiates catabolism of cholesterol rings C and D in actinobacteria , 2013, Molecular microbiology.

[64]  J. Livny,et al.  Identification of Novel sRNAs in Mycobacterial Species , 2013, PloS one.

[65]  Rick Lyons,et al.  The temporal expression profile of Mycobacterium tuberculosis infection in mice. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[66]  Jonathan B. Johnston,et al.  Cholesterol catabolism as a therapeutic target in Mycobacterium tuberculosis. , 2011, Trends in microbiology.

[67]  S. Fortune,et al.  Mycobacterial Esx-3 is required for mycobactin-mediated iron acquisition , 2009, Proceedings of the National Academy of Sciences.

[68]  E M Rubin,et al.  Mapping the Burkholderia cenocepacia niche response via high-throughput sequencing , 2009, Proceedings of the National Academy of Sciences.

[69]  John Chan,et al.  Responses of Mycobacterium tuberculosis to Growth in the Mouse Lung , 2005, Infection and Immunity.

[70]  M. Daffé,et al.  Mycolic acids: structures, biosynthesis, and beyond. , 2014, Chemistry & biology.

[71]  M. Lee,et al.  Development of a real‐time PCR‐based method for rapid differential identification of Mycobacterium species , 2007, Letters in applied microbiology.

[72]  N. Sampson,et al.  Mycobacterium tuberculosis utilizes a unique heterotetrameric structure for dehydrogenation of the cholesterol side chain. , 2013, Biochemistry.

[73]  Feng Chen,et al.  Statistical analysis of real-time PCR data , 2006, BMC Bioinformatics.

[74]  L. Dijkhuizen,et al.  The Steroid Catabolic Pathway of the Intracellular Pathogen Rhodococcus equi Is Important for Pathogenesis and a Target for Vaccine Development , 2011, PLoS pathogens.

[75]  W. Jacobs,et al.  Studies of a Ring-Cleaving Dioxygenase Illuminate the Role of Cholesterol Metabolism in the Pathogenesis of Mycobacterium tuberculosis , 2009, PLoS pathogens.

[76]  L. Ponnala,et al.  Deep RNA sequencing of L. monocytogenes reveals overlapping and extensive stationary phase and sigma B-dependent transcriptomes, including multiple highly transcribed noncoding RNAs , 2009, BMC Genomics.

[77]  C. Vilchèze,et al.  NAD+ auxotrophy is bacteriocidal for the tubercle bacilli , 2010, Molecular microbiology.

[78]  S. Sreevatsan,et al.  Host-Mycobacterium avium subsp. paratuberculosis interactome reveals a novel iron assimilation mechanism linked to nitric oxide stress during early infection , 2013, BMC Genomics.

[79]  I. Smith,et al.  Global Transcriptional Profile of Mycobacterium tuberculosis during THP-1 Human Macrophage Infection , 2007, Infection and Immunity.

[80]  Sandro Silva-Gomes,et al.  Heme Catabolism by Heme Oxygenase-1 Confers Host Resistance to Mycobacterium Infection , 2013, Infection and Immunity.

[81]  Samuel A. Assefa,et al.  A Strand-Specific RNA–Seq Analysis of the Transcriptome of the Typhoid Bacillus Salmonella Typhi , 2009, PLoS genetics.