Tatiana A. Semashko

Transcriptome analysis reveals novel regulatory mechanisms in a genome-reduced bacterium

The avian bacterial pathogen Mycoplasma gallisepticum is a good model for systems studies due to small genome and simplicity of regulatory pathways. In this study, we used RNA-Seq and MS-based proteomics to accurately map coding sequences, transcription start sites (TSSs) and transcript 3′-ends (T3Es). We used obtained data to investigate roles of TSSs and T3Es in stress-induced transcriptional responses. We identified 1061 TSSs at a false discovery rate of 10% and showed that almost all transcription in M. gallisepticum is initiated from classic TATAAT promoters surrounded by A/T-rich sequences. Our analysis revealed the pronounced operon structure complexity: on average, each coding operon has one internal TSS and T3Es in addition to the primary ones. Our transcriptomic approach based on the intervals between the two nearest transcript ends allowed us to identify two classes of T3Es: strong, unregulated, hairpin-containing T3Es and weak, heat shock-regulated, hairpinless T3Es. Comparing gene expression levels under different conditions revealed widespread and divergent transcription regulation in M. gallisepticum. Modeling suggested that the core promoter structure plays an important role in gene expression regulation. We have shown that the heat stress activation of cryptic promoters combined with the hairpinless T3Es suppression leads to widespread, seemingly non-functional transcription.

Specific pools of endogenous peptides are present in gametophore, protonema, and protoplast cells of the moss Physcomitrella patens.

BackgroundProtein degradation is a basic cell process that operates in general protein turnover or to produce bioactive peptides. However, very little is known about the qualitative and quantitative composition of a plant cell peptidome, the actual result of this degradation. In this study we comprehensively analyzed a plant cell peptidome and systematically analyzed the peptide generation process.ResultsWe thoroughly analyzed native peptide pools of Physcomitrella patens moss in two developmental stages as well as in protoplasts. Peptidomic analysis was supplemented by transcriptional profiling and quantitative analysis of precursor proteins. In total, over 20,000 unique endogenous peptides, ranging in size from 5 to 78 amino acid residues, were identified. We showed that in both the protonema and protoplast states, plastid proteins served as the main source of peptides and that their major fraction formed outside of chloroplasts. However, in general, the composition of peptide pools was very different between these cell types. In gametophores, stress-related proteins, e.g., late embryogenesis abundant proteins, were among the most productive precursors. The Driselase-mediated protonema conversion to protoplasts led to a peptide generation “burst”, with a several-fold increase in the number of components in the latter. Degradation of plastid proteins in protoplasts was accompanied by suppression of photosynthetic activity.ConclusionWe suggest that peptide pools in plant cells are not merely a product of waste protein degradation, but may serve as important functional components for plant metabolism. We assume that the peptide “burst” is a form of biotic stress response that might produce peptides with antimicrobial activity from originally functional proteins. Potential functions of peptides in different developmental stages are discussed.

RNA-Seq gene expression profiling of HepG2 cells: the influence of experimental factors and comparison with liver tissue.

BackgroundHuman hepatoma HepG2 cells are used as an in vitro model of the human liver. High-throughput transcriptomic sequencing is an advanced approach for assessing the functional state of a tissue or cell type. However, the influence of experimental factors, such as the sample preparation method and inter-laboratory variation, on the transcriptomic profile has not been evaluated.ResultsThe whole-transcriptome sequencing of HepG2 cells was performed using the SOLiD platform and validated using droplet digital PCR. The gene expression profile was compared to the results obtained with the same sequencing method in another laboratory and using another sample preparation method. We also compared the transcriptomic profile HepG2 cells with that of liver tissue. Comparison of the gene expression profiles between the HepG2 cell line and liver tissue revealed the highest variation, followed by HepG2 cells submitted to two different sample preparation protocols. The lowest variation was observed between HepG2 cells prepared by two different laboratories using the same protocol. The enrichment analysis of the genes that were differentially expressed between HepG2 cells and liver tissue mainly revealed the cancer-associated gene signature of HepG2 cells and the activation of the response to chemical stimuli in the liver tissue. The HepG2 transcriptome obtained with the SOLiD platform was highly correlated with the published transcriptome obtained with the Illumina and Helicos platforms, with moderate correspondence to microarrays.ConclusionsIn the present study, we assessed the influence of experimental factors on the HepG2 transcriptome and identified differences in gene expression between the HepG2 cell line and liver cells. These findings will facilitate robust experimental design in the fields of pharmacology and toxicology. Our results were supported by a comparative analysis with previous HepG2 gene expression studies.

Phase Transition of the Bacterium upon Invasion of a Host Cell as a Mechanism of Adaptation: a Mycoplasma gallisepticum Model

What strategies do bacteria employ for adaptation to their hosts and are these strategies different for varied hosts? To date, many studies on the interaction of the bacterium and its host have been published. However, global changes in the bacterial cell in the process of invasion and persistence, remain poorly understood. In this study, we demonstrated phase transition of the avian pathogen Mycoplasma gallisepticum upon invasion of the various types of eukaryotic cells (human, chicken, and mouse) which was stable during several passages after isolation of intracellular clones and recultivation in a culture medium. It was shown that this phase transition is manifested in changes at the proteomic, genomic and metabolomic levels. Eukaryotic cells induced similar proteome reorganization of M. gallisepticum during infection, despite different origins of the host cell lines. Proteomic changes affected a broad range of processes including metabolism, translation and oxidative stress response. We determined that the activation of glycerol utilization, overproduction of hydrogen peroxide and the upregulation of the SpxA regulatory protein occurred during intracellular infection. We propose SpxA as an important regulator for the adaptation of M. gallisepticum to an intracellular environment.

Draft Genome of the Nitrogen-Fixing Bacterium Pseudomonas stutzeri Strain KOS6 Isolated from Industrial Hydrocarbon Sludge

ABSTRACT Here we present a draft genome of Pseudomonas stutzeri strain KOS6. This strain was isolated from industrial hydrocarbon sludge as a diazotrophic microorganism. It represents one of the major parts of the culturable community of the waste and has potential importance for phytoremediation technology.

Application of sorting and next generation sequencing to study 5′-UTR influence on translation efficiency in Escherichia coli

Abstract Yield of protein per translated mRNA may vary by four orders of magnitude. Many studies analyzed the influence of mRNA features on the translation yield. However, a detailed understanding of how mRNA sequence determines its propensity to be translated is still missing. Here, we constructed a set of reporter plasmid libraries encoding CER fluorescent protein preceded by randomized 5΄ untranslated regions (5΄-UTR) and Red fluorescent protein (RFP) used as an internal control. Each library was transformed into Escherchia coli cells, separated by efficiency of CER mRNA translation by a cell sorter and subjected to next generation sequencing. We tested efficiency of translation of the CER gene preceded by each of 48 natural 5΄-UTR sequences and introduced random and designed mutations into natural and artificially selected 5΄-UTRs. Several distinct properties could be ascribed to a group of 5΄-UTRs most efficient in translation. In addition to known ones, several previously unrecognized features that contribute to the translation enhancement were found, such as low proportion of cytidine residues, multiple SD sequences and AG repeats. The latter could be identified as translation enhancer, albeit less efficient than SD sequence in several natural 5΄-UTRs.

Genome analysis of E. coli isolated from Crohn’s disease patients

BackgroundEscherichia coli (E. coli) has been increasingly implicated in the pathogenesis of Crohn’s disease (CD). The phylogeny of E. coli isolated from Crohn’s disease patients (CDEC) was controversial, and while genotyping results suggested heterogeneity, the sequenced strains of E. coli from CD patients were closely related.ResultsWe performed the shotgun genome sequencing of 28 E. coli isolates from ten CD patients and compared genomes from these isolates with already published genomes of CD strains and other pathogenic and non-pathogenic strains. CDEC was shown to belong to A, B1, B2 and D phylogenetic groups. The plasmid and several operons from the reference CD-associated E. coli strain LF82 were demonstrated to be more often present in CDEC genomes belonging to different phylogenetic groups than in genomes of commensal strains. The operons include carbon-source induced invasion GimA island, prophage I, iron uptake operons I and II, capsular assembly pathogenetic island IV and propanediol and galactitol utilization operons.ConclusionsOur findings suggest that CDEC are phylogenetically diverse. However, some strains isolated from independent sources possess highly similar chromosome or plasmids. Though no CD-specific genes or functional domains were present in all CD-associated strains, some genes and operons are more often found in the genomes of CDEC than in commensal E. coli. They are principally linked to gut colonization and utilization of propanediol and other sugar alcohols.

Reconstruction of Transcription Control Networks in Mollicutes by High-Throughput Identification of Promoters.

Bacteria of the class Mollicutes have significantly reduced genomes and gene expression control systems. They are also efficient pathogens that can colonize a broad range of hosts including plants and animals. Despite their simplicity, Mollicutes demonstrate complex transcriptional responses to various conditions, which contradicts their reduction in gene expression regulation mechanisms. We analyzed the conservation and distribution of transcription regulators across the 50 Mollicutes species. The majority of the transcription factors regulate transport and metabolism, and there are four transcription factors that demonstrate significant conservation across the analyzed bacteria. These factors include repressors of chaperone HrcA, cell cycle regulator MraZ and two regulators with unclear function from the WhiA and YebC/PmpR families. We then used three representative species of the major clades of Mollicutes (Acholeplasma laidlawii, Spiroplasma melliferum, and Mycoplasma gallisepticum) to perform promoter mapping and activity quantitation. We revealed that Mollicutes evolved towards a promoter architecture simplification that correlates with a diminishing role of transcription regulation and an increase in transcriptional noise. Using the identified operons structure and a comparative genomics approach, we reconstructed the transcription control networks for these three species. The organization of the networks reflects the adaptation of bacteria to specific conditions and hosts.

Transcription profiling data set of different states of Mycoplasma gallisepticum.

Mycoplasma gallisepticum belongs to class Mollicutes and causes chronic respiratory disease in birds. It has a reduced genome, lack of cell wall and many metabolic pathways, and also easy to culture and non-pathogenic to humans. Aforementioned made it is a convenient model for studying of systems biology of minimal cell. Studying the transcriptomic level of M. gallisepticum is interesting for both understanding of common principles of transcription regulation of minimal cell and response to definite influence for pathogen bacteria. For rapid investigation of gene expression we developed microarray design including 3366 probes for 678 genes. They included 665 protein coding sequences and 13 antisense RNAs from 816 genes and 17 ncRNAs present in Mycoplasma gallisepticum. The study was performed on Agilent one-color microarray with custom design and random-T7 polymerase primer for cDNA synthesis. Here we present the data for transcription profiling of M. gallisepticum under different types of exposures: genetic knock-out mutants, cell culture exposed to sublethal concentrations of antibiotics and well-characterized heat stress effect. Mutants have transposon insertion to hypothetical membrane protein, lactate dehydrogenase, helicase with unknown function, 1-deoxy-d-xylulose 5-phosphate reductoisomerase or potential sigma factor. For inhibition of important cell systems, treatment with carbonyl cyanide m-chlorophenylhydrazone (CCCP), novobiocin or tetracycline were chosen. Data are available via NCBI Gene Expression Omnibus (GEO) with the accession number GSE85777 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE85777)

Mycoplasma stress response: adaptive regulation or broken brakes?

The avian bacterial pathogen Mycoplasma gallisepticum is a good model for transcriptional regulation studies due to its small genome and relative simplicity. In this study, we used RNA-Seq experiments combined with MS-based proteomics to accurately map coding sequences (CDSs), transcription start sites (TSSs) and transcription terminators (TTs) and to decipher their roles in stress-induced transcriptional responses. We identified 1061 TSSs at an FDR (false discovery rate) of 10% and showed that almost all transcription in M. gallisepticum is initiated from classic TATAAT promoters, which are surrounded by A/T-rich sequences and rarely accompanied by a −35 element. Our analysis revealed the pronounced complexity of the operon structure: on average, each coding operon has one internal TSS and TT in addition to the primary ones. Our new transcriptomic approach based on the intervals between the two closest transcription initiators and/or terminators allowed us to identify two classes of TTs: strong, unregulated and hairpin-containing TTs and weak, heat shock-regulated and hairpinless TTs. Comparing the gene expression levels under different conditions (such as heat, osmotic and peroxide stresses) revealed widespread and divergent transcription regulation in M. gallisepticum. Modeling suggested that the structure of the core promoter plays a major role in gene expression regulation. We have shown that the heat stress activation of cryptic promoters combined with the suppression of hairpinless TTs leads to widespread, seemingly non-functional transcription.