Prevention of High Glucose-Mediated EMT by Inhibition of Hsp70 Chaperone.

Hyperglycemia may contribute to the progression of carcinomas by triggering epithelial-to-mesenchymal transition (EMT). Some proteostasis systems are involved in metastasis; in this paper, we sought to explore the mechanism of Hsp70 chaperone in EMT. We showed that knockdown of Hsp70 reduced cell migration capacity concomitantly with levels of mRNA of the Slug, Snail, and Twist markers of EMT, in colon cancer cells incubated in high glucose medium. Conversely, treatment of cells with Hsp70 inducer U-133 were found to elevate cell motility, along with the other EMT markers. To prove that inhibiting Hsp70 may reduce EMT efficiency, we treated cells with a CL-43 inhibitor of the HSF1 transcription factor, which lowered Hsp70 and HSF1 content in the control and induced EMT in carcinoma cells. Importantly, CL-43 reduced migration capacity, EMT-linked transcription factors, and increased content of epithelial marker E-cadherin in colon cancer cells of three lines, including one derived from a clinical sample. To prove that Hsp70 chaperone should be targeted when inhibiting the EMT pathway, we treated cancer cells with 2-phenylethynesulfonamide (PES) and demonstrated that the compound inhibited substrate-binding capacity of Hsp70. Furthermore, PES suppressed EMT features, cell motility, and expression of specific transcription factors. In conclusion, the Hsp70 chaperone machine efficiently protects mechanisms of the EMT, and the safe inhibitors of the chaperone are needed to hamper metastasis at its initial stage.

Development of qPCR platform with probes for quantifying prevalent and biomedically relevant human gut microbial taxa.

Nowadays the advent of innovative high-throughput sequencing allows obtaining high-quality microbiome profiling. However, PCR-based tests are still considered the "golden standard" for many clinical applications. Here, we designed a qPCR-based platform with fluorescent-labeled oligonucleotide probes for assessing human gut microbiome composition. The system allows conducting qualitative and semiquantitative analysis for 12 prokaryotic taxa that are prevalent in the human gut and associated with diseases, diet, age and other factors. The platform was validated by comparing microbiome profile data obtained with two different methods - the platform and high-throughput 16S rRNA sequencing - across 42 stool samples. The test can form the basis for precise and cost-efficient microbiome assay for large-scale surveys including clinical trials with interventions related to diet and disease risks.

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.

ResistoMap-online visualization of human gut microbiota antibiotic resistome.

Abstract: We created ResistoMap­a Web-based interactive visualization of the presence of genetic determinants conferring resistance to antibiotics, biocides and heavy metals in human gut microbiota. ResistoMap displays the data on more than 1500 published gut metagenomes of world populations including both healthy subjects and patients. Multiparameter display filters allow visual assessment of the associations between the meta-data and proportions of resistome. The geographic map navigation layer allows to state hypotheses regarding the global trends of antibiotic resistance and correlates the gut resistome variations with the national clinical guidelines on antibiotics application. Availability and Implementation: ResistoMap was implemented using AngularJS, CoffeeScript, D3.js and TopoJSON. The tool is publicly available at http://resistomap.rcpcm.org. Contact: yarygin@phystech.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

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.

Guidelines to Statistical Analysis of Microbial Composition Data Inferred from Metagenomic Sequencing.

Metagenomics, the application of high-throughput DNA sequencing for surveys of environmental samples, has revolutionized our view on the taxonomic and genetic composition of complex microbial communities. An enormous richness of microbiota keeps unfolding in the context of various fields ranging from biomedicine and food industry to geology. Primary analysis of metagenomic reads allows to infer semi-quantitative data describing the community structure. However, such compositional data possess statistical specific properties that are important to be considered during preprocessing, hypothesis testing and interpreting the results of statistical tests. Failure to account for these specifics may lead to essentially wrong conclusions as a result of the survey. Here we present a researcher introduced to the field of metagenomics with the basic properties of microbial compositional data including statistical power and proposed distribution models, perform a review of the publicly available software tools developed specifically for such data and outline the recommendations for the application of the methods.

Methods for The Metagenomic Data Visualization and Analysis.

Surveys of environmental microbial communities using metagenomic approach produce vast volumes of multidimensional data regarding the phylogenetic and functional composition of the microbiota. Faced with such complex data, a metagenomic researcher needs to select the means for data analysis properly. Data visualization became an indispensable part of the exploratory data analysis and serves a key to the discoveries. While the molecular-genetic analysis of even a single bacterium presents multiple layers of data to be properly displayed and perceived, the studies of microbiota are significantly more challenging. Here we present a review of the state-of-art methods for the visualization of metagenomic data in a multi-level manner: from the methods applicable to an in-depth analysis of a single metagenome to the techniques appropriate for large-scale studies containing hundreds of environmental samples.

MetaFast: fast reference-free graph-based comparison of shotgun metagenomic data.

MOTIVATION: High-throughput metagenomic sequencing has revolutionized our view on the structure and metabolic potential of microbial communities. However, analysis of metagenomic composition is often complicated by the high complexity of the community and the lack of related reference genomic sequences. As a start point for comparative metagenomic analysis, the researchers require efficient means for assessing pairwise similarity of the metagenomes (beta-diversity). A number of approaches were used to address this task, however, most of them have inherent disadvantages that limit their scope of applicability. For instance, the reference-based methods poorly perform on metagenomes from previously unstudied niches, while composition-based methods appear to be too abstract for straightforward interpretation and do not allow to identify the differentially abundant features. RESULTS: We developed MetaFast, an approach that allows to represent a shotgun metagenome from an arbitrary environment as a modified de Bruijn graph consisting of simplified components. For multiple metagenomes, the resulting representation is used to obtain a pairwise similarity matrix. The dimensional structure of the metagenomic components preserved in our algorithm reflects the inherent subspecies-level diversity of microbiota. The method is computationally efficient and especially promising for an analysis of metagenomes from novel environmental niches. AVAILABILITY AND IMPLEMENTATION: Source code and binaries are freely available for download at https://github.com/ctlab/metafast The code is written in Java and is platform independent (tested on Linux and Windows x86_64). CONTACT: ulyantsev@rain.ifmo.ru SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Assessment of k-mer spectrum applicability for metagenomic dissimilarity analysis.

BACKGROUND: A rapidly increasing flow of genomic data requires the development of efficient methods for obtaining its compact representation. Feature extraction facilitates classification, clustering and model analysis for testing and refining biological hypotheses. "Shotgun" metagenome is an analytically challenging type of genomic data - containing sequences of all genes from the totality of a complex microbial community. Recently, researchers started to analyze metagenomes using reference-free methods based on the analysis of oligonucleotides (k-mers) frequency spectrum previously applied to isolated genomes. However, little is known about their correlation with the existing approaches for metagenomic feature extraction, as well as the limits of applicability. Here we evaluated a metagenomic pairwise dissimilarity measure based on short k-mer spectrum using the example of human gut microbiota, a biomedically significant object of study. RESULTS: We developed a method for calculating pairwise dissimilarity (beta-diversity) of "shotgun" metagenomes based on short k-mer spectra (5 ≤ k ≤ 11). The method was validated on simulated metagenomes and further applied to a large collection of human gut metagenomes from the populations of the world (n=281). The k-mer spectrum-based measure was found to behave similarly to one based on mapping to a reference gene catalog, but different from one using a genome catalog. This difference turned out to be associated with a significant presence of viral reads in a number of metagenomes. Simulations showed limited impact of bacterial genetic variability as well as sequencing errors on k-mer spectra. Specific differences between the datasets from individual populations were identified. CONCLUSIONS: Our approach allows rapid estimation of pairwise dissimilarity between metagenomes. Though we applied this technique to gut microbiota, it should be useful for arbitrary metagenomes, even metagenomes with novel microbiota. Dissimilarity measure based on k-mer spectrum provides a wider perspective in comparison with the ones based on the alignment against reference sequence sets. It helps not to miss possible outstanding features of metagenomic composition, particularly related to the presence of an unknown bacteria, virus or eukaryote, as well as to technical artifacts (sample contamination, reads of non-biological origin, etc.) at the early stages of bioinformatic analysis. Our method is complementary to reference-based approaches and can be easily integrated into metagenomic analysis pipelines.

Large scale analysis of amino acid substitutions in bacterial proteomics.

BACKGROUND: Proteomics of bacterial pathogens is a developing field exploring microbial physiology, gene expression and the complex interactions between bacteria and their hosts. One of the complications in proteomic approach is micro- and macro-heterogeneity of bacterial species, which makes it impossible to build a comprehensive database of bacterial genomes for identification, while most of the existing algorithms rely largely on genomic data. RESULTS: Here we present a large scale study of identification of single amino acid polymorphisms between bacterial strains. An ad hoc method was developed based on MS/MS spectra comparison without the support of a genomic database. Whole-genome sequencing was used to validate the accuracy of polymorphism detection. Several approaches presented earlier to the proteomics community as useful for polymorphism detection were tested on isolates of Helicobacter pylori, Neisseria gonorrhoeae and Escherichia coli. CONCLUSION: The developed method represents a perspective approach in the field of bacterial proteomics allowing to identify hundreds of peptides with novel SAPs from a single proteome.

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.

Kapitza Resistance between Few-Layer Graphene and Water: Liquid Layering Effects.

The Kapitza resistance (RK) between few-layer graphene (FLG) and water was studied using molecular dynamics simulations. The RK was found to depend on the number of the layers in the FLG though, surprisingly, not on the water block thickness. This distinct size dependence is attributed to the large difference in the phonon mean free path between the FLG and water. Remarkably, RK is strongly dependent on the layering of water adjacent to the FLG, exhibiting an inverse proportionality relationship to the peak density of the first water layer, which is consistent with better acoustic phonon matching between FLG and water. These findings suggest novel ways to engineer the thermal transport properties of solid-liquid interfaces by controlling and regulating the liquid layering at the interface.

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

BACKGROUND: Protein 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. RESULTS: We 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. CONCLUSION: We 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.

Chromosome 18 transcriptoproteome of liver tissue and HepG2 cells and targeted proteome mapping in depleted plasma: update 2013.

We report the results obtained in 2012-2013 by the Russian Consortium for the Chromosome-centric Human Proteome Project (C-HPP). The main scope of this work was the transcriptome profiling of genes on human chromosome 18 (Chr 18), as well as their encoded proteome, from three types of biomaterials: liver tissue, the hepatocellular carcinoma-derived cell line HepG2, and blood plasma. The transcriptome profiling for liver tissue was independently performed using two RNaseq platforms (SOLiD and Illumina) and also by droplet digital PCR (ddPCR) and quantitative RT-PCR. The proteome profiling of Chr 18 was accomplished by quantitatively measuring protein copy numbers in the three types of biomaterial (the lowest protein concentration measured was 10(-13) M) using selected reaction monitoring (SRM). In total, protein copy numbers were estimated for 228 master proteins, including quantitative data on 164 proteins in plasma, 171 in the HepG2 cell line, and 186 in liver tissue. Most proteins were present in plasma at 10(8) copies/µL, while the median abundance was 10(4) and 10(5) protein copies per cell in HepG2 cells and liver tissue, respectively. In summary, for liver tissue and HepG2 cells a "transcriptoproteome" was produced that reflects the relationship between transcript and protein copy numbers of the genes on Chr 18. The quantitative data acquired by RNaseq, PCR, and SRM were uploaded into the "Update_2013" data set of our knowledgebase (www.kb18.ru) and investigated for linear correlations.

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

BACKGROUND: Human 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. RESULTS: The 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. CONCLUSIONS: In 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.

AGTR2 gene polymorphism is associated with muscle fibre composition, athletic status and aerobic performance.

Muscle fibre type is a heritable trait and can partly predict athletic success. It has been proposed that polymorphisms of genes involved in the regulation of muscle fibre characteristics may predispose the muscle precursor cells of a given individual to be predominantly fast or slow. In the present study, we examined the association between 15 candidate gene polymorphisms and muscle fibre type composition of the vastus lateralis muscle in 55 physically active, healthy men. We found that rs11091046 C allele carriers of the angiotensin II type 2 receptor gene (AGTR2; involved in skeletal muscle development, metabolism and circulatory homeostasis) had a significantly higher percentage of slow-twitch fibres than A allele carriers [54.2 (11.1) versus 45.2 (10.2)%; P = 0.003]. These data indicate that 15.2% of the variation in muscle fibre composition of the vastus lateralis muscle can be explained by the AGTR2 genotype. Next, we investigated the frequencies of the AGTR2 alleles in 2178 Caucasian athletes and 1220 control subjects. The frequency of the AGTR2 C allele was significantly higher in male and female endurance athletes compared with power athletes (males, 62.7 versus 51.7%, P = 0.0038; females, 56.6 versus 48.1%, P = 0.0169) and control subjects (males, 62.7 versus 51.0%, P = 0.0006; elite female athletes, 65.1 versus 55.2%, P = 0.0488). Furthermore, the frequency of the AGTR2 A allele was significantly over-represented in female power athletes (51.9%) in comparison to control subjects (44.8%, P = 0.0069). We also found that relative maximal oxygen consumption was significantly greater in male endurance athletes with the AGTR2 C allele compared with AGTR2 A allele carriers [n = 28; 62.3 (4.4) versus 57.4 (6.0) ml min(-1) kg(-1); P = 0.0197]. Taken together, these results demonstrate that the AGTR2 gene C allele is associated with an increased proportion of slow-twitch muscle fibres, endurance athlete status and aerobic performance, while the A allele is associated with a higher percentage of fast-twitch fibres and power-oriented disciplines.

Chromosome 18 transcriptome profiling and targeted proteome mapping in depleted plasma, liver tissue and HepG2 cells.

The final goal of the Russian part of the Chromosome-centric Human Proteome Project (C-HPP) was established as the analysis of the chromosome 18 (Chr 18) protein complement in plasma, liver tissue and HepG2 cells with the sensitivity of 10(-18) M. Using SRM, we have recently targeted 277 Chr 18 proteins in plasma, liver, and HepG2 cells. On the basis of the results of the survey, the SRM assays were drafted for 250 proteins: 41 proteins were found only in the liver tissue, 82 proteins were specifically detected in depleted plasma, and 127 proteins were mapped in both samples. The targeted analysis of HepG2 cells was carried out for 49 proteins; 41 of them were successfully registered using ordinary SRM and 5 additional proteins were registered using a combination of irreversible binding of proteins on CN-Br Sepharose 4B with SRM. Transcriptome profiling of HepG2 cells performed by RNAseq and RT-PCR has shown a significant correlation (r = 0.78) for 42 gene transcripts. A pilot affinity-based interactome analysis was performed for cytochrome b5 using analytical and preparative optical biosensor fishing followed by MS analysis of the fished proteins. All of the data on the proteome complement of the Chr 18 have been integrated into our gene-centric knowledgebase ( www.kb18.ru ).

DNA repair in Mycoplasma gallisepticum.

BACKGROUND: DNA repair is essential for the maintenance of genome stability in all living beings. Genome size as well as the repertoire and abundance of DNA repair components may vary among prokaryotic species. The bacteria of the Mollicutes class feature a small genome size, absence of a cell wall, and a parasitic lifestyle. A small number of genes make Mollicutes a good model for a "minimal cell" concept. RESULTS: In this work we studied the DNA repair system of Mycoplasma gallisepticum on genomic, transcriptional, and proteomic levels. We detected 18 out of 22 members of the DNA repair system on a protein level. We found that abundance of the respective mRNAs is less than one per cell. We studied transcriptional response of DNA repair genes of M. gallisepticum at stress conditions including heat, osmotic, peroxide stresses, tetracycline and ciprofloxacin treatment, stationary phase and heat stress in stationary phase. CONCLUSIONS: Based on comparative genomic study, we determined that the DNA repair system M. gallisepticum includes a sufficient set of proteins to provide a cell with functional nucleotide and base excision repair and mismatch repair. We identified SOS-response in M. gallisepticum on ciprofloxacin, which is a known SOS-inducer, tetracycline and heat stress in the absence of established regulators. Heat stress was found to be the strongest SOS-inducer. We found that upon transition to stationary phase of culture growth transcription of DNA repair genes decreases dramatically. Heat stress does not induce SOS-response in a stationary phase.