Scholarly article on topic 'Short RNA indicator sequences are not completely degraded by autoclaving'

Short RNA indicator sequences are not completely degraded by autoclaving Academic research paper on "Veterinary science"

0
0
Share paper
Academic journal
Sci. Rep.
OECD Field of science
Keywords
{""}

Academic research paper on topic "Short RNA indicator sequences are not completely degraded by autoclaving"

OPEN

SUBJECT AREAS:

MICROBIOLOGY TECHNIQUES

APPLIED MICROBIOLOGY

Received 20 September 2013

Accepted 23 January 2014

Published 12 February 2014

Short RNA indicator sequences are not completely degraded by autoclaving

Veena V. Unnithan1, Adrian Unc2,3, Valerisa Joe1 & Geoffrey B. Smith1

department of Biology, New Mexico State University, Las Cruces, NM, USA, 88003, 2Plant and Environmental Sciences, New Mexico State University, Las Cruces, NM, USA, 88003, 3Memorial University of Newfoundland, Boreal Ecosystems Research Initiative, Corner Brook, NL, Canada, A2H 6P9.

Short indicator RNA sequences (<100 bp) persist after autoclaving and are recovered intact by molecular amplification. Primers targeting longer sequences are most likely to produce false positives due to amplification errors easily verified by melting curves analyses. If short indicator RNA sequences are used for virus identification and quantification then post autoclave RNA degradation methodology should be employed, which may include further autoclaving.

Correspondence and requests for materials should be addressed to V.V.U. (vvu85@nmsu.

Quantitative PCR is a commonly used technique in microbiology as it employs a rapid, sensitive and specific detection of microbes that are difficult to culture or are non-culturable. One of the major limitations of this technique is the inability to differentiate between active and inactive viruses. Fully or partially intact genomes of inactivated viruses may produce positive PCR signals1-4. This can lead to over-estimation of the quantity of active viruses in samples.

Autoclaving at 121.1 °C induces cell and viral particle lysis and failure to recover viable organisms is considered as an evidence of its efficiency5. Degradation of genetic material is also assumed to occur and thus often PCR methods are used as verification in experimentation with autoclaved materials including autoclaved waters and wastewaters6. While many of these PCR methods target relatively larger indicator gene sequences, occasionally, most commonly for viruses, the target indicator genetic sequence may be a shorter sequence7-12. If degradation is incomplete, PCR signals due to non-degraded RNA sequences and false positive signals may confound subsequent experimentation and PCR based tests.

Previous reports have indicated that eukaryotic DNA13,14 or enveloped viruses DNA and RNA15,16 might survive standard autoclaving. These reports targeted sequences >100 bps and did not include sequencing verification of amplicons. However indicator viral RNA sequences, especially for non-enveloped viruses, are often short. If the longer sequences were more sensitive to degradation and fragmentation then it would be expected that shorter sequences might preferentially persist.

We have comprehensively verified the post-autoclave persistence of non-degraded RNA sequences of variable lengths for MS2 bacteriophage (ATCC® 15597-B1TM) using quantitative reverse transcriptase PCR (qRT PCR) followed by verification of the amplicons melting behavior, their sequence size and identity. For this we employed 16 primer sets that were expected to produce amplicons in the 70 bp to 2953 bp range and a series of three autoclaving cycles. Following quantitation by qRT PCR putative sequence identity was verified by melting curve analysis and gel electrophoresis amplicon size visualization. Amplicons were cloned into plasmids and their identity further confirmed by sequencing. MS2 is an icosahedral positive sense-, single stranded, non-enveloped RNA coliphages17 and is a common viral surrogate used in environmental water quality research8.

Results

After the 1st autoclave cycle plaque assays were negative, confirming inactivation of MS2. However three commonly used indicator primer sets7 (listed here as A, B and C, Figure 1 and Supplementary materials), producing amplicons in the 70 to 77 bp range, permitted qRT PCR quantification at 0.18%, 0.015% and 0.009% respectively, of the original concentration. After the 2nd and 3rd autoclave cycles the genetic material degraded such that qRT PCR quantification could only be obtained after more than 35 quantification cycles (Cq) (Figure 1.i.). Primer set A targeted the RNA replicase b chain while primer sets B and C targeted assembly proteins on the RNA genome of MS2 bacteriophage7.

After one autoclaving, amplicons in the 100 through 900 bps range (nine primer sets, D, E, F, G, H, I, J, L and O) could either not be quantified via qRT PCR or required over 35 Cq's (Supplementary information, Table 1).

Figure 1 | Quantification and quality control steps for three amplicons that were consistently recovered after repeated autoclaving.

(A), (B) and (C) indicate amplicons commonly employed as indicator sequences for the presence and quantification of MS2 coliphage7. Numeric subscripts indicate the number of autoclave cycles (e.g. Aj, Bj,Cj). If no numeric subscript is present then the capital letter indicates results obtained from non-autoclaved positive controls. j.i. The horizontal line represents the threshold Cq 40. The quantification cycles of A2 & A3, B2 & B3 and C2 & C3 are overlapping. j.iv. Plasmid extracts are indicated by the subscript "p" (e.g. pA1, pB1, pC1). j.v. is a paired comparison of Aj, Bj, Cj (bottom sequence in each pair) to A, B, C (top sequence in each pair). The first two pairs are sequenced in the 3'-5' direction and the third pair in the 5'-3' direction. The gels shown in this image have been cropped for convenience. The borders between separate gels have been marked with a black line.

www.nature.com/scientificreports \ > _

Table 1 | List of tested primer sets. Primer sets A, B and C were adapted from O'Connel et al., 20067. Primers E, G, I, J and O were adapted from Pecson et al., 20093. All other primers were designed by Primer-Blast19 at http://www.ncbi.nlm.nih.gov/tools/primer-blast/. The Tm and GC% were calculated by Oligoanalyzer 3.1 (Integrated DNA Technologies, Inc., Corralville IA). AG's were calculated on the mfold Web server20-21

Forward and reverse primers Amplicon

ID Sequence Tm (°C) GC% Length (bp) Tm (°C) GC% AG (kcal mol 1)

A 5 '-GCTCTGAGAGCGGCTCTATTG-3 ' 60.6 57.1 70 80 60 225.1 ±0)

3 '-CGTTATAGCGGACCGCGT-5 ' 59.97 61.1

B 5 '-TGTGGAGAGACAGGGCACTG-3 ' 61.48 60 77 79 55 226.9 ±0)

3 '-CAGTTGTTGGCCATACGGATT- 5 ' 58.91 47.6

C 5 '-GTCGCGGTAATTGGCGC-3 ' 59.62 64.7 77 79 56 221.9 ±0.6)

3 '-GGCCACGTGTTTTGATCGA- 5' 58.76 52.6

D 5 '-CCAGGTGCCTTCGATGTTCT-3 ' 60.04 55 157 86 62 264.1 ±0.7)

3 '-TCGAGCGATACGAGCAAGAC- 5' 59.97 55

E 5 '-GAAATCACCGACAGCATGAA-3 ' 56.43 45 244 89 64 2120.6 ±1.1)

3 '-AATCCCGGGTCCTCTCTTTA- 5' 57.44 50

F 5 '-TGCTCGTATCGCTCGAGAAC-3 ' 59.97 55 324 88 62 2160.1 ±1.3)

3'-TCCCGGGTCCTCTCTTTAGG- 5' 60.03 60

G 5 '-GGTCGGTGCTTTCATCAGA-3 ' 57.47 52.6 307 83 50 -118.4 ±1.4)

3 '-TGCCCAGAATATCATGGACTC 5' 57.22 47.6

H 5 '-GGCACTTGCCTACTACGGTT-3 ' 60.04 55 325 85 54 2133.7 ±1.4)

3 '-GGTGTATACCGAGACTGCCG- 5 ' 59.97 60

I 5 '-AAGGTGCCTACAAGCGAAGT-3 ' 59.6 50 335 87 58 2141.7 ±1.5)

3 '-TTCGTTTAGGGCAAGGTAGC- 5 ' 57.9 50

J 5 '-ATGTCAGATCCACGCCTCTA-3 ' 57.93 50 287 86 57 2128.1 ±1.5)

3 '-TTCATGCTGTCGGTGATTTC- 5 ' 56.43 45

K 5 '-GACTGGGGCCAAAACGAAAC-3 ' 59.97 55 2952 86 51 21191.2 ± 2.4)

3'-TCGAGCGATACGAGCAAGAC- 5' 59.97 55

L 5 '-ATAGAGCCCTCAACCGGAGT-3 ' 60.03 55 841 84 48 2288.5 ± 2.7)

3 '-AGAACGTGCATTGCCCAAAC 5' 59.97 50

M 5 '-CTGGGGCCAAAACGAAACAG-3 ' 59.97 55 1880 85 51 2722.8 ± 2.7)

3 '-AGAACGTGCATTGCCCAAAC 5' 59.97 50

N 5 '-GCTCCTATGGGGCACAAGTT-3 ' 60.03 55 993 85 51 2403.7 ±3.0)

3 '-GGTGTATACCGAGACTGCCG- 5 ' 59.97 60

O 5 '-GCATGGTTGTCGTCTCTAGGT-3 ' 59.8 52.38 301 84 53 -119.2 ±3.9)

3 '-ACTTTACGTACGCGCCAGTT- 5' 60.32 50

P 5 '-GTTTGGGCAATGCACGTTCT-3 ' 59.97 50 1021 85 52 2412.9 ±4.5)

3 '-GGTGTATACCGAGACTGCCG- 5 ' 59.97 60

Moreover the quantifiable amplicons were shorter than expected (<100 bps and as short as ~50 bps; see Supplementary Figure 3) and generally melted at lower temperatures. While for the non-auto-claved controls these primers produced amplicons with melting peaks (Tm) in the range of 80-90°C, after autoclaving they produced amplicons exhibiting melting peaks in the range of 65-80°C (see Supplementary Materials). The positive controls for the amplicons obtained with these nine primer sets were 100% identical to their expected sequence whereas all the autoclaved controls were unable to be sequenced because of their extremely small size. An exception was a medium sized amplicon (I, 335 bps) that, while barely recoverable after autoclaving (>35 Cq's) could still be recovered as intact sequence after autoclaving. Another four targeted sequences (primer sets, K, M, N and P), aimed at producing amplicons longer than 900 bps, while amplifiable in-vitro, could not be amplified in-vivo even from the positive control.

No non-specific amplification was observed after the 1st autoclave cycle for any of the A, B and C amplicons as described by the peak melting temperature alignment. The Tm of the positive controls and the quantified amplicons from these 3 primer sets were in the same range of 75-85°C (Figure 1. ii.). The amplicon sizes for these three primer sets were confirmed at <100 bps by agarose gel electrophoresis (Figure 1. iii.). After the 2nd and 3rd autoclave cycles recovered amplicons A produced only very diffuse wide bands on the electro-phoresis gels. Diffuse electrophoretic bands were not observed for amplicons obtained from the non-autoclaved control samples.

The three short amplicons A, B and C obtained after the 1st autoclave cycle as well as their respective positive controls were cloned and sequenced. Restriction digests of plasmids confirmed inserts shorter than 100 bp (Figure 1. iv.). Presence of the inserts in the plasmids was also verified by PCR (Figure 1. iv.). Subsequent sequencing confirmed the identity of the amplified sequences (Figure 1. v.). This indicates that short RNA regions were not completely degraded after the 1st cycle of autoclaving and could still be quantified by qRT PCR.

Discussion

These results suggests that short MS2 indicator gene sequences, of less than 100 bp, commonly used to quantify MS2 abundance7-10,12, may lead to false abundance estimates. A combination of culture based plaque assay and possibly targeting of larger indicator sequences in qRT PCR would therefore be recommended for an accurate quantitation of virions in environmental samples. Pecson et al.3 have shown that short heating events, 72°C for 3 min, do not eliminate MS2 amplifiable genomic fragments of 240-340 bp range. However simple PCR amplification does not guarantee an ampli-con's identity. Susceptibility to degradation through heat of different genomic regions varies due to both the particularities ofthe secondary structure and the level of genomic association with proteins3. The genome of MS2 has a complex secondary structure18. Melting temperatures (Tm) and the level of energy required to degrade secondary structures, expressed as changes in the Gibbs free energy (AG), pro-

vide a qualitative thermodynamic understanding of the structure of the amplicons. As expected amplicons'DG values had a significant positive correlation with their Tm (adjusted R2 = 66.8%, p = 0.001), and length (adjusted R2 = 95.4%, p < 0.001). Our results suggest that longer amplicon sequences have a tendency to produce non-specific amplifications evident on electrophoresis gels and through melting curve analyses (see Supplementary materials) but that shorter RNA fragments are not completely degraded by autoclaving and may be recovered intact by molecular amplification techniques. However while most of the tested > 100 bps amplicons were likely fragmented (see gel sizing in Supplementary Figure 3) at least one of the medium sized amplicon (I) was still recoverable at the expected size, supporting some of the earlier reports obtained for eukaryotic DNA13,14 or enveloped viruses DNA and RNA15,16, and with its sequence intact as confirmed by Sanger sequencing (data not shown). Nevertheless all > 100 bps amplicons required > 35 Cq's to be quantified.

Methods

Study design. Ten mL of the MS2 virus stock of 1010 pfu mL"1 was autoclaved three times with 3 mL of the sample removed after each autoclave cycle. Each autoclave cycle was verified to maintain a temperature of 1210C and pressure of 21 psi (144.8 kPa) for 20 minutes. Each autoclave cycle took a total of 50 minutes to complete. Non-autoclaved virus culture was used as positive control. Detection and enumeration of the virus was carried out by plaque assay and qRT PCR. The latter employed multiple primer sets that targeted 16 different RNA regions of the MS2 genome. The identity of the amplified PCR products was verified by cloning and subsequent sequencing. All experiments were carried out in triplicate.

MS2 bacteriophage. MS2 bacteriophage stock was obtained from the laboratory of Dr. C. Gerba at the University of Arizona and propagated using the double layer agar USEPA method 1601 (2001). Escherichia coliHS (pFamp)Rhost strain was purchased from American Type Culture Collection (ATCC-700891).

Primer design. A total of 16 different primer sets were selected based on amplicon length, GC%, and melting temperature to produce amplicons in the 69 to 2953 base pairs range (Methods Table 1). MS2 bacteriophage has a genome size of3569 bp18. Primers were synthesized by Eurofins MWG Operon (Huntsville, AL).

Quantitative reverse transcriptase polymerase chain reaction (qRT PCR). MS2

target copies were estimated by qRT PCR. RNA was extracted with the Invitrogen PureLink RNA/DNA viral extraction kit (Invitrogen, Carlsbad, CA) without the use of carrier RNA according to the manufacturer instructions. Extracts were stored at — 800C and qRT PCR was performed within 48 h of sampling. The iScript TM One-Step qRT PCR kit with SYBR Green (Bio-Rad, Hercules, CA) was used for the assay. Primers were added to a final concentration of 400 nM. The reaction was set up for 25 mL with 5 mL of template RNA. Non-template controls (NTC's) were included, where RNA template was replaced by nuclease free water. Assays were carried out on a Bio-Rad CFX Connect™ detection system (Biorad, Hercules CA). Thermocycling was carried out with an initial cycle of 500C for 10 min to synthesize cDNA and 95о C for 5 min followed by 40 cycles of 950C for 10 s and 600C for 30 s. This was followed by development of melting curves from 650C to 950C with 0.05 seconds incremental steps of 0.50C. All tests were run in triplicate.

RNA standards were prepared from total genomic RNA extracted from a high titer of purified MS2. The RNA extracted was quantified as 40 ng mL—1 with a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA). This was diluted into four 10-fold dilution standards ranging between 1 ng mL—1 and 1 pg mL"1 A genome size of 3,569 nucleotides with an average molecular mass of 330 Da per nucleotide7 was assumed. Thus a 1 ng mL"1 RNA was calculated to contain 5 X 108 copies mL"1 The quantification cycle (Cq) values of all 15 sets (4 set of standards and 1 NTC with 3 replicates each) were averaged using the Bio-Rad's CFX Connect software. The standards were run using primer sets A, B and C with efficiencies of 91.7%, 99.7% and 98.2% and R2 of 0.9928, 0.9991 and 0.9951 respectively. A linear regression master curve was created to correlate the Cq values to copy numbers.

Plaque assay. Infective MS2 virions were quantified by standard double agar layer plaque assay (USEPA Method 1601). Samples were vortexed for 60 s. An aliquot of 500 mL was mixed into 5 mL of 0.7% Tryptic Soy Agar, 0.015% Ampicillin sodium salt, and 0.015% Streptomycin sulphate. To this 100 mL of log-phase growth E. coli ATCC 700891 host were added. The mix, maintained at 460C, was vortexed for 5 s and then overlayed onto a 1.5% Tryptic Soy Agar layer with the same composition of antibiotics. After the agar solidified, it was incubated at 370C for 24 h. Infective MS2 was calculated in terms of plaque forming units (pfu) per mL. Four replicates were tested for each sample.

Cloning and sequencing. The amplified products from qRT PCR were analyzed on 3% agarose gel (50 V for 3 hrs), purified using QIAquick PCR purification kit (Qiagen), as per manufacturer's instructions, followed by elution in 100 mL of sterilized distilled water. Eluate concentrations were measured with a NanoDrop

ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA). Purified PCR products were ligated into pGEM T Easy vector (Promega). The ligation reaction mix contained 7.5 mL of 2X Rapid Ligation buffer, 1 mL of pGEM T easy vector, 5.5 mL of purified PCR product and 1 mL of T4 DNA Ligase with a final volume of 15 mL. Ligations were carried out at 160C in a water bath, overnight. Transformation was carried out with 15 mL ligation mix added to 200 mL of freshly prepared E. coli XL1 Blue competent cells and 3.4 mL of 10% (3-Mercaptoethanol solution. Transformation mixes were kept on ice for 30 min followed by heat shock at 420C for 90 sec and then placed immediately on ice for 5 min. Upon transformation 1000 mL of pre-warmed Luria Bertani (LB) broth was added and the mix was incubated at 370C for 1 K hrs on a horizontal shaker set at 180 rpm. The transformed cells were centrifuged at 4000 rpm for 10 min. The pellet was resuspended in 100 mL of the supernatant and plated on LB agar amended with carbenicillin (100 mg mL_1), 20% Isopropyl b-D-1-thiogalactopyranoside (IPTG), and 2% X-gal. Plates were incubated at 370C for 14-18 h. Colonies picked from plates were used for plasmid extraction (QIAprep Spin Miniprep Kit, Qiagen). Plasmids were quantified on a NanoDrop ND-1000 spectrophotometer. Positive clones were confirmed by visualization of EcoR1 restriction digests of the extracted plasmids on a 3% agarose gel (50 V run for 3 h). Plasmids with a confirmed insert were sequenced at Functional Biosciences (Madison, WI).

1. Katayama, H., Nakamura, M. & Ohgaki, S. [PCR determination of inactivated RNA coliphage Qß] Advances in Water and Wastewater Technology [Matsuo, T., Hanaki, K., Takizawa, S. & Satoh, H. (eds)] [101-108] (Elsevier, Amsterdam, 2001).

2. Blackmer, F., Reynolds, K. A., Gerba, C. P. & Pepper, I. L. Use of integrated cell culture-PCR to evaluate the effectiveness of poliovirus inactivation by chlorine. Appl Environ Microb 66, 2267-2268 (2000).

3. Pecson, B. M., Martin, L. V. & Kohn, T. Quantitative PCR for determining the infectivity of bacteriophage MS2 upon inactivation by heat, UV-B radiation, and singlet oxygen: advantages and limitations of an enzymatic treatment to reduce false-positive results. Appl Environ Microb 75, 5544-5554 (2009).

4. Ma, J. F., Straub, T. M., Pepper, I. L. & Gerba, C. P. Cell culture and PCR determination of poliovirus inactivation by disinfectants. Appl Environ Microb 60, 4203-4206 (1994).

5. Joslyn, L. J. [Sterilization by heat] Disinfection, Sterilisation, and Preservation [Block, S. S. (ed)] [695-728] (Lippincott Williams & Wilkins, Philadelphia, 2001).

6. Limsawat, S. & Ohgaki, S. Fate ofliberated viral RNA in wastewater determined by PCR. Appl Environ Microb 63, 2932-2933 (1997).

7. O'Connell, K. P. et al Real-time fluorogenic reverse transcription-PCR assays for detection of bacteriophage MS2. Appl Environ Microb 72, 478-483 (2006).

8. Dreier, J., Stormer, M. & Kleesiek, K. Use of bacteriophage MS2 as an internal control in viral reverse transcription-PCR assays. J Clin Microbiol 43, 4551-4557 (2005).

9. Julian, T. R., Tamayo, F. J., Leckie, J. O. & Boehm, A. B. Comparison of surface sampling methods for virus recovery from fomites. Appl Environ Microb 77, 6918-6925 (2011).

10. Wolf, S., Hewitt, J., Rivera-Aban, M. & Greening, G. E. Detection and characterization of F1 RNA bacteriophages in water and shellfish: application of a multiplex real-time reverse transcription PCR. J Virol Methods 149, 123-128 (2008).

11. Porter-Jordan, K. et al. Nested polymerase chain reaction assay for the detection of cytomegalovirus overcomes false positives caused by contamination with fragmented DNA. J Med Virol 30, 85-91 (1990).

12. Venkatesan, D. Isolation and Detection of Indicator MS2 Coliphage in different environment and sea foods by PEG Precipitation and GAC-UAPB-RT-PCR method. Advanced BioTech 7, 26-32 (2008).

13. Gefrides, L. A., Powell, M. C., Donley, M. A. & Kahn, R. UV irradiation and autoclave treatment for elimination of contaminating DNA from laboratory consumables. Forensic Sci Int-Gen 4, 89-94 (2010).

14. Dwyer, D. E. & Saksena, N. Failure of ultra-violet irradiation and autoclaving to eliminate PCR contamination. Mol Cell Probe 6, 87-88 (1992).

15. Espy, M. J., Uhl, J. R., Sloan, L. M. & Rosenblatt, J. E. Detection of vaccinia virus, herpes simplex virus, varicella-zoster virus, and Bacillus anthracis DNA by LightCyclerNested polymerase chain reaction after autoclaving: Implications for biosafety of bioterrorism agents in. Mayo Clin Proc 77, 624-628 (2002).

16. Elhafi, G., Naylor, C., Savage, C. & Jones, R. Microwave or autoclave treatments destroy the infectivity of infectious bronchitis virus and avian pneumovirus but allow detection by reverse transcriptase-polymerase chain reaction. Avian Pathol 33, 303-306 (2004).

17. Van Duin, J. & Tsareva, N. [Single-stranded RNA phages] The Bacteriophages [Calendar, R. L. (ed)] [175-196] (Oxford University Press, Oxford, 2006).

18. Fiers, W. R. et al. Complete nucleotide sequence of bacteriophage MS2 RNA: primary and secondary structure of the replicase gene. Nature 260, 500-507 (1976).

19. Ye, J. etal. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC bioinformatics 13, 134 (2012).

20. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic acids res 31, 3406-3415 (2003).

21. SantaLucia, J. A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. P Natl Acad Sci USA 95.4,1460-1465 (1998).

Acknowledgments

This work was supported by the New Mexico State University Agricultural Experiment Station, by the U.S. Department of Energy under contract DE-EE0003046 awarded to the National Alliance for Advanced Biofuels and Bioproducts and by NIH through grant R25GM061222.

Author contributions

V.V.U. and A.U. designed the experiment, V.V.U. and V.J. carried out the experiment, V.V.U. and A.U. analyzed the data, V.V.U., A.U. and G.S. wrote the manuscript. All authors reviewed the manuscript.

Additional information

Supplementary information accompanies this paper at http://www.nature.com/scientificreports

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Unnithan, V.V., Unc, A., Joe, V. & Smith, G.B. Short RNA indicator sequences are not completely degraded by autoclaving. Sci. Rep. 4, 4070; D0I:10.1038/srep04070 (2014).

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported license. To view a copy of this license,

visit http://creativecommons.Org/licenses/by-nc-sa/3.0