Skip to main content
  • Research article
  • Open access
  • Published:

Intrinsic terminators in Mycoplasma hyopneumoniae transcription

Abstract

Background

Mycoplasma hyopneumoniae, an important pathogen of swine, exhibits a low guanine and cytosine (GC) content genome. M. hyopneumoniae genome is organised in long transcriptional units and promoter sequences have been mapped upstream of all transcription units. These analysis provided insights into the gene organisation and transcription initiation at the genome scale. However, the presence of transcriptional terminator sequences in the M. hyopneumoniae genome is poorly understood.

Results

In silico analyses demonstrated the presence of putative terminators in 82% of the 33 monocistronic units (mCs) and in 74% of the 116 polycistronic units (pCs) considering different classes of terminators. The functional activity of 23 intrinsic terminators was confirmed by RT-PCR and qPCR. Analysis of all terminators found by three software algorithms, combined with experimental results, allowed us to propose a pattern of RNA hairpin formation during the termination process and to predict the location of terminators in the M. hyopneumoniae genome sequence.

Conclusions

The stem-loop structures of intrinsic terminators of mycoplasma diverge from the pattern of terminators found in other bacteria due the low content of guanine and cytosine. In M. hyopneumoniae, transcription can end after a transcriptional unit and before its terminator sequence and can also continue past the terminator sequence with RNA polymerases gradually releasing the RNA.

Background

Transcription is a highly regulated multi-step process roughly divided into initiation, elongation and termination. Prokaryotic transcription termination can occur via factor-dependent or factor-independent termination and is known to play key roles in regulating genetic systems. Factor-independent termination occurs at defined sequence regions known as intrinsic terminators and factor-dependent termination relies on the destabilisation of transcription complexes by a Rho regulatory protein [1]. Sequence features within intrinsic terminators have been studied in different organisms resulting in models for how these motifs contribute to overall termination efficiency, and have been used to develop computational methods to identify terminator elements within genome sequences [2-4].

Mycoplasmas are cell wall-less prokaryotes of the class Mollicutes, phylogenetically related to Gram-positive bacteria and characterised by having small genomes with a low GC content. Mycoplasma hyopneumoniae is considered the etiological agent of porcine enzootic pneumonia, a disease with global distribution and causing significant economic losses in the pig farming industry [5].

In recent years, the genomes of many Mycoplasma species have been completely sequenced, including those of some M. hyopneumoniae strains [6-8]. The sequencing of several mycoplasma genomes has provided an extensive comparative analysis of gene content among different species and information related to gene organisation in large transcriptional units (TUs) [8-10]. Moreover, promoter sequences have been defined and mapped for the M. hyopneumoniae genome [11,12] and recently, transcriptome analysis has validated the organization of M. hyopneumoniae genome in long transcriptional units [10].

However, despite the genome-scale sequencing efforts, prediction and recognition of mycoplasma terminator elements is poorly known. The unusually low GC content of intergenic regions (IRs) in the mycoplasma genome [13] may be the main reason for the lack of success in application of existing computational methods developed to predict putative intrinsic terminators. Moreover, gene coding for Rho regulatory protein was not found in mycoplasma genomes suggesting the absence of factor-dependent transcription termination in these organisms [4].

Studies investigating the presence of intrinsic terminators in mycoplasma were initially contradictory. Analysis of average RNA folding energy near stop codons demonstrated that no stem-loops were formed in Mycoplasma genitalium or Mycoplasma pneumoniae coding sequences (CDS), indicating the existence of qualitatively different and uncharacterised mechanisms for transcription termination [14]. However, the development of algorithms to predict terminators in low GC content genomes allowed the identification of stem-loop structures in some positions of mycoplasma genomes [3,15-17]. Furthermore, several studies showed that intrinsic terminators have a functional role in mycoplasmas [18-20] suggesting that the intrinsic terminators could be the main mode of termination.

Previously, we have described the genome organisation profile and mapped promoter sequences upstream of all transcription units of the M. hyopneumoniae genome [9,11,12]. These results provided insights into the gene organisation and transcription initiation at the genome scale. To further understand the mechanism of transcription in M. hyopneumoniae genome, in the current study we have analysed and predicted the presence of terminator sequences downstream of the transcription units. We have also determined the role of these predicted sequences, revealing that intrinsic terminators are the main mechanism of transcription termination in M. hyopneumoniae.

Methods

In silico analysis of terminators

The prediction of terminator sequences was performed in M. hyopneumoniae 7448 (NC_007332) using three software algorithms: TransTermHP [16], WebGesTer [17] and ARNold [21]. The software WebGesTer restricts the search for palindrome sequences in a region ranging from −20 to +270 bp from the CDS stop codon and allows the search of non-canonical terminators (without the U-tract). TransTermHP algorithm restricts the search of palindrome sequences to regions in which occurs the presence of at least three thymines in sequence and, therefore, unable to find non-canonical terminators. ARNold combines two algorithms, Erpin and RNAmotif, which are based on terminators of Escherichia coli and Bacillus subtilis and do not consider the gene context in which the terminators occur. Figure 1 shows the workflow of our terminator prediction classification. Class 1 terminators (tc1) were defined using genome localisation parameters as follows: terminator sequences located at the 3’ end of the last gene of the polycistronic unit (pC); terminator sequences located at the 3’ end of the monocistronic unit (mC); terminator sequences located at the 3’ end of genes within the pCs. For subsequent analyses steps we have selected all terminator sequences located at the 3’ end of pCs or mCs and established new criteria for terminator classification: i) class 2 terminator (tc2) sequences were considered when positioned at the end of pC or mC and predicted by at least two different algorithms; ii) class 3 terminators (tc3) classification considered all tc1 with two additional features: first the distance between the terminator first nucleotide and the stop codon of the target gene should range from −11 to 200 base pairs and second, the values of Gibbs free energy (ΔG) should be less than −4 kcal/mol; and iii) class 4 terminators (tc4) were all tc1 that contained only one of the features of class 3 terminators. The terminators in the M. hyopneumoniae genome were mapped using Artemis software [22].

Figure 1
figure 1

In silico design workflow utilised to classify terminators predicted in M. hyopneumoniae 7448. Predicted terminators were analysed by the software ARNold, TransTermHP and WebGesTer.

Primer design

Specific RT-PCR and qPCR primers were designed (based on the M. hyopneumoniae genome sequence — GenBank access: NC_007332) to target transcription before and after the selected predicted position of the terminators. Figure 2 illustrates an example (M. hyopneumoniae 7448 TU_10) of the position of primer pairs to validate the putative terminator. Figure 2A shows a global view of the putative terminator w_014 and Figure 2B demonstrates a detailed position of the primers with the product length expected. As shown in Figure 2B we designed two pairs of primers for each terminator analysis for both RT-PCR and qPCR methodologies. Briefly, the forward primers were positioned inside the coding region of the selected CDS (named Te.F and Ex.F or Ge.F when the primer forward was the same for the two reactions). One reverse primer was positioned between the stop codon of the CDS and the start of the putative terminator (upstream terminator – Te.R). The second reverse primer was positioned immediately after the putative terminator (downstream terminator – Ex.R). A total of 60 primer pairs were prepared for confirmation of 28 putative terminators (Additional file 1). Primers were designed in Vector NTI Advance 10 (Invitrogen, USA).

Figure 2
figure 2

Schematic representation of primers design for experimental analyses by RT-PCR or qPCR. A: Representation of polycistronic unit 10 with the terminator w_014 (hairpin). B: Enlargement of highlighted region in A, showing the position of primers for terminator validation. Black arrows Te.F, Te.R, Ex.F and Ex.R represent the primers U.10.Te.F, U.10.Te.R, U.10.Ex.F and U.10.Ex.R, respectively. (U.10 – polycistronic unit 10, Te – terminator, Ex – external, F – forward and R – reverse).

Culture conditions and RNA isolation

M. hyopneumoniae strain 7448 isolated from an infected swine (Embrapa, Santa Catarina, Brazil) [6] was grown in 25 ml of Friis broth [23] at 37°C for 24 h with gentle agitation in a roller drum. Total RNA was isolated with RNeasy Mini Kit (Qiagen, USA). For cell lyses, 0.7 ml of RNeasy Lysis Buffer (RLT buffer) in the presence of 0.134 M of β-mercaptoethanol was used per cultivation flask. The purification was performed according to the manufacturer’s instructions, with on-column DNaseI digestion using the RNase-Free DNase Set (Qiagen, Germany) and a second round of treatment with DNase I (Fermentas, USA). DNA absence was monitored to below PCR-detectable levels. Extracted RNA was analysed by gel electrophoresis and quantified in the Qubit system (Invitrogen, USA). Finally, RNA quality and integrity were determined by the evaluation of the RNA Integrity Number (RIN) using the Agilent 2100 Bioanalyzer (Agilent, USA). Values equal to or greater than 9.5 indicated sufficient quality.

Reverse transcriptase PCR (RT-PCR) assay

For first-strand cDNA synthesis we used 1 μg of total RNA, 10 pmol of specific reverse primer (Additional file 1) and 10 mM of deoxynucleotide triphosphates. The mixture was heated to 70°C for 5 min and then incubated on ice for 5 min. First-strand buffer (Promega, USA), 0.1 M dithiothreitol, 40 U RNase inhibitor (Promega, USA) and 200 U M-MLV RT (Moloney Murine Leukemia Virus Reverse Transcriptase—Promega, USA) were then added to a total volume of 20 μl. The reaction was incubated at 37°C for 60 min followed by 15 min at 70°C for enzyme inactivation. A negative control was prepared in parallel, differing only by the absence of the RT enzyme.

PCRs included 1 U GoTaq DNA polymerase (Promega, USA), 5X of GoTaq buffer, 1 mM of each deoxynucleotide triphosphate, 10 pmol of each primer (Additional file 1) and 1 μl of the first-strand cDNA reaction in a final volume of 25 μl. A negative control of RT-PCR was prepared in parallel, which differed only by the absence of cDNA, and no genomic DNA was added to the reaction mixture for the PCR negative control. A PCR positive control was prepared using the genomic DNA of M. hyopneumoniae 7448 as the template. The PCR conditions were: 1 cycle at 94°C for 5 min followed by 30 cycles of 94°C for 30 s; denaturation and extension temperature and time varied according to each primer pair (Additional file 1). The final extension step was at 72°C for 10 min. Reaction products were analysed in 1.2% agarose gels.

qPCR experimental design

Quantitative RT-PCR was performed using 1:50 cDNA prepared as described for RT-PCR as template and Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen, USA) on the StepOne Real-Time PCR System (Applied Biosystems, USA). The qPCR reactions were carried out at 90°C for 2 min and 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 55°C for 1 min each. The specificity of the synthesised products and the absence of primer dimers were visualised using a melting curve analysis for each reaction. Amplification efficiency for each primer pair was calculated using the LinRegPCR software application [24] and the mean efficiency values for each primer were added to Additional file 1. This efficiency value was used for the quantification analysis.

Relative expression of mRNA was calculated by the 2-ΔCt method [25]. To control for all the experiments the threshold cycle (CT) values were normalised to the reference gene MHP7448_0333 [12,26]. The CT of each test target represents the average of three reactions. Three independent biological replicates were done for each target gene. We performed statistical analysis using GraphPad Prism 6 software. A two-tailed unpaired t-test was used to test for differences in the relative expression values between the regions before and after the terminator structure (P < 0.05).

Results

Prediction of terminator sequences

Computational analysis of M. hyopneumoniae genome using three different algorithms predicted 1068 terminators: 439 terminators using ARNold software, 334 using TransTermHP software and 295 using WebGesTer software (Table 1 and Additional file 2). These terminators were named using the initial letter of the software (a = ARNold, t = TransTermHP and w = WebGesTer) followed by a number. For example terminator t_001 is the first terminator predicted by TransTermHP software. Terminators named “wf” represent terminators classified as unstable terminators in the WebGesTer software (Additional files 2, 3 and 4).

Table 1 Terminators predicted in Mycoplasma hyopneumoniae

The classification criteria defined in Figure 1 were applied to validate the predicted terminators. Class 1 terminators, named tc1, should be positioned at the 3’ end of the CDS. Using this criterion, 556 terminators of the 1068 predicted terminators were selected (52%) (Table 1). The remaining predicted terminators were localised within coding regions or on the antisense strand and therefore, were not further analysed in this work.

The 556 tc1 terminators were analysed according to the genomic context of M. hyopneumoniae. The M. hyopneumoniae genome is organised in 33 monocistronic units (mCs) and 116 polycistronic units (pCs) containing two or more ORFs [9]. Therefore, systematic terminator localisation was performed in the downstream regions of all annotated ORFs in the genome. This analysis revealed the presence of 65 tc1 terminators at the end of mCs, 255 tc1 terminators at the end of pCs and 236 tc1 terminators at the 3’ end of some ORFs that were located internal in some polycistronic units (identified as internal genes) (Table 1). Moreover, more than one predicted terminator sequence can be found in some transcriptional units. Class 1 terminators were identified in all 33 mCs and in 106 of the 116 pCs (Additional files 3 and 4) suggesting that transcription termination occurs preferentially at the ends of mCs and pCs.

The class 1 terminators (tc1) were further divided into class 2, class 3 and class 4 terminators. To ensure the reliability of the generated data, the mC and pC tc1 terminators located at the same position by at least two algorithms were considered putative terminators and named class 2 terminators (tc2; Additional file 5). Using this approach we were able to identify class tc2 terminators in seven of the 33 mCs (21%) and in 39 of the 116 pCs (34%) (Table 2). Detailed analysis of the characteristics of these terminators revealed the presence of two distinguishing features: i) the distance between the start of the terminator and the stop codon of the target gene showed a range of −11 bp to 200 bp and ii) the values of Gibbs free energy (ΔG) were less than −4 kcal/mol (Figure 3). The presence of two class tc2 terminators at a distance greater than 230 bp from the stop codon (see Figure 3A) was identified at 3’ end of TU_09 (fusA) and TU_114 (MHP7448_r2). However, further analysis revealed the presence of other tc2 terminators at the 3’ end of these two pCs at position 39 or 41 for TU_09 and at positions 167, 185 or 186 for TU_114 (Additional file 2).

Table 2 Terminator distribution according to the genome organisation of M. hyopneumoniae
Figure 3
figure 3

Characteristics of class 2 terminators. The graphics show the distribution of the class 2 terminators according to distance to stop codon (A) and the values of Gibbs free energy determined by the algorithm that predicted the terminator (B). Negative distances indicate the overlap with the stop codon of the gene.

To increase the number of mCs and pCs containing putative terminators (similar to class 2 terminators) all the remaining tc1 terminators were screened using the features defined in class 2 terminators (see Figure 3). Terminators that satisfied both criteria (distance to stop codon and ΔG) were reclassified as class 3 terminators (tc3) and those that displayed only one of the two criteria were named class 4 terminators (tc4). It is important to point out that, in contrast to class 2 terminators, the class 3 terminators were defined by only one algorithm. Class 3 terminators were identified in 14 of 33 mCs (43%) and in 36 of 116 pCs (31%) increasing the number of mCs and pCs with putative terminators to 21 (64%) and 75 (65%), respectively (Table 2). Class 4 terminators (tc4) were found in six mCs (18%) and eleven pCs (9%) of the M. hyopneumoniae genome.

Experimental validation of M. hyopneumoniae terminators

RT-PCR and qPCR were used to test whether the predicted terminators were related to transcription termination in M. hyopneumoniae. The terminator validation hypothesis was based on the following premises: i) mC or pC transcripts ended at the predicted terminator sequences; and ii) no mC or pC transcripts were present after the predicted terminator sequences. Primer pairs (Figure 2) were positioned to validate this hypothesis, similar to the terminator validation previously used by Arrebola et al. [27].

RT-PCR was performed for 15 putative terminator sequences. Among the putative terminator sequences analysed, four belonged to class tc2, one to class tc3 and 10 to class tc4 (Table 3). The Figure 4 shows the validation of a tc2 terminator (w_014) in TU_10, used as example, where two primer pairs were designed for RT-PCR analysis. As shown in Figure 2, primer pairs were positioned to analyse the presence of transcripts upstream and immediately downstream of the predicted terminator. TU_10 is composed of genes MHP7448_0081 and nox. At the end of the nox gene, three tc2 terminators were identified (Additional file 2). The experimental analysis considered the w_014 terminator located at 39 bp of the nox gene stop codon. This terminator is formed of a stem of 11 nucleotides with only one mismatch, one loop of three nucleotides, the U-tract and has a ΔG of −7.08 kcal/mol (Figure 4A). As demonstrated in Figure 4B (lines 3 and 7) amplification was observed only when primers were positioned upstream of the terminator sequence indicating that the w_014 terminator is responsible for transcription termination in M. hyopneumoniae TU_10.

Table 3 Predicted terminators analysed by RT-PCR
Figure 4
figure 4

Experimental analyses of polycistronic unit 10 (TU_10) by RT-PCR. A: Structure and characteristics of class 2 terminator w_014, predicted by the software WebGesTer. B: agarose gel electrophoresis of RT-PCR products obtained from primers used to analyse expression in TU_10. The reactions were performed in the presence (+) and absence (−) of reverse transcriptase. Lanes 2, 3, 4 and 5 represent the reaction with the primer positioned upstream of the terminator. Lanes 6, 7, 8 and 9 represent the reaction with the primer positioned downstream of the terminator. Results of the reaction with positive cDNA are presented in lanes 3 and 7. Control reactions (DNA) using genomic DNA as the template are presented in lanes 2 and 6. Lanes 4 and 8 represent the control reaction of cDNA synthesis while lanes 5 and 9 are the control of the PCR reaction. Lane 1 contains a 100 bp ladder.

Transcription termination was also demonstrated in the presence of tc2 terminators t_264 and w_102 located at TU_89 and TU_113, respectively (Table 3). In six tc4 terminators (wf_491, wf_513, wf_375, wf_137, wf_202 and wf_215), amplification product was detected only when using primers located upstream of the putative terminator sequence, suggesting that the presence of tc4 terminators at the 3’ end of mC_28, mC_33, TU_60, TU_66, TU_98 and TU_103 were responsible for transcription termination (Table 3). These data shows that the putative intrinsic terminators are functional in mycoplasmas.

In the other class 2 terminator (w_076), class 3 terminator (w_053) and four class 4 terminators (wf_204, w_202, wf_075, w_147), amplification products were detected when primers located downstream of the terminator were used and, therefore, could not be validated by RT-PCR. However, it is possible that some terminators, although incapable of efficient transcription termination, could cause a decrease in the level of transcription. To confirm this hypothesis, RNA expression level was quantified using qPCR.

Real-time quantitative PCR was performed for selected class 2, class 3 and class 4 terminators based on the feasibility of primer design, as the M. hyopneumoniae genome has an intergenic low GC content. A total of 16 terminators were analysed by qPCR (Table 4): eight tc2 (Figure 5), seven tc3 (Figure 6) and one tc4 terminators (Figure 7). In order to compare RT-PCR and qPCR methodologies, terminator w_102 of TU_113, a tc2 terminator previously confirmed by RT-PCR (Table 3), was also analysed by qPCR (Figure 5).

Table 4 Predicted terminators analysed by qPCR
Figure 5
figure 5

Analysis of relative expression by qPCR of class 2 terminators. Graphics represent the analysis of the terminator present in each of the pCs (terminator w_109 in TU_09, terminator t_022 in TU_11, terminator t_055 in TU_20, terminator a_207 in TU_58, terminator t_173 in TU_62, terminator a_310 in TU_83, terminator w_093 in TU_105 and terminator w_102 in TU_113). Column “Terminator” shows the level of expression of the primer positioned upstream of the terminator and the column “Extern” represents the data for the primer positioned downstream of the terminator. Data are presented as mean ± standard deviation of three independent experiments. Asterisks indicate statistically significant differences in levels of expression downstream of the terminator; *0.01 < P < 0.05; **0.001 < P <0.01; ***P < 0.001.

Figure 6
figure 6

Analysis of relative expression by qPCR of class 3 terminators. Graphic represent the analysis of terminator present in mCs and pCs (w_016 in mC_04, w_053 in mC_15, w_082 in mC_24, w_140 in TU_31, w_055 in TU_47, wf_140 in TU_68 and w_181 in TU_77). Column “Terminator” shows the level of expression of the primer positioned upstream of the terminator and the column “Extern” represents the data of the primer positioned downstream of the terminator. Data are presented as mean ± standard deviation of three independent experiments. Asterisks indicate statistically significant differences in expression levels downstream of the terminator; *0.01 < P < 0.05; **0.001 < P < 0.01.

Figure 7
figure 7

Analysis of relative expression by qPCR of class 4 terminator. Graphic represents the analysis of the terminator wf_204 present in mC_30. Column “Terminator” shows the level of expression of the primer positioned upstream of the terminator and the column “Extern” represent the data of the primer positioned downstream of the terminator. Data are presented as means ± standard deviation of three independent experiments. Asterisk indicates a statistically significant difference in expression levels downstream of the terminator wf_204 in mC_30; *P < 0.05.

For the majority of tested terminators primer efficiency was greater than 80%. Only the primers used to analyse the terminator in TU_62 showed, in both primer pairs, a lower efficiency of 65% (Additional file 1). Moreover, the primer pairs designed for terminator sequences w_102 in TU_113 and t_055 in TU_20 revealed differences in efficiency between upstream and downstream primers of around 50% (Additional file 1). This result is probably related to the absence of detected transcripts in the qPCR (Figure 5).

The eight class 2 transcriptional terminators analysed by qPCR were localised at 3’ end of the last gene in TU_09, TU_11, TU_20, TU_58, TU_62, TU_83, TU_115 and TU_113 (Figure 5). Real-time quantitative PCR demonstrated that only transcripts found for TU_09, TU_11, TU_20, TU_58 and TU_113 had similar expression levels to the gene used for normalization (MHP7448_333) in these experiments. The transcriptional terminators were characterised by comparison between transcripts level using primers upstream of the putative terminator sequence (Terminator) and downstream of the terminator sequence (External). Terminators w_109 in TU_09, t_055 in TU_20 and w_102 in TU_113 can be considered valid terminators as no transcripts were detected downstream of the terminator sequence (Figure 5 see External). The other tc2 class terminators, t_022 in TU_11, a_207 in TU_58, t_173 in TU_62, a_310 in TU_83, and w_093 in TU_105 displayed a significant (0.001 < P < 0.05) decrease in transcription when downstream primers were used (see Figure 5 Terminator versus External) and also could be related to transcriptional termination.

Transcriptional terminators classified as class 3 (w_016 in mC_04, w_053 in mC_15, w_082 in mC_24, w_140 in TU_31, w_055 in TU_47, wf_140 in TU_68 and w_181 in TU_77) and class 4 (wf_204 in mC_30) terminators were also analysed by qPCR (Figure 6 and Figure 7). Similar to some tc2 class terminators, all tc3 and tc4 class terminators displayed a significant (0.01 < P < 0.05) decrease in transcription when downstream primers were used, with exception of w_082 in mC_24 (Figure 6 and Figure 7; see Terminator versus External). Interestingly, a second tc3 (w_083) putative terminator was localised at 74 bp downstream of the first terminator w_082 (Additional file 2). The presence of two putative terminators in this monocistronic unit suggests that they could have complementary functions.

Aiming to summarize our results Figure 8 shows the localisation of both promoter sequences and intrinsic terminators according the genome organisation of M. hyopneumoniae. Previous works of our group have demonstrated that the genome is organized preferably in long transcriptional units containing promoter sequences upstream the start codon of the first gene of the units [9,11]. Among the putative class 2 and class 3 terminators, represented in Figure 8, the functional terminators at the end of TU_09, TU_10, TU_11 and mC_04, were also experimentally demonstrated by RT-PCR and qPCR analyses (Tables 3 and 4).

Figure 8
figure 8

M. hyopneumoniae genome organisation containing promoter and terminator transcription sites. Figure shows the transcription unit organisation in the genome region from 87,839 to 114,900 base pair. Genome organisation and promoter identification were previously defined [9,11]. Black arrows represent promoter sequence and black circle represent terminator sequence.

Discussion

The importance of the characterisation of Rho-independent termination in mycoplasma can be related with the accurate prediction of transcription units in these bacteria. In general, the properties of intrinsic terminators, such as sequence and the structural features of the hairpin and, in some cases, the U-tract of the nascent transcript, could be involved in transcription termination. Previous works have demonstrated the predominance and conservation of Rho-independent or intrinsic termination among the Mollicutes class [3,4,16].

We identified transcription terminators for M. hyopneumoniae in silico using three algorithms that demonstrated differences in efficiency of prediction. This method allowed the prediction of 1068 terminators in M. hyopneumoniae genome. Moreover, detailed analysis revealed the prediction of at least one terminator sequence in 253 of the 705 M. hyopneumoniae ORFs (35.9%).

In comparison to E. coli, the intrinsic terminators found in the M. hyopneumoniae genome display relatively high values of Gibbs free energy density of stem-loop formation correlated to the low content of guanine and cytosine in the Mollicutes genome [13]. Similar values were also found in intrinsic terminators of M. genitalium and M. pneumoniae and are probably responsible for previous difficulties identifying intrinsic terminators in M. pneumoniae based on analysis of average RNA folding energy near stop codons [3,14]. The presence of a U-tract following the stem-loop structure was detected for the majority of the predicted M. hyopneumoniae terminators. However, approximately 20% of M. hyopneumoniae stem-loop formation revealed absence of the U-tract tail of the nascent transcript. These results are similar to those found in other mycoplasma species, suggesting that the U-tract tail is not an essential requisite for intrinsic terminators in this genus [4]. Nevertheless, the efficiency of transcription termination is maintained in intrinsic terminators without the U-tract [28,29].

Detailed localisation analysis of the 1068 predicted terminators allowed us to select only those found at the 3’ end of pCs or mCs. Using this approach we were able to distinguish inter- from intracistronic terminator-like structures and select 320 putative Rho-independent terminators at the end of mCs and pCs. The prediction efficiency was different among the software, reflecting the specific algorithm used for each one. The best results were obtained with the software WebGesTer, as is the only one of the three software packages that takes into consideration the gene context, considering a region ranging from −20 to +270 bp from the stop codon of CDS, and allows the detection of terminators without a U-tract [17]. Of the 295 terminators predicted by this program 279 were classified as class tc1. The software TransTermHP identified 191 terminators class tc1. The lower number of predicted terminators found by the TransTermHP software is related to the restrict search of palindromic sequences in regions with the presence of at least three thymines in sequence. Therefore, the TransTermHP program is able to found only canonical sequences [16]. The ARNold software use the patterns defined for E. coli and B.subtilis and do not take in consideration the gene context in which the terminator occur [21]. This probably explains why only 86 of the 439 (19%) predicted terminators, classified as class 1 terminators, were found by this program. The remaining terminators predicted by ARNold program were localised within coding regions or on the antisense strand and therefore were not analysed in this paper.

Class tc1 terminators were identified in all 33 mCs and in 106 of 116 pCs. This method was unable to localise terminators in only 10 (8.6%) pCs. Interestingly, in seven of them (TU_15, TU_16, TU_52, TU_72, TU_78, TU_104 and TU_106) the last ORF is classified as hypothetical. Moreover, terminator sequences were found upstream of the last ORF in TU_15, TU_16, TU_72 and TU_106 suggesting that transcription could be terminated before the last ORF in these polycistronic units. In three pCs (TU_72, TU_104 and TU_115) the presence of a gene similar to one that codes for transposase was detected. Therefore, we suggest that Rho-independent termination occurs at specific sequences in M. hyopneumoniae genome.

Aiming to distinguish M. hyopneumoniae Rho-independent terminators from random stem-loop structure sequences we established rules to discover their distinguishing properties. We used class 2 terminator classification (tc2) to propose a model of intrinsic terminators and as a decision rule to define as putative terminators in M. hyopneumoniae genome. The decision rule had four parameters: i) an inverted repeat in the primary DNA sequence positioned at the 3’ end of pC or mC followed or not by a short stretch of thymine residues; ii) predicted by at least two different algorithms; iii) distance between the terminator first nucleotide and the stop codon of the target gene range from −11 to 200 bp; and iv) the values of Gibbs free energy (ΔG) should be less than −4 kcal/mol. However, the application of these four rules validated terminators in silico in only 21% of the mCs (seven of the 33) and in 34% of the pCs (39 of the 116).

Therefore, to increase the number of putative terminators of the 320 class 1 terminators we lowered the specificity of the class 2 decision rule as follows: class tc3 — class 1 intrinsic terminators validated by only one algorithm, maintaining rules i, iii and iv, described above; and class tc4 — class 1 intrinsic terminators validated by only one algorithm, maintaining rules i and iii or iv. Using these new rules the numbers of transcriptional units containing putative terminators analysed in silico increased to 27 mCs and 86 pCs (tc2, tc3 and tc4). Therefore, putative intrinsic terminators have been found, by in silico analysis, in 76% of transcription units of the M. hyopneumoniae genome (82% of the 33 mCs and in 74% of the 116 pCs).

Previous studies have demonstrated the presence of intrinsic terminators in genomes of M. genitalium (20%) and M. pneumoniae (24%) [16], and also in specific genes such as the 16S-23S rRNA operon of M. hyopneumoniae [30], the P1 operon of M. genitalium [31], the MgPa and P65 operons, and the ldh gene of M. pneumoniae [32-34]. Moreover, transcriptional terminators in bacterial genomes belonging to the Firmicutes phylum have also been predicted revealing a high level of conservation [3].

Little information on transcription termination in M. hyopneumoniae is available. To validate the in silico putative terminators, experiments were performed for analysis of the transcription end guided by putative terminators. The experimental analyses evaluated 28 of the 113 putative intrinsic terminators belonging to class 2, class 3 and class 4 terminators. The number of terminators experimentally analysed (24%) was limited due to difficulties in designing primers for some locations as these intergenic regions are characterised by a high adenine and thymine content. In eleven putative terminator sequences, located at the 3’ end of two mCs and nine pCs, RT-PCR or qPCR demonstrated transcription termination. However, for some intrinsic terminators (12 terminators) transcription continued through the terminator sequence but the level of transcription decreased. The degree of transcriptional change that occurs during quantitative PCR suggests that RNA polymerase continues past the terminator sequence and gradually release the RNA, as demonstrated by transcriptome analysis [10]. As the termination efficiency of an intrinsic terminator is directly related to the stability of the stem-loop structure [4] the relatively high ΔG values in intrinsic terminators could be responsible for the read-through of some M. hyopneumoniae terminator sequences. Supporting our findings, the presence of functional intrinsic terminators was previously demonstrated for genes in other mycoplasma species such as hmw of M. pneumoniae, vmpaU of Mycoplasma agalactiae and for the operon ftsZ of M. genitalium [18-20].

Conclusion

In silico predictions, combined with experimental analysis, confirmed the presence of intrinsic terminators in M. hyopneumoniae genome. The Figure 8 summarizes the current knowledge related to transcription in M. hyopneumoniae. The localisation of terminators sequences at the 3’ end of mCs and pCs supports previous findings that the M. hyopneumoniae genome is organised preferentially in TUs containing two or more genes [9,10]. Moreover, the experimental results suggest that at least some terminators have a functional role in mycoplasma. Therefore, the presence of terminator sequences associated with the identification of promoter sequences in M. hyopneumoniae transcriptional units (mC and pCs) [12] suggest that transcription runs on from an upstream promoter and terminates at specific stem-loop regions.

Abbreviations

IR:

Intergenic region

CDS:

Coding sequences

ORF:

Open reading frame

tc1 :

class 1 terminator

tc2 :

class 2 terminator

tc3 :

class 3 terminator

tc4 :

class 4 terminator

ΔG:

Gibbs free energy

RT-PCR:

Reverse transcriptase PCR

qPCR:

quantitative real time PCR

TU:

Transcription unit

pC:

polycistronic unit

mC:

monocistronic unit

CT:

Threshold cycle

References

  1. Peters JM, Vangeloff AD, Landick R. Bacterial transcription terminators: the RNA 3'-end chronicles. J Mol Biol. 2011;412(5):793–813.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. d'Aubenton Carafa Y, Brody E, Thermes C. Prediction of rho-independent Escherichia coli transcription terminators. A statistical analysis of their RNA stem-loop structures. J Mol Biol. 1990;216(4):835–58.

    Article  PubMed  Google Scholar 

  3. de Hoon MJ, Makita Y, Nakai K, Miyano S. Prediction of transcriptional terminators in Bacillus subtilis and related species. PLoS Comput Biol. 2005;1(3):e25.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Mitra A, Angamuthu K, Jayashree HV, Nagaraja V. Occurrence, divergence and evolution of intrinsic terminators across eubacteria. Genomics. 2009;94(2):110–6.

    Article  CAS  PubMed  Google Scholar 

  5. Thacker EL. Diagnosis of Mycoplasma hyopneumoniae. Anim Health Res Rev. 2004;5(2):317–20.

    Article  PubMed  Google Scholar 

  6. Vasconcelos AT, Ferreira HB, Bizarro CV, Bonatto SL, Carvalho MO, Pinto PM, et al. Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae. J Bacteriol. 2005;187(16):5568–77.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Liu W, Feng Z, Fang L, Zhou Z, Li Q, Li S, et al. Complete genome sequence of Mycoplasma hyopneumoniae strain 168. J Bacteriol. 2011;193(4):1016–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Siqueira FM, Thompson CE, Virginio VG, Gonchoroski T, Reolon L, Almeida LG, et al. New insights on the biology of swine respiratory tract mycoplasmas from a comparative genome analysis. BMC Genomics. 2013;14:175.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Siqueira FM, Schrank A, Schrank IS. Mycoplasma hyopneumoniae transcription unit organization: genome survey and prediction. DNA Res. 2011;18(6):413–22.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Siqueira FM, Gerber AL, Guedes RL, Almeida LG, Schrank IS, Vasconcelos AT, et al. Unravelling the transcriptome profile of the Swine respiratory tract mycoplasmas. PLoS One. 2014;9(10):e110327.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Weber SS, Sant'Anna FH, Schrank IS. Unveiling Mycoplasma hyopneumoniae promoters: sequence definition and genomic distribution. DNA Research. 2012;19(2):103–15.

    Article  PubMed Central  CAS  Google Scholar 

  12. Siqueira FM, de Souto WS, Cattani AM, Schrank IS. Genome organization in Mycoplasma hyopneumoniae: identification of promoter-like sequences. Mol Biol Rep. 2014;41(8):5395–402.

    Article  CAS  PubMed  Google Scholar 

  13. Muto A, Osawa S. The guanine and cytosine content of genomic DNA and bacterial evolution. Proc Natl Acad Sci U S A. 1987;84(1):166–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Washio T, Sasayama J, Tomita M. Analysis of complete genomes suggests that many prokaryotes do not rely on hairpin formation in transcription termination. Nucleic Acids Res. 1998;26(23):5456–63.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Musatovova O, Dhandayuthapani S, Baseman JB. Transcriptional starts for cytadherence-related operons of Mycoplasma genitalium. FEMS Microbiol Lett. 2003;229(1):73–81.

    Article  CAS  PubMed  Google Scholar 

  16. Kingsford CL, Ayanbule K, Salzberg SL. Rapid, accurate, computational discovery of Rho-independent transcription terminators illuminates their relationship to DNA uptake. Genome Biol. 2007;8(2):R22.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Mitra A, Kesarwani AK, Pal D, Nagaraja V. WebGeSTer DB–a transcription terminator database. Nucleic Acids Res. 2011;39(Database issue):D129–35.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Waldo RH, Popham PL, Romero-Arroyo CE, Mothershed EA, Lee KK, Krause DC. Transcriptional analysis of the hmw gene cluster of Mycoplasma pneumoniae. J Bacteriol. 1999;181(16):4978–85.

    PubMed Central  CAS  PubMed  Google Scholar 

  19. Benders GA, Powell BC, Hutchison CA. Transcriptional analysis of the conserved ftsZ gene cluster in Mycoplasma genitalium and Mycoplasma pneumoniae. J Bacteriol. 2005;187(13):4542–51.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Chopra-Dewasthaly R, Citti C, Glew MD, Zimmermann M, Rosengarten R, Jechlinger W. Phase-locked mutants of Mycoplasma agalactiae: defining the molecular switch of high-frequency Vpma antigenic variation. Mol Microbiol. 2008;67(6):1196–210.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Naville M, Ghuillot-Gaudeffroy A, Marchais A, Gautheret D. ARNold: a web tool for the prediction of Rho-independent transcription terminators. RNA Biol. 2011;8(1):11–3.

    Article  CAS  PubMed  Google Scholar 

  22. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.

    Article  CAS  PubMed  Google Scholar 

  23. Friis NF. Some recommendations concerning primary isolation of Mycoplasma suipneumoniae and Mycoplasma flocculare a survey. Nordisk Veterinaer Medicin. 1975;27(6):337–9.

    CAS  Google Scholar 

  24. Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37(6):e45.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  26. Madsen ML, Puttamreddy S, Thacker EL, Carruthers MD, Minion FC. Transcriptome changes in Mycoplasma hyopneumoniae during infection. Infect Immun. 2008;76(2):658–63.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Arrebola E, Carrión VJ, Cazorla FM, Pérez-García A, Murillo J, de Vicente A. Characterisation of the mgo operon in Pseudomonas syringae pv. syringae UMAF0158 that is required for mangotoxin production. BMC Microbiol. 2012;12:10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Unniraman S, Prakash R, Nagaraja V. Alternate paradigm for intrinsic transcription termination in eubacteria. J Biol Chem. 2001;276(45):41850–5.

    Article  CAS  PubMed  Google Scholar 

  29. Ingham CJ, Hunter IS, Smith MC. Rho-independent terminators without 3' poly-U tails from the early region of actinophage øC31. Nucleic Acids Res. 1995;23(3):370–6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Taschke C, Herrmann R. Analysis of transcription and processing signals of the 16S-23S rRNA operon of Mycoplasma hyopneumoniae. Mol Gen Genet. 1986;205(3):434–41.

    Article  CAS  PubMed  Google Scholar 

  31. Inamine JM, Loechel S, Hu PC. Analysis of the nucleotide sequence of the P1 operon of Mycoplasma pneumoniae. Gene. 1988;73(1):175–83.

    Article  CAS  PubMed  Google Scholar 

  32. Inamine JM, Loechel S, Collier AM, Barile MF, Hu PC. Nucleotide sequence of the MgPa (mgp) operon of Mycoplasma genitalium and comparison to the P1 (mpp) operon of Mycoplasma pneumoniae. Gene. 1989;82(2):259–67.

    Article  CAS  PubMed  Google Scholar 

  33. Krause DC, Proft T, Hedreyda CT, Hilbert H, Plagens H, Herrmann R. Transposon mutagenesis reinforces the correlation between Mycoplasma pneumoniae cytoskeletal protein HMW2 and cytadherence. J Bacteriol. 1997;179(8):2668–77.

    PubMed Central  CAS  PubMed  Google Scholar 

  34. Halbedel S, Stülke J. Tools for the genetic analysis of Mycoplasma. Int J Med Microbiol. 2007;297(1):37–44.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from the Brazilian National Research Council (CNPq) and the Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS). TEF and FMS received scholarships from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Irene Silveira Schrank.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TEF, FMS and ISS conceived and designed the experiments. In silico experiment: TEF. RT-PCR and qPCR experiments: TEF and FMS. Data analyse and writing the manuscript: TEF, FMS and ISS. All authors read and approved the final manuscript.

Additional files

Additional file 1:

List of primers used in RT-PCR and qPCR, showing the sequence, melting point and product size of the primers pairs used in RT-PCR and qPCR. Amplification efficiency was calculated by LinRegPCR for primers used in qPCR.

Additional file 2:

Classes 1, 2, 3 and 4 putative terminators predicted in M. hyopneumoniae . Characteristics of predicted terminators in the region 3’ end after the coding sequence. These terminators were classified in classes 1, 2, 3 and 4 according to the criteria described in the text.

Additional file 3:

Mapping of terminators predicted in monocistronic units (mC) of M. hyopneumoniae . Presence of terminators class 1 in monocistronic units according the genome organisation of M. hyopneumoniae.

Additional file 4:

Mapping of terminators at the end of polycistronic units and in internal genes of M. hyopneumoniae genome. Presence of terminators class 1 in polycistronic units according the genome organisation of M. hyopneumoniae.

Additional file 5:

Definition of class 2 terminators. The scheme represents a region of 101600 at 111600 bp in M. hyopneumoniae genome with the polycistronic units 10 (nox) and 11 (deaA). TU_10 showed three predicted terminators in the same position; two were predicted by TransTermHP (red) and one by WebGesTer (yellow). TU_11 has three predicted terminators in the same position that were predicted by all three software algorithms. The terminator in orange was predicted by ARNold.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fritsch, T.E., Siqueira, F.M. & Schrank, I.S. Intrinsic terminators in Mycoplasma hyopneumoniae transcription. BMC Genomics 16, 273 (2015). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-015-1468-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-015-1468-6

Keywords