Analysis of the noncoding RNA regulatory networks of H37Rv- and H37Rv△1759c-infected macrophages

Abstract

Noncoding RNAs regulate the process of Mycobacterium tuberculosis (M. tb) infecting the host, but there is no simultaneous transcriptional information of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) and the global regulatory networks of non-coding RNA. Rv1759c, a virulence factor, is a member of protein family containing the proline-glutamic acid (PE) in M. tb, which can increase M. tb survival. To reveal the noncoding RNA regulatory networks and the effect of Rv1759c on non-coding RNA expression during M. tb infection, we collected samples of H37Rv- and H37Rv△1759c-infected macrophages and explored the full transcriptome expression profile. We found 356 mRNAs, 433 lncRNAs, 168 circRNAs, and 12 miRNAs differentially expressed during H37Rv infection, 356 mRNAs, 433 lncRNAs, 168 circRNAs, and 12 miRNAs differentially expressed during H37Rv△1759c infection. We constructed lncRNA/circRNA-miRNA-mRNA regulatory networks during H37Rv and H37Rv△1759c infection. We demonstrated the role of one of the hubs of the networks, hsa-miR-181b-3p, for H37Rv survival in macrophages. We discovered that the expression changes of 68 mRNAs, 92 lncRNAs, 26 circRNAs, and 3 miRNAs were only related to the deletion of Rv1759c by comparing the transcription profiles of H37Rv and H37Rv△1759c. Here, our study comprehensively characterizes the transcriptional profiles in THP1-derived-macrophages infected with H37Rv and H37Rv△1759c, which provides support and new directions for in-depth exploration of noncoding RNA and PE/PPE family functions during the infection process.