Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2004, p. 1686-1693, Vol. 42, No. 4
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.4.1686-1693.2004
and Kevin L. Russell1*
Department of Defense Center for Deployment Health Research, Naval Health Research Center, San Diego, California 92186-5122,1 Brooks City-Base, San Antonio, Texas 78235,2 Division of Biological Sciences, University of California, San Diego, La Jolla, California 920933
Received 25 September 2003/ Returned for modification 10 October 2003/ Accepted 25 November 2003
|
|
|---|
|
|
|---|
The classical method for subtyping adenoviral isolates is whole-genome digestion via a stepwise, systematic restriction enzyme analysis (REA) process. Restriction enzyme classification methods use a numbering and lettering system appended to the serotype number to distinguish unique strains (16). The letters "a" through "k" represent restriction enzyme whole-genome electrophoretic banding patterns by the referent enzyme BamHI. An Arabic numeral is added when additional enzymes are used to further distinguish whole genomes.
Intensively studied in the 1950s and 1960s, adenoviruses were found to infect up to ca. 80% of military recruits and lead to hospitalization in up to 20% (9). Adenovirus type 4 (Ad4) and Ad7 were the primary serotypes responsible for this morbidity and together constituted 60% of all hospitalized cases of acute respiratory disease (ARD) among military recruits (9). Live, enteric-coated oral vaccines, which induce immune responses through selective infection of the gastrointestinal tract, were developed first for Ad4 and later for Ad7a (3, 28). These vaccines were shown to be safe and highly effective in the immunization of military trainees (28), and routine administration to recruits began in 1971 (10). The Ad4 and Ad7 vaccines together lowered AV morbidity by 95 to 99% and total ARD morbidity by 50 to 60% during the period of vaccine use (10). However, the production of the vaccines was discontinued in 1996, and the remaining lots were rationed until supplies were exhausted in early 1999 (10). Recently, the U.S. Department of Defense awarded a contract to Barr Laboratories, Inc., for resumed production of the vaccines.
Population-based febrile respiratory illness surveillance was initiated by the Naval Health Research Center in 1996 to document the epidemiology of adenoviruses during and after the period of vaccine loss (10). The surveillance program was originally established at four recruit training camps in the United States to define the burden of adenoviruses (10) but was later expanded to eight sites and included testing for other viral agents (24). This surveillance documented large increases in adenovirus morbidity and several fatal cases after the vaccine was exhausted (2, 10, 25), suggesting that the initial vaccine was efficacious for the majority of circulating pathogenic strains. However, several cases of Ad4 and Ad7 infection were discovered among previously vaccinated individuals, raising the possibility that newly emergent strains of adenovirus had appeared.
Recent research on the evolution of circulating adenoviruses has engendered concern about the efficacy of the old vaccine against current strains. In order to determine the suitability of the original vaccine strains for a new vaccine, a study of strain variation among the circulating Ad4 and Ad7 serotypes was conducted in 1999 (5). The antigenic regions of the hexon gene from prototype, vaccine, community-acquired, and military wild-type strains collected from 1953 to 1997 were sequenced and compared. Whereas the hexon antigens of Ad7 were generally conserved over time, an Ad4 variant (strain Z-G 95-873) with nine amino acid changes in the hexon antigens was found to have been circulating since 1995 (5). These changes were noted to confer decreased neutralization. The ability of this strain to cause infection among vaccinated individuals, however, was not investigated. Another recent study analyzed Ad7 isolates from the United States from the years 1966 to 2000 by whole-genome REA (8). The study noted the appearance of two genome types previously undocumented in North America, Ad7d2 and Ad7h, which indicated a shift in the predominant Ad7 genotype circulating in the United States. The hexon protein of Ad7d2 contains only one unique amino acid substitution, but it has possible antigenic implications.
The primary objective of the present study was to determine whether the adenovirus infections among previously vaccinated individuals who became ill between 1996 and 2000 resulted from prototype strains, recently discovered adenovirus variants, or a completely new variant. A 1,500-bp region of the adenovirus hexon capsid protein gene contains seven discrete hypervariable regions (HVR1 to HVR7) that account for >99% of serotype-specific variation and code for the type-specific epitopes on the protein (4, 5). This region of the hexon gene from adenoviruses isolated from vaccinated individuals and unvaccinated individuals was sequenced and compared to the vaccine strains, the Ad4 variant strain, and the other existing Ad4 and Ad7 hexon sequences in GenBank. The results of the present study will help guide future vaccine development initiatives.
|
|
|---|
38°C and either a cough or sore throat, are recorded. Recruits meeting the case definition are asked to permit the collection of a throat swab specimen and to answer a brief questionnaire. Specimens are stored in viral transport medium (Remel, Lenexa, Kans.) and sent in batches with associated data to the Naval Health Research Center Respiratory Disease Laboratory, San Diego, Calif., for analysis. Classic viral culture and isolation are performed using R-Mix shell vials (Diagnostic Hybrids, Athens, Ohio) and A549 (American Type Culture Collection, Manassas, Va.) and RMK (Diagnostic Hybrids) cells. Cultures exhibiting viral growth are processed for identification of the infecting viral pathogen by using immunofluorescence. A random selection of isolates positive for adenovirus are further analyzed by microneutralization with type-specific antisera to determine serotype, as previously described (6, 17). Viral isolate selection. Adenovirus-positive samples, negative for other respiratory pathogens, were considered for inclusion in the study. In this manner, adenovirus was likely the causative agent of illness in the symptomatic recruit at the time of collection. The vaccination status of the host (vaccinated or unvaccinated) and serotype (Ad4 or Ad7), as previously determined by microneutralization, were then considered. Vaccinated individuals were defined as those with febrile respiratory illness beginning >10 days after documented receipt of the Ad4 and Ad7 vaccines. The Ad4 sample set consisted of 10 isolates from vaccinated individuals and 30 from unvaccinated individuals and represented six of the eight recruit camps participating in surveillance. The Ad7 sample set consisted of 11 isolates from vaccinated individuals and 1 from an unvaccinated individual. A low incidence of Ad7 in recent years, with the exception of an outbreak at Great Lakes in the Fall of 1997 (25), limited the availability of clinical isolates from unvaccinated individuals for this set. However, additional sequences from GenBank supplemented the unvaccinated Ad7 sequences used in this analysis.
PCR amplification and sequencing. All selected samples were regrown in A549 cells from original patient specimens. Viral DNA was extracted with the QIAamp DNA Blood Mini Kit (Qiagen, Valencia, Calif.). The early conserved region and serotype-defining hypervariable regions of the adenovirus hexon, corresponding approximately to the first 1,500 bp of the gene, was amplified by PCR with the primers UP and 5L, as previously described (5). The reaction mixture contained 1x PCR Buffer II (Applied Biosystems, Foster City, Calif.), 3 mM MgCl2, 200 µM concentrations of each deoxynucleoside triphosphate, 0.5 µM concentrations of each primer, 2.5 U of Taq polymerase (Applied Biosystems), and 5 µl of DNA template in a total volume of 100 µl. The cycling consisted of an initial denaturation at 94°C for 2 min, followed by 30 cycles of 94°C for 1 min, 45°C for 1 min, and 72°C for 3 min, with a final extension of 72°C for 7 min. The amplicons were purified with QIAquick PCR purification kit (Qiagen).
Purified adenovirus hexon amplicons were cycle sequenced with the Perkin-Elmer 9700 thermocycler (PE Biosystems, Foster City, Calif.), using the Big Dye terminator cycle sequencing kit V3.0 (Applied Biosystems). The sequencing consisted of 25 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min. All products were purified by using CentriSep columns (Princeton Separations, Adelphia, N.J.). Direct nucleotide sequencing was performed with an ABI 310 genetic analyzer (Applied Biosystems) according to the manufacturer's protocols. The primers used for sequencing included two PCR primers (UP and 5L), four internal primers (1Lm, 1Rm, 3L, and 3R) as previously described (5), and four novel internal primers (R630, R822, R440, and R674) (Table 1). Contig assembly was performed using the DNAstar software package (DNA Star, Inc., Madison, Wis.). BioEdit (v5.0.9) (11) was used for multiple sequence alignments, DNA translations, and additional sequence analysis.
|
View this table: [in a new window] |
TABLE 1. Adenovirus hexon and fiber gene primers
|
Analysis. Phylogenetic trees were calculated with PAUP version 4.0 b9 (26) for all nucleotide and predicted amino acid sequences by neighbor-joining and parsimony methods. Sequences were analyzed with reference to the trees to reveal character states relevant to phylogenetic branching. Synonymous substitution rates, defined as the number of synonymous substitutions per 1,000 synonymous sites, and nonsynonymous substitution rates, defined as the number of nonsynonymous substitutions per nonsynonymous site, were computed for all pairwise sequence comparisons within the hypervariable and conserved regions of each serotype by using MEGA, version 2.1 (13). These values were used to detect differential evolutionary selective pressure between hypervariable and conserved regions of the hexon.
GenBank accession numbers. Novel nucleotide sequences were submitted to the GenBank database. The accession numbers for submitted Ad4 sequences included AY337237, AY337238, AY337239, AY337240, AY337241, AY337242, AY337243, AY337244, AY337245, AY337246, AY337247, AY337248, AY337249, AY337250, AY337251, and AY337252 (for Ad4 isolates) and AY337253, AY337254, AY337255, AY337256, AY337257, and AY337258 (for Ad7 isolates).
|
|
|---|
![]() View larger version (54K): [in a new window] |
FIG. 1. Proportional distribution of adenovirus serotypes from febrile respiratory surveillance in military recruit camps from 1996 to 2002.
|
![]() View larger version (21K): [in a new window] |
FIG. 2. Ad4 and Ad7 hexon maximum-parsimony protein tree with canine Ad1 root. Isolates sequenced in the present study are labeled according to the following format: last two digits of accession number of identical sequence.location of isolation.year of isolation.vaccination status. FB, Fort Benning; FJ, Fort Jackson; FLW, Fort Leonard Wood; LACK, Lackland; GL, Great Lakes; MCRD, Marine Corps Recruit Depot; V, vaccinated; UV, unvaccinated. Full accession numbers appear once in parentheses. Accession numbers in boldface correspond to those in Fig. 3. Multiple accession numbers among identical protein sequences are due to DNA sequences with silent mutations. Bootstrap values greater than 50 (of 100 total bootstraps) are shown. Rescaled consistency index = 0.97.
|
![]() View larger version (41K): [in a new window] |
FIG. 3. Amino acid polymorphisms within aligned Ad7 (a) and Ad4 (b) sequence sets. Conserved positions are omitted, and residues identical to first sequences are represented by a period. Numbers refer to the amino acid position according to alignment with the Gomen prototype (AF065065) (a) and the Ad4 vaccine strain (AF065063) (b). Isolates correspond to those in boldface in Fig. 2.
|
The variant group was distinguished from the Ad4 vaccine-strain group by 10 amino acid substitutions (Fig. 3b). Of the 33 military samples in the variant group, 22 had identical amino acid sequences (Fig. 3b, sequence AY337244). Two exceptions (49.FLW.97.V and 51.LACK.00.NV) contained a nonconservative K216N substitution, as defined by differences in size and charge (4), which they shared only with the Z-G 95-873 variant and the vaccine strain (Fig. 3b).
Synonymous versus nonsynonymous substitutions. Pairwise comparison of conserved regions within each serotype yielded significantly greater mean synonymous substitution rates than mean nonsynonymous substitution rates (Table 2). In contrast, pairwise comparison of hypervariable regions within each serotype yielded greater mean nonsynonymous substitution rates than mean synonymous substitution rates. The differences between the means in the hypervariable regions were not significant; however, 103 of 231 independent pairwise comparisons among the Ad7 sequences had significantly greater nonsynonymous substitution rates, whereas only 10 had significantly greater synonymous substitution rates. Likewise, with the Ad4 sequences 319 of 780 comparisons had significantly greater nonsynonymous substitution rates, and 26 had significantly greater synonymous substitution rates. The large number of nonsignificant pairwise comparisons was primarily the result of a high level of sequence redundancy in our data set.
|
View this table: [in a new window] |
TABLE 2. Mean synonymous and nonsynonymous substitution rates in conserved and hypervariable regions of Ad7 and Ad4
|
|
|
|---|
The major phylogenetic bifurcation within Ad7 hexon sequences was previously shown (5, 8, 14) and provided genetic evidence for a distinction originally seen only at a phenotypic level by cross-neutralization and REA (5). An eightfold reduction in neutralization of the Gomen Ad7 prototype compared to the S-1050 Ad7a strain with Ad7a antiserum was also previously shown (5). This neutralizing difference may have contributed to a breakthrough of the sole strain from a vaccinated patient in the prototype group (58.GL.96.V; Fig. 2), but the low number of prototype-like strains isolated in the present study indicated that they are not currently a major epidemiological threat in recruit camps.
A recent REA study on the 1997 Great Lakes outbreak showed >70% of the samples were of the Ad7d2 genome type (8). It was also found in the same study, as elsewhere (14), that an L443Q hexon substitution in HVR7 is consistently linked with Ad7d and Ad7d2, as distinct from Ad7a, Ad7b, Ad7c, Ad7g, and Ad7h. Of 10 samples in our study from the same outbreak but from vaccinated individuals, 7 also contained the L443Q substitution, indicating that these samples are most likely of the Ad7d2 genome type. The substitution is nonconservative and dramatically affects the hydropathicand probably structuralcharacteristics of the protein by transforming the nine-residue surrounding region in HVR7 from a secondary-structure conformation that exposes primarily hydrophobic amino acid residues to a conformation that exposes primarily hydrophilic residues (14). It has therefore been suggested that this substitution imparts antigenic implications (8), which may explain the high number of Ad7s of this type isolated from vaccinated individuals in our study. In contrast, however, evidence showing a rapid decrease in AV-associated illness after introduction of the vaccine in this outbreak (25) strongly suggests that the strain is at least partially susceptible to vaccine-induced immunity. Cross-neutralization data between Ad7a and a strain containing the L443Q substitution would provide more evidence of the antigenic implications.
Several other substitutions occurring among the Ad7s may have contributed to infection in vaccinated individuals. The nonconservative S432Y substitution in samples 56.GL.97.V and 57.GL.97.V (Fig. 3a) resided in HVR7 and may have directly altered the epitope. In addition, isolated nonconservative substitutions in nearby conserved regions may have indirectly affected the antigens. The hexon epitopes are known to be conformational, since only the protein in its native trimeric form and not linear hexon peptides or heat-denatured monomeric proteins neutralize in vitro (4). Therefore, a change in a structural region, such as the P480S substitution (Fig. 3a, 55.GL.97.V) in the conserved region at the base of the loop containing HVR7 (1, 4), for example, may drastically affect protein folding in the antigenic regions.
Given the predominance of Ad4 among cases of respiratory illness in recent years, as well as the predominance of the variant genotype within Ad4 (Fig. 1), it appears that this single strain was responsible for nearly all adenovirus-associated respiratory illness in military recruit camps since 1999. The continued increase in the dominance of this strain after vaccine cessation suggests that its proliferation was not primarily a result of vaccine use. It is difficult to say how well the prior vaccine strain will clinically protect against this variant, since the majority of the variant isolates were isolated after vaccine use was drastically reduced. However, it was recently shown in an in vitro cross-neutralization analysis that the Z-G 95-873 variant has a fourfold reduction in neutralization in comparison with the Ad4 RI-67 prototype strain (5), a strain similar in genotype to the vaccine strain. The majority (32 of 35) of genotypes in the variant group differed from the Z-G 95-873 variant, as well as the vaccine strain, by the N216K amino acid substitution (Fig. 3b, sequence AY337247). If this substitution changes the epitope at all, it would render it further unrecognizable by vaccine-induced antibodies; however, it is unknown whether the substitution would have such an effect, given its nonconservative character in a conserved region. Nevertheless, due to the current prevalence of the Ad4 variant with the N216K mutation, any new vaccine should be optimized for effectiveness against this strain.
As with Ad7, several isolated substitutions in the Ad4 set have the potential to alter the epitope. Three nonconservative amino acid changes (Fig. 3b, P240T in sample 39.GL.97.V, N246K, and S420F) resided in antigenic regions. Nonidentical nucleotide mutations code for the N246K substitution, which suggests convergent evolution at this site, perhaps driven by selection. However, given that only one sample from a vaccinated individual was found with each mutation, our data do not conclusively show their clinical effects.
Nucleotide sites are classified as either nonsynonymous or synonymous depending on whether or not, respectively, a mutation at the site will lead to an amino acid substitution. The first and second nucleotide sites in a codon are largely nonsynonymous, and the third site, due to redundancy in the genetic code, is largely synonymous. The degree to which a site is classified as either synonymous or nonsynonymous is determined by accounting for the coding imparted by all possible mutations. Within a gene region, the number of nonsynonymous substitutions per nonsynonymous site relative to synonymous substitutions per synonymous site indicates the type of evolutionary selective pressure on the region. Equal rates of nonsynonymous and synonymous mutations indicate neutral drift, an excess of synonymous substitution rates indicates conservative evolutionary pressure (purifying selection), and an excess of nonsynonymous substitution rates indicates selection for diversity (diversifying selection). Evidence of diversifying selection is often found in protein regions that benefit from variation such as surface antigens of pathogens (7) and in immunoglobulin variable regions (27). In our study, the nonsynonymous substitution rates in the antigenic (hypervariable) regions of the sequences suggest that these regions are affected by diversifying selection. In contrast, the nonantigenic (conserved) regions show strong purifying selection, as is expected from regions with a structural function. Given the propensity for variation in the hypervariable regions, one would expect that the selective pressure of vaccination efforts might eventually allow the virus to evolve resistance to vaccine-induced immune responses. Antigenic cross-neutralization analysis between new vaccine strains and current circulating strains should be closely monitored.
A recent study of the adenovirus hexon (23) proposed a reassignment of the hypervariable regions based on stretches of diversity in an alignment of 40 partial hexon amino acid sequences from both human and nonhuman adenoviruses. The hypervariable regions as presented in the alignment in the study generally overlap with the assignments as previously defined for Ad2 (4) but not with those previously defined for Ad4 and Ad7 (5). An analysis of mutation rates in the reassigned hypervariable regions of the sequences used in our study shows higher synonymous than nonsynonymous substitution rates in the hypervariable regions (data not shown), which is the opposite of what was found with the previous assignments. Additional support for the previous region assignments is provided by a visual inspection of the alignment used to define the new assignments. This reveals that definite motifs in a subset of sequences, including that of Ad2, are shifted out of phase with those of Ad4, Ad7, and several others. Such a shift would affect identity scores assigned to each region and hence the placement of hypervariable regions in the alignment. Our data therefore support the hypervariable positions as previously defined for Ad4 and Ad7.
Initial use of the Ad4 vaccine in the 1960s was associated with Ad7 and Ad21 emerging as the chief etiological agents of ARD in military recruit camps. At the time, this constituted the first outbreak of Ad21 in the Western hemisphere (21). The isolation of three of the four Ad5s from vaccinated individuals indicates that this serotype may have similarly filled an evolutionary niche left open during the period of Ad4 and Ad7 vaccine use, but to a much smaller extent. Recent evidence shows that Ad5 manifests a lower binding affinity and replicates more slowly in A549 cells relative to other adenoviruses in subgenus A, B (includes Ad7), D, and E (includes Ad4) (18). The A549-mediated growth typically used in the serological diagnosis of AV may fail to produce an Ad5 titer sufficient for detection by microneutralization. The epidemiological impact of Ad5 on military recruit populations during the period of vaccine use may therefore have been underestimated. Molecular surveillance should be used to monitor for the possible emergence of Ad21, Ad5, or other serotypes upon reinstitution of the Ad4 and Ad7 vaccines.
Although the present study cannot definitively explain the reason for past vaccine failures, it shows the distribution of strains with possible partial resistance to the previous vaccines. Given the predominance of the Ad4 variant in recent years, an effective vaccine against it may significantly reduce adenovirus-associated respiratory infection in the U.S. military. Although Ad7 is not circulating at present, protection against the Ad7d2 strain would prevent its opportunistic proliferation upon controlling Ad4. Other serotypes have not shown Ad7's capacity to fill the Ad4 niche but, with diversifying selection acting on the hexon epitopes, continued surveillance for newly emergent strains should be conducted.
The views expressed here are those of the authors and do not reflect the official policy or position of the Department of the Navy, the Department of Defense, or the U.S. Government. The work is approved for public release (distribution unlimited). This research has been conducted in compliance with all applicable federal regulations governing the protection of human subjects in research under protocol 32267.
We thank Anthony W. Hawksworth of the Department of Defense Center for Deployment Health Research for organization and analysis of febrile illness surveillance data. We also thank Leta K. Crawford-Miksza and David P. Schnurr of the Viral and Rickettsial Disease Laboratory, Division of Communicable Disease Control, California State Department of Health Services, for adenovirus antisera and expert consultation.
Present address: College of Public Health, University of Iowa, Iowa City, IA 52242. ![]()
|
|
|---|
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»