Tuberculosis (TB) remains one of the most common infectious diseases caused by Mycobacterium tuberculosis complex (MTBC). To panoramically analyze MTBC's genomic methylation, we completed the genomes of 12 MTBC strains (Mycobacterium bovis; M. bovis BCG; M. microti; M. africanum; M. tuberculosis H37Rv; H37Ra; and 6 M. tuberculosis clinical isolates) belonging to different lineages and characterized their methylomes using single-molecule real-time (SMRT) technology. We identified three m6A sequence motifs and their corresponding methyltransferase (MTase) genes, including the reported mamA, hsdM and a newly discovered mamB. We also experimentally verified the methylated motifs and functions of HsdM and MamB. Our analysis indicated the MTase activities varied between 12 strains due to mutations/deletions. Furthermore, through measuring ‘the methylated-motif-site ratio’ and ‘the methylated-read ratio’, we explored the methylation status of each modified site and sequence-read to obtain the ‘precision methylome’ of the MTBC strains, which enabled intricate analysis of MTase activity at whole-genome scale. Most unmodified sites overlapped with transcription-factor binding-regions, which might protect these sites from methylation. Overall, our findings show enormous potential for the SMRT platform to investigate the precise character of methylome, and significantly enhance our understanding of the function of DNA MTase.

Tuberculosis (TB) has been the leading lethal infectious disease worldwide since 2014, and about one third of the world’s population has a latent TB infection (LTBI). This is largely attributed to the difficulties in diagnosis and treatment of TB and LTBI patients. Exosomes offer a new perspective on investigation of the process of TB infection. In this study, we performed small RNA sequencing to explore small RNA profiles of serum exosomes derived from LTBI and TB patients and healthy controls (HC). Our results revealed distinct miRNA profile of the exosomes in the three groups. We screened 250 differentially expressed miRNAs including 130 specifically expressed miRNAs. Some miRNAs were further validated to be specifically expressed in LTBI (hsa-let-7e-5p, hsa-let-7d-5p, hsa-miR-450a-5p, and hsa-miR-140-5p) and TB samples (hsa-miR-1246, hsa-miR-2110, hsa-miR-370-3P, hsa-miR-28-3p, and hsa-miR-193b-5p). Additionally, we demonstrated four expression panels in LTBI and TB groups, and six expression patterns among the three groups. These specifically expressed miRNAs and differentially expressed miRNAs in different panels and patterns provide potential biomarkers for detection/diagnosis of latent and active TB using exosomal miRNAs. Additionally, we also discovered plenty of small RNAs derived from genomic repetitive sequences, which might play roles in host immune responses along with Mtb infection progresses. Overall, our findings provide important reference and an improved understanding about miRNAs and repetitive region-derived small RNAs in exosomes during the Mtb infectious process, and facilitate the development of potential molecular targets for detection/diagnosis of latent and active tuberculosis.

Increasing evidence shows that gut microbiota dysbiosis may play important roles in the occurrence and progression of Parkinson's disease (PD), but the findings are inconsistent. Besides, the effect of family environment on gut microbiota dysbiosis remains unclear. Methods: We characterized the gut microbial compositions of 63 PD patients, 63 healthy spouses (HS) and 74 healthy people (HP) using 16S rRNA sequencing. Clinical phenotypes and microbial composition were analyzed comprehensively. Results: There were markedly different microbial compositions among PD, HS and HP samples by alpha/beta diversity. We also found differential microbial compositions among Hoehn & Yahr stage/disease duration. Eight inflammation-associated microbial genera shared a continuously increase trend with increased Hoehn & Yahr stage and disease duration, indicating characteristic bacteria associated with deterioration in PD. Additionally, seven bacterial markers were identified for accurately differentiating PD patients from the controls (area under the curve [AUC]: 0.856). Conclusions: Our study shows altered gut microbiota in PD patients. Importantly, inflammation-associated microbial genera may play roles in PD progression. Differential microbial compositions in HS and HP samples demonstrate that the gut microbiota are also affected by family environment. Disease-associated metagenomics studies should consider the family environmental factor. Our research provides an important reference and improves the understanding of gut microbiota in PD patients.

Here, we report a multidrug-resistant hypervirulent Klebsiella pneumoniae (MDR-HvKP) strain of sequence type 23 (ST23) with a rare hybrid plasmid harboring virulence genes and blaCTX-M-24, and we analyze the genetic basis for relationship between genotypes and MDR-hypervirulence phenotypes. Further analysis indicates that the hybrid plasmid is formed by IS903D-mediated intermolecular transposition of the blaCTX-M-24 gene into the virulence plasmid. The emergence of MDR-HvKP strains, especially those carrying drug-resistant virulent plasmids, poses unprecedented threats/challenges to public health. This is a dangerous trend and should be closely monitored.

Tuberculosis (TB) has surpassed HIV as the leading infectious disease killer worldwide since 2014. The main pathogen, Mycobacterium tuberculosis (Mtb), contains ~4,000 genes that account for ~90% of the genome. However, it is still unclear which of these genes are primary/secondary, which are responsible for generality/individuality, and which interconvert during evolution. Here we utilized a pan-genomic analysis of 36 Mtb genomes to address these questions. We identified 3,679 Mtb core (i.e., primary) genes, determining their phenotypic generality (e.g., virulence, slow growth, dormancy). We also observed 1,122 dispensable and 964 strain-specific secondary genes, reflecting partially shared and lineage-/strain-specific individualities. Among which, five L2 lineage-specific genes might be related to the increased virulence of the L2 lineage. Notably, we discovered 28 Mtb “Super Core Genes” (SCGs: more than a copy in at least 90% strains), which might be of increased importance, and reflected the “super phenotype generality.” Most SCGs encode PE/PPE, virulence factors, antigens, and transposases, and have been verified as playing crucial roles in Mtb pathogenicity. Further investigation of the 28 SCGs demonstrated the interconversion among SCGs, single-copy core, dispensable, and strain-specific genes through copy number variations (CNVs) during evolution; different mutations on different copies highlight the delicate adaptive-evolution regulation amongst Mtb lineages. This reflects that the importance of genes varied through CNVs, which might be driven by selective pressure from environment/host-adaptation. In addition, compared with Mycobacterium bovis (Mbo), Mtb possesses 48 specific single core genes that partially reflect the differences between Mtb and Mbo individuality.

Circular RNA (circRNA) is a novel RNA molecule that has become a research focus recently. Although some research indicated that the circRNAs in different subcellular compartments could execute different regulatory functions, a panoramic analysis of the subcellular distribution and the transport mechanism of circRNA is still required. In this study, we comprehensively analyzed the subcellular distribution/ characteristics and the transport mechanism, through systemically investigating the circRNA profiles among the subcellular fractions of HepG2 cell (nucleus, cytoplasm, mitochondria, ribosome, cytosol and exosome). CircRNAs were widely distributed among the subcellular fractions except in the mitochondria, with differences in the subcellular distribution/characteristics in terms of classification, length, GC content, alternative circularization and parental gene function. Further analysis indicated this might be due to the selective transportation mediated by the transport-related RNA binding proteins (RBPs). The circRNAs may follow the same transportation mechanism of linear RNAs, in which the RBPs specially recognize/transport the RNAs with the corresponding binding motifs. Interestingly, we found that the exosome could selectively package the circRNAs containing the purine-rich 5ʹ-GMWGVWGRAG-3ʹ motif, with the characteristic of 'garbage dumping' and 'intercellular signaling' functions. Besides, although we observed numerous circRNAs enriched in the ribosome, we did not reliably identify any unique-peptides from circRNAs using 3D-LC-MS/MS strategy. This suggests that circRNAs rarely function as translation templates in vivo like lincRNA. Our findings not only indicates the differential distributions/characteristics among the subcellular fractions, but also reveals the possible transportation mechanism. This provides an improved understanding of the life history and molecular behavior of circRNA in cells.

Tuberculosis (TB) has exceeded HIV as the most lethal infectious disease globally for two consecutive years. Moreover, one third of the world’s population is estimated to have latent tuberculosis infection (LTBI). This is mainly because of difficulties associated with diagnosis and treatment for both TB and LTBI patients. Exosomes provide a promising research tool for TB diagnosis and treatment because they are released from various cells containing valuable biochemical information related to disease. In this study, we performed RNA-sequencing analysis on exosomes derived from clinical specimens of healthy controls (HC), active tuberculosis (ATB), and LTBI patients. Our results revealed the distinct gene expression profiles of the exosomes from LTBI and ATB patients. (1) We identified many distinct up-regulated and down-regulated differentially expressed genes (DEGs) in LTBI and ATB samples, and further screened the top-20 DEGs which might provide a potential panel for differentiation of HC, LTBI, and ATB. (2) We classified all the DEGs into six expression patterns, screened the top-20 genes in each pattern, and mainly focused on those highly expressed in LTBI and ATB. (3) Some Mycobacterium tuberculosis (Mtb) RNAs were only enriched in the exosomes of LTBI samples. (4) Pathway and function analysis further indicated down-regulated signaling pathways/immune response and up-regulated apoptosis/necrosis. Our findings indicate the selective packaging of RNA cargoes into exosomes under different stages of Mtb infection, while facilitating the development of potential targets for the diagnosis, prevention and treatment of tuberculosis.

Background The increase in multidrug-resistant tuberculosis (MDR-TB) severely hampers tuberculosis prevention and control in China, a country with the second highest MDR-TB burden globally. The first nationwide drug-resistant tuberculosis surveillance program provides an opportunity to comprehensively investigate the epidemiological/drug-resistance characteristics, potential drug-resistance mutations, and effective population changes of Chinese MDR-TB. Methods We sequenced 357 MDR strains from 4600 representative tuberculosis-positive sputum samples collected during the survey (70 counties in 31 provinces). Drug-susceptibility testing was performed using 18 anti-tuberculosis drugs, representing the most comprehensive drug-resistance profile to date. We used 3 statistical and 1 machine-learning methods to identify drug-resistance genes/single-nucleotide polymorphisms (SNPs). We used Bayesian skyline analysis to investigate changes in effective population size. Results Epidemiological/drug-resistance characteristics showed different MDR profiles, co-resistance patterns, preferred drug combination/use, and recommended regimens among 7 Chinese administrative regions. These factors not only reflected the serious multidrug co-resistance and drug misuse but they were also potentially significant in facilitating the development of appropriate regimens for MDR-TB treatment in China. Further investigation identified 86 drug-resistance genes/intergenic regions/SNPs (58 new), providing potential targets for MDR-TB diagnosis and treatment. In addition, the effective population of Chinese MDR-TB displayed a strong expansion during 1993–2000, reflecting socioeconomic transition within the country. The phenomenon of expansion was restrained after 2000, likely attributable to the advances in diagnosis/treatment technologies and government support. Conclusions Our findings provide an important reference and improved understanding of MDR-TB in China, which are potentially significant in achieving the goal of precision medicine with respect to MDR-TB prevention and treatment.