Open Access

The katG mRNA of Mycobacterium tuberculosis and Mycobacterium smegmatis is processed at its 5' end and is stabilized by both a polypurine sequence and translation initiation

BMC Molecular Biology20089:33

DOI: 10.1186/1471-2199-9-33

Received: 06 November 2007

Accepted: 04 April 2008

Published: 04 April 2008

Abstract

Background

In Mycobacterium tuberculosis and in Mycobacterium smegmatis the furA-katG loci, encoding the FurA regulatory protein and the KatG catalase-peroxidase, are highly conserved. In M. tuberculosis furA-katG constitute a single operon, whereas in M. smegmatis a single mRNA covering both genes could not be found. In both species, specific 5' ends have been identified: the first one, located upstream of the furA gene, corresponds to transcription initiation from the furA promoter; the second one is the katG mRNA 5' end, located in the terminal part of furA.

Results

In this work we demonstrate by in vitro transcription and by RNA polymerase Chromatin immunoprecipitation that no promoter is present in the M. smegmatis region covering the latter 5' end, suggesting that it is produced by specific processing of longer transcripts. Several DNA fragments of M. tuberculosis and M. smegmatis were inserted in a plasmid between the sigA promoter and the lacZ reporter gene, and expression of the reporter gene was measured. A polypurine sequence, located four bp upstream of the katG translation start codon, increased beta-galactosidase activity and stabilized the lacZ transcript. Mutagenesis of this sequence led to destabilization of the mRNA. Analysis of constructs, in which the polypurine sequence of M. smegmatis was followed by an increasing number of katG codons, demonstrated that mRNA stability requires translation of at least 20 amino acids. In order to define the requirements for the 5' processing of the katG transcript, we created several mutations in this region and analyzed the 5' ends of the transcripts: the distance from the polypurine sequence does not seem to influence the processing, neither the sequence around the cutting point. Only mutations which create a double stranded region around the processing site prevented RNA processing.

Conclusion

This is the first reported case in mycobacteria, in which both a polypurine sequence and translation initiation are shown to contribute to mRNA stability. The furA-katG mRNA is transcribed from the furA promoter and immediately processed; this processing is prevented by a double stranded RNA at the cutting site, suggesting that the endoribonuclease responsible for the cleavage cuts single stranded RNA.

Background

mRNA decay is known to play an important role in the post-transcriptional regulation of gene expression in bacteria. Studies in Escherichia coli indicated a possible model, in which the endoribonuclease RNase E binds and cuts the mRNA, followed by RNase II and/or PNPase that degrade the fragments generated by RNase E cleavage [15]. The accessibility of the 5' end of an mRNA to RNase E has been shown to be an important mediator of stability in E. coli [6, 7].

Much less is known about mRNA decay in other bacteria, particularly in gram positive species. Several studies, performed in Bacillus subtilis and Bacillus thuringiensis, indicated that the 5' terminal leader sequence has a relevant role in mRNA stabilization, either by containing a polypurine rich sequence or by the presence of a predicted 5'-terminal stem-loop structure [812]. Extensive mutagenesis was used to define how these two motifs influence mRNA stability and the presence of an RNase E-like ribonuclease, acting on the 5' end, has been suggested. Its action may be blocked by either the stem-loop structure and/or ribosomes located in proximity of the 5' terminus of the transcript [8, 12].

In mycobacteria few examples of mRNAs stability control have been described [1315]. The most interesting is the stability of the DNA gyrase mRNA in Mycobacterium smegmatis, which is 5' protected by a stem-loop structure, followed by a Shine-Dalgarno sequence for translation initiation. Disruption of the stem-loop caused loss of transcript stability [15].

In this work we report that the mycobacterial katG mRNA originates from specific processing of the longer bicistronic transcript, covering furA and katG. The katG mRNA encodes the catalase-peroxidase, responsible for the degradation of toxic oxygen compounds. In all studied mycobacterial species, the DNA region immediately upstream of katG encodes FurA [1619]. Fur-like proteins are ubiquitous in bacteria, and act as transcriptional repressors that exhibit an iron-dependent DNA binding activity and regulate several genes involved in iron metabolism [20, 21]. In Mycobacteria this role is not achieved by FurA, but by the iron dependent transcription regulator (IdeR) which plays a critical role in maintaining the intracellular iron homeostasis [22].

In recent years, genes directly or indirectly controlled by Fur-like proteins have been discovered, leading to the hypothesis that Fur proteins have a wider role in bacterial gene expression [23, 24].

In mycobacteria, the organization of the furA-katG region is conserved, suggesting that it may constitute a single operon [16, 17, 19]. This hypothesis was confirmed by identification of a single furA-katG transcript in Mycobacterium bovis BCG [18, 25]. In apparent contrast, two different 5' ends were identified by S1 mapping: the first coincides with the initiation codon of furA and is generated by transcription initiation at the furA promoter (Fig. 1A; [18, 25, 26]); the second is 54 bp upstream of katG translation start codon. Master et al. [25] proposed that a katG promoter was present in this region (Fig. 1A and 1B).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig1_HTML.jpg
Figure 1

The furA - katG region in M. tuberculosis and M. smegmatis. Panel A. Schematic representation of the furA-katG region. Genes are indicated by open boxes, the polypurine sequences (PPS) by a black box. The coordinates of the M. tuberculosis sequence are arbitrary, with coordinate +1 corresponding to the first codon of furA (coordinate 2,156,590 in the reverse strand of the M. tuberculosis complete genome sequence deposited in the GenBank:AL123456). Coordinates of the M. smegmatis sequence are arbitrary, with coordinate +1 corresponding to one bp immediately upstream of the furA start codon (GenBank:NC_008596, [42]). The position of the 5' ends of katG mRNA (+428 in M. tuberculosis, [25]; 443 and 478 in M. smegmatis, [27]) are indicated by vertical arrows below the map. Panel B. Sequence of part of the furA-katG region of M. tuberculosis and M. smegmatis. The 5' ends of the katG mRNA are highlighted in grey, the stop codon of furA and the start codon of katG are in bold. The polypurine sequence upstream of katG is underlined. Panel C. Complementarity between the 16S rRNA of M. tuberculosis and M. smegmatis and the PPS upstream of katG initiation codon. The PPS is in bold and the translation initiation codon is boxed. The accession number of 16S rRNA of M. tuberculosis is GenBank:NC_000962.2 [16] and of M. smegmatis is GenBank:NC_008596 [42].

In M. smegmatis specific transcripts for furA and katG were identified by Northern blotting, but a single transcript encompassing both furA and katG genes was not detected [27]. Specific 5' ends were mapped by primer extension analysis [26, 27]: the first (coordinate 1 in Fig. 1A) was due to transcription initiation from the furA promoter; two additional 5' ends were mapped upstream of katG, the first one, corresponding to coordinate 443, falls in the terminal part of the furA gene; the second is less intense and corresponds to coordinate 478, located in the furA-katG intergenic region. No specific promoter sequences were identified upstream of either 443 or 478.

In this work we addressed the question whether the katG 5' ends were generated by transcription initiation or by mRNA processing. Our results support the latter alternative. Furthermore, a polypurine sequence and translation initiation are both involved in katG mRNA stability.

Results

The 5' end of katG mRNA is not generated by initiation of transcription

In a previous work, Milano et al. [27] reported the presence of a promoter region for katG expression in M. smegmatis, located in the terminal part of the upstream furA gene. In order to demonstrate the presence of this promoter, we amplified the M. smegmatis region that overlaps the 5' ends upstream of katG (coordinate 437–552) and cloned it in plasmid pMV261, substituting the hsp60 promoter region (pMYS648). In vitro transcription was performed in the presence of [α32P]CTP and M. smegmatis cell extract, comparing the transcripts synthesized by pMV261, from the constitutive hsp60 promoter, and pMYS648, from the potential katG promoter (Fig. 2). pMV261 produced an RNA of about 300 nt that may represent the RNA starting from the hsp60 promoter and terminating at the transcription terminator downstream. On the contrary, no transcript of the expected size for an mRNA starting at the hypothetical pkatG and terminating at the transcription terminator (about 200 nt) was detected when pMYS648 was used, suggesting that no promoter was present in this region. This data is in agreement with the absence of any promoter consensus sequences upstream of the 5' ends.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig2_HTML.jpg
Figure 2

In vitro transcription. In vitro transcription was performed with 10 μg of crude M. smegmatis extract in the presence of [α-32P]-CTP. The templates were plasmid pMV261, which carries the M. bovis BCG hsp60 promoter, and plasmid pMYS648, in which the 437–552 M. smegmatis region has been substituted to the hsp60 promoter region. The samples were run on a 5% denaturing polyacrylamide gel. Molecular weight markers (MWM) run in the same gel are indicated on the right. The about 300 bp transcript synthesized by pMV261 is indicated by an arrow.

In order to exclude that a katG promoter was present in this region, we performed an RNA polymerase Chromatin Immunoprecipitation assay (RNApol-ChIP) on M. smegmatis mc2155. After rifampicin treatment in order to block RNA polymerase on promoter sequences [28], formaldehyde cross-linking, and sonication, the DNA-RNA polymerase complexes were immunoprecipitated with mouse monoclonal antibodies to the β subunit of E. coli RNA polymerase, as described in the Methods section. As rifampicin inhibits transcript elongation, we expected RNA polymerase to be trapped at any potential promoter. The precipitated DNA was then analyzed by quantitative RT-PCR to determine the amount of recovered DNA.

Specific oligonucleotides were used to amplify different 150–200 bp DNA fragments of the furA-katG region (fragments A, B, C, and D in Fig. 3). The amplified regions with the indication of the ratio relative to the amount of the fragment containing the pfurA promoter (fragment B) are reported in Fig. 3. Compared to fragment B, all the other regions were recovered with low frequency: in particular, the region where the potential katG promoter had previously been located (fragment C) was recovered with a relative value of 0.10. Thus, we concluded that a promoter was not present in this region and that the katG mRNA 5' end was produced by processing of a longer transcript.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig3_HTML.jpg
Figure 3

Results of RNAPol-ChIP experiments. Schematic representation of the oxidoreductase-furA-katG region of M. smegmatis. The genes are indicated by boxes and the arrow corresponds to the furA promoter. The lines below the map correspond to different regions amplified by a pair of oligonucleotides. Fragment A: -689/-492; B: -81/91 (furA promoter); C: 374/557; D: 1155/1344. Below the lines, the values of the ratio of each fragment relative to the pfurA DNA.

The polypurine sequence preceding the katG start codon stabilizes the downstream transcript

In Bacillus subtilis and Bacillus thuringiensis polypurine sequences (PPS) in the 5' terminal part of several transcripts were reported to increase mRNA stability [911]. By sequence analysis of the region immediately upstream of katG in M. tuberculosis and M. smegmatis we identified a PPS four bases upstream of the katG translation start codon (Fig. 1B; GGAAGGAA at coordinates 470–477 in M. tuberculosis; GAAAGGAAA at coordinates 487–495 in M. smegmatis). We tested whether the PPS in the mycobacterial katG mRNA was also able to stabilize the transcript.

An integrative plasmid was constructed, pMYS694, that carries the M. smegmatis sigA promoter [29] upstream of the lacZ reporter (Fig. 4A). Then, either the 437–531 region of M. smegmatis (pMYS690) or the 405–531 region of M. tuberculosis (pMYT733) were cloned between the sigA promoter and the lacZ gene. These regions cover the terminal part of the furA gene, the furA-katG intergenic region with the wild type PPS and the first codons of katG (Fig. 4A). Two mutant plasmids were also constructed, in which the PPS sequence was mutagenized (pMYS692 and pMYT735 carry CTCTGGAGG and CCTCCCTC, respectively; see Fig. 4A).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig4_HTML.jpg
Figure 4

Effect of the polypurine sequence on RNA stability. Panel A. pMYS694 carries the M. smegmatis -195/+27 DNA region (where +1 is the transcription initiation site from psigA [29]), containing the sigA promoter, upstream of lacZ in pSM128. The lacZ gene is preceded by its ribosome binding site (44 nt from the cloning site). pMYS690 carries the 437–531 region of M. smegmatis downstream of psigA, containing the wild type PPS (GAAAGGAAA); pMYS692 carries the same region with the indicated mutations in the PPS (CTCTGGAGG). pMYT733 and pMYT735 carry the 405–531 region of M. tuberculosis downstream of psigA. In pMYT733 carries the wild type PPS (GGAAGGAA), whereas in pMYT735 the PPS was substituted by CCTCCCTC. Beta-galactosidase activity and the half life of the lacZ mRNA were determined as indicated in the Methods section. Average beta-galactosidase activity in Miller Units of three to eight different replicas ± standard deviation is reported. The average relative increment in beta-galactosidase activity and in lacZ mRNA half life are also indicated. N.T. = not tested. Panel B. RNA was extracted from M. smegmatis strains mc2155(pMYS690) and mc2155(pMYS692), and primer extension performed with oligonucleotide 809, internal to lacZ, as described in Methods. The DNA sequence obtained with the same oligonucleotide on plasmid pMYS690 was run in the same gel.

M. smegmatis mc2155 was transformed, and beta-galactosidase activity measured. As reported in Fig. 4A, beta-galactosidase activity expressed by pMYS694 was increased 3.5-fold in pMYS690 and 3.9-fold in pMYT733, which carry the wild type PPS. On the contrary, beta-galactosidase expressed from pMYS692 and pMYT735, carrying the mutated PPS, did not change significantly compared to the control pMYS694 (0.7 and 1.6 fold, respectively).

Furthermore, we measured the half-life of lacZ mRNA from the plasmids carrying M. smegmatis fragments. The mc2155 strains carrying pMYS694/690/692 were treated with rifampicin (see Methods) to block transcription and the RNAs were extracted from the cultures at 0, 2, 4, 10, and 30 minutes after the addition of the antibiotic. Quantitative RT-PCR was performed in order to measure the amount of RNA at the different time points and to evaluate the half lives in the three strains. The half life of the lacZ mRNA was increased from 0.4 minutes in pMYS694 to 2.2 minutes in pMYS690 (5.5 fold increase), whereas no substantial increase was observed in pMYS692, which carries the mutated PPS. These results suggest that the identified PPS is able to stabilize the katG-lacZ mRNA. Moreover, we observed a good correlation between the beta-galactosidase activity and the lacZ mRNA stability.

Primer extension analysis was performed on the RNA extracted from mc2155(pMYS690) and (pMYS692) (Fig. 4B). Two bands were detected in pMYS690: one corresponding to the previously mapped 5' end of sigA [29], and the other corresponding to the 5' end of katG (443 [27]) in M. smegmatis. The band at 443 appears very strong, in accordance with the stabilizing effect of the downstream PPS. The two signals were also present in pMYS692, but in this case the 443 band was visibly fainter, according to the decreased stability caused by the mutated PPS.

Relevance of translation in katG mRNA stability

In other gram positive bacteria it has been shown that ribosome binding and ternary complex formation are responsible for the PPS mediated stability, whereas translation appears not to be involved [11, 31, 32]. On the other hand, in many E. coli RNAs translation seems to contribute directly to transcript stability [3335].

In order to evaluate the role of translation in increasing mRNA stability in our system, we constructed several plasmids in which DNA fragments of M. smegmatis were inserted in pMYS694 between the sigA promoter and the lacZ reporter gene (Fig. 5). pMYS690 carries an M. smegmatis region that starts from coordinate 437 and extends to the first 23 codons of katG, followed by a translation stop codon (TGA). pMYS727 carries the same DNA fragment, but the katG translation start codon was mutagenized to CCG.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig5_HTML.jpg
Figure 5

Role of translation in mRNA stability. Beta-galactosidase expression of M. smegmatis mc2155 transformed with the different plasmids, derived from pMYS694 (in the legend to Fig. 4A), by inserting M. smegmatis katG DNA regions of different length, downstream of the sigA promoter. All the plasmids carry the same region from 437 to the start codon of katG (coordinate 499, see Fig. 1), and differ for the length of the following open reading frame: pMYS690 and pMYS727 carry the same DNA region, but the katG translation start codon TTG has been substituted in pMYS727 with a CCG codon; pMYS720 is mutated in the start codon of katG (TTG to TAA codon); pMYS723 and pMYS744 carry 11 and 19 codons of katG, respectively. Strain mc2155 was transformed with the plasmids and beta-galactosidase expression was measured, as described in Methods. Average beta-galactosidase activity in Miller Units of three to eight different replicas ± standard deviation is reported. The average relative increment in beta-galactosidase activity is also indicated. The amount of lacZ transcript was evaluated by quantitative RT-PCR. N.T. = not tested.

In all these constructs translation of the downstream lacZ gene proceeds independently, by ribosome binding to the lacZ ribosome binding site and starting translation from its ATG start codon. Thus, beta-galactosidase activity and the stability of the lacZ transcript can be compared in the different constructs.

Both the amount of beta-galactosidase and of lacZ transcript were increased more than 3 fold in pMYS690 compared to pMYS694, whereas in pMYS727 relative expression was even lower than in the control. These results suggest that, beside the PPS, translation of the katG open reading frame is also directly involved in transcript stability.

Other plasmids were constructed in which the length of the katG open reading frame was varied: in pMYS720 a stop codon substitutes the katG start codon (TTG to TAA); in pMYS723 and pMYS744 the translation stop codon is located after 11 and 19 codons, respectively. We measured beta-galactosidase activity as an indication of the lacZ mRNA stability. pMYS720 expressed beta-galactosidase activity to a level similar to the pMYS694 control, confirming the relevance of translation for transcript stability. An increased expression of beta-galactosidase could be observed in pMYS723 and in pMYS744, in which translation can proceed for 11 and 19 codons, respectively. These data indicate that translation of about twenty or more codons is necessary for achieving full transcript stability.

A single stranded mRNA region is necessary for the correct processing of the katG mRNA

To characterize the requirements for the processing of the katG mRNA 5' end, several plasmids were constructed by insertion of M. smegmatis DNA regions from 437 to 552, downstream of psigA: pMYS739 carries the wild type sequence, pMYS740 carries a 30 bp internal deletion (from 450 to 480), in pMYS743 five bp flanking the processing site have been substituted, and in pMYS749 five bp were substituted in order to allow the formation of a double stranded region corresponding to the processing site (see Fig. 6). For each plasmid the predicted structure of the RNA from sigA is reported in Fig. 6. RNAs were extracted from mc2155 strain transformed with the plasmids and primer extension was performed with oligo 809, internal to the downstream lacZ gene. The results indicated that in the wild type pMYS739 two signals were detectable: one corresponding to the 5' end due to sigA, the second corresponding to the processing site at 443. The only small effect of the 30 bp deletion in pMYS740 was the presence of a double signal at 442 and 443. The transcript of pMYS743, which carries a different sequence around the processing site, is normally processed at 443, whereas the transcript of pMYS749, in which a double stranded region is created, was no more correctly processed. In order to confirm that a double stranded RNA around 443 could not be processed, we constructed an other plasmid in which different bases were substituted in order to create a double stranded RNA, and also in this case processing at 443 did not occur (data not shown). These data indicate that the ribonuclease responsible for the processing needs a single stranded RNA region for cutting.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2199-9-33/MediaObjects/12867_2007_Article_290_Fig6_HTML.jpg
Figure 6

Structure prediction and primer extension on mRNAs synthesized by different plasmids. The structure of the transcripts starting at the sigA promoter in pMYS739, pMYS740, pMYS743 and pMYS749 were obtained with the m-fold structure prediction program of Zuker [41] and primer extension analyses were performed on RNAs extracted from mc2155 transformed with the different plasmids, with oligonucleotide 809. Sequence reactions obtained with the same oligonucleotide were run in the same gel. The position of the sigA promoter, the 443 processing site and the PPS region are indicated on the right. The arrows indicate the position of the cut on the sequence structure. The mutations are circled.

Discussion

The 5' end of the katG transcript in M. tuberculosis and M. smegmatis is a processing site

In a previous work [27] we suggested the presence of a promoter for M. smegmatis katG expression in the terminal part of the furA gene and we identified two 5' ends at coordinate 443 and 478. However, in this work we tried without success to use the cloned fragment for in vitro transcription: no RNA could be produced by pMYS648, in which the M. smegmatis 437–552 region has been cloned. Moreover, neither promoter consensus sequences were present upstream of 443 nor this region could be recovered by RNA-polymerase ChIP assays. Thus, we concluded that the 5' end is not due to transcription initiation rather to processing of a longer transcript.

In M. tuberculosis a longer transcript, starting from pfurA, has been observed [18, 25]. In M. smegmatis we were unable to detect a polycistronic furA-katG transcript by RT-PCR experiments [27]. This suggests that immediate processing may follow furA-katG transcript synthesis in M. smegmatis.

In the constructs in which the sigA promoter region has been cloned immediately upstream of the M. smegmatis fragment (pMYS690 and pMYS692), we could demonstrate the presence of a longer RNA, starting at psigA. Primer extension analysis also revealed the presence of a 5' end at 443, that is the same position in which the M. smegmatis furA-katG mRNA is cut. However, no primer extension signal was observed at 478. In fact, this latter 5' end was no longer detected in many primer extension experiments performed with RNAs extracted from mc2155 (unpublished data). We suggest that this signal, always less intense than the 5' end at 443, may arise by specific pausing of the reverse transcriptase.

Role of a polypurine sequence in katG mRNA stability

By cloning downstream of the sigA promoter the M. tuberculosis or the M. smegmatis regions that cover the katG gene translation start point, beta-galactosidase expression and lacZ mRNA half life were increased more than threefold. These regions contain a PPS four bp upstream of the katG translation start codon. Mutagenesis of the PPS abolished both effects, thus suggesting the central role of this sequence in transcript stability.

The complementarity of the PPS to the 16S rRNA sequence (Fig. 1C) indicates that it could act as the katG ribosome binding site. Mutations altering the complementarity to the 16S rRNA decreased the amount of the lacZ transcript, suggesting that ribosome binding protects the RNA from degradation.

The role of ribosome binding sites in mRNA stabilization has been demonstrated in many gram-positive and gram-negative bacteria. It was shown that many mRNAs, like lacZ and ompA of E. coli [33, 34, 36], aprE, gsiB and ermC of B. subtilis [31, 32, 37, 38], cryIIIA of B. thuringiensis [11], can be significantly destabilized by mutations in the ribosome binding site that interfere with ribosome binding by reducing its complementarity to 16S rRNA.

This is the first example in mycobacteria, in which a polypurine rich sequence is able to confer stability to the downstream RNA likely by ribosome binding.

Translation is required for stabilization

Most examples of mRNA stability mechanisms derived from gram-positive bacteria seem to depend on ribosome binding, but not on translation [912, 31, 32, 38]. On the other hand, several examples in E. coli indicate that translation has also a role [3336, 39, 40].

We could demonstrate that in our mycobacterial system translation is required for stabilization, and more than 20 codons are needed for conferring full stability to the transcript. Thus, we suggest that the ribonuclease that processes the transcript is impaired both by ribosome bound to the PPS and by translation of the first part of katG.

Regulation of KatG expression and processing of the katG transcript

In previous works [25, 27], by means of two different reporter systems, katG expression was suggested to be originated from specific transcription initiation. Considering our new findings, indicating the absence of any promoter and the presence of a stabilizing sequence upstream of the katG start codon, a new explanation should be given to previous results. We hypothesize that the reporter activity could be explained by the presence of an unknown promoter upstream of the cloning site in the vector sequence. Transcription starting from this promoter may produce a transcript that, upon processing at 443, is stabilized by both the polypurine sequence and translation initiation. This led to the expression of the downstream reporter gene. Thus, a reporter gene assay, usually used to identify promoter sequences, led us to detect an RNA stabilizing sequence.

Therefore, in agreement with what has been suggested by Pym et al. [18], in M. tuberculosis and in M. smegmatis transcription of the katG gene may originate from the upstream pfurA promoter. A single transcript containing furA and katG is synthesized and processed to originate two different transcripts, one covering furA, the second covering katG. The latter is stabilized by the PPS and translation of the first codons of katG. In conclusion, expression of katG appears to be under the direct control of the FurA protein that controls transcription initiation at pfurA [26].

Although the transcript is processed both in M. tuberculosis and in M. smegmatis, the efficiency of the processing event differs. This may be the consequence of a better substrate in M. smegmatis and/or a more active RNAse. It would be possible to replace the M. smegmatis furA-katG intergenic region by the corresponding M. tuberculosis region and test whether a single mRNA for both genes can be found. If this occurs, one can conclude that the substrate makes the difference. On the contrary, if the transcript is completely processed, M. smegmatis RNAse is more active.

The secondary structure of the transcripts covering the M. smegmatis furA-katG region were obtained by mfold-structure prediction program of Zuker [41]. The structure of the region surrounding the processing sites is single stranded both in the transcripts starting at the furA promoter (not shown) and in the RNAs starting at psigA (Fig. 6; pMYS739). Moreover, mutagenesis of this region, either by deletion of about 30 bp or by substitution of the bases surrounding the processing site did not prevent processing at 443. On the contrary, a double stranded RNA region impaired the processing, suggesting that the RNase involved in the processing of the furA-katG transcript requires single stranded RNA.

Several endoribonucleases have been identified in M. tuberculosis and in M. smegmatis complete genome sequence by homology with known enzymes in other bacteria [16, 42]. Experiments are in progress to identify which mycobacterial ribonuclease is responsible for this processing.

Conclusion

This is the first reported case in mycobacteria, in which both a polypurine sequence and translation initiation are shown to contribute to mRNA stability. The furA-katG mRNA is transcribed from the furA promoter and likely processed by a single stranded RNAse.

Methods

Bacterial strains and culture conditions

M. smegmatis mc2155 [43] was grown in LD medium [44] containing 0.2% (vol/vol) glycerol and 0.05% (vol/vol) Tween 80 and supplemented when necessary with spectinomycin (100 μg/ml). Escherichia coli DH10B [45] was used as host strain for cloning and plasmid propagation.

Plasmids

pMV261 is a multicopy mycobacterial vector, carrying the promoter region of the M. bovis BCG hsp60 gene about 300 bp upstream of the rrnABt1 transcriptional terminator [46]. pMYS648 was constructed by inserting the PCR amplified M. smegmatis 437–552 region into pMV261, digested Xba I – Msc I, replacing the hsp60 regulatory region.

pSM128 is an integrative vector carrying the lacZ reporter gene [47]. Plasmids derived from pSM128, containing DNA fragments from either M. tuberculosis or M. smegmatis, were obtained by PCR amplification of the regions, digestion and cloning in the Sca I site, upstream of the lacZ reporter gene. Some of the mutated DNA fragments were obtained by the overlap extension technique [48]. The sequence of the oligonucleotides used for plasmid constructions will be provided upon request. All the inserted fragments were sequenced.

pMYS739, pMYS740, pMYS743 and pMYS749 carry the psigA promoter [29] and the 437–552 region of M. smegmatis upstream of the lacZ reporter gene. pMYS740 carries the 437–552 region of M. smegmatis with a deletion of the 450–480 region; pMYS743 carries the 437–552 region of M. smegmatis with mutagenesis of the nucleotides in position 442–446 (from AACTT to CCAAG) and pMYS749 carries the 437–552 region of M. smegmatis with mutagenesis of the nucleotides in position 452–456 (from TCACA to AGTTC).

Preparation of crude extracts

Crude extracts of M. smegmatis were prepared from 100 ml log phase cultures (OD600 = 0.8). The cells were pelleted, washed in 20 mM TrisHCl pH 8, 10 mM MgCl2, 0.1 mM EDTA, 0.1 mM DTT, 50 mM KCl, and resuspended in 500 μl of 20 mM TrisHCl pH 8, 10 mM MgCl2, 0.1 mM EDTA, 0.1 mM DTT, 50 mM KCl, 1 mM PMSF, 10% glycerol. The cells were lysed by sonication (2 × 20 seconds at 40% amplitude). The cell lysate was spun for 10 min at 15000 × g and the supernatant (about 10 mg/ml) was used as the crude extract.

In vitro transcription assay

Approximately 10 μg of crude extract and 4 μg of CsCl-banded plasmid DNA were incubated for 10 min at 37°C in assay buffer (40 mM Tris pH 7.9, 6 mM MgCl2, 2 mM spermidine, 10 mM NaCl, 100 μg/ml bsa, 20 units of rnasin), followed by addition of UTP, ATP and GTP at 0.4 mM, CTP at 0.2 μM, and 5 μCi [α-32P]-CTP per 20 μl reaction. After incubation for 30 min at 37°C, 80 μl of 200 mM EDTA and 20 μg of yeast tRNA were added and the reaction mixture was extracted with phenol-chloroform and ethanol precipitated. The pellet containing the labelled RNA was resuspended in sample buffer and analysed on a 5% denaturing polyacrylamide gel.

Rifampicin treatment

For RNA half-life measurements, cultures were grown up to OD600 = 0.8, treated with rifampicin (500 μg/ml), and, at different time points, RNA was isolated as described below.

RNA extraction

Fifteen ml of M. smegmatis mc2155 cultures carrying the different plasmids (OD600 = 0.8) were pelleted and resuspended in 100 μl of Tris-EDTA buffer; 75 μl of RNA lysis buffer (4 M guanidinium thiocyanate, 0.01 M Tris pH 7.5, 0.97% β-mercaptoethanol) were added and the suspension was sonicated (20 seconds at 40% amplitude). The RNA was then purified using the SV Total RNA Isolation System according to the manufacturer's protocol (Promega).

Primer extension

Twenty μg of RNA extracted as described above were used in Primer Extension analysis [49] with oligonucleotide 809 (5'-TGCGTTTGTTTGCACGAACC-3') internal to lacZ.

Beta-Galactosidase assay

Independent cultures of M. smegmatis mc2155 carrying the indicated plasmids were grown at 37°C to OD600 = 0.8. Cells were collected by centrifugation, resuspended in 500 μl of TEDP (0.1 M Tris-HCl, 1 mM EDTA, 1 mM DTT and 1 mM PMSF), and disrupted by sonication (2 pulses of 20 seconds and 40% amplitude). Beta-galactosidase activity of the extracts was measured as described in Miller [50]. The enzyme activity was expressed as nanomoles of o-nitrophenol-beta-galactopyranoside converted to o-nitrophenol minute-1 milligram of protein-1. Each experiment was performed three times.

Quantitative RT-PCR

Ten μg of RNA extracted as described above were retrotranscribed with Superscript II Reverse Transcriptase (Invitrogen) using random primers. Double stranded DNA binding dye iQ SYBR Green Supermix (Bio-Rad) in an iCycler iQ real time PCR detection system from Bio-Rad was used to quantify the number of target cDNA molecules in the different samples. Each reaction was run in triplicate and the melting curves were constructed. The relative expression of lacZ mRNA was determined using the Δ (Δ Ct) method [51], with either the 16S rRNA or the sigA mRNA as standard.

RNA pol-ChIP

Ten ml M. smegmatis culture was grown up to OD600 = 0.8, rifampicin was added to a final concentration of 500 μg/ml and the culture incubated for 5 minutes at 37°C with shaking. The culture was then fixed with 1% formaldehyde at room temperature for two minutes under gentle shaking and the reaction stopped by addition of glycine (final concentration 0.125 M) for 10 minutes with gentle shaking. Cross-linked cells were pelleted by centrifugation, washed twice in PBS and once in STET (150 mM NaCl, 10 mM Tris-HCl pH 8, 10 mM EDTA, 0.25% Triton X-100), and finally resuspended in 2 ml of TE (10 mM Tris-HCl pH 8, 1 mM EDTA). The sample was sonicated on ice with 20 pulses of 20 seconds and 40% amplitude. The average DNA fragment size obtained was approximately 0.5 kb. After sonication insoluble debris was removed by centrifugation and the extract precleared with 40 μl 50% protein G-agarose slurry (Sigma, previously blocked with 100 μg/ml lambda DNA, 500 μg/ml tRNA, 1 mg/ml BSA) for 3 h at 4°C. 100 μl of the precleared extract were saved as input fraction. 900 μl were incubated under rotation overnight at 4°C, with 10 μl of 8RB13 mouse monoclonal antibodies to β subunit of E. coli RNA polymerase (Neoclone, cross-reacts with RNA polymerases of many others prokaryotes) in 1× RIPA buffer (140 mM NaCl, 10 mM Tris HCl pH 8, 1 mM EDTA, 1% Triton X-100, 0.1% SDS, 0.1% Na deoxycholate) and RNA polymerase-DNA complexes were immunoprecipitated with 50 μl 50% protein G-agarose slurry for 4 h under rotation at 4°C. The immunocomplex was recovered by centrifugation, washed four times with 1× RIPA buffer, once with LiCl buffer (250 mM LiCl, 10 mM Tris-HCl pH 8, 1 mM EDTA, 0.5% NP-40, 0.5% Na deoxycholate), twice with TE, and resuspended with 100 μl TE (IP fraction). 900 μl of the precleared extract were treated as above but in the absence of the antibodies (no ab fraction). The DNA of the input, IP and no ab fractions was recovered by treatment for 30 min at 37°C with RNase A (20 μg/ml) and overnight digestion with 50 μg/ml proteinase K in 0.5% SDS at 37°C. Crosslinking was reversed for 6 h at 65°C. DNA was extracted once with phenol/chlorophorm, ethanol precipitated and resuspended in 50 μl TE.

The DNA of the samples (1:10 dilution) was analyzed by quantitative PCR, as described above. Four primer pairs, useful to amplify different fragments in the furA/katG region of M. smegmatis, were used. The coordinates of the DNA fragments are reported in the legend to Figure 3.

The amount of the amplified DNA was calculated using the standard curve of each primer pair. The results of the no ab sample were subtracted from those of the IP sample, the resulting figures divided by their corresponding input and finally the result of each primer pair divided with that of the pfurA region pair.

Folding and hybridization prediction

The RNA sequence of the transcripts synthesized by pMYS739, pMYS740, pMYS743 and pMYS749 have been analyzed by the mfold program of Zucker [41].

Declarations

Acknowledgements

We are grateful to Dr. Alberto Danielli for useful suggestions for ChIP experiments, and to Prof. Gianni Dehò for discussing the results. Marcella Paris, Chiara Camisaschi, Elisabetta Prette and Elisabetta Di Florio contributed to some experiments. This work was supported by EC-VI Framework Contract n. LSHP-CT-2005-018923, and by a grant n. 2003059340_002 from MIUR (Italy).

Authors’ Affiliations

(1)
Dipartimento di Scienze Biomolecolari e Biotecnologie, Università degli Studi di Milano
(2)
EPFL-SV-GHI-UPCOL-Station

References

  1. Donovan WP, Kushner SR: Polynucleotide phosphorylase and ribonuclease II are required for cell viability and mRNA turnover in Escherichia coli K-12. Proc Natl Acad Sci USA 1986, 83: 120-124. 10.1073/pnas.83.1.120PubMed CentralView ArticlePubMedGoogle Scholar
  2. Coburn GA, Mackie GA: Degradation of mRNA in Escherichia coli : an old problem with some new twists. Prog Nucleic Acid Res Mol Biol 1999, 62: 55-108.View ArticlePubMedGoogle Scholar
  3. Grünberg-Manago M: Messenger RNA stability and its role in control of gene expression in bacteria and phages. Annu Rev Genet 1999, 33: 193-227. 10.1146/annurev.genet.33.1.193View ArticlePubMedGoogle Scholar
  4. Régnier P, Arraiano CM: Degradation of mRNA in bacteria: emergence of ubiquitous features. Bioessays 2000, 22: 235-244. 10.1002/(SICI)1521-1878(200003)22:3<235::AID-BIES5>3.0.CO;2-2View ArticlePubMedGoogle Scholar
  5. Kushner SR: mRNA decay in Escherichia coli comes of age. J Bacteriol 2002, 184: 4658-4665. 10.1128/JB.184.17.4658-4665.2002PubMed CentralView ArticlePubMedGoogle Scholar
  6. Emory SA, Bouvet P, Belasco JG: A 5'-terminal stem-loop structure can stabilize mRNA in Escherichia coli . Genes Dev 1992, 6: 135-148. 10.1101/gad.6.1.135View ArticlePubMedGoogle Scholar
  7. Bricker AL, Belasco JG: Importance of a 5' stem-loop for longevity of papA mRNA in Escherichia coli . J Bacteriol 1999, 181: 3587-3590.PubMed CentralPubMedGoogle Scholar
  8. Snapper SB, Melton RE, Mustafa S, Kieser T, Jacobs WR Jr: Isolation and characterization of efficient plasmid transformation mutants of Mycobacterium smegmatis . Mol Microbiol 1990, 4: 1911-1919. 10.1111/j.1365-2958.1990.tb02040.xView ArticlePubMedGoogle Scholar
  9. Agaisse H, Lereclus D: Structural and functional analysis of the promoter region involved in full expression of the cryIIIA toxin gene of Bacillus thuringiensis . Mol Microbiol 1994, 13: 97-107. 10.1111/j.1365-2958.1994.tb00405.xView ArticlePubMedGoogle Scholar
  10. Hue KK, Cohen SD, Bechhofer DH: A polypurine sequence that acts as a 5' mRNA stabilizer in Bacillus subtilis . J Bacterio 1995, 177: 3465-3471.Google Scholar
  11. Agaisse H, Lereclus D: STAB-SD: a Shine-Dalgarno sequence in the 5' untranslated region is a determinant of mRNA stability. Mol Microbiol 1996, 20: 633-643. 10.1046/j.1365-2958.1996.5401046.xView ArticlePubMedGoogle Scholar
  12. Sharp JS, Bechhofer DH: Effect of 5'-proximal elements on decay of a model mRNA in Bacillus subtilis . Mol Microbiol 2005, 57: 484-495. 10.1111/j.1365-2958.2005.04683.xView ArticlePubMedGoogle Scholar
  13. Hu Y, Coates AR:Transcription of two sigma 70 homologue genes,sigA and sigB ,in stationary-phase Mycobacterium tuberculosis . J Bacteriol1999,181: 469-476.PubMed CentralPubMedGoogle Scholar
  14. Hu Y, Butcher PD, Mangan JA, Rajandream MA, Coates AR: Regulation of hmp gene transcription in Mycobacterium tuberculosis : effects of oxygen limitation and nitrosative and oxidative stress. J Bacteriol 1999, 181: 3486-3493.PubMed CentralPubMedGoogle Scholar
  15. Unniraman S, Chatterji M, Nagaraja V: A hairpin near the 5' end stabilises the DNA gyrase mRNA in Mycobacterium smegmatis . Nucleic Acids Res 2002, 30: 5376-5381. 10.1093/nar/gkf697PubMed CentralView ArticlePubMedGoogle Scholar
  16. Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, Gordon SV, Eiglmeier K, Gas S, Barry CE III, Tekaia F, Badcock K, Basham D, Brown D, Chillingworth T, Connor R, Davies R, Devlin K, Feltwell T, Gentles S, Hamlin N, Holroyd S, Hornsby T, Jagels K, Krogh A, McLean J, Moule S, Murphy L, Oliver K, Osborne J, Quail MA, Rajandream MA, Rogers J, Rutter S, Seeger K, Skelton J, Squares R, Squares S, Sulston JE, Taylor K, Whitehead S, Barrell BG: Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature 1998, 393: 537-544. 10.1038/31159View ArticlePubMedGoogle Scholar
  17. Pagán-Ramos E, Song J, McFalone M, Mudd MH, Deretic V: Oxidative stress response and characterization of the oxyR-ahpC and furA-katG loci in Mycobacterium marinum . J Bacteriol 1998, 180: 4856-4864.PubMed CentralPubMedGoogle Scholar
  18. Pym AS, Domenech P, Honoré N, Song J, Deretic V, Cole ST: Regulation of catalase-peroxidase (KatG) expression, isoniazid sensitivity and virulence by furA of Mycobacterium tuberculosis . Mol Microbiol 2001, 40: 879-889. 10.1046/j.1365-2958.2001.02427.xView ArticlePubMedGoogle Scholar
  19. Zahrt TC, Song J, Siple J, Deretic V: Mycobacterial FurA is a negative regulator of catalase-peroxidase gene katG . Mol Microbiol 2001, 39: 1174-1185. 10.1111/j.1365-2958.2001.02321.xView ArticlePubMedGoogle Scholar
  20. de Lorenzo V, Wee LVS, Herrero M, Neilands JB: Operator sequences of the aerobactin operon of plasmid ColV-K30 binding the ferric uptake regulation (fur) repressor. J Bacteriol 1987, 169: 2624-2630.PubMed CentralPubMedGoogle Scholar
  21. Escolar L, Pérez-Martín J, de Lorenzo V: Opening the iron box: transcriptional metalloregulation by the Fur protein. J Bacteriol 1999, 181: 6223-6229.PubMed CentralPubMedGoogle Scholar
  22. Rodriguez GM: Control of iron metabolism in Mycobacterium tuberculosis . Trends Microbiol 2006, 14: 320-327. 10.1016/j.tim.2006.05.006View ArticlePubMedGoogle Scholar
  23. Hall HK, Foster JW: The role of fur in the acid tolerance response of Salmonella typhimurium is physiologically and genetically separable from its role in iron acquisition. J Bacteriol 1996, 178: 5683-5691.PubMed CentralPubMedGoogle Scholar
  24. Vasil ML, Ochsner UA: The response of Pseudomonas aeruginosa to iron: genetics, biochemistry and virulence. Mol Microbiol 1999, 34: 399-413. 10.1046/j.1365-2958.1999.01586.xView ArticlePubMedGoogle Scholar
  25. Master S, Zahrt TC, Song J, Deretic V: Mapping of Mycobacterium tuberculosis katG promoters and their differential expression in infected macrophages. J Bacteriol 2001, 183: 4033-4039. 10.1128/JB.183.13.4033-4039.2001PubMed CentralView ArticlePubMedGoogle Scholar
  26. Sala C, Forti F, Di Florio E, Canneva F, Milano A, Riccardi G, Ghisotti D: Mycobacterium tuberculosis FurA autoregulates its own expression. J Bacteriol 2003, 185: 5357-5362. 10.1128/JB.185.18.5357-5362.2003PubMed CentralView ArticlePubMedGoogle Scholar
  27. Milano A, Forti F, Sala C, Riccardi G, Ghisotti D: Transcriptional regulation of furA and katG upon oxidative stress in Mycobacterium smegmatis . J Bacteriol 2001, 183: 6801-6806. 10.1128/JB.183.23.6801-6806.2001PubMed CentralView ArticlePubMedGoogle Scholar
  28. Campbell EA, Korzheva N, Mustaev A, Murakami K, Nair S, Goldfarb A, Darst SA: Structural mechanism for rifampicin inhibition of bacterial RNA polymerase. Cell 2001, 104: 901-912. 10.1016/S0092-8674(01)00286-0View ArticlePubMedGoogle Scholar
  29. Gomez M, Doukhan L, Nair G, Smith I: sigA is an essential gene in Mycobacterium smegmatis . Mol Microbiol 1998, 29: 617-628. 10.1046/j.1365-2958.1998.00960.xView ArticlePubMedGoogle Scholar
  30. Jürgen B, Schweder T, Hecker M: The stability of mRNA from the gsiB gene of Bacillus subtilis is dependent on the presence of a strong ribosome binding site. Mol Gen Genet 1998, 258: 538-545. 10.1007/s004380050765View ArticlePubMedGoogle Scholar
  31. Hambraeus G, Karhumaa K, Rutberg B: A 5' stem-loop and ribosome binding but not translation are important for the stability of Bacillus subtilis aprE leader mRNA. Microbiology 2002, 148: 1795-1803.View ArticlePubMedGoogle Scholar
  32. Sharp JS, Bechhofer DH: Effect of translational signals on mRNA decay in Bacillus subtilis . J Bacteriol 2003, 185: 5372-5379. 10.1128/JB.185.18.5372-5379.2003PubMed CentralView ArticlePubMedGoogle Scholar
  33. Hansen MJ, Chen LH, Fejzo ML, Belasco JG: The ompA 5' untranslated region impedes a major pathway for mRNA degradation in Escherichia coli . Mol Microbiol 1994, 12: 707-716. 10.1111/j.1365-2958.1994.tb01058.xView ArticlePubMedGoogle Scholar
  34. Arnold TE, Yu J, Belasco JG: mRNA stabilization by the ompA 5' untranslated region: two protective elements hinder distinct pathways for mRNA degradation. RNA 1998, 4: 319-330.PubMed CentralPubMedGoogle Scholar
  35. Nogueira T, de Smit M, Graffe M, Springer M: The relationship between translational control and mRNA degradation for the Escherichia coli threonyl-tRNA synthetase gene. J Mol Biol 2001, 310: 709-722. 10.1006/jmbi.2001.4796View ArticlePubMedGoogle Scholar
  36. Joyce SA, Dreyfus M: In the absence of translation, RNase E can bypass 5' mRNA stabilizers in Escherichia coli . J Mol Biol 1998, 282: 241-254. 10.1006/jmbi.1998.2027View ArticlePubMedGoogle Scholar
  37. DiMari JF, Bechhofer DH: Initiation of mRNA decay in Bacillus subtilis . Mol Microbiol 1993, 7: 705-717. 10.1111/j.1365-2958.1993.tb01161.xView ArticlePubMedGoogle Scholar
  38. Drider D, DiChiara JM, Wei J, Sharp JS, Bechhofer DH: Endonuclease cleavage of messenger RNA in Bacillus subtilis . Mol Microbiol 2002, 43: 1319-1329. 10.1046/j.1365-2958.2002.02830.xView ArticlePubMedGoogle Scholar
  39. Deana A, Ehrlich R, Reiss C: Silent mutations in the Escherichia coli ompA leader peptide region strongly affect transcription and translation in vivo. Nucleic Acids Res 1998, 26: 4778-4782. 10.1093/nar/26.20.4778PubMed CentralView ArticlePubMedGoogle Scholar
  40. Deana A, Belasco JG: Lost in translation: the influence of ribosomes on bacterial mRNA decay. Genes Dev 2005, 19: 2526-2533. 10.1101/gad.1348805View ArticlePubMedGoogle Scholar
  41. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 2003, 31: 3406-3415. 10.1093/nar/gkg595PubMed CentralView ArticlePubMedGoogle Scholar
  42. Fleischmann RD, Dodson RJ, Haft DH, Merkel JS, Nelson WC, Fraser CM: Mycobacterium smegmatis str. mc 2 155. National Center for Biotechnology, NIH, Bethesda, MD 20894 USA; 2006.Google Scholar
  43. Snapper SB, Melton RE, Mustafa S, Kieser T, Jacobs WR Jr: Isolation and characterization of efficient plasmid transformation mutants of Mycobacterium smegmatis . Mol Microbiol 1990, 4: 1911-1919. 10.1111/j.1365-2958.1990.tb02040.xView ArticlePubMedGoogle Scholar
  44. Sabbattini P, Forti F, Ghisotti D, Dehò G: Control of transcription termination by an RNA factor in bacteriophage P4 immunity: identification of the target sites. J Bacteriol 1995, 177: 1425-1434.PubMed CentralPubMedGoogle Scholar
  45. Grant SG, Jessee J, Bloom FR, Hanahan D: Differential plasmid rescue from transgenic mouse DNAs into Escherichia coli methylation-restriction mutants. Proc Natl Acad Sci USA 1990, 87: 4645-4649. 10.1073/pnas.87.12.4645PubMed CentralView ArticlePubMedGoogle Scholar
  46. Stover CK, de la Cruz VF, Fuerst TR, Burlein JE, Benson LA, Bennett LT, Bansal GP, Young JF, Lee MH, Hatfull GF, Snapper SB, Barletta RG, Jacobs WR Jr, Bloom BR: New use of BCG for recombinant vaccines. Nature 1991, 351: 456-460. 10.1038/351456a0View ArticlePubMedGoogle Scholar
  47. Timm J, Lim EM, Gicquel B: Escherichia coli -mycobacteria shuttle vectors for operon and gene fusions to lacZ : the pJEM series. J Bacteriol 1994, 176: 6749-6753.PubMed CentralPubMedGoogle Scholar
  48. Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR: Site-directed mutagenesis by overlap extension using the polymerase chain reaction. Gene 1989, 77: 51-59. 10.1016/0378-1119(89)90358-2View ArticlePubMedGoogle Scholar
  49. Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual. 2nd edition. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York, USA; 1989.Google Scholar
  50. Miller JM: Experiments in molecular genetics. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York, USA; 1972.Google Scholar
  51. Ririe KM, Rasmussen RP, Wittwer CT: Product differentiation by analysis of DNA melting curves during the polymerase chain reaction. Anal Biochem 1997, 245: 154-160. 10.1006/abio.1996.9916View ArticlePubMedGoogle Scholar

Copyright

© Sala et al; licensee BioMed Central Ltd. 2008

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement