Non-human peptides revealed in blood reflect the composition of intestinal microbiota.

BACKGROUND: The previously underestimated effects of commensal gut microbiota on the human body are increasingly being investigated using omics. The discovery of active molecules of interaction between the microbiota and the host may be an important step towards elucidating the mechanisms of symbiosis. RESULTS: Here, we show that in the bloodstream of healthy people, there are over 900 peptides that are fragments of proteins from microorganisms which naturally inhabit human biotopes, including the intestinal microbiota. Absolute quantitation by multiple reaction monitoring has confirmed the presence of bacterial peptides in the blood plasma and serum in the range of approximately 0.1 nM to 1 µM. The abundance of microbiota peptides reaches its maximum about 5 h after a meal. Most of the peptides correlate with the bacterial composition of the small intestine and are likely obtained by hydrolysis of membrane proteins with trypsin, chymotrypsin and pepsin - the main proteases of the gastrointestinal tract. The peptides have physicochemical properties that likely allow them to selectively pass the intestinal mucosal barrier and resist fibrinolysis. CONCLUSIONS: The proposed approach to the identification of microbiota peptides in the blood, after additional validation, may be useful for determining the microbiota composition of hard-to-reach intestinal areas and monitoring the permeability of the intestinal mucosal barrier.

Snapper: high-sensitive detection of methylation motifs based on Oxford Nanopore reads.

MOTIVATION: The Oxford Nanopore technology has a great potential for the analysis of methylated motifs in genomes, including whole-genome methylome profiling. However, we found that there are no methylation motifs detection algorithms, which would be sensitive enough and return deterministic results. Thus, the MEME suit does not extract all Helicobacter pylori methylation sites de novo even using the iterative approach implemented in the most up-to-date methylation analysis tool Nanodisco. RESULTS: We present Snapper, a new highly sensitive approach, to extract methylation motif sequences based on a greedy motif selection algorithm. Snapper does not require manual control during the enrichment process and has enrichment sensitivity higher than MEME coupled with Tombo or Nanodisco instruments that was demonstrated on H.pylori strain J99 studied earlier by the PacBio technology and on four external datasets representing different bacterial species. We used Snapper to characterize the total methylome of a new H.pylori strain A45. At least four methylation sites that have not been described for H.pylori earlier were revealed. We experimentally confirmed the presence of a new CCAG-specific methyltransferase and inferred a gene encoding a new CCAAK-specific methyltransferase. AVAILABILITY AND IMPLEMENTATION: Snapper is implemented using Python and is freely available as a pip package named "snapper-ont." Also, Snapper and the demo dataset are available in Zenodo (10.5281/zenodo.10117651).

Phigaro: high-throughput prophage sequence annotation.

SUMMARY: Phigaro is a standalone command-line application that is able to detect prophage regions taking raw genome and metagenome assemblies as an input. It also produces dynamic annotated 'prophage genome maps' and marks possible transposon insertion spots inside prophages. It is applicable for mining prophage regions from large metagenomic datasets. AVAILABILITY AND IMPLEMENTATION: Source code for Phigaro is freely available for download at https://github.com/bobeobibo/phigaro along with test data. The code is written in Python. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Correction to: Draft genome sequences of Hirudo medicinalis and salivary transcriptome of three closely related medicinal leeches.

An amendment to this paper has been published and can be accessed via the original article.

Draft genome sequences of Hirudo medicinalis and salivary transcriptome of three closely related medicinal leeches.

BACKGROUND: Salivary cell secretion (SCS) plays a critical role in blood feeding by medicinal leeches, making them of use for certain medical purposes even today. RESULTS: We annotated the Hirudo medicinalis genome and performed RNA-seq on salivary cells isolated from three closely related leech species, H. medicinalis, Hirudo orientalis, and Hirudo verbana. Differential expression analysis verified by proteomics identified salivary cell-specific gene expression, many of which encode previously unknown salivary components. However, the genes encoding known anticoagulants have been found to be expressed not only in salivary cells. The function-related analysis of the unique salivary cell genes enabled an update of the concept of interactions between salivary proteins and components of haemostasis. CONCLUSIONS: Here we report a genome draft of Hirudo medicinalis and describe identification of novel salivary proteins and new homologs of genes encoding known anticoagulants in transcriptomes of three medicinal leech species. Our data provide new insights in genetics of blood-feeding lifestyle in leeches.

New Insights into Therapy-Induced Progression of Cancer.

The malignant tumor is a complex heterogeneous set of cells functioning in a no less heterogeneous microenvironment. Like any dynamic system, cancerous tumors evolve and undergo changes in response to external influences, including therapy. Initially, most tumors are susceptible to treatment. However, remaining cancer cells may rapidly reestablish the tumor after a temporary remission. These new populations of malignant cells usually have increased resistance not only to the first-line agent, but also to the second- and third-line drugs, leading to a significant decrease in patient survival. Multiple studies describe the mechanism of acquired therapy resistance. In past decades, it became clear that, in addition to the simple selection of pre-existing resistant clones, therapy induces a highly complicated and tightly regulated molecular response that allows tumors to adapt to current and even subsequent therapeutic interventions. This review summarizes mechanisms of acquired resistance, such as secondary genetic alterations, impaired function of drug transporters, and autophagy. Moreover, we describe less obvious molecular aspects of therapy resistance in cancers, including epithelial-to-mesenchymal transition, cell cycle alterations, and the role of intercellular communication. Understanding these molecular mechanisms will be beneficial in finding novel therapeutic approaches for cancer therapy.

Challenges of the Human Proteome Project: 10-Year Experience of the Russian Consortium.

This manuscript collects all the efforts of the Russian Consortium, bottlenecks revealed in the course of the C-HPP realization, and ways of their overcoming. One of the main bottlenecks in the C-HPP is the insufficient sensitivity of proteomic technologies, hampering the detection of low- and ultralow-copy number proteins forming the "dark part" of the human proteome. In the frame of MP-Challenge, to increase proteome coverage we suggest an experimental workflow based on a combination of shotgun technology and selected reaction monitoring with two-dimensional alkaline fractionation. Further, to detect proteins that cannot be identified by such technologies, nanotechnologies such as combined atomic force microscopy with molecular fishing and/or nanowire detection may be useful. These technologies provide a powerful tool for single molecule analysis, by analogy with nanopore sequencing during genome analysis. To systematically analyze the functional features of some proteins (CP50 Challenge), we created a mathematical model that predicts the number of proteins differing in amino acid sequence: proteoforms. According to our data, we should expect about 100 000 different proteoforms in the liver tissue and a little more in the HepG2 cell line. The variety of proteins forming the whole human proteome significantly exceeds these results due to post-translational modifications (PTMs). As PTMs determine the functional specificity of the protein, we propose using a combination of gene-centric transcriptome-proteomic analysis with preliminary fractionation by two-dimensional electrophoresis to identify chemically modified proteoforms. Despite the complexity of the proposed solutions, such integrative approaches could be fruitful for MP50 and CP50 Challenges in the framework of the C-HPP.

Long-term impact of fecal transplantation in healthy volunteers.

BACKGROUND: Fecal microbiota transplantation (FMT) has been recently approved by FDA for the treatment of refractory recurrent clostridial colitis (rCDI). Success of FTM in treatment of rCDI led to a number of studies investigating the effectiveness of its application in the other gastrointestinal diseases. However, in the majority of studies the effects of FMT were evaluated on the patients with initially altered microbiota. The aim of our study was to estimate effects of FMT on the gut microbiota composition in healthy volunteers and to monitor its long-term outcomes. RESULTS: We have performed a combined analysis of three healthy volunteers before and after capsule FMT by evaluating their general condition, adverse clinical effects, changes of basic laboratory parameters, and several immune markers. Intestinal microbiota samples were evaluated by 16S rRNA gene and shotgun sequencing. The data analysis demonstrated profound shift towards the donor microbiota taxonomic composition in all volunteers. Following FMT, all the volunteers exhibited gut colonization with donor gut bacteria and persistence of this effect for almost ∼1 year of observation. Transient changes of immune parameters were consistent with suppression of T-cell cytotoxicity. FMT was well tolerated with mild gastrointestinal adverse events, however, one volunteer developed a systemic inflammatory response syndrome. CONCLUSIONS: The FMT leads to significant long-term changes of the gut microbiota in healthy volunteers with the shift towards donor microbiota composition and represents a relatively safe procedure to the recipients without long-term adverse events.

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

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.

The A-allele of the FTO Gene rs9939609 Polymorphism Is Associated With Decreased Proportion of Slow Oxidative Muscle Fibers and Over-represented in Heavier Athletes.

Guilherme, JPLF, Egorova, ES, Semenova, EA, Kostryukova, ES, Kulemin, NA, Borisov, OV, Khabibova, SA, Larin, AK, Ospanova, EA, Pavlenko, AV, Lyubaeva, EV, Popov, DV, Lysenko, EA, Vepkhvadze, TF, Lednev, EM, Govorun, VM, Generozov, EV, Ahmetov, II, and Lancha Junior, AH. The A-allele of the FTO gene rs9939609 polymorphism is associated with decreased proportion of slow oxidative muscle fibers and over-represented in heavier athletes. J Strength Cond Res 33(3): 691-700, 2019-The purpose of this study was to explore the frequency of the FTO T > A (rs9939609) polymorphism in elite athletes from 2 cohorts (Brazil and Russia), as well as to find a relationship between FTO genotypes and muscle fiber composition. A total of 677 athletes and 652 nonathletes were evaluated in the Brazilian cohort, whereas a total of 920 athletes and 754 nonathletes were evaluated in the Russian cohort. It was found a trend for a lower frequency of A/A genotype in long-distance athletes compared with nonathletes (odds ratio [OR]: 0.65; p = 0.054). By contrast, it was found an increased frequency of the A-allele in Russian power athletes. The presence of the T/A + A/A genotypes rather than T/T increased the OR of being a Russian power athlete compared with matched nonathletes (OR: 1.45; p = 0.002). Different from that observed in combat sports athletes of lighter weight categories, the A-allele was also over-represented in combat sports athletes of heavier weight categories. The presence of the T/A + A/A genotypes rather than T/T increased the OR of being a combat sports athlete of heavier weight categories compared with nonathletes (OR: 1.79; p = 0.018). Regarding the muscle fibers, we found that carriers of the A/A genotype had less slow-twitch muscle fibers than T-allele carriers (p = 0.029). In conclusion, the A/A genotype of the FTO T > A polymorphism is under-represented in athletes more reliant on a lean phenotype and associated with decreased proportion of slow-twitch muscle fibers, while is over-represented in strength and heavier athletes.

Genetic diversity of Escherichia coli in gut microbiota of patients with Crohn's disease discovered using metagenomic and genomic analyses.

BACKGROUND: Crohn's disease is associated with gut dysbiosis. Independent studies have shown an increase in the abundance of certain bacterial species, particularly Escherichia coli with the adherent-invasive pathotype, in the gut. The role of these species in this disease needs to be elucidated. METHODS: We performed a metagenomic study investigating the gut microbiota of patients with Crohn's disease. A metagenomic reconstruction of the consensus genome content of the species was used to assess the genetic variability. RESULTS: The abnormal shifts in the microbial community structures in Crohn's disease were heterogeneous among the patients. The metagenomic data suggested the existence of multiple E. coli strains within individual patients. We discovered that the genetic diversity of the species was high and that only a few samples manifested similarity to the adherent-invasive varieties. The other species demonstrated genetic diversity comparable to that observed in the healthy subjects. Our results were supported by a comparison of the sequenced genomes of isolates from the same microbiota samples and a meta-analysis of published gut metagenomes. CONCLUSIONS: The genomic diversity of Crohn's disease-associated E. coli within and among the patients paves the way towards an understanding of the microbial mechanisms underlying the onset and progression of the Crohn's disease and the development of new strategies for the prevention and treatment of this disease.

ACI-1 beta-lactamase is widespread across human gut microbiomes in Negativicutes due to transposons harboured by tailed prophages.

Antibiotic resistance is increasing among pathogens, and the human microbiome contains a reservoir of antibiotic resistance genes. Acidaminococcus intestini is the first Negativicute bacterium (Gram-negative Firmicute) shown to be resistant to beta-lactam antibiotics. Resistance is conferred by the aci1 gene, but its evolutionary history and prevalence remain obscure. We discovered that ACI-1 proteins are phylogenetically distinct from beta-lactamases of Gram-positive Firmicutes and that aci1 occurs in bacteria scattered across the Negativicute clade, suggesting lateral gene transfer. In the reference A. intestini RyC-MR95 genome, we found transposons residing within a tailed prophage context are likely vehicles for aci1's mobility. We found aci1 in 56 (4.4%) of 1,267 human gut metagenomes, mostly hosted within A. intestini, and, where could be determined, mostly within a consistent mobile element constellation. These samples are from Europe, China and the USA, showing that aci1 is distributed globally. We found that for most Negativicute assemblies with aci1, the prophage observed in A. instestini is absent, but in all cases aci1 is flanked by varying transposons. The chimeric mobile elements we identify here likely have a complex evolutionary history and potentially provide multiple complementary mechanisms for antibiotic resistance gene transfer both within and between cells.

The Pathogenesis of the Demyelinating Form of Guillain-Barre Syndrome (GBS): Proteo-peptidomic and Immunological Profiling of Physiological Fluids.

Acute inflammatory demyelinating polyneuropathy (AIDP) - the main form of Guillain-Barre syndrome-is a rare and severe disorder of the peripheral nervous system with an unknown etiology. One of the hallmarks of the AIDP pathogenesis is a significantly elevated cerebrospinal fluid (CSF) protein level. In this paper CSF peptidome and proteome in AIDP were analyzed and compared with multiple sclerosis and control patients. A total protein concentration increase was shown to be because of even changes in all proteins rather than some specific response, supporting the hypothesis of protein leakage from blood through the blood-nerve barrier. The elevated CSF protein level in AIDP was complemented by activization of protein degradation and much higher peptidome diversity. Because of the studies of the acute motor axonal form, Guillain-Barre syndrome as a whole is thought to be associated with autoimmune response against neurospecific molecules. Thus, in AIDP, autoantibodies against cell adhesion proteins localized at Ranvier's nodes were suggested as possible targets in AIDP. Indeed, AIDP CSF peptidome analysis revealed cell adhesion proteins degradation, however no reliable dependence on the corresponding autoantibodies levels was found. Proteome analysis revealed overrepresentation of Gene Ontology groups related to responses to bacteria and virus infections, which were earlier suggested as possible AIDP triggers. Immunoglobulin blood serum analysis against most common neuronal viruses did not reveal any specific pathogen; however, AIDP patients were more immunopositive in average and often had polyinfections. Cytokine analysis of both AIDP CSF and blood did not show a systemic adaptive immune response or general inflammation, whereas innate immunity cytokines were up-regulated. To supplement the widely-accepted though still unproven autoimmunity-based AIDP mechanism we propose a hypothesis of the primary peripheral nervous system damaging initiated as an innate immunity-associated local inflammation following neurotropic viruses egress, whereas the autoantibody production might be an optional complementary secondary process.

AGTR2 and sprint/power performance: a case-control replication study for rs11091046 polymorphism in two ethnicities.

We aimed to replicate, in a specific athletic event cohort (only track and field) and in two different ethnicities (Japanese and East European, i.e. Russian and Polish), original findings showing the association of the angiotensin-II receptor type-2 gene (AGTR2) rs11091046 A>C polymorphism with athlete status. We compared genotypic frequencies of the AGTR2 rs11091046 polymorphism among 282 track and field sprint/power athletes (200 men and 82 women), including several national record holders and Olympic medallists (214 Japanese, 68 Russian and Polish), and 2024 control subjects (842 men and 1182 women) (804 Japanese, 1220 Russian and Polish). In men, a meta-analysis from the two combined cohorts showed a significantly higher frequency of the C allele in athletes than in controls (odds ratio: 1.62, P=0.008, heterogeneity index I 2 =0%). With regard to respective cohorts, C allele frequency was higher in Japanese male athletes than in controls (67.7% vs. 55.9%, P=0.022), but not in Russian/Polish male athletes (61.9% vs. 51.0%, P=0.172). In women, no significant results were obtained by meta-analysis for the two cohorts combination (P=0.850). The AC genotype frequency was significantly higher in Russian/Polish women athletes than in controls (69.2% vs. 42.1%, P=0.022), but not in Japanese women athletes (P=0.226). Our results, in contrast to previous findings, suggested by meta-analysis that the C allele of the AGTR2 rs11091046 polymorphism is associated with sprint/power track and field athlete status in men, but not in women.

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

BACKGROUND: Escherichia 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. RESULTS: We 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. CONCLUSIONS: Our 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.

Binding site of restriction-modification system controller protein in Mollicutes.

BACKGROUND: Bacteria of the class Mollicutes underwent extreme reduction of genomes and gene expression control systems. Only a few regulators are known to date. In this work, we describe a novel group of transcriptional regulators that are distributed within different Mollicutes and control the expression of restriction-modification systems (RM-systems). RESULTS: We performed cross-species search of putative regulators of RM-systems (C-proteins) and respective binding sites in Mollicutes. We identified a set of novel putative C-protein binding motifs distributed within Mollicutes. We studied the most frequent motif and respective C-protein on the model of Mycoplasma gallisepticum S6. We confirmed our prediction and identified key nucleotides important for C-protein binding. Further we identified novel target promoters of C-protein in M. gallisepticum. CONCLUSIONS: We found that C-protein of M. gallisepticum binds predicted conserved direct repeats of the (GTGTTAN5)2 motif. Apart from its own operon promoter, HsdC can bind to the promoters of the clpB chaperone gene and a tRNA cluster.

State of the Art of Chromosome 18-Centric HPP in 2016: Transcriptome and Proteome Profiling of Liver Tissue and HepG2 Cells.

A gene-centric approach was applied for a large-scale study of expression products of a single chromosome. Transcriptome profiling of liver tissue and HepG2 cell line was independently performed using two RNA-Seq platforms (SOLiD and Illumina) and also by Droplet Digital PCR (ddPCR) and quantitative RT-PCR. Proteome profiling was performed using shotgun LC-MS/MS as well as selected reaction monitoring with stable isotope-labeled standards (SRM/SIS) for liver tissue and HepG2 cells. On the basis of SRM/SIS measurements, protein copy numbers were estimated for the Chromosome 18 (Chr 18) encoded proteins in the selected types of biological material. These values were compared with expression levels of corresponding mRNA. As a result, we obtained information about 158 and 142 transcripts for HepG2 cell line and liver tissue, respectively. SRM/SIS measurements and shotgun LC-MS/MS allowed us to detect 91 Chr 18-encoded proteins in total, while an intersection between the HepG2 cell line and liver tissue proteomes was ∼66%. In total, there were 16 proteins specifically observed in HepG2 cell line, while 15 proteins were found solely in the liver tissue. Comparison between proteome and transcriptome revealed a poor correlation (R2 ≈ 0.1) between corresponding mRNA and protein expression levels. The SRM and shotgun data sets (obtained during 2015-2016) are available in PASSEL (PASS00697) and ProteomeExchange/PRIDE (PXD004407). All measurements were also uploaded into the in-house Chr 18 Knowledgebase at http://kb18.ru/protein/matrix/416126 .

Comprehensive analysis of draft genomes of two closely related pseudomonas syringae phylogroup 2b strains infecting mono- and dicotyledon host plants.

BACKGROUND: In recent years, the damage caused by bacterial pathogens to major crops has been increasing worldwide. Pseudomonas syringae is a widespread bacterial species that infects almost all major crops. Different P. syringae strains use a wide range of biochemical mechanisms, including phytotoxins and effectors of the type III and type IV secretion systems, which determine the specific nature of the pathogen virulence. RESULTS: Strains 1845 (isolated from dicots) and 2507 (isolated from monocots) were selected for sequencing because they specialize on different groups of plants. We compared virulence factors in these and other available genomes of phylogroup 2 to find genes responsible for the specialization of bacteria. We showed that strain 1845 belongs to the clonal group that has been infecting monocots in Russia and USA for a long time (at least 50 years). Strain 1845 has relatively recently changed its host plant to dicots. CONCLUSIONS: The results obtained by comparing the strain 1845 genome with the genomes of bacteria infecting monocots can help to identify the genes that define specific nature of the virulence of P. syringae strains.

Proteome-metabolome profiling of ovarian cancer ascites reveals novel components involved in intercellular communication.

Ovarian cancer ascites is a native medium for cancer cells that allows investigation of their secretome in a natural environment. This medium is of interest as a promising source of potential biomarkers, and also as a medium for cell-cell communication. The aim of this study was to elucidate specific features of the malignant ascites metabolome and proteome. In order to omit components of the systemic response to ascites formation, we compared malignant ascites with cirrhosis ascites. Metabolome analysis revealed 41 components that differed significantly between malignant and cirrhosis ascites. Most of the identified cancer-specific metabolites are known to be important signaling molecules. Proteomic analysis identified 2096 and 1855 proteins in the ovarian cancer and cirrhosis ascites, respectively; 424 proteins were specific for the malignant ascites. Functional analysis of the proteome demonstrated that the major differences between cirrhosis and malignant ascites were observed for the cluster of spliceosomal proteins. Additionally, we demonstrate that several splicing RNAs were exclusively detected in malignant ascites, where they probably existed within protein complexes. This result was confirmed in vitro using an ovarian cancer cell line. Identification of spliceosomal proteins and RNAs in an extracellular medium is of particular interest; the finding suggests that they might play a role in the communication between cancer cells. In addition, malignant ascites contains a high number of exosomes that are known to play an important role in signal transduction. Thus our study reveals the specific features of malignant ascites that are associated with its function as a medium of intercellular communication.

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.