Revealing of Mycobacterium marinum transcriptome by RNA-seq

Rok from B. subtilis: Bridging genome structure and transcription regulation

Bacterial genomes are folded and organized into compact yet dynamic structures, called nucleoids. Nucleoid orchestration involves many factors at multiple length scales, such as nucleoid-associated proteins and liquid-liquid phase separation, and has to be compatible with replication and transcription. Possibly, genome organization plays an intrinsic role in transcription regulation, in addition to classical transcription factors. In this review, we provide arguments supporting this view using the Gram-positive bacterium Bacillus subtilis as a model. Proteins BsSMC, HBsu and Rok all impact the structure of the B. subtilis chromosome. Particularly for Rok, there is compelling evidence that it combines its structural function with a role as global gene regulator. Many studies describe either function of Rok, but rarely both are addressed at the same time. Here, we review both sides of the coin and integrate them into one model. Rok forms unusually stable DNA-DNA bridges and this ability likely underlies its repressive effect on transcription by either preventing RNA polymerase from binding to DNA or trapping it inside DNA loops. Partner proteins are needed to change or relieve Rok-mediated gene repression. Lastly, we investigate which features characterize H-NS-like proteins, a family that, at present, lacks a clear definition.

Targeted RNA-Seq Reveals the M. tuberculosis Transcriptome from an In Vivo Infection Model

Simple Summary

High-throughput sequencing techniques such as RNA-seq allow a more detailed characterization of the gene expression profile during in vivo infections. However, using this strategy for intracellular pathogens such as Mycobacterium tuberculosis (Mtb) entails technical limitations. Some authors have resorted to flow cytometers to separate infected cells or significantly increase sequencing depth to obtain pathogens’ gene expression. However, these options carry additional expenses in specialized equipment. We propose an experimental protocol based on differential cell lysis and a probe-based ribosomal depletion to determine the gene expression of Mtb and its host during in vivo infection. This method allowed us to increase the number of observed expressed genes from 13 using a traditional RNA-seq approach to 702. In addition, we observed the expression of genes essential for establishing the infection, codifying proteins such as PE-PGRS, lipoproteins lppN and LpqH, and three ncRNAs (small RNA MTS2823, transfer-messenger RNA RF00023, and ribozyme RF00010). We believe our method represents a valuable alternative to current RNA-seq approaches to study host–pathogen interactions and will help explore host–pathogen mechanisms in tuberculosis and other similar models of intracellular infections.

Abstract

The study of host-pathogen interactions using in vivo models with intracellular pathogens like Mycobacterium tuberculosis (Mtb) entails technical limitations, such as: (i) Selecting an efficient differential lysis system to enrich the pathogen cells; (ii) obtaining sufficient high-quality RNA; and (iii) achieving an efficient rRNA depletion. Thus, some authors had used flow cytometers to separate infected cells or significantly increase the sequencing depth of host–pathogen RNA libraries to observe the pathogens’ gene expression. However, these options carry additional expenses in specialized equipment typically not available for all laboratories. Here, we propose an experimental protocol involving differential cell lysis and a probe-based ribosomal depletion to determine the gene expression of Mtb and its host during in vivo infection. This method increased the number of observed pathogen-expressed genes from 13 using the traditional RNA-seq approach to 702. After eliminating rRNA reads, we observed that 61.59% of Mtb sequences represented 702 genes, while 38.41% represented intergenic regions. Some of the most expressed genes codified for IS1081 (Rv2512c) transposase and eight PE-PGRS members, such as PGRS49 and PGRS50. As expected, a critical percent of the expressed genes codified for secreted proteins essential for infection, such as PE68, lppN, and LpqH. Moreover, three Mtb ncRNAs were highly expressed (small RNA MTS2823, transfer-messenger RNA RF00023, and ribozyme RF00010). Many of the host-expressed genes were related to the inflammation process and the expression of surfactant proteins such as the Sftpa and Sftpc, known to bind Mtb to alveolar macrophages and mi638, a microRNA with no previous associations with pulmonary diseases. The main objective of this study is to present the method, and a general catalog of the Mtb expressed genes at one point of the in vivo infection. We believe our method represents a different approach to the existing ones to study host–pathogen interactions in tuberculosis and other similar intracellular infections, without the necessity of specialized equipment.

Comparative transcriptome analysis reveals key insights into male sterility in Salvia miltiorrhiza Bunge

Background

Large-scale heterosis breeding depends upon stable, inherited male sterility lines. We accidentally discovered a male sterility line (SW-S) in the F₁progeny of a Salvia miltiorrhiza Bunge from Shandong, China (purple flowers) crossed with a S. miltiorrhiza f. alba from Sichuan, China (white flowers). We sought to provide insights into the pollen development for male sterility in S. miltiorrhiza.

Methods

The phenotypic and cytological features of the SW-S and fertile control SW-F were observed using scanning electron microscopy and paraffin sections to identify the key stage of male sterility. Transcriptome profiles were recorded for anthers at the tetrad stage of SW-S and SW-F using Illumina RNA-Seq.

Results

The paraffin sections showed that sterility mainly occurred at the tetrad stage of microspore development, during which the tapetum cells in the anther compartment completely fell off and gradually degraded in the sterile line. There was little-to-no callose deposited around the microspore cells. The tetrad microspore was shriveled and had abnormal morphology. Therefore, anthers at the tetrad stage of SW-S and fertile control SW-F were selected for comparative transcriptome analysis. In total, 266,722,270 clean reads were obtained from SW-S and SW-F, which contained 36,534 genes. There were 2,571 differentially expressed genes (DEGs) in SW-S and SW-F, of which 63.5% were downregulated. Gene Ontology (GO) enrichment analysis indicated that the differentially expressed genes were enriched in 56 functional groups (GO terms); of these, all DEGs involved in microgametogenesis and developmental maturation were downregulated in SW-S. These results were confirmed by quantitative RT-PCR. The two GO terms contained 18 DEGs, among which eight DEGs (namely: GPAT3, RHF1A, phosphatidylinositol, PFAS, MYB96, MYB78, Cals5, and LAT52) were related to gamete development. There were 10 DEGs related to development and maturation, among which three genes were directly related to pollen development (namely: ACT3, RPK2, and DRP1C). Therefore, we believe that these genes are directly or indirectly involved in the pollen abortion of SW-S. Our study provides insight into key genes related to sterility traits in S. miltiorrhiza, and the results can be further exploited in functional and mechanism studies.

Insight into the biology of Mycobacterium mucogenicum and Mycobacterium neoaurum clade members

Nontuberculous mycobacteria, NTM, are of growing concern and among these members of the Mycobacterium mucogenicum (Mmuc) and Mycobacterium neoaurum (Mneo) clades can cause infections in humans and they are resistant to first-line anti-tuberculosis drugs. They can be isolated from different ecological niches such as soil, tap water and ground water. Mycobacteria, such as Mmuc and Mneo, are classified as rapid growing mycobacteria, RGM, while the most familiar, Mycobacterium tuberculosis, belongs to the slow growing mycobacteria, SGM. Modern “omics” approaches have provided new insights into our understanding of the biology and evolution of this group of bacteria. Here we present comparative genomics data for seventeen NTM of which sixteen belong to the Mmuc- and Mneo-clades. Focusing on virulence genes, including genes encoding sigma/anti-sigma factors, serine threonine protein kinases (STPK), type VII (ESX genes) secretion systems and mammalian cell entry (Mce) factors we provide insight into their presence as well as phylogenetic relationship in the case of the sigma/anti-sigma factors and STPKs. Our data further suggest that these NTM lack ESX-5 and Mce2 genes, which are known to affect virulence. In this context, Mmuc- and Mneo-clade members lack several of the genes in the glycopeptidolipid (GLP) locus, which have roles in colony morphotype appearance and virulence. For the M. mucogenicum type strain, Mmuc^T, we provide RNASeq data focusing on mRNA levels for sigma factors, STPK, ESX proteins and Mce proteins. These data are discussed and compared to in particular the SGM and fish pathogen Mycobacterium marinum. Finally, we provide insight into as to why members of the Mmuc- and Mneo-clades show resistance to rifampin and isoniazid, and why Mmuc^T forms a rough colony morphotype.

Comparative Genome Analysis of 2 Mycobacterium Tuberculosis Strains from Pakistan: Insights Globally Into Drug Resistance, Virulence, and Niche Adaptation

Multidrug-resistant Mycobacterium tuberculosis is a global threat particularly in developing countries like Pakistan. In this study, we identified 2 M tuberculosis strains, mnpk and swlpk, by 16S RNA genes, sequenced their draft genome, and compared the 2 genomes with reference strain H37Rv and gene expression analysis of selected virulent genes. Phylogenetic analysis of M tuberculosis strains, mnpk and swlpk, using 16S RNA genes revealed that the strains are closely related with reference strain H37Rv. The draft genome sequence of mnpk and swlpk contains 4305 and 4295 protein-coding genes, respectively, having 99.9% with high collinearity when compared with H37Rv. Although some important drug-resistant genes such as fabG, faDE24, and iniA were missing, genome mining also revealed key drug-resistant genes such as katG, inhA, rpoA, rpoB, and rpoC against first-line isoniazid and rifampicin drug. The strain mnpk and swlpk encodes 257 putative and 86 verified virulent genes including type 7 secretion system (T7SS) key genes. The variation in the expression profile of selected T7SS genes, particularly low expression level of EspK, raised concern that the mechanism of virulence of mnpk and swlpk might be different from H37Rv strains as espK is associated with ATPase EccC1a and EccC1b which showed high expression level. Briefly, this study shows that the strains mnpk and swlpk are linked with H37Rv having 99% similarity in genomes, but the absence of drug-resistant genes and variation in key genes’ expression profile espK, EccE1, PPE41, and espC provide a rationale for the future investigation of M tuberculosis mnpk and swlpk pathogenesis via RNA sequencing, single-nucleotide polymorphisms, as well as gene manipulation.

When Dicty Met Myco, a (Not So) Romantic Story about One Amoeba and Its Intracellular Pathogen

In recent years, Dictyostelium discoideum has become an important model organism to study the cell biology of professional phagocytes. This amoeba not only shares many molecular features with mammalian macrophages, but most of its fundamental signal transduction pathways are conserved in humans. The broad range of existing genetic and biochemical tools, together with its suitability for cell culture and live microscopy, make D. discoideum an ideal and versatile laboratory organism. In this review, we focus on the use of D. discoideum as a phagocyte model for the study of mycobacterial infections, in particular Mycobacterium marinum. We look in detail at the intracellular cycle of M. marinum, from its uptake by D. discoideum to its active or passive egress into the extracellular medium. In addition, we describe the molecular mechanisms that both the mycobacterial invader and the amoeboid host have developed to fight against each other, and compare and contrast with those developed by mammalian phagocytes. Finally, we introduce the methods and specific tools that have been used so far to monitor the D. discoideum-M. marinum interaction.

Transcriptome Landscape of Mycobacterium smegmatis

The non-pathogenic bacterium Mycobacterium smegmatis mc²155 has been widely used as a model organism in mycobacterial research, yet a detailed study about its transcription landscape remains to be established. Here we report the transcriptome, expression profiles and transcriptional structures through growth-phase-dependent RNA sequencing (RNA-seq) as well as other related experiments. We found: (1) 2,139 transcriptional start sites (TSSs) in the genome-wide scale, of which eight samples were randomly selected and further verified by 5′-RACE; (2) 2,233 independent monocistronic or polycistronic mRNAs in the transcriptome within the operon/sub-operon structures which are classified into five groups; (3) 47.50% (1016/2139) genes were transcribed into leaderless mRNAs, with the TSSs of 41.3% (883/2139) mRNAs overlapping with the first base of the annotated start codon. Initial amino acids of MSMEG_4921 and MSMEG_6422 proteins were identified by Edman degradation, indicating the presence of distinctive widespread leaderless features in M. smegmatis mc²155. (4) 150 genes with potentially wrong structural annotation, of which 124 proposed genes have been corrected; (5) eight highly active promoters, with their activities further determined by β-galactosidase assays. These data integrated the transcriptional landscape to genome information of model organism mc²155 and lay a solid foundation for further works in Mycobacterium.

Transcriptome Analysis of the Intracellular Facultative Pathogen Piscirickettsia salmonis: Expression of Putative Groups of Genes Associated with Virulence and Iron Metabolism

The intracellular facultative bacteria Piscirickettsia salmonis is one of the most important pathogens of the Chilean aquaculture. However, there is a lack of information regarding the whole genomic transcriptional response according to different extracellular environments. We used next generation sequencing (NGS) of RNA (RNA-seq) to study the whole transcriptome of an isolate of P. salmonis (FAVET-INBIOGEN) using a cell line culture and a modified cell-free liquid medium, with or without iron supplementation. This was done in order to obtain information about the factors there are involved in virulence and iron acquisition. First, the isolate was grown in the Sf21 cell line; then, the bacteria were cultured into a cell-free liquid medium supplemented or not with iron. We identified in the transcriptome, genes associated with type IV secretion systems, genes related to flagellar structure assembly, several proteases and sigma factors, and genes related to the development of drug resistance. Additionally, we identified for the first time several iron-metabolism associated genes including at least two iron uptake pathways (ferrous iron and ferric iron uptake) that are actually expressed in the different conditions analyzed. We further describe putative genes that are related with the use and storage of iron in the bacteria, which have not been previously described. Several sets of genes related to virulence were expressed in both the cell line and cell-free culture media (for example those related to flagellar structure; such as basal body, MS-ring, C-ring, proximal and distal rod, and filament), which may play roles in other basic processes rather than been restricted to virulence.

DgcA, a diguanylate cyclase from Xanthomonas oryzae pv. oryzae regulates bacterial pathogenicity on rice

Xanthomonas oryzae pv. oryzae (Xoo) is the causal agent of rice blight disease as well as a serious phytopathogen worldwide. It is also one of the model organisms for studying bacteria-plant interactions. Current progress in bacterial signal transduction pathways has identified cyclic di-GMP as a major second messenger molecule in controlling Xanthomonas pathogenicity. However, it still remains largely unclear how c-di-GMP regulates the secretion of bacterial virulence factors in Xoo. In this study, we focused on the important roles played by DgcA (XOO3988), one of our previously identified diguanylate cyclases in Xoo, through further investigating the phenotypes of several dgcA-related mutants, namely, the dgcA-knockout mutant ΔdgcA, the dgcA overexpression strain OdgcA, the dgcA complemented strain CdgcA and the wild-type strain. The results showed that dgcA negatively affected virulence, EPS production, bacterial autoaggregation and motility, but positively triggered biofilm formation via modulating the intracellular c-di-GMP levels. RNA-seq data further identified 349 differentially expressed genes controlled by DgcA, providing a foundation for a more solid understanding of the signal transduction pathways in Xoo. Collectively, the present study highlights DgcA as a major regulator of Xoo virulence, and can serve as a potential target for preventing rice blight diseases.

REMap: Operon map of M. tuberculosis based on RNA sequence data

A map of the transcriptional organization of genes of an organism is a basic tool that is necessary to understand and facilitate a more accurate genetic manipulation of the organism. Operon maps are largely generated by computational prediction programs that rely on gene conservation and genome architecture and may not be physiologically relevant. With the widespread use of RNA sequencing (RNAseq), the prediction of operons based on actual transcriptome sequencing rather than computational genomics alone is much needed. Here, we report a validated operon map of Mycobacterium tuberculosis, developed using RNAseq data from both the exponential and stationary phases of growth. At least 58.4% of M. tuberculosis genes are organized into 749 operons. Our prediction algorithm, REMap (RNA Expression Mapping of operons), considers the many cases of transcription coverage of intergenic regions, and avoids dependencies on functional annotation and arbitrary assumptions about gene structure. As a result, we demonstrate that REMap is able to more accurately predict operons, especially those that contain long intergenic regions or functionally unrelated genes, than previous operon prediction programs. The REMap algorithm is publicly available as a user-friendly tool that can be readily modified to predict operons in other bacteria.

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

Comparative Sigma Factor-mRNA Levels in Mycobacterium marinum under Stress Conditions and during Host Infection

We have used RNASeq and qRT-PCR to study mRNA levels for all σ-factors in different Mycobacterium marinum strains under various growth and stress conditions. We also studied their levels in M. marinum from infected fish and mosquito larvae. The annotated σ-factors were expressed and transcripts varied in relation to growth and stress conditions. Some were highly abundant such as sigA, sigB, sigC, sigD, sigE and sigH while others were not. The σ-factor mRNA profiles were similar after heat stress, during infection of fish and mosquito larvae. The similarity also applies to some of the known heat shock genes such as the α-crystallin gene. Therefore, it seems probable that the physiological state of M. marinum is similar when exposed to these different conditions. Moreover, the mosquito larvae data suggest that this is the state that the fish encounter when infected, at least with respect to σ-factor mRNA levels. Comparative genomic analysis of σ-factor gene localizations in three M. marinum strains and Mycobacterium tuberculosis H37Rv revealed chromosomal rearrangements that changed the localization of especially sigA, sigB, sigD, sigE, sigF and sigJ after the divergence of these two species. This may explain the variation in species-specific expression upon exposure to different growth conditions.

Phylogenetic analysis of vitamin B12-related metabolism in Mycobacterium tuberculosis

Comparison of genome sequences from clinical isolates of Mycobacterium tuberculosis with phylogenetically-related pathogens Mycobacterium marinum, Mycobacterium kansasii, and Mycobacterium leprae reveals diversity amongst genes associated with vitamin B12-related metabolism. Diversity is generated by gene deletion events, differential acquisition of genes by horizontal transfer, and single nucleotide polymorphisms (SNPs) with predicted impact on protein function and transcriptional regulation. Differences in the B12 synthesis pathway, methionine biosynthesis, fatty acid catabolism, and DNA repair and replication are consistent with adaptations to different environmental niches and pathogenic lifestyles. While there is no evidence of further gene acquisition during expansion of the M. tuberculosis complex, the emergence of other forms of genetic diversity provides insights into continuing host-pathogen co-evolution and has the potential to identify novel targets for disease intervention.

Small RNAs in mycobacteria: an unfolding story

Mycobacteria represent a class of powerful pathogens, including those causing tuberculosis and leprosy, which continue to be worldwide health challenges. In the last 20 years, an abundance of non-coding, small RNAs (sRNAs) have been discovered in model bacteria and gained significant attention as regulators of cellular responses, including pathogenesis. Naturally, a search in mycobacteria followed, revealing over 200 sRNAs thus far. Characterization of these sRNAs is only beginning, but differential expression under environmental stresses suggests relevance to mycobacterial pathogenesis. This review provides a comprehensive overview of the current knowledge of sRNAs in mycobacteria, including historical perspective and techniques used for identification and characterization.