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

Faculty of Veterinary Science, The University of Melbourne, Werribee, Victoria 3030, Australia
Received 12 August 2008/ Returned for modification 16 October 2008/ Accepted 18 November 2008
|
|
|---|
|
|
|---|
Diagnosis of FAdV infections can be made from the observation of gross and histopathological changes in the liver and the use of electron microscopy (3, 11) and various serological tests, such as enzyme-linked immunosorbent assay (3, 18, 45), agar gel immunodiffusion, counterimmunoelectrophoresis, indirect hemagglutination, immunofluorescence (3), and Southern hybridization (10, 12). Identification of the serotype(s) involved is very useful for epidemiological tracing and is of critical importance where vaccination is to be used for the control of the disease (3, 20). Typing of the virus conventionally requires isolation in cell culture, followed by a virus neutralization assay (16); however, the implementation of this method is lengthy and labor intensive, and cross-reactivity between serotypes can sometimes render results inconclusive. Tests using PCR together with DNA sequencing (22, 59) and/or restriction enzyme analysis (29, 48, 50) have recently been used for comparatively faster FAdV typing. However, these methods are still time consuming, often require extensive interpretation, and can be relatively expensive to use as a routine typing tool.
Recent research has found that the technique of combining PCR with high-resolution melting (HRM)-curve analysis provides a useful and cost-effective alternative for the direct analysis of genetic variation, particularly when large numbers of samples are to be analyzed. HRM-curve analysis of PCR amplicons has been shown to aid in clinical (52) and/or epidemiological (5, 36) investigations, by providing means for rapid and effective genotyping, variation scanning, microbial detection, and species determination (13, 21, 26, 53). When melted in the presence of a saturating intercalating fluorescent dye, such as LCGreen I or Styo 9 (17, 31, 46), and acquiring fluorescence data over small temperature increments, amplicons containing different sequences can be discriminated based on the melting transition of the PCR product and the resulting melt curve shape.
Commonly, amplicons of fewer than 200 bp have been used for detection of single nucleotide polymorphisms, mutation scanning, and genotyping (26), as slight shifts in the melting domain are easily detected. The use of larger amplicons has been reported to reduce sensitivity and specificity of HRM-curve analysis (6, 24, 25, 60). However, some studies have reported that the melting of larger amplicons, up to 1,000 bp, may produce multiple melting domains and can be used for genotyping (8, 21, 46, 47). The presence of more than one melting domain, as depicted by multiple peaks in the conventional melt curve profile, adds further variation to the normalized HRM-curve profile, thus increasing the power of the HRM-curve analysis technique for sequence differentiation and genotyping (21, 26, 58).
The adenovirus hexon is the major capsid protein of FAdV and contains type-, group-, and subgroup-specific antigenic determinants (18, 28, 35). Several reports have used PCR of different regions of the FAdV viral gene for detection, differentiation, and phylogenetic analysis of FAdVs (29, 48, 50, 59). The hexon loop 1 (L1) region, flanked by pedestal 1 (P1), represents the most variable region and has, when used in PCR coupled with DNA sequencing or restriction enzyme analysis, proved most successful in identifying and differentiating some or all of the 12 FAdV serotypes (2, 29, 30, 37, 48).
The aim of the present study was to investigate the application of HRM-curve analysis of PCR amplicons of different sizes from the hexon gene L1 and P1 region of representative strains of the 12 FAdV serotypes, using the DNA intercalating dye Syto 9, to establish a single closed-tube test method for FAdV serotype differentiation. Both conventional and normalized dissociation plots were generated, and genotyping was applied to differentiate the serotypes, and the results were compared with those of existing FAdV typing methods of PCR and restriction enzyme analysis of the hexon gene and DNA sequencing of PCR amplicons.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. FAdV species, serotypes (according to the current ICTV nomenclature system and the previous European Union and United States systems), strains used in this study, and other known strains within each serotype (19, 28, 48)
|
Extraction of viral DNA. Viral DNA was extracted from virus propagated in cultures and from Australian FAdV vaccine. A 10-µl volume of each virus sample was solubilized in 300 µl of RLT lysis buffer (Qiagen, Valencia, CA) containing 3 µl of β-mercaptoethanol (1.12 g/ml) and incubated at 4°C overnight. A volume of 15 µl of Qiaex II suspension (Qiagen) and 300 µl of 70% ethanol was combined and then added to the sample mixture. This mixture was loaded into a multispin MSK-100 column (Axygen, Inc., Hayward, CA) and placed in a wash tube. Columns were centrifuged for 30 s at 13,500 x g in a microcentrifuge, and the flowthrough was discarded. The matrix was washed once with 600 µl buffer RW1 (Qiagen) and twice with 500 µl buffer RPE (Qiagen). For each wash, columns were centrifuged for 30 s at 13,500 x g and the flowthrough was discarded, and then the columns were centrifuged for 90 s at 19,000 x g to dry the matrix. The matrix was overlaid with 30 µl of diethyl pyrocarbonate-treated water and incubated at room temperature for 1 min, and DNA was eluted by centrifugation of the column for 60 s at 13,500 x g. The DNA elution step was repeated to obtain a final elution volume of 60 µl. The elutions containing extracted viral DNA were used immediately or stored at –20°C prior to use in PCR.
Real-time PCR for amplification of three regions of FAdV hexon gene. Three sets of forward and reverse degenerate oligonucleotide primers (Table 2) binding to the conserved P1 region of the FAdV hexon gene were used to amplify different regions of P1 and L1 (Fig. 1). The HEX-S forward and reverse primers were designed by comparison of the published hexon gene nucleotide sequences of all 12 FAdV serotypes (National Center for Biotechnology Information, GenBank [4], sequence accession numbers AF339914, AF508946, AF508948, AF508951, AF508952, AF508954, AF508955, AF508956, AF508958, AF508959, AF339920, and AF339925) to amplify a region of P1 containing as small as 2-bp changes between serotypes. All primers were manufactured by Invitrogen Pty Ltd.
|
View this table: [in a new window] |
TABLE 2. Primers used in this study
|
![]() View larger version (6K): [in a new window] |
FIG. 1. Schematic presentation of Pl and L1 regions of the FAdV-1 CELO hexon gene and the primers used in this study. The base-pair numbering of the primer positions is based on the FAdV-1 CELO genes for hexon proteins, GenBank accession number Z67970, nt 1380 to 4208 (1).
|
The reaction mixtures were subjected to 94°C for 2 min, and then, for the Hexon A/B and Hex L1 primer sets, 40 cycles of 94°C for 20 s, 62°C (56°C for Hex L1 primer set) for 20 s, and 72°C for 30 s, and for the HEX-S primer set, 50 cycles of 94°C for 15 s, 60°C for 15 s, and 72°C for 20 s. Optical measurements in green channel (excitation at 470 nm and detection at 510 nm) were recorded during the extension step. After completion of PCR cycles, a final extension of 72°C for 2 min was performed.
The PCR products were analyzed by electrophoresis in 1.5% agarose gels stained with GelRed (Biotium, Hayward, CA), and visualized by UV transillumination.
Restriction enzyme analysis. Restriction enzyme analysis of Hexon A/B PCR products was performed as described previously (29), with some modifications. Briefly, 4 µl of PCR product from each reference serotype (except for FAdV-5, which yielded very low-level or no amplification with Hexon A/B primers) were digested for 3 h with 10 U of each restriction enzyme (New England BioLabs) in a total volume of 20 µl. Restriction fragments were separated in 2% agarose gel stained with GelRed and visualized by UV transillumination.
HRM-curve acquisition and analysis. HRM-curve analysis was performed on the Rotor-Gene 6000 using identical volumes of PCR products subjected to a temperature increasing from 80 to 94°C at intervals (ramps) of 0.1°C/s or 0.15°C/s. HRM-curve analysis was performed using Rotor-Gene 6000 1.7.87 software and the HRM algorithm provided. All amplicons were tested in triplicate to detect variations induced by technical errors.
The conventional melt curves were generated automatically and contained one or more peaks, depending on the amplicon tested. To generate the normalized HRM curves, normalization regions were applied to the curves as follows: HEX-S, 84.97 to 85.47 and 90.48 to 91.00; Hex L1, 84.50 to 85.00 and 90.50 to 91.00; and Hexon A/B, 81.00 to 81.40 and 90.50 to 90.90. Genotypes were defined by selecting a representative sample from each serotype. The software then auto-called the genotype of each sample and provided confidence percentages (C%) as an integrity check. A confidence threshold was not applied to the analysis in order to view the genotype and C% assigned to each sample replicate. Profiles giving a C% of less than 95% to any of the existing profiles were considered distinctly different profiles. The means of the C% of the sample replicates correctly assigned to a representative genotype, together with the standard error (SE), were calculated using Microsoft Office Excel (2003).
DNA sequence analysis. Amplicons were gel purified using the QIAquick gel purification kit (Qiagen) according to the manufacturer's instructions, eluted in 30 µl of buffer E, and then subjected to automated sequencing (BigDye Terminator version 3.1, Applied Biosystems) in both directions using the same primers used for PCR. The sequences in the approximately 590-bp region (amplified by primers Hex L1-s and Hex L1-as) were analyzed using ClustalW and DNAdist in BioManager (Australian National Genomic Information Service, Sydney Bioinformatics) and BioEdit Sequence Alignment Editor (version 6.0.9.0).
Nucleotide sequence accession number. GenBank accession numbers were assigned to the nucleotide sequences of the FAdV reference strains as follows: FAdV-1 CELO, EU979367; FAdV-2 SR48, EU979368; FAdV-3 SR49, EU979369; FAdV-4 KR5, EU979370; FAdV-5 340, EU979371; FAdV-6 CR119, EU979372; FAdV-7 YR36, EU979373; FAdV-8a TR59, EU979374; FAdV-8b 764, EU979375; FAdV-9 A02, EU979376; FAdV-10 C2B, EU979377; and FAdV-11 UF71, EU979378.
|
|
|---|
The three primer sets all generated detectable amplicons (by agarose gel electrophoresis) from each of the 12 reference serotypes, with the exception of FAdV-5, which showed low-level or no amplification with the Hexon A/B primer set (Fig. 2).
![]() View larger version (12K): [in a new window] |
FIG. 2. Agarose gel electrophoresis of Hexon A/B, Hex L1, and HEX-S PCR products of the hexon gene from FAdV serotypes 1, 5, and 8a, plus a nontemplate control (NTC). MW, molecular weight marker (PCR marker, Sigma).
|
Restriction enzyme analysis of Hexon A/B PCR products failed to classify some FAdV reference serotypes. The region between nucleotides (nt) 144 and 1040 (amplified by primers Hexon A and B) of the FAdV reference serotypes used in this study were subjected to restriction enzyme analysis using the method described previously (29). Serotypes 1, 3, 6, 7, 8a, and 9 were differentiated by successive use of the restriction endonucleases BsiWI, StyI, and MluI. An additional enzyme (BglI) was used to differentiate serotypes 4 and 10. Table 3 shows the different restriction enzyme digestion patterns exhibited, denoted by profiles A to J.
|
View this table: [in a new window] |
TABLE 3. Results of restriction enzyme analysis, denoted by profiles A to J, and normalized HRM-curve analysis, indicating genotype and mean C% of sample replicates for 12 FAdV reference serotypes
|
Due to low-level or no amplification by the Hexon A/B primers, FAdV-5 was not included in the restriction enzyme analysis.
HRM-curve analysis of the Hex L1 PCR product differentiated all FAdV serotypes from each other. PCR products of the region of the hexon P1 gene amplified by the HEX-S primers (as defined in Table 2) from FAdV reference serotypes were subjected to HRM-curve analysis. In the conventional melt curve analysis, all serotypes generated one peak between 86°C and 88°C, with FAdV-8a generating two additional shoulder peaks between 88.5°C and 90.5°C (Fig. 3A). Examination of the FAdV-8a PCR product by agarose gel electrophoresis resulted in a single stained band (Fig. 2), indicating that the additional peaks were not due to other amplifications in the PCR. The melt curve profiles were consistently produced from PCR products generated on different days and from templates belonging to different DNA extractions. Some serotypes exhibited very similar melt curve profiles and peak temperatures, making it difficult to distinguish them from each other (Fig. 3A). These included FAdV-4 and -10, -8b and -9, and -6 and -11. Analysis of the normalized HRM curves (Fig. 3B) revealed that while some serotypes could easily be distinguished, for example FAdV-1 and -8a, the curve profiles of other serotypes were very similar and could therefore not be visually differentiated. When genotyping was applied to the normalized HRM curves (Table 3), a single replicate of some serotypes was auto-called as a different genotype at a C% between 98.9 and 99.4% (e.g., one replicate of FAdV-4 genotyped as FAdV-10 with a C% of 99.3). Therefore the normalized HRM-curve profiles of these serotypes could not be considered to represent distinctly different genotypes.
![]() View larger version (34K): [in a new window] |
FIG. 3. (A) Conventional melt curves of HEX-S PCR products of the P1 region of the hexon gene from 12 FAdV reference serotypes. (B) Normalized HRM curves of HEX-S PCR products of the P1 region of the hexon gene from 12 FAdV reference serotypes.
|
![]() View larger version (41K): [in a new window] |
FIG. 4. (A) Conventional melt curves of Hexon PCR products of the P1 and L1 regions of the hexon gene from 12 FAdV reference serotypes. (B) Normalized HRM curves of Hexon PCR products of the P1 and L1 regions of the hexon gene from 12 FAdV reference serotypes.
|
![]() View larger version (42K): [in a new window] |
FIG. 5. (A) Conventional melt curves of Hex L1 PCR products of the P1 and L1 regions of the hexon gene from 12 FAdV reference serotypes. (B) Normalized HRM curves of Hex L1 PCR products of the P1 and L1 regions of the hexon gene from 12 FAdV reference serotypes.
|
Differences in melting curve profiles were not necessarily related to nucleotide sequence identities but mainly to distribution of GC contents throughout the amplicons. In order to determine if differences in melting profiles correlated with sequence identities, the homologies of the nucleotide sequence of the Hex L1 PCR products from each FAdV serotype were compared (Table 4). Sequence identities ranged from 28% to 99%. More than 94% identity was found between FAdV-3 and -9 and FAdV-4 and -10, with an identity of 99% found between serotypes 2 and 11. The lowest identity (28%) was between FAdV-10 and -11.
|
View this table: [in a new window] |
TABLE 4. Sequence identity matrix and percentages of GC content corresponding to the region of nt 301 to 890 of the FAdV-1 CELO genes for hexon proteinsa
|
Similarly, the pair-wise alignment of FAdV-4 and -10 showed 23 nt differences, of which 7 were inverted. The overall GC compositions of the amplicons were almost identical, with 58.8% for FAdV-4 and 58.7% for FAdV-10 (Table 4). FAdV-4 and -10 shared similar GC contents of 58.6% and 58.0%, respectively, for the region spanning nt 1 to 193 of the 578-bp amplicon. However FAdV-10 had a higher GC content than FAdV-4 in the middle region (nt 194 to 386) with 62.2% and 60.6%, respectively. In the last third of the amplicon (nt 387 to 578), FAdV-10 had a lower GC content than FAdV-4, with 55.7% and 57.3%, respectively.
The serotypes sharing the highest identity, FAdV-2 and -11, were shown to differ by only six nucleotides, with no inverted differences. The overall GC contents of the 593-bp amplicons differed slightly, with 54.0% for FAdV-2 and 53.6% for FAdV-11 (Table 4). The GC contents of the region spanning nt 1 to 300 of the amplicons were similar with 55.0% and 55.3% for FAdV-2 and -11, respectively. The following 150 bp (nt 301 to 450), however, contained four of the six nucleotide differences, all of which contained a G or C nucleotide for FAdV-2 and an A or T nt for FAdV-11; therefore, FAdV-2 had a higher GC content (54.7%) than FAdV-11 (52.0%) in this region. The remaining 143-bp region of the amplicon (nt 451 to 593) contained similar low GC contents for FAdV-2 (51.1%) and FAdV-11 (51.8%).
Hex-L1 PCR HRM-curve analysis genotyped a vaccine and a field strain according to their serotype. In order to assess the capacity of the newly developed Hex-L1 PCR HRM-curve analysis for genotyping additional FAdV isolates, two viruses of known serotype (8b) were subjected to HRM-curve analyses. The viruses used were the live attenuated Australian FAdV vaccine and FAdV field isolate FI-W888-05. Analysis of the conventional melt curve revealed that these viruses generated identical patterns, two peaks at 86.2 and 88.5°C, which were most similar to reference strain FAdV-8b. Examination of the normalized HRM curves also revealed that the curves generated from the Australian vaccine and field strains had resemblance to FAdV-8b. When genotyping was applied in the presence of all reference strains, replicates of both the vaccine and field isolate were genotyped as FAdV-8b, with mean confidence percentages of 68.9% ± 1.89 SE and 57.6% ± 0.81 SE, respectively.
|
|
|---|
In previous studies the use of techniques such as microtiter virus neutralization tests (11, 16, 33, 38-40), PCR alone (49, 50, 59) or in combination with either sequencing (2, 14, 15, 30, 32, 48) or restriction fragment polymorphism (RFLP) (29, 44, 54), and a combination of these techniques (7, 9, 10, 20, 37, 41, 42, 57) have been used to identify FAdV serotypes. Although each of these techniques is readily available to most laboratories, they have limitations, including the lengthy processes involved, need for extensive interpretation, and potential for indeterminate results, for example, the cross-neutralization of serum samples (7, 10, 16, 27).
The results of this study indicate that restriction enzyme analysis, or RFLP, failed to classify FAdV serotypes 2 and 11 (Table 3) using enzymes previously shown to distinguish 12 FAdV serotypes (29, 30). Therefore, an additional restriction enzyme would be required to differentiate these two serotypes by this method. It should be noted that the representative strain of each serotype used in this study did not conform with those used previously (29). The ambiguity of the nomenclature of FAdV serotypes under two different systems, European Union and U.S. (Table 1), has made the correlation of data relating to serotypes sourced from different geographical locations difficult. The International Committee on Taxonomy of Viruses recently implemented a single international nomenclature system for FAdVs (19). The results of the current study suggest that RFLP may not be able to distinguish serotypes from the same species and/or strains from the same serotype, which may be of importance in epidemiological studies of disease outbreaks. In addition, RFLP requires a high quantity of PCR product, which was shown not to be feasible in this study for FAdV-5 amplification by primers Hexon A and B (Fig. 2).
In this study HRM-curve analysis was applied to amplicons of three different regions of the FAdV hexon gene, using conventional PCR reagents, oligonucleotide primers, and Syto 9. Amplicons of the 191-bp HEX-S region from FAdV serotypes could not easily be discriminated from each other by analysis of their HRM curves (Table 3). This may have resulted from small amplicons containing limited sequence variations.
The amplicons of the approximately 897-bp region generated by primers Hexon A and B produced conventional melt curves with two or more major and minor peaks, with or without shoulders (Fig. 4A), with most serotypes distinguishable from others using the normalized HRM curves (Table 3). However, as these primers yielded low-level or no amplification of FAdV-5, further analysis was not feasible for this serotype. Therefore, HRM-curve analysis of the Hexon A/B PCR product was deemed not suitable for genotyping of FAdV reference serotypes.
In contrast, HRM-curve analysis of Hex L1 PCR products proved a highly sensitive and specific means of genotyping all 12 FAdV serotypes. The complexity of the conventional melt curves of the approximately 590-bp amplicons increased the power of the normalized HRM-curve analysis to differentiate serotypes. All serotypes generated one or more major peaks and were visually distinct from each other in conventional melt curve profiles (Fig. 5A), and replicates of each serotype were correctly genotyped using the normalized HRM curves with C% of more than 99% (Table 3). Even where two serotypes exhibited a single major peak with an identical melting temperature (FAdV-3 and -9 in Fig. 5A), in the normalized HRM graph these serotypes could be clearly differentiated (Fig. 5B) and replicates were correctly genotyped with more than 99% confidence (Table 3). The capacity of the newly developed Hex L1 PCR HRM-curve analysis to genotype nonreference strains was also demonstrated by correctly genotyping two further strains (an Australian FAdV vaccine and a field strain) according to their known serotypes in the presence of all 12 reference strains. Even though the design of the Hex L1 primers was based on the GenBank sequences of 27 different FAdV strains (48), further proof of the universality of these primers would become evident if they were utilized in further research involving a number of other FAdV strains.
In this study, FAdV-2 and -11 sequences shared 99.0% identity in the region amplified by the Hex L1 primers (Table 4) but nevertheless produced different conventional melt curve profiles (Fig. 5A), and replicates in the normalized HRM graph were correctly genotyped with close to 100% confidence (Table 3). This demonstrates the sensitivity of the technique developed here and its capacity for detecting very small sequence variations.
The melting profile of a PCR product is dependent upon length, sequence, GC content, and heterozygosity (17, 51, 53). The GC composition of a DNA fragment is not necessarily a good predictor of the melting temperature; however, the actual DNA melting analysis can detect finer-scale differences in sequences, such as localized domains of different GC composition (31). It has been suggested that melting domains are usually 50- to 300-bp in length; therefore, amplicons up to 800-bp may have multiple melting domains (58). Thus, it can be assumed that multiple and/or shoulder peaks in the conventional melt curve profile of amplicons of the L1 region of the hexon gene of some FAdV serotypes, as seen in Fig. 5A, may be due to the presence of multiple melting domains differing in GC content. For example, the melting temperatures of the three peaks for FAdV-8a, 86.7°C, 88.5°C, and 90.2°C (Fig. 5A), may be attributable to the different GC compositions of each third of the amplified sequence (51.5%, 57.3%, and 62.8%, respectively). Interestingly, the peak with the lowest melting temperature closely matches that of FAdV-3 and -9, which have overall GC contents (51.8% and 51.4%, respectively) close to that of the FAdV-8a region (51.5%) from nt 399 to 596. Even where sequence differences were small and localized, as seen between FAdV-2 and -11, the conventional melt curve profiles (Fig. 5A) were distinctly different as a result of contrasting GC compositions within a confined region of each amplicon. This supports the notion that the composition and distribution of GC nucleotides in an amplicon directly influences the shape and complexity of the conventional melt curve, and in turn, the power of the normalized HRM curve as a tool for genotyping.
The melt curve profiles generated in this study were highly reproducible and in agreement with those predicted by the web-based DNA melting simulation program POLAND (Heinrich-Heine University in Dusseldorf, Germany, Institute of Biophysics [http://www.biophys.uni-duesseldorf.de/local/POLAND/poland.html]) (55) (results not shown), a program that has been used to aid in the design of primers to amplify a region of DNA that will result in certain melting-curve characteristics (8, 23). This also confirmed that the complex melting profile of PCR amplicons was solely the result of DNA dissociation (47). The variability between different strains within each FAdV serotype could potentially be assessed using sequences from the GenBank database in the POLAND program. The Hex L1 nucleotide sequences from some FAdVs with known serotype identities were obtained from GenBank (AF508957, FAdV-8 strain 58; DQ323986, FAdV-9 isolate Stanford; U26221, FAdV serotype 10; and AF339925, FAdV-12 strain 380) and analyzed using POLAND as part of this study. The resultant melt curve shapes were highly similar, with regard to the number, height, and temperature of the peaks, to the respective reference strains used in this study (results not shown). Combining this technique with HRM-curve analysis could therefore provide serotype identification and potentially strain identification in laboratories that can perform HRM-curve analysis but do not have access to the reference strains of FAdV serotypes. However, as genotyping is performed using the normalized HRM data and further evaluated by the resultant confidence percentages in the RotorGene software, the use of the POLAND program to assess and compare sequences from GenBank is limited to visual comparisons of conventional melt curve shapes and is therefore not definitive.
This study presents for the first time a robust genotyping technique based on the closed-tube method of real-time PCR and HRM-curve analysis for FAdVs. The method is an accurate, rapid, and cost-effective alternative to existing serotype identification methods, and it can be used as a model for developing genotyping techniques for other pathogens. In addition, during this study, the basis for differences in melt curve profiles resulting from reference FAdVs were elucidated.
We thank Intervet International (Boxmeer, Netherlands) for providing the FAdV reference serotypes and Tom Grimes for his assistance and input.
Published ahead of print on 26 November 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»