Co-Expression Networks for Causal Gene Identification Based on RNA-Seq Data of Corynebacterium pseudotuberculosis

RIP-seq reveals RNAs that interact with RNA polymerase and primary sigma factors in bacteria

Bacteria have evolved structured RNAs that can associate with RNA polymerase (RNAP). Two of them have been known so far-6S RNA and Ms1 RNA but it is unclear if any other types of RNAs binding to RNAP exist in bacteria. To identify all RNAs interacting with RNAP and the primary σ factors, we have established and performed native RIP-seq in Bacillus subtilis, Corynebacterium glutamicum, Streptomyces coelicolor, Mycobacterium smegmatis and the pathogenic Mycobacterium tuberculosis. Besides known 6S RNAs in B. subtilis and Ms1 in M. smegmatis, we detected MTS2823, a homologue of Ms1, on RNAP in M. tuberculosis. In C. glutamicum, we discovered novel types of structured RNAs that associate with RNAP. Furthermore, we identified other species-specific RNAs including full-length mRNAs, revealing a previously unknown landscape of RNAs interacting with the bacterial transcription machinery.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

OBJECTIVE

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

METHODS

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

RESULTS

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

CONCLUSION

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

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.