Previous Article | Next Article ![]()
Journal of Clinical Microbiology, February 2009, p. 435-440, Vol. 47, No. 2
0095-1137/09/$08.00+0 doi:10.1128/JCM.01247-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Second Department of Virology, National Institute of Infectious Diseases, 4-7-1 Gakuen, Musashimurayama-shi, Tokyo 208-0012,1 Department of Virology, Kitasato Research Center for Environmental Science, 1-15-1 Kitasato, Sagamihara-shi, Kanagawa 228-8555,2 Roche Diagnostics K.K., 2-6-1 Shiba, Tokyo 105-0014, Japan3
Received 2 July 2008/ Returned for modification 8 August 2008/ Accepted 2 December 2008
|
|
|---|
|
|
|---|
Recently, the noroviruses were divided into five genetically distinct genogroups, but the majority of human noroviruses can be divided into two genetically distinct genogroups, genogroup I and genogroup II (GI and GII), which can be further subdivided into at least 14 GI and 17 GII genotypes (7). Norovirus genotype identities are generally maintained across the open reading frames (ORFs). However, a number of norovirus strains have failed to maintain their sequence identities for RNA-dependent RNA polymerase and VP1, and they have been shown to be recombinant (9, 10, 11, 18). Evidence suggests that the recombination site is at the conserved polymerase and capsid junction between ORF1 and ORF2.
The purpose of this study was to develop a more rapid, technically simpler detection and typing method for noroviruses in clinical samples. Data on the sensitivity and specificity of the method are reported, and the applicability of this technology to the clinical diagnosis of noroviruses is discussed.
|
|
|---|
Cloning. Because we had only small quantities of the patient stool samples used in this study, we constructed plasmids for each genotype and used these to optimize the HRM assay. We selected six patient stool samples to use as the reference strains and constructed plasmids.
DNAs containing the target sequences amplified with primer set COG1F-G1-SKR (381 bp) for GI typing and primer set COG2F-G2-SKR (387 bp) for GII typing were cloned into the pCR 2.1 TOPO vector (Invitrogen, Tokyo, Japan) according to the manufacturer's instructions.
Conventional RT-PCR. A conventional RT-PCR was performed so that the findings could be compared with those of the HRM assay. A 25-µl aliquot of purified RNA was added to 3 µl of a reaction mixture containing 1x DNase I buffer and 2 U of DNase (Promega, Tokyo, Japan). This reaction mixture was incubated for 30 min at 37°C, followed by 5 min at 75°C to inactivate the DNase. The DNase-treated RNA (30 µl) was added to a reaction mixture containing 1x SuperScript II RT buffer (Invitrogen, Carlsbad, CA), 2.5 mM dithiothreitol (Invitrogen), 0.4 mM each deoxynucleoside triphosphate (Roche Diagnostics GmbH, Mannheim, Germany), 1 U RNase inhibitor (Toyobo, Tokyo, Japan), and 10 U SuperScript RT II (Invitrogen) in a final reaction volume of 50 µl. Reverse transcription was performed at 42°C for 1 h, followed by inactivation at 72°C for 15 min.
PCR was carried out according to the method described by Kojima et al. (10). For GI, we used sense primer COG1F and antisense primer G1SKR and produced a 381-bp amplicon. For GII, we used sense primer COG2F and antisense primer G2SKR to produce a 387-bp amplicon. PCR was performed with 5 µl of cDNA in a 50-µl reaction mixture containing 33 pmol of each sense and antisense primer, 1x Taq polymerase buffer (Toyobo), 0.2 mM each deoxynucleoside triphosphate, and 2.5 U of Taq polymerase (Toyobo). After initial denaturation at 94°C for 3 min, 40 cycles of amplification were performed using the GeneAmp PCR system 9600 (Perkin-Elmer Applied Biosystems, Foster City, CA). Each cycle consisted of denaturation at 94°C for 1 min, primer annealing at 55°C for 1 min, and extension at 72°C for 2 min, followed by a final extension at 72°C for 15 min. Products were visualized under UV light in a 2% agarose gel stained with ethidium bromide.
Sequence analysis. PCR-generated amplicons were extracted from the gel and purified using the QIAquick gel extraction kit (Qiagen). The nucleotide sequences of the amplified fragments were directly determined with the Terminator cycle sequencing kit (version 1.1) and the ABI 3100-Avant sequencer (Perkin-Elmer ABI, Boston, MA). Sequence analysis was performed using GENETYX (for MAC, version 10.1.4; Genetyx Co., Tokyo, Japan). The sequences were compared, by use of Cluster-X, with those of the reference strains of noroviruses obtained from GenBank, and distances from each of the tested genotypes were calculated using Kimura's two-parameter method (9). Positive-control plasmids were developed by using cDNA, as described above.
HRM curve acquisition and analysis. HRM curve acquisition and analysis were performed on the LightCycler 480 instrument (Roche). Plasmids (1 x 107 copies/well) or cDNA (2 µl/well) was added to 20 µl of a reaction mixture containing 500 nM each primer (G1-SKF-G1-SKR [330 bp] and G2-SKF-G2-SKR [344 bp]), 3.5 mM MgCl2, and 1x LightCycler 480 HRM master mix, which includes ResoLight HRM dye (Roche). Both the GI genotype 8 (GI/8) plasmid and the GII/10 plasmid were used as reference genotypes for the HRM analysis. All samples were tested in duplicate in a 96-well plate at every examination.
PCR cycling for HRM curve acquisition was run under the following conditions: 1 cycle of 95°C for 10 min; 45 cycles of 95°C for 10 s, 62°C to 58°C in 0.5°C/cycle increments for 10 s, and 72°C for 15 s; and 1 cycle of 95°C for 60 s and 40°C for 60 s. Then the fragment was melted by raising the temperature from 60°C to 97°C, with an increment of 0.11°C/s, in order to obtain information on melting profiles.
Melting-curve analysis was performed using the LightCycler 480 gene-scanning software module as reported by Wittwer et al. (20). This software analyzes the HRM curve data to identify changes in the shape of the curve that indicate sequence polymorphisms. Software programs employ a three-step analysis (Fig. 1) as follows. (i) The first step is to normalize the raw melting-curve data by setting the premelt (initial fluorescence) and postmelt (final fluorescence) signals of all samples to uniform values. Premelt signals are uniformly set to a relative value of 100%, while postmelt signals are set to a relative value of 0%, as shown in Fig. 1A. (ii) The second step is to shift the temperature along the x axis of the normalized melting curves at the point where the entire double-stranded DNA is completely denatured (Fig. 1B). (iii) The final step is to further analyze the differences in melting-curve shape by subtracting the curves from a reference curve (also called the "base curve"), thus generating a difference plot curve, which helps to cluster samples into groups that have similar melting curves (i.e., those with the same genotype), as shown in Fig. 1C.
![]() View larger version (33K): [in a new window] |
FIG. 1. Principle of high-resolution amplicon melting-curve analysis to produce fluorescence difference plot curves. Melting-curve acquisitions and analyses were performed on the LightCycler 480 instrument (Roche). PCR products were denatured at 95°C for 60 s, renatured at 40°C for 60 s, and then melted at 60°C to 97°C with 25 signal acquisitions per degree. Melting-curve analysis was then performed using the LightCycler 480 gene-scanning software module. This software analyzes the HRM curve data to identify changes in the shape of the curve (difference plot curve) that indicate sequence polymorphisms. The amplicon melting curve is automatically calculated by the software in a three-step process, shown in the three panels. (A) Normalization. The raw data (left) are normalized to both 100% (at the beginning of the melting) and 0% (at the end of the melting). (Left) Raw data; (right) the same curve after normalization. (B) Temperature shift. (Left) The point (ellipse) at which the entire double-stranded DNA was completely denatured; (right) the same curve after the temperature has been shifted along the x axis. The temperature shifting with the curve allows comparisons between samples. (C) Fluorescence difference plot. (Left) Three points that show the difference between the base curve and the samples; (right) fluorescence difference plot curve obtained by subtracting each curve from the base curve. Because of the shift in the temperature of the curve, the temperature axis no longer reflects absolute temperature but rather reflects temperature differences relative to the superimposed segment of the curves. Data for three samples (a, b, and c) are shown.
|
Sensitivity of HRM analysis. To monitor the sensitivity of the HRM assay, serial 10-fold dilutions (1 x 107 copies/well to 1 copy/well) of the positive-control plasmids were prepared. The HRM reaction was performed in duplicate for each dilution.
|
|
|---|
The difference plot curves of GI/8 for genogroup I and GII/10 for genogroup II were tentatively used as base curves in this study (Fig. 2 and 3). We identified other GI genotypes by comparing their unique subtracted difference plot curves from the HRM analysis with that of GI/8. The subtracted difference plot curves of GI/4 and GI/9 were distinguishable from that of GI/8, as shown in Fig. 2A. For GII genotyping, we selected GII/10 as the base curve and examined two additional samples. We distinguished GII/3 and GII/4 from GII/10 (base curve) (Fig. 2B).
![]() View larger version (31K): [in a new window] |
FIG. 2. Genotyping of noroviruses within a 380-bp fragment using high-resolution amplicon melting-curve analysis. (A) Difference plots from HRM analysis of the norovirus GI genotypes. GI/8 was the reference genotype (base curve) for the HRM analysis. Samples contained 107 copies/reaction. Each sample was measured in duplicate in a 96-well plate at every examination. (B) Difference plots from HRM analysis of norovirus GII genotypes. GII/10 was the reference genotype (base curve) for the HRM analysis. Samples contained 107 copies/reaction.
|
![]() View larger version (42K): [in a new window] |
FIG. 3. Alignment of the nucleotide sequences and difference plots of four different clones of the GII/4 plasmid samples. Each clone of the four GII/4 plasmid samples was measured. GII/10 was the baseline sample for the HRM analysis. (Bottom) Alignment of the nucleotide sequences of four GII/4 plasmids. The positions of primers are shown. (Top) Difference plots from HRM analysis of four GII/4 plasmids compared with GII/10.
|
In order to investigate the detection limits of the HRM assay, the six sample plasmids were examined four to nine times, as shown in Table 1. We found that 10 copies of sample RNA per reaction mixture was the detection limit for all genotypes tested. Based on these observations, we generated an HRM database and analysis schema for rapid identification of the norovirus genotypes used.
|
View this table: [in a new window] |
TABLE 1. End points of detection levels of the HRM assaya
|
|
View this table: [in a new window] |
TABLE 2. Comparative study of HRM with RT-PCR, sequencing, and real-time PCRa
|
![]() View larger version (29K): [in a new window] |
FIG. 4. Results of HRM analysis of the sample from patient 5. Patient 5's sample was identified as GII/4 by RT-PCR and sequencing, as shown in Table 2. To monitor the sensitivity of the HRM assay, serial 10-fold dilutions (1 x 107 copies/well to 1 copy/well) of the GII/4 positive-control plasmids were prepared. The HRM reaction was performed in duplicate for each dilution. In this study, the RNA from patient 5 was identified as GII/4 by the HRM analysis. The examination gave the same pattern as the reference GII/4 plasmid.
|
|
|
|---|
Most of the methods currently available for norovirus genotyping use genotype-specific primers to amplify genes, allowing for the detection of one specific pathogen at a time. For sample screening and genotype identification in an outbreak, multiplexing or multiple reactions followed by post-PCR electrophoresis to separate amplicons of different sizes are usually required. Although real-time RT-PCR has been used successfully for the detection of noroviruses, TaqMan or hybridization probes are usually required, increasing the expense on a cost-per-sample basis. This type of assay is based on a phylogenetic analysis that reveals the preferential association between norovirus genes and coevolution of these genes.
Because the quantity of patient stool sample that we stored was small, we tried to construct plasmids and examined them to clarify the conditions of the HRM assay. We selected six strains from patient stools for reference and constructed plasmids. The HRM assay could directly detect the genotypes of noroviruses with difference plot curve patterns in patient stool samples.
The detection limit of this method is equivalent to 1 to 10 copies of viral RNA in the reaction mixture, and it has a dynamic linear range for quantification across at least 7 log units of RNA concentration.
We demonstrated that noroviruses can be differentiated simply through the use of one PCR primer pair for real-time RT-PCR, followed by HRM of the ORF1-ORF2 junction. PCR products likely contain information for at least partial phylogenetic characterization and genotyping of clinically important noroviruses. Compared to other genotyping methods, HRM-based genotyping gives superior resolving power. When combined with rapid-cycle PCR, HRM analysis requires minimal time, and the cost of materials is typically less than that for standard methods. The time required for the differentiation of norovirus genotypes is considerably shorter when PCR and sequencing are performed directly on clinical samples.
The HRM assay demonstrated good sensitivity, detecting 107 to around 1 copy per reaction mixture of RNA collected directly from patient stool samples for both GI and GII noroviruses. A comparative study with patient feces revealed good correlation with the reference plasmids between the sequencing and HRM results. When we obtained normalized HRM data for a reference plasmid, it was simple to classify the genotype in less than 60 min. Our approach is also a very convenient method to classify noroviruses.
The circulating norovirus genotypes in Japan from 2000 have been reported by the Infectious Agents Surveillance Report to be GI/4 (48%), GI/8 (20%), and GI/9 (12%) for GI and GII/4 (90%), GII/3 (4%), and GII/6 (2%) for GII. In this study, we confirmed that those typed norovirus strains can be identified by HRM analysis. However, untyped noroviruses could not be identified, because the primer pairs commonly used were not suitable.
The conventional RT-PCR and real-time RT-PCR methods are commonly used for the detection of typed noroviruses in outbreaks and sporadic cases in Japan. The primer pairs commonly used for GI and/or GII can be used only for the detection of well-known and typed strains of noroviruses. TEM is still detecting untyped, circulating SRSVs (Table 2). In the present study, we also found five samples that were negative by conventional RT-PCR, real-time PCR, and the HRM assay; nevertheless, they contained SRSVs detected by TEM.
RT-PCR, real-time PCR, and the HRM assay with the commonly used primer pairs are currently unable to detect untyped SRSVs. Further study will need to focus on selecting and identifying untyped caliciviruses. In addition, the broad reactivity of this assay will be validated using a panel of stool samples containing typed GI and/or GII genotypes in the near future.
Published ahead of print on 10 December 2008. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»