Characterisation of genes differentially expressed in macrophages by virulent and attenuated Mycobacterium tuberculosis through RNA-Seq analysis

Comparative transcriptomics reveals common and strain-specific responses of human macrophages to infection with Mycobacterium tuberculosis and Mycobacterium bovis BCG

Mycobacterium tuberculosis (MTB) and Mycobacterium bovis (M. bovis) are closely related pathogenic mycobacteria known to cause chronic pulmonary infections in both humans and animals. Despite sharing nearly identical genomes and virulence factors, these two bacteria display variations in host tropism, epidemiology, and clinical presentations. M. bovis Bacillus Calmette-Guérin (BCG) is an attenuated strain of M. bovis commonly utilized as a vaccine for tuberculosis (TB). Nevertheless, the molecular underpinnings of these distinctions and the intricacies of host-pathogen interactions remain areas of ongoing research. In this study, a comparative transcriptomic analysis was conducted on human leukemia macrophages (THP-1) infected with either MTB H37Rv or M. bovis BCG (Tokyo strain) to elucidate common and strain-specific responses at the transcriptional level. RNA sequencing was utilized to characterize the transcriptomes of human primary macrophages infected with MTB or BCG at 6 and 24 h post-infection. The findings indicate that both MTB and BCG induce substantial and dynamic alterations in the transcriptomes of THP-1, with a notable overlap in the quantity and extent of differentially expressed genes (DEGs). Moreover, gene ontology (GO) enrichment analysis unveiled shared pathways related to immune response, cytokine signaling, and apoptosis. The immune response of macrophages to bacterial infections at 6 h exhibited significantly greater intensity compared to that at 24 h. Furthermore, distinct gene sets displaying notable variances between MTB and BCG infections were identified. The profound impact of MTB infection on macrophage gene expression, particularly within the initial 6 h, was evident. Additionally, downregulation of pathways such as Focal adhesion, Rap1 signaling pathway, and Regulation of actin cytoskeleton was observed. The pathways associated with inflammation reactions and cell apoptosis exhibited significant differences, with BCG triggering macrophage apoptosis and MTB enhancing the survival of intracellular bacteria. Our findings reveal that MTB and BCG provoke similar yet distinct transcriptional responses in human macrophages, indicating variations in their pathogenesis and ability to adapt to host environments. These results offer novel insights into the molecular mechanisms governing host-pathogen interactions and may contribute to a deeper understanding of TB pathogenesis.

Elevated SLC7A2 expression is associated with an abnormal neuroinflammatory response and nitrosative stress in Huntington's disease

We previously identified solute carrier family 7 member 2 (SLC7A2) as one of the top upregulated genes when normal Huntingtin was deleted. SLC7A2 has a high affinity for L-arginine. Arginine is implicated in inflammatory responses, and SLC7A2 is an important regulator of innate and adaptive immunity in macrophages. Although neuroinflammation is clearly demonstrated in animal models and patients with Huntington's disease (HD), the question of whether neuroinflammation actively participates in HD pathogenesis is a topic of ongoing research and debate. Here, we studied the role of SLC7A2 in mediating the neuroinflammatory stress response in HD cells. RNA sequencing (RNA-seq), quantitative RT-PCR and data mining of publicly available RNA-seq datasets of human patients were performed to assess the levels of SLC7A2 mRNA in different HD cellular models and patients. Biochemical studies were then conducted on cell lines and primary mouse astrocytes to investigate arginine metabolism and nitrosative stress in response to neuroinflammation. The CRISPR-Cas9 system was used to knock out SLC7A2 in STHdhQ7 and Q111 cells to investigate its role in mediating the neuroinflammatory response. Live-cell imaging was used to measure mitochondrial dynamics. Finally, exploratory studies were performed using the Enroll-HD periodic human patient dataset to analyze the effect of arginine supplements on HD progression. We found that SLC7A2 is selectively upregulated in HD cellular models and patients. HD cells exhibit an overactive response to neuroinflammatory challenges, as demonstrated by abnormally high iNOS induction and NO production, leading to increased protein nitrosylation. Depleting extracellular Arg or knocking out SLC7A2 blocked iNOS induction and NO production in STHdhQ111 cells. We further examined the functional impact of protein nitrosylation on a well-documented protein target, DRP-1, and found that more mitochondria were fragmented in challenged STHdhQ111 cells. Last, analysis of Enroll-HD datasets suggested that HD patients taking arginine supplements progressed more rapidly than others. Our data suggest a novel pathway that links arginine uptake to nitrosative stress via upregulation of SLC7A2 in the pathogenesis and progression of HD. This further implies that arginine supplements may potentially pose a greater risk to HD patients.

Pathogenic mycoplasmas of humans regulate the long noncoding RNAs in epithelial cells

Non-coding RNAs (ncRNAs), specifically long ncRNAs (lncRNAs), regulate cellular processes by affecting gene expression at the transcriptional, post-transcriptional, and epigenetic levels. Emerging evidence indicates that pathogenic microbes dysregulate the expression of host lncRNAs to suppress cellular defense mechanisms and promote survival. To understand whether the pathogenic human mycoplasmas dysregulate host lncRNAs, we infected HeLa cells with Mycoplasma genitalium (Mg) and Mycoplasma penumoniae (Mp) and assessed the expression of lncRNAs by directional RNA-seq analysis. HeLa cells infected with these species showed up-and-down regulation of lncRNAs expression, indicating that both species can modulate host lncRNAs. However, the number of upregulated (200 for Mg and 112 for Mp) and downregulated lncRNAs (30 for Mg and 62 for Mp) differ widely between these two species. GREAT analysis of the noncoding regions associated with differentially expressed lncRNAs showed that Mg and Mp regulate a discrete set of lncRNA plausibly related to transcription, metabolism, and inflammation. Further, signaling network analysis of the differentially regulated lncRNAs exhibited diverse pathways such as neurodegeneration, NOD-like receptor signaling, MAPK signaling, p53 signaling, and PI3K signaling, suggesting that both species primarily target signaling mechanisms. Overall, the study's results suggest that Mg and Mp modulate lncRNAs to promote their survival within the host but in distinct manners.

Comparative interleukins and chemokines analysis of mice mesenchymal stromal cells infected with Mycobacterium tuberculosis H37Rv and H37Ra

Inflammatory pathways involving Mesenchymal stromal cells (MSCs) play an important role in Mycobacterium tuberculosis (Mtb) infection. H37Rv (Rv) is a standard virulent strain, however, H37Ra (Ra) is a strain with reduced virulence. Interleukins and chemokines production are known to promote inflammation resistance in mammalian cells and is recently reported to regulate mycobacterial immunopathogenesis via inflammatory responses. MSCs are very important cells during Mtb infection. However, the different expressions of interleukins and chemokines in the process of Mtb-infected MSCs between Ra and Rv remain unclear. We used the techniques of RNA-Seq, Q-RT-PCR, ELISA, and Western Blotting. We have shown that Rv infection significantly increased mRNA expressions of Mndal, Gdap10, Bmp2, and Lif, thereby increasing more differentiation of MSCs compared with Ra infection in MSCs. Further investigation into the possible mechanisms, we found that Rv infection enhanced more inflammatory response (Mmp10, Mmp3, and Ptgs2) through more activation of the TLR2-MAP3K1-JNK pathway than did Ra infection in MSCs. Further action showed that Rv infection enhanced more Il1α, Il6, Il33, Cxcl2, Ccl3, and Ackr3 production than did Ra infection. Rv infection showed more expressions of Mmp10, Mmp3, Ptgs2, Il1α, Il6, Il33, Cxcl2, Ccl3, and Ackr3 possibly through more active TLR2-MAP3K1-JNK pathway than did Ra infection in MSCs. MSCs may therefore be a new candidate for the prevention and treatment of tuberculosis.

Advances in computational frameworks in the fight against TB: The way forward

Around 1.6 million people lost their life to Tuberculosis in 2021 according to WHO estimates. Although an intensive treatment plan exists against the causal agent, Mycobacterium Tuberculosis, evolution of multi-drug resistant strains of the pathogen puts a large number of global populations at risk. Vaccine which can induce long-term protection is still in the making with many candidates currently in different phases of clinical trials. The COVID-19 pandemic has further aggravated the adversities by affecting early TB diagnosis and treatment. Yet, WHO remains adamant on its "End TB" strategy and aims to substantially reduce TB incidence and deaths by the year 2035. Such an ambitious goal would require a multi-sectoral approach which would greatly benefit from the latest computational advancements. To highlight the progress of these tools against TB, through this review, we summarize recent studies which have used advanced computational tools and algorithms for-early TB diagnosis, anti-mycobacterium drug discovery and in the designing of the next-generation of TB vaccines. At the end, we give an insight on other computational tools and Machine Learning approaches which have successfully been applied in biomedical research and discuss their prospects and applications against TB.

Involvement of 2′-5′ oligoadenylate synthetase-like protein in the survival of Mycobacterium tuberculosis avirulent strain in macrophages

Mycobacterium tuberculosis (M. tuberculosis) can replicate in the macrophage by interfering with many host protein functions. While it is far from known these host proteins for controlling M. tuberculosis infection. Herein, we infected macrophages including THP-1 and Raw264.7 cells with M. tuberculosis and identified the differentially expressed genes (DEGs) in the interferon signaling pathway. Among them, 2′-5′ oligoadenylate synthetase-like (OASL) underwent the greatest upregulation in M. tuberculosis-infected macrophages. Knockdown of the expression of OASL attenuated M. tuberculosis survival in macrophages. Further, bioinformatics analysis revealed the potential interaction axis of OASL-TAB3- Rv0127, which was further validated by the yeast-two-hybrid (Y2H) assay and Co-IP. This interaction axis might regulate the M. tuberculosis survival and proliferation in macrophages. The study reveals a possible role of OASL during M. tuberculosis infection as a target to control its propagation.

Macrophage NFATC2 mediates angiogenic signaling during mycobacterial infection

During mycobacterial infections, pathogenic mycobacteria manipulate both host immune and stromal cells to establish and maintain a productive infection. In humans, non-human primates, and zebrafish models of infection, pathogenic mycobacteria produce and modify the specialized lipid trehalose 6,6'-dimycolate (TDM) in the bacterial cell envelope to drive host angiogenesis toward the site of forming granulomas, leading to enhanced bacterial growth. Here, we use the zebrafish-Mycobacterium marinum infection model to define the signaling basis of the host angiogenic response. Through intravital imaging and cell-restricted peptide-mediated inhibition, we identify macrophage-specific activation of NFAT signaling as essential to TDM-mediated angiogenesis in vivo. Exposure of cultured human cells to Mycobacterium tuberculosis results in robust induction of VEGFA, which is dependent on a signaling pathway downstream of host TDM detection and culminates in NFATC2 activation. As granuloma-associated angiogenesis is known to serve bacterial-beneficial roles, these findings identify potential host targets to improve tuberculosis disease outcomes.

Untargeted metabolomics analysis reveals Mycobacterium tuberculosis strain H37Rv specifically induces tryptophan metabolism in human macrophages

Background

Tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tb) remains a global health issue. The characterized virulent M. tb H37Rv, avirulent M. tb H37Ra and BCG strains are widely used as reference strains to investigate the mechanism of TB pathogenicity. Here, we attempted to determine metabolomic signatures associated with the Mycobacterial virulence in human macrophages through comparison of metabolite profile in THP-1-derived macrophages following exposure to the M. tb H37Rv, M. tb H37Ra and BCG strains.

Results

Our findings revealed remarkably changed metabolites in infected macrophages compared to uninfected macrophages. H37Rv infection specifically induced 247 differentially changed metabolites compared to H37Ra or BCG infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed H37Rv specifically induces tryptophan metabolism. Moreover, quantitative PCR (qPCR) results showed that indoleamine 2,3-dioxygenase 1 (IDO1) and tryptophan 2,3-dioxygenase 2 (TDO2) which converts the tryptophan to a series of biologically second metabolites were up-regulated in H37Rv-infected macrophages compared to H37Ra- or BCG-infected macrophages, confirming the result of enhanced tryptophan metabolism induced by H37Rv infection. These findings indicated that targeting tryptophan (Trp) metabolism may be a potential therapeutic strategy for pulmonary TB.

Conclusions

We identified a number of differentially changed metabolites that specifically induced in H37Rv infected macrophages. These signatures may be associated with the Mycobacterial virulence in human macrophages. The present findings provide a better understanding of the host response associated with the virulence of the Mtb strain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02659-y.

Diagnosis of pulmonary tuberculosis via identification of core genes and pathways utilizing blood transcriptional signatures: a multicohort analysis

Background

Blood transcriptomics can be used for confirmation of tuberculosis diagnosis or sputumless triage, and a comparison of their practical diagnostic accuracy is needed to assess their usefulness. In this study, we investigated potential biomarkers to improve our understanding of the pathogenesis of active pulmonary tuberculosis (PTB) using bioinformatics methods.

Methods

Differentially expressed genes (DEGs) were analyzed between PTB and healthy controls (HCs) based on two microarray datasets. Pathways and functional annotation of DEGs were identified and ten hub genes were selected. They were further analyzed and selected, then verified with an independent sample set. Finally, their diagnostic power was further evaluated between PTB and HCs or other diseases.

Results

62 DEGs mostly related to type I IFN pathway, IFN-γ-mediated pathway, etc. in GO term and immune process, and especially RIG-I-like receptor pathway were acquired. Among them, OAS1, IFIT1 and IFIT3 were upregulated and were the main risk factors for predicting PTB, with adjusted risk ratios of 1.36, 3.10, and 1.32, respectively. These results further verified that peripheral blood mRNA expression levels of OAS1, IFIT1 and IFIT3 were significantly higher in PTB patients than HCs (all P < 0.01). The performance of a combination of these three genes (three-gene set) had exceeded that of all pairwise combinations of them in discriminating TB from HCs, with mean AUC reaching as high as 0.975 with a sensitivity of 94.4% and a specificity of 100%. The good discernibility capacity was evaluated d via 7 independent datasets with an AUC of 0.902, as well as mean sensitivity of 87.9% and mean specificity of 90.2%. In regards to discriminating PTB from other diseases (i.e., initially considered to be possible TB, but rejected in differential diagnosis), the three-gene set equally exhibited an overall strong ability to separate PTB from other diseases with an AUC of 0.999 (sensitivity: 99.0%; specificity: 100%) in the training set, and 0.974 with a sensitivity of 96.4% and a specificity of 98.6% in the test set.

Conclusion

The described commonalities and unique signatures in the blood profiles of PTB and the other control samples have considerable implications for PTB biosignature design and future diagnosis, and provide insights into the biological processes underlying PTB.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-022-02035-4.

Amino Acid Solute Carrier Transporters in Inflammation and Autoimmunity

The past decade exposed the importance of many homeostasis and metabolism related proteins in autoimmunity disease and inflammation. Solute carriers (SLCs) are a group of membrane channels that can transport amino acids, the building blocks of proteins, nutrients, and neurotransmitters. This review summarizes the role of SLCs amino acid transporters in inflammation and autoimmunity disease. In detail, the importance of Glutamate transporters SLC1A1, SLC1A2, and SLC1A3, mainly expressed in the brain where they help prevent glutamate excitotoxicity, is discussed in the context of central nervous system disorders such as multiple sclerosis. Similarly, the cationic amino acid transporter SLC7A1 (CAT1), which is an important arginine transporter for T cells, and SLC7A2 (CAT2), essential for innate immunity. SLC3 family proteins, which bind with light chains from the SLC7 family (SLC7A5, SLC7A7 and SLC7A11) to form heteromeric amino acid transporters, are also explored to describe their roles in T cells, NK cells, macrophages and tumor immunotherapies. Altogether, the link between SLC amino acid transporters with inflammation and autoimmunity may contribute to a better understanding of underlying mechanism of disease and provide novel potential therapeutic avenues. Significance Statement SIGNIFICANCE STATEMENT In this review, we summarize the link between SLC amino acid transporters and inflammation and immune responses, specially SLC1 family members and SLC7 members. Studying the link may contribute to a better understanding of related diseases and provide potential therapeutic targets and useful to the researchers who have interest in the involvement of amino acids in immunity.

Commonalities of Mycobacterium tuberculosis Transcriptomes in Response to Defined Persisting Macrophage Stresses

As the goal of a bacterium is to become bacteria, evolution has imposed continued selections for gene expression. The intracellular pathogen Mycobacterium tuberculosis, the causative agent of tuberculosis, has adopted a fine-tuned response to survive its host's methods to aggressively eradicate invaders. The development of microarrays and later RNA sequencing has led to a better understanding of biological processes controlling the relationship between host and pathogens. In this study, RNA-seq was performed to detail the transcriptomes of M. tuberculosis grown in various conditions related to stresses endured by M. tuberculosis during host infection and to delineate a general stress response incurring during persisting macrophage stresses. M. tuberculosis was subjected to long-term growth, nutrient starvation, hypoxic and acidic environments. The commonalities between these stresses point to M. tuberculosis maneuvering to exploit propionate metabolism for lipid synthesis or to withstand propionate toxicity whilst in the intracellular environment. While nearly all stresses led to a general shutdown of most biological processes, up-regulation of pathways involved in the synthesis of amino acids, cofactors, and lipids were observed only in hypoxic M. tuberculosis. This data reveals genes and gene cohorts that are specifically or exclusively induced during all of these persisting stresses. Such knowledge could be used to design novel drug targets or to define possible M. tuberculosis vulnerabilities for vaccine development. Furthermore, the disruption of specific functions from this gene set will enhance our understanding of the evolutionary forces that have caused the tubercle bacillus to be a highly successful pathogen.

Mycoplasma genitalium and M. pneumoniae Regulate a Distinct Set of Protein-Coding Genes in Epithelial Cells

Mycoplasma genitalium and M. pneumoniae are two significant mycoplasmas that infect the urogenital and respiratory tracts of humans. Despite distinct tissue tropisms, they both have similar pathogenic mechanisms and infect/invade epithelial cells in the respective regions and persist within these cells. However, the pathogenic mechanisms of these species in terms of bacterium-host interactions are poorly understood. To gain insights on this, we infected HeLa cells independently with M. genitalium and M. pneumoniae and assessed gene expression by whole transcriptome sequencing (RNA-seq) approach. The results revealed that HeLa cells respond to M. genitalium and M. pneumoniae differently by regulating various protein-coding genes. Though there is a significant overlap between the genes regulated by these species, many of the differentially expressed genes were specific to each species. KEGG pathway and signaling network analyses revealed that the genes specific to M. genitalium are more related to cellular processes. In contrast, the genes specific to M. pneumoniae infection are correlated with immune response and inflammation, possibly suggesting that M. pneumoniae has some inherent ability to modulate host immune pathways.

Small Noncoding RNAs and Their Role in the Pathogenesis of Mycobacterium tuberculosis Infection

Mycobacterium tuberculosis possesses a significant arsenal of strategies to combat immune defense of the host organism. Small noncoding RNAs, which constitute the largest group of regulatory RNAs, play an important role in the host–pathogen interactions and represent one of the levels of the regulation of interactions of microbial cells with their environment. The regulatory role of small RNAs in pathogenic bacteria is essential when rapid adaptation to the changing environmental conditions with further synchronization of metabolic reactions are required to ensure microbial survival and infection progression. During the past few years, eight small RNAs from M. tuberculosis have been functionally characterized, and targets for four of them have been identified. Small RNAs from M. tuberculosis and other pathogenic microorganisms were found to be one of the most important functional factors in the adaptive response to changing environmental conditions.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

Background

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne’s disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (−); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip.

Results

The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/−) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10^− 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids.

Conclusion

In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07487-4.

The Transcriptional Regulatory Network of Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis is a Gram-positive, facultative intracellular, pathogenic bacterium that infects several different hosts, yielding serious economic losses in livestock farming. It causes several diseases including oedematous skin disease (OSD) in buffaloes, ulcerative lymphangitis (UL) in horses, and caseous lymphadenitis (CLA) in sheep, goats and humans. Despite its economic and medical-veterinary importance, our understanding concerning this organism’s transcriptional regulatory mechanisms is still limited. Here, we review the state of the art knowledge on transcriptional regulatory mechanisms of this pathogenic species, covering regulatory interactions mediated by two-component systems, transcription factors and sigma factors. Key transcriptional regulatory players involved in virulence and pathogenicity of C. pseudotuberculosis, such as the PhoPR system and DtxR, are in the focus of this review, as these regulators are promising targets for future vaccine design and drug development. We conclude that more experimental studies are needed to further understand the regulatory repertoire of this important zoonotic pathogen, and that regulators are promising targets for future vaccine design and drug development.

New Conjugated Compound T5 Epidioxy-Sterol-ANB Inhibits the Growth of Mycobacterium tuberculosis Affecting the Cholesterol and Folate Pathways

The upsurge and persistence of drug resistant strains of Mycobacterium tuberculosis (Mtb) is an important limitant to the battery of drugs available for the elimination of tuberculosis (TB). To avoid future scarcity of antibiotics against Mtb, it is important to discover new effective anti-mycobacterial agents. In this study, we present data from a series of experiments to determine in vitro and in vivo anti-mycobacterial activity of a library of epidioxy-sterol analogs. We test 15 compounds for their ability to reduce the viability of Mtb. We found that one compound called T5 epidioxy-sterol-ANB display significant potency against Mtb in vitro specifically inside macrophages but without effectivity in axenic cultures. A viability assay confirms that this T5 compound is less toxic for macrophages in vitro as compared to the current Mtb drug Rifampicin at higher concentrations. We use a transcriptomic analysis of Mtb inside macrophages after T5 epidioxy-sterol-ANB treatment, and we found a significant down-regulation of enzymes involved in the cholesterol and folic acid pathways. In vivo, significant differences were found in the lungs and spleen CFUs of Mtb infected mice treated with the T5 epidioxy-sterol-ANB as compared with the untreated control group, which provides additional evidence of the effectivity of the T5 compound. Altogether these results confirm the potential of this T5 epidioxy-sterol-ANB compound against Mtb.

Transcriptome Changes of Mycobacterium marinum in the Process of Resuscitation From Hypoxia-Induced Dormancy

Nearly one-third of the world's population is latently infected with Mycobacterium tuberculosis (M. tb), which represents a huge disease reservoir for reactivation and a major obstacle for effective control of tuberculosis. During latent infection, M. tb is thought to enter nonreplicative dormant states by virtue of its response to hypoxia and nutrient-deprived conditions. Knowledge of the genetic programs used to facilitate entry into and exit from the nonreplicative dormant states remains incomplete. In this study, we examined the transcriptional changes of Mycobacterium marinum (M. marinum), a pathogenic mycobacterial species closely related to M. tb, at different stages of resuscitation from hypoxia-induced dormancy. RNA-seq analyses were performed on M. marinum cultures recovered at multiple time points after resuscitation. Differentially expressed genes (DEGs) at each time period were identified and analyzed. Co-expression networks of transcription factors and DEGs in each period were constructed. In addition, we performed a weighted gene co-expression network analysis (WGCNA) on all genes and obtained 12 distinct gene modules. Collectively, these data provided valuable insight into the transcriptome changes of M. marinum upon resuscitation as well as gene module function of the bacteria during active metabolism and growth.