- Research article
- Open Access
DNA repair in Mycoplasma gallisepticum
BMC Genomics volume 14, Article number: 726 (2013)
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.
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.
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.
Genomic DNA is constantly subject to damage. This includes the misincorporation of nucleotides by DNA-polymerases and chemical modification either by endogenous metabolites or external compounds. DNA damage is a frequent event, and cell division is possible because most of the lesions are effectively repaired. The DNA repair mechanisms of bacteria are numerous and universal. The inactivation of single genes in the DNA repair system is generally not lethal because of a significant amount of redundancy among DNA repair pathways .
Bacteria of the class Mollicutes features a small though sufficient genome to grow on a cell-free medium and a low GC content (31% for Mycoplasma gallisepticum). Their genome size typically ranges from 580 thousand to 1.4 million bp. Most mollicute species are obligate parasites [2–4]. As a result, they are constantly exposed to stress conditions caused by host defense, including elevated temperature and reactive oxygen species . It is interesting how organisms with reduced genomes retain genome stability under adverse conditions.
DNA repair systems are divided into groups based on mechanism of action. Mismatch repair (MMR) fix non-pairing bases, which originated from the misincorporation of nucleotides by DNA-polymerase . Chemically modified nucleotides are restored by base excision repair (BER). It recognizes specific modifications by dedicated proteins and replaces modified nucleotides with the correct ones . Nucleotide excision repair (NER) deals with a broad spectrum of DNA damage. It removes a section of the defective DNA chain and fills it based on a complementary chain . DNA lesions that are substrate of the BER system can be restored by the NER system but in a less efficient manner . Double strand breaks or crosslinks are removed by homologous recombination. However, this requires information from homologous DNA molecules (not from a complementary chain, as for BER and NER), which is time consuming and involves massive DNA synthesis [1, 8]. Therefore, recombination activates as a part of an SOS-response only in cases of severe DNA damage . SOS-response in Mollicutes is of particular interest due to absence of LexA repressor in all Mollicutes. An SOS-response may involve a number of repair enzymes, but key players are recombination protein and error-prone polymerases . In bacteria, an SOS-response is activated by massive DNA damage . Activation of an SOS-response inevitably leads to elevation of a mutation rate as well , which is unwanted for a cell in normal conditions. Absence of LexA repressor in Mollicutes was considered evidence that an SOS-response is not functional in these bacteria [2, 10].
Mismatch repair was considered most degenerated in mycoplasmas. Homologs of MutH, MutS and MutL proteins were not found in all Mollicutes. This led to the conclusion that mutation rate is elevated in mycoplasmas ; however, it was later shown that it does not differ from the one in E. coli[15, 16]. It was shown recently that histone-like protein HU (Hup2) in M. gallisepticum is able to recognize mispaired bases in DNA in contrast to HU from Acholeplasma laidlawii with functional MMR [17, 18]. It can also compensate for inactivation of hupAB proteins in E. coli. This may indicate its role in mismatch repair.
In the current study, we attempted to characterize the DNA repair system of Mycopalsma gallisepticum. We identified members of different DNA repair pathways based on genomic data and publications. We quantified mRNA of respective proteins per cell in a set of adverse conditions and identified most of them on the protein level. We showed induction of SOS-response in adverse conditions on a transcriptional level in the absence of a known regulator.
In silico DNA repair system reconstruction
Nucleotide excision repair
Nucleotide excision repair of M. gallisepticum encompasses all required members on the protein level (Tables 1 and 2). This includes UvrA (lesion binding), UvrB (local melt of DNA), UvrC (excision of damaged site), and UvrD-helicase (removal of damaged site), as well as DNA-ligase and DNA-polymerase . The interesting question here is which polymerase is recruited for DNA repair because DNA-polymerase I is absent in M. gallisepticum as well as in most of Mollicutes. It seems that DNA-polymerase III serves both DNA replication and repair. Alternatively, the DNA repair function can be taken by DNA-polymerase IV.
Base excision repair
Base excision repair starts for glycosylases, which recognize specific lesions in DNA. M. gallisepticum have two enzymes of this type, including uracil- (Ung) and formamidopyrimidine-DNA-glycosylases (MutM). The next player is AP-endonuclease, which nicks DNA at 5′ position from the AP-site. The only AP-endonuclease in Mollicutes is Nfo, which belongs to the endonuclease IV family. This family of endonucleases features 3′-phosphatase activity but lacks 5′ -phosphatase and 3′-5′ exonuclease activity in contrast to the endonuclease III family . In well-studied bacteria like E. coli, the AP-site is removed and repaired by DNA-polymerase I (PolA), which has 5′-3′ exonuclease activity. In M. gallisepticum and most of Mollicutes, this mechanism is likely substituted with DNA-polymerase III and Exo protein, which is an exonuclease, homologous to the 5′-3′ exonuclease domain of PolA. Exo seems to originate from PolA as a result of truncation with the loss of the polymerase domain (Figure 1). All components of BER except uracil-glycosylase were identified on the protein level [20–22].
We found two additional proteins that can play a role in mismatch repair. They are MGA_0195 and MGA_0793. The first one contains an endonuclease domain that belongs to the same superfamily as the endonuclease domain of MutH. The second one contains a vsr domain, which takes part in the repair of mispaired guanines. MGA_0195 and MGA_0793 were identified on the protein level (Table 2). A hallmark of M. gallisepticum is the absence of known exonucleases (Exo I, Exo VII, RecJ, Exo X ) of the MMR pathway. However, missing activity can be compensated by exonuclease Exo, a truncated PolA.
We identified one DNA methylation enzyme: HsdM. This is site-specific methylase that methylate adenine in position six.
Recombinational repair system
The genome of M. gallisepticum carries genes that encode key proteins required for recombination - recombinase RecA, RecR, and RecO (MGA_0016), as well as genes ruvA and ruvB, encoding a DNA helicase RuvAB (involved in the migration of the DNA chains), as well as two genes encoding enzymes that allow Holliday junction resolution - DNA resolvase RecU and MGA_0836. In addition, there is a gene smc, encoding Smc-cohesin capable of implementing cohesion after the replication of chromosomes, thus participating in the recombinational repair .
SOS response system
The following participants of SOS-system were found in the genome of M. gallisepticum: recombinase A (recA), recombinase R (recR), helicase complex (ruvA, ruvB), nuclease-helicase complex UvrABC, and DNA-dependent DNA polymerase IV (dinB) which homologous to the mutagenic DNA pol IV of E. coli. In bacteria, this polymerase (along with DNA polymerase V, which is absent in all members of the class Mollicutes) is able to use damaged DNA template . At the same time, we didn’t find any homologs of known regulators of the bacterial SOS-system (LexA , HdiR ) in any of the analyzed genomes of Mollicutes. We found all of the prospective members of the SOS-response system at the mRNA level; while at the protein level, we were unable to identify only helicase RuvAB (Table 2).
Presence of DNA-repair transcripts in the cell
Comparative analysis based only on genomic data does not allow us to draw conclusions about genes expression and the functional activity of the proteins encoded by them. To test the activity of the annotated genes, we carried out a quantitative analysis of transcription levels for the genes of the repair system by quantitative reverse transcription PCR (droplet-digital PCR and real-time PCR - see Methods). We also counted the number of transcripts copies per single bacterial genome. The results of the expression assay for the studied genes at the mRNA and protein levels are shown in Table 1. Transcription was detected for all of analyzed genes.
It should be noted that the number of copies of mRNA varies from one per one hundred (recA, ruvA, ung, dinB) to one per 2.5 copies of genomic DNA (hup2). If we assume that each cell contains, on average, one copy of genomic DNA, the presence of transcripts is rather low, and in most of the cells there are no transcripts of the genes’ encoding repair system proteins.
Presence of DNA-repair proteins in the cell
In order to assess the presence of repair systems members at the protein level, we used the methods of liquid chromatography-mass spectrometry for the comprehensive identification of all the proteins in the cells of M. gallisepticum. We have identified 561 proteins in total (data not shown) by using the following selection criteria: the number of unique peptides for protein - at least two - and a global FDR of 1% (based on the analysis in PSPEP, with threshold unused score for the protein equaled 0.4). Thus 17,221 unique peptides have been identified (global FDR 1% PSPEP). It should be noted that due to the fact that the algorithm Paragon, used for searching, works on the basis of probability factors, the total number of peptides includes not only non-modified tryptic peptides but semi-tryptic and non-tryptic peptides as well, and peptides containing all possible amino acid modifications considered in the default search algorithm of Paragon.
Analysis of the presence of the DNA repair proteins (Table 2) showed the presence of most of them in the mycoplasma’s proteome - 18 proteins for 22 genes are listed in Table 2. All proteins in Table 2 were unique according to the groups of spectrums using the Pro Group algorithm.
Repair system genes response on stress (T, NaCl, H2O2, stationary phase of growth, antibiotics)
According to the literature data, the SOS-repair system of Mollicutes have undergone reduction . In order to understand whether an SOS-response takes place in M. gallisepticum, we decided to test changes in the levels of mRNAs of the genes’ encoding proteins of the repair system in response to stress influences. We exposed the cells to sub-lethal stress (see Methods) - temperature, osmotic and peroxide stress, and antibiotics. Such effects were selected because they are physiologically common for parasites, as the interaction with the host organism cells encounter reactions of inflammation: heat, immune system peroxide attack, and antibiotics treatment. Figure 2A shows the color map, representing the change of the transcription profiles for the genes’ encoding DNA repair proteins. Raw data are presented in Additional file 1.
Ciprofloxacin is known as an SOS-inducer in different bacteria [13, 27]. It was used as a reference condition to identify members of SOS-response in M. gallisepticum. Ciprofloxacin treatment resulted in the upregulation of recA, recR, ung, ruvA, ligA, parC, and dinB. Remarkably, the 15-fold induction of error-prone polymerase (dinB) was the strongest effect of the ciprofloxacin treatment. Tetracycline treatment induced upregulation of the SOS-response participants recA and recR but not dinB. DNA gyrase (gyrA and gyrB), which is a target of ciprofloxacin, was induced by tetracycline but not ciprofloxacin. Osmotic stress induced the upregulation of dinB and ung but not recA and recR. Peroxide stress had little effect on the transcription of genes involved in DNA repair. Heat stress invoked the strongest response among all types of stress. Response to heat stess includes genes involved in responses to ciprofloxacin, uvrABCD, nfo glycosylase, and DNA gyrase (gyrA and gyrB). It may be concluded that the stress response system of M. gallisepticum consists of several regulons.
The stationary growth phase differs remarkably from the exponential phase in terms of gene expression as well as stress response. The stationary phase demonstrates a significant downregulation of most of the studied genes. Heat stress in the stationary phase did not result in the change of expression of any studied gene in contrast to its exponential phase. It may indicate that the transcription process in the stationary phase is modulated by some global regulator, like an alternative sigma factor.
The transcription profile for the repair system genes in M. gallisepticum under heat shock conditions (Figure 2B) suggests three groups of genes. The first group is characterized by a gradual growth of the transcription level and includes the majority of the genes. The second group includes the genes that respond to the heat stress by dramatical increasing of the transcription level following by decreasing or stagnation of the last one. The third group includes two genes of histone-like proteins hup1 and hup2. They are characterized by a steady decline of transcription (Figure 2B).
The study of the mycoplasma repair system, as bacteria with a reduced genome, gives an indication of the minimum number of genes required to maintain genomic stability. Comparative analysis shows that the repair system of M. gallisepticum includes fewer genes than the one of E. coli. However, if we give preference to functional rather than numerical parameters, we can see the presence of the key elements of all major repair systems in the absence of overlapping units.
The results show that the possible composition of the M. gallisepticum repair system may be greater than it was assumed up until now. In particular, we found several previously unknown mycoplasma proteins. One of them - a protein MGA_0793 - was annotated as a hypothetical protein of unknown function. Results of the alignment revealed homology with the Vsr protein, which was involved in the mismatch repair of DNA containing an unpaired guanine in E. coli. Interestingly, the HU (Hup2) protein identified earlier in M. gallisepticum, capable of binding DNA mismatches, does not bind unpaired T-G, A-G, and G-G base pairs . The Figure 3 shows a hypothetical model of the MMR system in M. gallisepticum.
We identified one DNA methylation enzyme: HsdM. This is site-specific methylase that methylate adenine in position six. DNA methylation is required for MMR pathways to distinguish the correct strand to use as a template . It was shown that methylation occurs in mycoplasmas [29–31]. However dam-methyltransferase required for the methylation of GATC sites is absent in M. gallisepticum, the discrimination of the chains may be due to the interaction of the repair complex with the B-subunit of DNA polymerase III in replication - such a mechanism was previously shown for a number of microorganisms (Gram-positive bacteria) . Thus we can at least say that M. gallisepticum has an enzyme capable of recognizing and binding DNA mismatches, as well as an enzyme hypothetically capable of the excision of the damaged DNA fragment. But the question remains, which enzyme of mycoplasma is a functional analog of MutH protein, which is required to make a single-strand break in the DNA being repaired.
We found that M. gallisepticum has the full pathway of nucleotide excision repair and recombination of DNA. From all of DNA recombination participants, only RecR have been previously found in Mycoplasma gallisepticum at the protein level [20–22]. In this study, we identified RecA and Smc proteins as well (Table 2).
Of particular interest for the analysis of gene excision is the base excision repair system, as it has always been represented by a large number of proteins, each of which recognizes a different type of damage. There are three known types of damage: oxidation, alkylation, and deamination. DNA glycosylase MutM in E. coli, whose homolog is present in the genome of M. gallisepticum, corrects the most frequently occurring DNA damage caused by endogenous oxidative stress . Uracil-glycosylase removes uracil, which arises spontaneously by the deamination of cytosine or by erroneous inclusion during replication. Mollicutes have only one of two AP-endonucleases - endonuclease IV. It’s interesting that endonuclease IV (Nfo) has an additional activity - it recognizes oxidized bases (hydroxycytosine, dihydroxyuracil, and dihydroxythymine) and makes a single-stranded gap upstream of the damage, which is used as a primer for the repair by polymerase and ligase . The Figure 3 shows a hypothetical model of the BER system in M. gallisepticum.
Assuming that Mollicutes are bacteria with a minimal genome and are capable of self-reproduction, we come to the conclusion that they possess the number of repair system proteins that are necessary and sufficient for life in the cell-free medium.
The results of the gene transcription assay indicate that the mRNA molecules are not present in every cell in the population (see Table 1). These data are consistent with the literature and may be associated with the long lifetime of the functional protein in contrast to the mRNA . At the proteomic level, we were able to identify a large portion (80%) of repair proteins, including DinB and RecA, members of SOS-response, which is an additional indication of the functionality of the repair system in M. gallisepticum.
An interesting result of the transcription-profiling assay was the induction of SOS-response genes in different shocks. This fact interested us for several reasons. First of all, such a response has not been shown previously in Mollicutes, including large-scale studies of transcription responses in M. pneumoniae, which is one of the closest relatives of M. gallisepticum according to phylogenetic studies . The other interest is that the genomes of all members of the class Mollicutes don’t have any known regulator of the SOS-response system, and therefore several authors consider the SOS-system as non-functional in mycoplasmas [2, 10]. However, our results are consistent with published data obtained by the transcriptional analysis of non-relatives to mycoplasmas bacteria, where the regulatory system of the SOS-response was described [1, 36, 37]. These observations may indicate the functionality of the SOS-response on the one hand and the presence of an unknown regulator on the other. In favor of the hypothesis of the presence of such a regulator, this may also indicate the presence of a number of genes in the genome, whose protein products, according to the annotation, have sequence-specific DNA-binding domains and could potentially act as transcription factors (data not shown). In addition rapid induction at the level of gene transcription in the second group (Figure 2B) may be indicative of the presence of a repressor that acts similar to the previously described LexA-repressor in the E. coli. The subsequent decline of the mRNA of parC and dinB genes may be due to the presence of a negative feedback regulation. This is particularly likely in the case of dinB, which encodes an alternative DNA polymerase as its activation can lead to dangerously high levels of mutagenesis. The increase in mRNA levels of most of the repair system genes may be a consequence of the global regulatory mechanism rather than the result of a single transcription factor.
Demonstrated here, the induction of DNA polymerase IV in different types of stress is consistent with the literature and may be a mechanism of adaptation to stress by increasing endogenous levels of mutagenesis [9, 13, 38, 39].
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 unduce SOS-response in a stationary phase.
Strains and conditions
Mycoplasma gallisepticum S6 was cultivated on a liquid medium  at 37°C for 12 and 24 hours for exponential and stationary phases, respectively (Additional file 2). Cells were passaged twice for 24 hours, starting from frozen culture prior to the experiment. Cells were passaged in 1:10 dilution. Cells were harvested by centrifugation at 8 000 g and 4°C for 10 min.
Cells were harvested as described above and lysed with CTAB buffer (2% CTAB, 100 mM Tris–HCl, pH = 8.0, 20 mM EDTA, 1.4 M NaCl) at 60°C for 30 min with subsequent chloroform extraction (1:1) and isopropanol precipitation (1:1) with the addition of 10% v/v sodium acetate .
RNA extraction and cDNA synthesis
Total RNA was extracted directly from the cell culture using TrizolLS (Invitrogen) reagent (1:3 cell culture: TrizolLS) with subsequent chloroform extraction (1:5 chloroform: TrizolLS) and isopropanol precipitation (1:1). RNA was treated by DNase I (Thermo Scientific) and used for cDNA synthesis with H-minus Mu-MLV reverse transcriptase (Thermo Scientific).
Real-time PCR and droplet-digital PCR
Real-time PCR was performed using iQ SYBR Green Supermix (Bio-Rad) and CFX96™ Real-Time PCR Detection System (Bio-Rad) PCR machine. Droplet digital PCR allows direct quantification of DNA molecules in a sample . It was performed using ddPCR™ Supermix for Probes (Bio-Rad) and QX100 system (droplet generator and droplet reader) along with a DNA Engine Tetrad 2 (Bio-Rad) PCR machine. Real-time and ddPCR data was analyzed with CFX Manager and QuantaSoft (Bio-Rad) software, respectively. Primers and probes are listed in Additional file 3. All PCR experiments were carried out in three biological and two technical replicates.
Quantification of RNA copy number per cell
The copy number of RNA and DNA molecules was measured using ddPCR. Final data was normalized with the respect of cell culture volume which results in a copy number of RNA and DNA per unit of culture volume. The copy number of RNA per cell was estimated as a ratio of RNA per DNA copy number based on the assumption that cells have one copy of genomic DNA.
Determination of sub-lethal conditions
Sub-lethal conditions were determined as conditions when stressful conditions are maximal but most of the cells are still viable. Working under such conditions ensures the maximal response in the absence of massive cell death, which can significantly affect results. Cell viability was estimated by the determination of colony forming units that are formed by cells after stress. Sub-lethal conditions for different stresses were found to be the following: 46°C for 1 hour for heat stress, 1.2 M NaCl for 1 hour, 0.02% H2O2 for 1 hour, 2 μg of ciprofloxacin for 4 hours, and 8 μg of tetracycline for 1 hour.
Protein extraction and 1D electrophoresis
Cells harvested as described were washed twice in a wash buffer (150 mM NaCl, 50 mM Tris–HCl, 2 mM MgCl2, pH = 7.4). Cells were lysed in 20 μl of 1% SDS in 100 mM NH4HCO3 and incubated in an ultrasonic bath for 15 min with subsequent centrifugation at 10000 g at 4°С for 5 min. Supernatant was taken, and protein concentration was determined by Bicinchoninic acid protein assay kit (Sigma). 20 μl of 2x Laemmli reagents were then added, and samples were incubated at 95°С for 5 min. Then 50 μg of protein was loaded to polyacrylamide gel (10×0.1 cm, 12% polyacrylamide), and electrophoresis was performed according to Laemmli  (10 mA current). Electrophoresis was stopped when the front dye reached 1.5 cm in separating gel.
Trypsinolysis in polyacrylamide gel
The polyacrylamide gel was fixed in a fixation buffer (20% CH3OH and 10% CH3COOH) for 30 min and washed twice in H20. The gel was cut in pieces 1×1 mm, transferred into tubes, and treated with 10 mM DTT and 100 mM NH4HCO3 for 30 min at 56°C. Then proteins were alkylated with 55 mM iodoacetamide in 100 mM NH4HCO3 for 20 min in the dark. Then water was removed from gel pieces by the addition of 100% acetonitrile.
Dehydrated samples were treated with a 150 μl trypsin solution (40 mM NH4HCO3, 10% acetonitrile, 20 ng/μl Trypsin Gold, mass spectrometry grade, Promega). Samples were incubated for 60 min at 40°C and for 16–18 h. at 37°C. Peptides were extracted one time by 5% formic acid and two times by 50% acetonitrile with 5% formic acid. Extracts were joined and dried on a vacuum centrifuge at 45°C. Precipitate was diluted in 50 μl of 5% acetonitrile with 0.1% formic acid.
Peptides were analyzed using a TripleTOF 5600+ (ABSciex) mass-spectrometer with NanoSpray III ion source and NanoLC Ultra 2D + chromatograph (Eksigent). Chromatographic separation was carried out in the gradient of acetonitrile in water (5 to 40% of acetonitrile in 120 min) with 0.1% formic acid on 75×150 μm columns with Phenomenex Luna C18 3 μm sorbent and a flow rate of 300 nL/min.
The IDA mode of a mass-spectrometer was used to analyze peptides. Based on the first MS1 spectrum (mass range for analysis and subsequent ion selection for MS2 analysis was 300–1250 m/z, signal accumulation was 250 ms), 50 parent ions with maximum intensity in the current spectrum were chosen for subsequent MS/MS analysis (resolution of quadrupole UNIT was 0.7 Da, measurement mass range was 200–1800 m/z, optimization of ion beam focus was to obtain maximal sensitivity, signal accumulation was 50 ms for each parent ion). Nitrogen was used for collision dissociation with fixed average energy of 40 V. Collision energy was linearly increased from 25 to 55 V during signal accumulation time (50 ms). Parental ions that had already been analyzed were excluded from analysis for 15 sec.
Analysis of mass-spectrometry data
Raw data was analyzed with ProteinPilot 4.5 revision 1656 (ABSciex) using search algorithm Paragon 126.96.36.199, revision 1654 (ABSciex), and standard search settings to search against a database of all proteins of M. gallisepticum S6 (genbank id: AFFR01000000). The following parameters were used for this search: alkylation of cysteine by iodoacetamide, trypsin digestion, TripleTOF 5600 equipment, and a deep search with additional statistical analysis of results reliability. Specters were grouped with default settings by a ProGroup algorithm build-in to ProteinPilot. The statistical analysis of results reliability (and identification of threshold value of unused score) was carried out by a ProteomicS Performance Evaluation Pipeline Software (PSPEP) algorithm build-in to ProteinPilot.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository  with the dataset identifier PXD000249 and DOI 10.6019/PXD000249.
Comparative analysis and in-silico reconstruction of DNA repair system
We have used Gene Ontology database (http://www.geneontology.org/) and literature analysis, we compiled the list of all the genes involved in DNA repair in E. coli and (or) B. subtilis. For all of the selected genes we found homologues in the genome of M. gallisepticum using blastn and blastp algorithms (e-value <1e-25). We aligned selected proteins of M. gallisepticum with its respective homologues of E. coli and (or) B. subtilis in order to analyze the amino acid substitutions in the active center with the help of ClustalW2 algorithm (http://www.ebi.ac.uk/Tools/msa/clustalw2/) and PDB database (http://www.rcsb.org/pdb/home/home.do). The results are shown in the Table 1 and the Additional file 4.
Statistical analysis of mRNA expression levels by quantitative RT-PCR (qRT-PCR)
Data presented are the average of three individual biological experiments with calculated standard deviation (Additional file 1); within each experiment, technical duplicates were performed. The level of each mRNA (log2) identified in each condition was compared with the control, which was exponential growth phase for all conditions and stationary phase for stationary phase heat stress. To identify the significance of a mRNA level change, we used a t-test with multiple hypothesis testing correction by the Benjamini-Hochberg procedure .
Identification of genes expression patterns during heat stress
The Gene expression levels (log2) during 15 and 30 min heat stress were averaged. After that, a rank of each mRNA was calculated for control, 5 min and an average 15–30 min stress. The condition with the maximal expression gained a rank of “1”, and the condition with the minimal expression gained a rank of “3”.
Lenhart JS, Schroeder JW, Walsh BW, Simmons LA: DNA repair and genome maintenance in Bacillus subtilis. Microbiol Mol Biol Rev. 2012, 76: 530-564. 10.1128/MMBR.05020-11.
Razin S, Yogev D, Naot Y: Molecular biology and pathogenicity of mycoplasmas. Microbiol Mol Biol Rev. 1998, 62: 1094-1156.
Alexeev D, Kostrjukova E, Aliper A, Popenko A, Bazaleev N, Tyakht A, Selezneva O, Akopian T, Prichodko E, Kondratov I, Chukin M, Demina I, Galyamina M, Kamashev D, Vanyushkina A, Ladygina V, Levitskii S, Lazarev V, Govorun V: Application of Spiroplasma melliferum proteogenomic profiling for the discovery of virulence factors and pathogenicity mechanisms in host-associated spiroplasmas. J Proteome Res. 2012, 11: 224-236. 10.1021/pr2008626.
Vanyushkina AA, Kamashev DE, Altukhov IA, Govorun VM: Identification of intracellular Spiroplasma melliferum metabolites by the HPLC-MS method. Biochemistry (Mosc). 2012, 77: 864-877. 10.1134/S000629791208007X.
Iyer RR, Pluciennik A, Burdett V, Modrich PL: DNA mismatch repair: functions and mechanisms. Chem Rev. 2006, 106: 302-323. 10.1021/cr0404794.
Dalhus B, Laerdahl JK, Backe PH, Bjørås M: DNA base repair–recognition and initiation of catalysis. FEMS Microbiol Rev. 2009, 33: 1044-1078. 10.1111/j.1574-6976.2009.00188.x.
Ischenko AA, Saparbaev MK: Alternative nucleotide incision repair pathway for oxidative DNA damage. Nature. 2002, 415: 183-187. 10.1038/415183a.
Dillingham MS, Kowalczykowski SC: RecBCD enzyme and the repair of double-stranded DNA breaks. Microbiol Mol Biol Rev. 2008, 72: 642-671. 10.1128/MMBR.00020-08. Table of Contents
Janion C: Inducible SOS response system of DNA repair and mutagenesis in Escherichia coli. Int J Biol Sci. 2008, 4: 338-344.
Carvalho M, Fonseca MM, Scortecci C, Alfredo C, Blaha G, Agnez-lima LF: DNA repair in reduced genome: The Mycoplasma model. DNA Repair. 2005, 360: 111-119.
Schlacher K, Leslie K, Wyman C, Woodgate R: DNA polymerase V and RecA protein, a minimal mutasome. Mol Cell. 2005, 17: 561-572. 10.1016/j.molcel.2005.01.006.
Foster PL: Stress-Induced Mutagenesis in Bacteria. Crit Rev Biochem. 2007, 42 (5): 373-397. 10.1080/10409230701648494.
Dörr T, Lewis K, VuliÄ M: SOS response induces persistence to fluoroquinolones in Escherichia coli. PLoS Genet. 2009, 5: e1000760-10.1371/journal.pgen.1000760.
Woese C: Bacterial evolution. Microbiol Rev. 1987, 51 (2): 221-271.
Gruson D, Pereyre S, Renaudin H, Charron A, Be C: In Vitro Development of Resistance to Six and Four Fluoroquinolones in Mycoplasma pneumoniae and Mycoplasma hominis. Respectively. 2005, 49: 1190-1193.
Curti E, McDonald JP, Mead S, Woodgate R: DNA polymerase switching: effects on spontaneous mutagenesis in Escherichia coli. Mol Microbiol. 2009, 71: 315-331. 10.1111/j.1365-2958.2008.06526.x.
Kamashev D, Oberto J, Serebryakova M, Gorbachev A, Zhukova Y, Levitskii S, Mazur AK, Govorun V: Mycoplasma gallisepticum Produces a Histone-like Protein That Recognizes Base Mismatches in DNA. Biochemistry. 2011, 50 (40): 8692-8702. 10.1021/bi2009097.
Levitskiy SA, Sycheva AM, Kharlampieva DD, Oberto J, Kamashev DE, Serebryakova MV, Moshkovskii SA, Lazarev VN, Govorun VM: Purification and functional analysis of recombinant Acholeplasma laidlawii histone-like HU protein. Biochimie. 2011, 93: 1102-1109. 10.1016/j.biochi.2011.03.005.
Lage C, De Pádula M, De Alencar TAM, Gonçalves SRDF, Vidal LDS, Cabral-Neto J, Leitão AC: New insights on how nucleotide excision repair could remove DNA adducts induced by chemotherapeutic agents and psoralens plus UV-A (PUVA) in Escherichia coli cells. Mutation Research/Reviews in Mutation Research. 2003, 544: 143-157. 10.1016/j.mrrev.2003.06.003.
Fisunov GY, Alexeev DG, Bazaleev NA, Ladygina VG, Galyamina MA, Kondratov IG, Zhukova NA, Serebryakova MV, Demina IA, Govorun VM: Core proteome of the minimal cell: comparative proteomics of three mollicute species. PLoS One. 2011, 6: e21964-10.1371/journal.pone.0021964.
Demina IA, Serebriakova MV, Ladygina VG, Galiamina MA, Zhukova NA, Alekseev DG, Fisunov GI, Govorun VM: Comparative proteomic characteristic of Mycoplasmas (Mollycutes). Bioorg Khim. 2011, 37: 70-80.
Demina IA, Serebryakova MV, Ladygina VG, Rogova MA, Zgoda VG, Korzhenevskyi DA, Govorun VM: Proteome of the bacterium Mycoplasma gallisepticum. Biochemistry (Mosc). 2009, 74: 165-174. 10.1134/S0006297909020072.
Graumann PL, Knust T: Dynamics of the bacterial SMC complex and SMC-like proteins involved in DNA repair. Chromosome Res. 2009, 17: 265-275. 10.1007/s10577-008-9014-x.
Wolff E, Kim M, Hu K, Yang H, Miller JH: Polymerases Leave Fingerprints : Analysis of the Mutational Spectrum in Escherichia coli rpoB To Assess the Role of Polymerase IV in Spontaneous Mutation. J Bacteriol. 2004, 186: 2004-2900.
Sciences ML: Review The bacterial LexA transcriptional repressor. Cell Mol Life Sci. 2009, 66: 82-93. 10.1007/s00018-008-8378-6.
Savijoki K, Ingmer H, Frees D, Vogensen FK, Palva A, Varmanen P: Heat and DNA damage induction of the LexA-like regulator HdiR from Lactococcus lactis is mediated by RecA and ClpP. Mol Microbiol. 2003, 50: 609-621. 10.1046/j.1365-2958.2003.03713.x.
Cirz RT, O’Neill BM, Hammond JA, Head SR, Romesberg FE: Defining the Pseudomonas aeruginosa SOS response and its role in the global response to the antibiotic ciprofloxacin. J Bacteriol. 2006, 188: 7101-7110. 10.1128/JB.00807-06.
Robertson AB, Matson SW: Reconstitution of the very short patch repair pathway from Escherichia coli. J Biol Chem. 2012, 287: 32953-32966. 10.1074/jbc.M112.384321.
Razin A, Razin S: Methylated bases in mycoplasmal DNA. Nucleic Acids Res. 1980, 8: 1383-1390. 10.1093/nar/8.6.1383.
Nur I, Szyf M, Razin A, Glaser G, Rottem S, Razin S: Procaryotic and eucaryotic traits of DNA methylation in spiroplasmas (mycoplasmas). J Bacteriol. 1985, 164: 19-24.
Lluch-Senar M, Luong K, Lloréns-Rico V, Delgado J, Fang G, Spittle K, Clark TA, Schadt E, Turner SW, Korlach J, Serrano L: Comprehensive Methylome Characterization of Mycoplasma genitalium and Mycoplasma pneumoniae at Single-Base Resolution. PLoS Genet. 2013, 9: e1003191-10.1371/journal.pgen.1003191.
Wang D, Kreutzer DA, Essigmann JM: Mutagenicity and repair of oxidative DNA damage: insights from studies using defined lesions. Mutat Res. 1998, 400: 99-115. 10.1016/S0027-5107(98)00066-9.
Maier T, Schmidt A, Güell M: Quantification of mRNA and protein and integration with protein turnover in a bacterium. Mol Syst. 2011, 7: 1-12.
Güell M, Van Noort V, Yus E, Chen W-H, Leigh-Bell J, Michalodimitrakis K, Yamada T, Arumugam M, Doerks T, Kühner S, Rode M, Suyama M, Schmidt S, Gavin A-C, Bork P, Serrano L: Transcriptome complexity in a genome-reduced bacterium. Science. 2009, 326: 1268-1271. 10.1126/science.1176951.
Mandelco L, Sechrest J, Lawrence TG, Van EJ, Weisburg WG, Tully JG, Rose DL, Petzel JP, Oyaizu H, Yang D: A Phylogenetic Analysis of the Mycoplasmas : Basis for Their Classification. J Bacteriol. 1989, 171: 6445-6467.
Sioud M, Boudabous A, Cekaite L: Transcriptional responses of Bacillus subtillis and thuringiensis to antibiotics and anti-tumour drugs. Int J Mol Med. 2009, 23: 33-39.
Au N, Kuester-schoeck E, Mandava V, Bothwell LE, Canny SP, Chachu K, Colavito SA, Fuller SN, Groban ES, Hensley LA, Brien TCO, Shah A, Tierney JT, Tomm LL, Gara TMO, Goranov AI, Grossman AD, Lovett CM: Genetic Composition of the Bacillus subtilis SOS System. J Bacteriol. 2005, 187: 7655-7666. 10.1128/JB.187.22.7655-7666.2005.
Giuliodori AM, Gualerzi CO, Soto S, Vila J, Tavío MM: Review on bacterial stress topics. Ann N Y Acad Sci. 2007, 1113: 95-104. 10.1196/annals.1391.008.
Layton JC, Foster PL: Error-Prone DNA Polymerase IV Is Regulated by the Heat Shock Chaperone GroE in Escherichia coli. J Bacteriol. 2005, 187: 449-457. 10.1128/JB.187.2.449-457.2005.
Morowitz H, Maniloff J: Analysis of the life cycle of Mycoplasma gallisepticum. J Bacteriol. 1966, 91: 1638-1644.
Stewart CN, Via LE: A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques. 1993, 14: 748-750.
Hindson BJ, Ness KD, Masquelier DA, Belgrader P, Heredia NJ, Makarewicz AJ, Bright : High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal Chem. 2011, 83: 8604-8610. 10.1021/ac202028g.
Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage. Nature. 1970, 15: 680-685.
Vizcaíno JA, Côté RG, Csordas A, Dianes JA, Fabregat A, Foster JM, Griss J, Alpi E, Birim M, Contell J, O’Kelly G, Schoenegger A, Ovelleiro D, Pérez-Riverol Y, Reisinger F, Ríos D, Wang R, Hermjakob H: The PRoteomics IDEntifications (PRIDE) database and associated tools: status in 2013. Nucleic Acids Res. 2013, 41 (Database issue): D1063-D1069.
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Stat Soc. 1995, 57: 289-300.
This work was supported by a Grant (12-04-31696) from Russian Foundation for Basic Research (RFBR).
The authors declare that they have no competing interests.
AG designed the project, performed experimental and bioinformatics data and wrote the paper; GF designed the project, performed bioinformatics, qRT-PCR data and wrote the paper; MI performed qRT-PCR data and corrected the paper; DE performed functional protein domain search; PM performed bioinformatics and statistical calculations and wrote the paper; DA performed calculation for mass-spectrometry data and corrected the paper; OP performed protein extraction and sample preparation for chromato-mass-spectrometry; TG performed trypsinolysis for chromato-mass-spectrometry; SK performed mass-spectrometry data and wrote the paper; DK designed the project and wrote the paper, VG designed and supervised the project and wrote the paper. All authors read and approved the final manuscript.
Alexey Y Gorbachev, Gleb Y Fisunov contributed equally to this work.