ORIGINAL RESEARCH article

Front. Microbiol., 05 June 2018

Sec. Virology

Volume 9 - 2018 | https://doi.org/10.3389/fmicb.2018.01180

Comparative Genomic Analysis of Re-emergent Human Adenovirus Type 55 Pathogens Associated With Adult Severe Community-Acquired Pneumonia Reveals Conserved Genomes and Capsid Proteins

  • 1. Guangzhou Key Laboratory of Drug Research for Emerging Virus Prevention and Treatment, Guangdong Provincial Key Laboratory of Tropical Disease Research, School of Public Health, Southern Medical University, Guangzhou, China

  • 2. Treatment and Research Center for Infectious Diseases, 302 Military Hospital of China, Beijing, China

  • 3. Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Beijing, China

  • 4. Dermatology Hospital, Southern Medical University, Guangzhou, China

Abstract

Human adenovirus type 55 (HAdV-B55) is a recently identified acute respiratory disease (ARD) pathogen in HAdV species B with a recombinant genome between renal HAdV-B11 and respiratory HAdV-B14. Since HAdV-B55 first appeared in China school in 2006, no more ARD cases associated with it had been reported until 2011, when there was an outbreak of adult severe community-acquired pneumonia (CAP) in Beijing, China. Reported here is the bioinformatics analysis of the re-emergent HAdV-B55 responsible for this outbreak. Recombination and protein sequence analysis re-confirmed that this isolate (BJ01) was a recombinant virus with the capsid hexon gene from HAdV-B11. The selection pressures for the three capsid proteins, i.e., hexon, penton base, and fiber genes, were all negative, along with very low non-synonymous (dN) and synonymous (dS) substitutions/site (<0.0007). Phylogenetic analyses of the whole genome and the three major capsid genes of HAdV-B55 revealed the close phylogenetic relationship among all HAdV-B55 strains. Comparative genomic analysis of this re-emergent HAdV-B55 strain (BJ01; 2011) with the first HAdV-B55 strain (QS-DLL; 2006) showed the high genome identity (99.87%), including 10 single-nucleotide non-synonymous substitutions, 11 synonymous substitutions, 3 insertions, and one deletion in non-coding regions. The major non-synonymous substitutions (6 of 10) occurred in the protein pVI in its L3 region, which protein has different functions at various stages of an adenovirus infection, and may be associated with the population distribution of HAdV-B55 infection. No non-synonymous substitutions were found in the three major capsid proteins, which proteins are responsible for type-specific neutralizing antibodies. Comparative genomic analysis of the re-emergent HAdV-B55 strains associated with adult severe CAP revealed conserved genome and capsid proteins, providing the foundation for the development of effective vaccines against this pathogen. This study also facilitates the further investigation of HAdV-B55 epidemiology, molecular evolution, patterns of pathogen emergence and re-emergence, and the predication of genome recombination between adenoviruses.

Introduction

Human adenoviruses (HAdV) are associated in humans with a broad spectrum of clinical diseases, such as acute respiratory disease (ARD), community-acquired pneumonia (CAP), as well as conjunctivitis, gastrointestinal, opportunistic infections in immune-deficient individuals, and possibly, obesity (Dingle and Langmuir, 1968; Chu and Pavan-Langston, 1979; Wood, 1988; Dhurandhar et al., 1992; Kojaoghlanian et al., 2003; Hakim and Tleyjeh, 2008). Since the first HAdV isolate was found in 1953 (Hilleman et al., 1953; Rowe et al., 1953), more than 80 HAdV genotypes have been identified and classified into seven species (A to G) using the whole-genome sequence as the typing standard1 (Jones et al., 2007; Ishiko et al., 2008; Walsh et al., 2009, 2010, 2011; Liu et al., 2011, 2012; Matsushima et al., 2011, 2012, 2013; Robinson et al., 2011, 2013b; Seto et al., 2011; Dehghan et al., 2012; Singh et al., 2012, 2013; Zhou et al., 2012; Alissa Alkhalaf et al., 2014; Hage et al., 2015; Hashimoto et al., 2018). HAdV species B are further divided into subspecies B1 (HAdV-B3, -B7, -B16, -B21, and -B50) and B2 (HAdV-B11, -B14, -B34, -B35, and -B55) according to their restriction enzyme digestion patterns (Wadell et al., 1980). Subspecies B1 cause respiratory disease except for HAdV-B50. The members in subspecies B2 cause renal or urinary tract disease, except for HAdV-B14 and -B55, which are associated with ARD (Echavarria, 2008).

Human adenoviruses are usually typed according to serum neutralization and hemagglutination-inhibition tests, both of which are associated with the three major capsid proteins, i.e., hexon, fiber, and penton base. Based on partial characterization of its hexon and fiber epitopes, as well as the restriction enzyme patterns (Zhang Q. et al., 2016), HAdV-B55 was first identified as HAdV-B11a (Hierholzer et al., 1974). This virus may have occurred earlier from an ARD outbreak in a military training base (Spain, 1969), and subsequently re-emerged in several outbreaks of ARD in Turkey (2004), Singapore (2005) discontinuously (Hierholzer et al., 1974; Chmielewicz et al., 2005; Zhu et al., 2009; Kajon et al., 2010; Li et al., 2014; Lu et al., 2014). Complete genome analysis suggested that this new type, HAdV-B55, is a Trojan horse, which contains a recombinant genome from both a renal pathogen, HAdV-B11, which provides the antigenic epitope, and a respiratory pathogen, HAdV-B14, which confers the cell tropism, biological and pathogenicity properties (Walsh et al., 2010). The strain QS-DLL is the first completely characterized HAdV-B55 based on the whole genomic analysis (Walsh et al., 2010). It caused an ARD outbreak in China, from March to April 2006, with 254 patients affected, of which 247 were students in a senior high school and one died during this outbreak (Zhu et al., 2009).

Since this first isolate QS-DLL of HAdV-B55 was identified in 2006 (Yang et al., 2009; Walsh et al., 2010), no more ARD cases associated with this virus had been reported in China until 2011, when 16 HAdV-B55-like strains were identified based on partial gene sequencing in a CAP outbreak in Beijing, China in February and March (Cao et al., 2014). Of these, strain BJ01 was isolated from a throat swab of a 29-year-old adult who presented with severe CAP (Gu et al., 2012). The whole genome (34,773 bp) of strain BJ01 (Human adenovirus B human/CHN/BJ01/2011/55[P14H11F14]) was sequenced using the Sanger method, following PCR amplification of targeted overlapping regions. (GenBank accession number JX491639). (Zhang et al., 2012a). In this study, comparative genomic analysis of this re-emergent HAdV-B55 strain (BJ01) with the first HAdV-B55 strain (QS-DLL), as well as the bioinformatic analysis of HAdV-B55 strains circulating recently were performed. The findings will contribute to the studies of the evolution and epidemiology of HAdV-B55, as well as provide both the foundation for the development of effective vaccines and public health strategies against this pathogen.

Materials and Methods

Cells, Virus Stock, and Genomic DNA Extraction

This specimen was originally isolated from a throat swab of an adult diagnosed with severe CAP at Beijing Chao-Yang Hospital of China (March 2011) (Gu et al., 2012; Zhang et al., 2012a). Written informed consent was obtained from the patient. The protocol was approved by the Ethics Committee of Beijing Chao-Yang Hospital, Capital Medical University in accordance with the Declaration of Helsinki. The throat swab was inoculated and cultured in A549 cells (ATCC), and grown in Dulbecco’s minimum essential medium supplemented with 100 IU/ml penicillin, 100 μg/ml streptomycin, and 2% (v/v) fetal calf serum, at an atmosphere of 5% (v/v) carbon dioxide. Viral genomic DNA was extracted from infected A549 cells, as described earlier (Le et al., 2001). PCR amplification and molecular typing was performed according to the previous study (Han et al., 2013).

Genome Sequencing and Annotation

The genome of strain BJ01 was sequenced using a Sanger primer-walking method, following PCR amplification of targeted overlapping regions (Zhang et al., 2012a,b). Both 5- and 3-ends (including both inverted terminal repeats) were sequenced directly by 5 primer Ad14-LTRS1 (5-ACAGAAGACTTTCACACGGT-3) and 3 primer Ad14-LTRS2 (5-GGTCCCTCTAAATACACATACA-3), respectively, using extracted genomic DNA as sequencing template. The sequence data provided average three- to fivefold coverage of the genome. Ambiguous sequences and gaps were re-amplified and re-sequenced. All the sequencing ladders were assembled using the SEQMAN (DNAStar; Madison, WI). The genome was further annotated according to the previous report (Walsh et al., 2010).

Genome Recombination Analysis

Genome sequence recombination analysis was done using Simplot 3.5.12, including the options of Simplot and Bootscan analysis (Lole et al., 1999). The software MAFFT was used to align the HAdV-B sequences3. The parameters of gap opening penalty (1.53), Offset value (0), and the Scoring matrix for nucleotide sequences (200 PAM/k = 2) were set by default. The LAGAN (Limited Area Global Alignment of Nucleotides) program of mVISTA was used for the pair-wise comparisons of genomes4 (Brudno et al., 2003).

Phylogenetic Analysis

Phylogenetic analysis was performed by Molecular Evolutionary Genetics Analysis (MEGA) 5.0.55 (Tamura et al., 2007). Phylogenetic trees with 1,000 bootstrap replications were constructed using a maximum-composite-likelihood method. All the other parameters were set by default.

Substitution Rate Analysis of HAdV-B55 Hexon, Penton Base, and Fiber Genes

The numbers of non-synonymous (dN) and synonymous (dS) substitutions per site between sequences and dN/dS ratios were calculated. The analysis was performed with MEGA 5.05 (Tamura et al., 2011) using the Nei–Gojobori model (Nei and Gojobori, 1986). The nucleotide sequences from 23 hexon, 43 fiber, and 17 penton base genes available in GenBank were included.

All the genome sequences included for substitution rate analysis were summarized in Table 1. For HAdV-B55, additional genome typing details were also presented. The sequences of the hexon, penton base, and fiber genes were either extracted from the genome files or deposited as single gene entries.

Table 1

TypeStrainYear isolatedCountrySequenceGenBank accession no.
HAdV-A12Huie1954USAGenomeAC_000005
HAdV-B3GB1953USAGenomeAY599834
HAdV-B7Gomen1952USAGenomeAY594255
HAdV-B11Slobitski1956USAGenomeNC_011202
HAdV-B11aak372001EGYGenomeJX423385
12222005SGNGenomeFJ597732
AK362005ARGGenomeJX423384
HAdV-B14de Wit1955NetherlandsGenomeAY803294
HAdV-B55BJ01#2011CHNGenomeJX491639
CQ-8142010CHNGenomeJX123027
BD67282013CHNGenomeKJ883520
TJ-2013-902013CHNGenomeKF908851
CQ-29032012CHNGenomeJX123029
QS-DLL2006CHNGenomeFJ643676
BD67292013CHNGenomeKJ883521
QZ012011CHNGenomeKJ883522
CQ-16572011CHNGenomeJX123028
XZ2012-4922012CHNGenomeKC857701
Hebei/BD012012CHNGenomeKP896478
Liaoning/LS012013CHNGenomeKP896483
Tianjin/TJ012013CHNGenomeKP896484
JS2015012015CHNGenomeKX289874
60-GD-20162016CHNGenomeKY070248
AFMC 16-00112016S. KoreaGenomeKX494979
AH-CHN/CZ-TC82012CHNHexonKC551973
HAdV11-QS2010CHNHexonDQ874353
SHX-P012011CHNHexonKC999882
Shanxi-Y162011CHNHexonKF911353
BJ102013CHNHexonKM458628
872016CHNFiberKY070249
602016CHNFiberKY070250
592016CHNFiberKY070251
812016CHNFiberKY070252
1112016CHNFiberKY070253
1612016CHNFiberKY070254
1692016CHNFiberKY070255
732016CHNFiberKY070256
1232016CHNFiberKY070257
802016CHNFiberKY070258
BJ102013CHNFiberKP270921
YT2011_12-442011CHNFiberKC510748
YT2011_12-1152011CHNFiberKC510749
BJ2011_6-802011CHNFiberKC510750
YT2011_12-312011CHNFiberKC510751
BJ2011_10-542011CHNFiberKC510752
BJ2011_7-352011CHNFiberKC510753
BJ2011_10-252011CHNFiberKC510754
BJ2011_1-322011CHNFiberKC510755
BJ2011_4-282011CHNFiberKC510756
BJ2011_10-392011CHNFiberKC510757
BJ2011_1-312011CHNFiberKC510758
BJ2011_1-292011CHNFiberKC510759
BJ2011_1-1002011CHNFiberKC510760
BJ2011_1-982011CHNFiberKC510761
YT2011_12-422011CHNFiberKC510762
HAdV-B16ch.791955USAGenomeAY601636
HAdV-B21AV-16451956Saudi ArabiaGenomeAY601633
SAdV-21Bertha1954USAGenomeAC_000010
HAdV-B34Compton1972USAGenomeAY737797
HAdV-B35Holden1973USAGenomeAY128640
HAdV-B50Wan1988USAGenomeAY737798
HAdV-B6687-9221987ARGGenomeJN860676
HAdV-B68Arg 827/042004ARGGenomeJN860678
HAdV-C1Adenoid 711953USAGenomeAF534906
HAdV-D9Hicks1954USAGenomeAJ854486
HAdV-E4RI-671952USAGenomeAY594253
HAdV-F40Dugan1979NetherlandsGenomeNC_001454
HAdV-G52T03-22442003USAGenomeDQ923122

The genome, hexon, penton base and fiber sequences of adenovirus species A-G used in this study.

They are either deposited as single gene entries in GenBank or extracted from the genome files. #The HAdV-B55 genome analyzed in the present study. They are likely typed as HAdV-B55 according to the genome sequences.

Sequences Used in the Study

The HAdV sequences of the hexon, penton base, fiber genes, and the genomes used for phylogenic analyses were retrieved from GenBank and summarized in Table 1, with additional sequence origin details if available (strain names, countries, years and GenBank accession numbers). The prototype strains were also used.

Results

Nucleotide Sequence Analysis of HAdV-B55 Strain BJ01

The genome data of strain BJ01 was deposited into GenBank (accession number JX491639), and named formally as “Human adenovirus 55 isolate HAdV-B/CHN/BJ01/2011/55[P14H11F14]”, further referred to as “BJ01”.

Figure 1 presented the genomic organization and transcription map of strain BJ01. The genome contained four early, two intermediate, and five late transcription units, and 34,773 bp in length, with a base composition of 26.16% A, 25.08% T, 24.38% G, 24.39% C, and a GC content of 48.76%, which was similar with the other members of subspecies B2 (mean of 49%) (Seto et al., 2010a). A total of 39 coding sequences were identified (Table 2). The 137-bp ITR of strain BJ01 was identical to that of the first HAdV-B55 strain (QS-DLL; 2006) (Walsh et al., 2010).

Table 2

RegionProductLocation
E1A29.1 kDa protein58–1166,1251–1459
E1A25.7 kDa protein587–1073,1251–1459
E1A6.5 kDa protein587–658,1251–1355
E1B20 kDa protein1629–2171
E1B54.9 kDa protein1934–3418
pIXpIX protein3501–3920
IVa2IVa2 protein3989–5322,5601–5613
E2BDNA polymerase5092–8664,13645–13653
E2BpTP8463–10424,13645–13653
L143 kDa protein10671–11831
L1pIIIa11857–13620
L2penton base protein13701–15374
L2protein VII15379–15957
L2protein V precursor16000–17055
L2protein X17084–17314
L3protein VI17395–18135
L3hexon protein18251–21091
L323 kDa protein21131–21760
E2ADNA binding protein21838–23394
L4100 kDa hexon-assembly associated protein23425–25863
L433 kDa protein25595–25913,26083–26444
L422 kDa protein25595–26170
L4protein VIII26494–27177
E311.7 kDa protein27177–27494
E314.6 kDa protein27448–27843
E318.4 kDa protein27828–28328
E320.1 kDa protein28348–28893
E320.8 kDa protein28911–29462
E310.1 kDa protein29506–29781
E314.9 kDa protein29786–30190
E315 kDa protein30183–30590
UU protein30614..30778
L5fiber protein30793–31770
E4Orf6/7 protein31806–32057,32780–32953
E4Orf6 protein32054–32953
E4Orf4 protein32856–33224
E4Orf3 protein33233–33586
E4Orf2 protein33583–33972
E4Orf1 protein34017–34394

The genome sequence annotation of HAdV-B55 strain BJ01.

FIGURE 1

Protein Sequence Analysis of HAdV-B55 Strain BJ01

The amino acid sequence percent identities of the proteins of HAdV-B55 strain BJ01 were presented in Table 3. Selected proteins of the representative prototypes in each HAdV species, including all species B adenoviruses, were used for comparison. The amino acid sequence analysis revealed that the proteins of strain BJ01 had the highest amino acid percent identities with that of HAdV-B14 (strain de Wit) except for the hexon protein, which is responsible for virus characterization by serum neutralization. The hexon protein of strain BJ01 had the highest amino acid identity with HAdV-B11 (98.4%), rather than HAdV-B14 (92.2%). The amino acid sequence identities indicated that strain BJ01 is a recombinant between HAdV-B14 and HAdV-B11, with the HAdV-B11 hexon gene inserted chimeric into the HAdV-B14 genome backbone (Table 3). Interestingly, E3 11.7 kDa protein of HAdV-B55 strain BJ01 had the same amino acid sequence identity with both HAdV-B14 and HAdV-B35, indicating that these three types share the conserved and similar E3 11.7 kDa protein. Among all the proteins analyzed, the highest sequence diversity was found in the fiber gene of all the seven species.

Table 3

ProteinE1AE1BE2BL1L2L3E2AL4E3L5E4
29.1 kDa20 kDaDNA polymerasepIIIaPenton baseHexonDBPpVIII11.7 kDaFiberORF6

HAdV-B3595.898.393.299.598.694.462.999.199.162.197.7
HAdV-B3497.799.498.799.395.391.399.499.198.162.297.7
HAdV-B1497.7100.099.799.899.592.299.899.699.199.199.3
HAdV-B1196.698.393.099.398.498.499.299.198.192.398.0
HAdV-B5079.488.990.891.291.890.783.494.789.661.596.0
HAdV-B2179.489.090.690.392.092.483.494.389.661.596.0
HAdV-B1679.888.991.092.985.285.383.494.789.651.489.0
HAdV-B779.487.889.992.785.386.682.694.390.691.196.7
HAdV-B379.087.289.993.085.586.383.894.389.656.797.3
SAdV-B2179.887.991.592.091.590.589.095.289.655.284.1
HAdV-B6679.487.890.792.985.586.382.494.790.656.197.0
HAdV-B6869.588.990.792.985.285.783.294.789.651.489.0
HAdV-A1242.442.270.271.072.276.847.177.461.016.849.5
HAdV-C136.349.274.074.568.476.655.379.751.417.59.3
HAdV-D941.653.074.379.076.381.259.880.665.126.364.6
HAdV-E454.760.084.786.782.782.372.789.974.523.070.0
HAdV-F4038.146.170.170.472.578.645.079.014.820.645.3
HAdV-G5237.746.672.470.572.778.748.577.757.124.648.0

Amino acid percent identities of select HAdV-55 strain BJ01 proteins with representative HAdVs from all species, including all species B adenoviruses.

The proteins above span the whole genome. The highest values from each selected protein are in bold.

Nucleotide Substitution Rates and Selection Pressures for HAdV-B55 Major Capsid Protein Genes

The synonymous and non-synonymous substitutions as well as the selection pressures for the three major capsid protein genes of HAdV-B55 were examined and compared. All the dN/dS ratios of the three genes were less than 1 (Table 4). This is to be expected because the evolution is dominated by negative selection, which removes mutations harmful to fitness (Zhao et al., 2014). Specifically, the hexon gene had least non-synonymous substitutions per site and the lowest ratio of dN/dS (0.06254). The non-synonymous substitutions and dN/dS ratio of the fiber gene was also low. However, it was relatively higher than that of the penton base and hexon genes. Overall, the major capsid proteins of HAdV-B55, hexon, penton base, and fiber, are highly conserved when compared with other types, such as HAdV-B7 (Zhao et al., 2014).

Table 4

GeneLength(bp)No. of sequencesMean Non-synonymous Substations/siteMean Synonymous Substitutions/sitedN/dS
Hexon2839230.000040.000640.06254
Penton base1674170.000480.000590.81874
Fiber978430.000370.000390.96781

Nucleotide substitution rate and selection pressures for the HAdV-55 major capsid protein genes.

The mean numbers of non-synonymous (dN) and synonymous(dS) substitutions per site from between sequences were noted and dN/dS were calculated. This selection pressure analysis was conducted using the maximum-composite-likelihood model, and included nucleotide sequences from 23 hexon, 43 fiber, and 17 penton base genes of HAdV-55 available from GenBank. All positions containing gaps were eliminated automatically. Evolutionary analyses were performed with MEGA 5.05 version.

Genome Recombination Analysis of HAdV-B55 Strain BJ01

The software mVISTA-LAGAN was used to compare the genome sequence identity between HAdV-B55 strain BJ01, HAdV-B11, and -B14. The visualization alignment showed that HAdV-B55 BJ01 had higher genome sequence identity with HAdV-B14 (98.879%) globally than with HAdV-B11 (97.605%), except for the L3 region (the hexon gene), where strain BJ01 had higher sequence identity with HAdV-B11 than HAdV-B14 (Figure 2A). This indicates that HAdV-B55 BJ01 was a recombinant composed of the HAdV-B14 genome backbone and the chimeric HAdV-B11 hexon.

FIGURE 2

Recombination analysis of both the genome and the hexon gene of HAdV-B55 strain BJ01 along with other HAdV-B viruses were shown in Figures 2BD. Simplot analysis strongly indicates that there was a recombination event between HAdV-B14 and HAdV-B11 (Figure 2B). Apparently, except for the hexon region, BJ01 had the highest similarity with HAdV-B14 instead of HAdV-B11 (Figure 2B). Bootscan analysis of the genome further confirmed the recombination (Figure 2C). The fine-resolution Bootscan analysis of the divergent hexon genes showed that the recombination event was within the hexon gene. Loop 1 (HVR1–6; nt 407 to 1,034) and Loop 2 (HVR7; nt 1,363 to 1,484) were derived from HAdV-B11 (the green regions). The other portions were derived from HAdV-B14, so was the rest of the genome (Figure 2D).

Phylogenetic Analysis of the Whole Genome, the Hexon, Penton Base, and Fiber Genes of the HAdV-B55

The phylogenetic analysis of the whole genomic sequences confirmed that the genomes of HAdV-B members were closely related to each other, forming a clade (Figure 3A). HAdV-B55 was closer to HAdV-B14 than HAdV-B11 genome (Figure 3A), which was consistent with the global pairwise alignment result in Figure 2A and the genome percent identities between the two types (98.879% vs 97.605%). The phylogenetic analysis of the fiber base and penton base genes drew the same conclusion as the whole genomes: they were closer to HAdV-B14 than HAdV-B11 (Figures 3C,D). However, the phylogenetic tree of the hexon genes showed a different distribution: HAdV-B55 was closer to HAdV-B11 than HAdV-B14 (Figure 3B), which was consistent with the Bootscan analysis in Figures 2C,D as well as with the amino acid percent identities of the hexon protein (the highest identity with HAdV-B11: 98.4%) (Table 3).

FIGURE 3

Comparative Genomic Analysis of the First Two HAdV-B55 Strains Emerged in China: QS-DLL and BJ01

Comparative genomics analysis showed that strain BJ01 had a high genome identity with strain QS-DLL (99.87%). However, when compared between the two genomes in high resolution detail, 22 single-nucleotide substitutions, 3 insertions and one deletion were identified between the two genomes. (Table 5).

Table 5

RegionProductPositionMutation in DNASubstitution in protein
Non-synonymousSynonymous
E1ANCR1463TCCATATCCGTGT
E1B20kDa1721A→TI→L
1885T→CA→A
55kDa2170T→CS→P
2823A→GK→K
IXNCR3910AAAAAA
E2BDNA polymerase8352A→G|→|
pTP10309G→AP→S
L143kDa11293G→AA→A
L2pV16892C→TR→R
NCR17310-17311AA
NCR17328A
NCR17374G→A
L3pVI17383G→AD→N
17573CGA→G
17579GAR→K
17598GAQ→Q
17617GTV→F
17643TGN→K
17676GCQ→H
Hexon19861T→CL→L
E2ADNA Binding Protein22846A→CL→L
L4Hexon-assembly Associated Protein25509A→GL→L
E311.7kDa27198T→CH→H
20.1kDa28827G→AK→K
14.9kDa30017G→TA→S

Comparative genomics analysis of two HAdV-B55 strains: QS-DLL and BJ01.

Nucleic acid and amino acid sequence changes are noted, along with their genome locations and coding consequences. The reference genome is HAdV-B55 strain QS-DLL. NCR, non-coding region; , insertion; , deletion; –, no change or not applicable. The major non-synonymous substitutions in pVI are italicized.

Twenty-one of twenty-two single-nucleotide substitutions were located in coding sequences from eight regions (E1B, E2B, L1, L2, L3, E2A, L4, E3). Ten single-nucleotide substitutions were non-synonymous substitutions, while 11 were synonymous substitutions (Table 5). Of the 10 non-synonymous substitutions, six were located in the pVI protein in the L3 region (Figure 4), which is one of the four cement proteins of HAdVs (IIIa, VI, VIII, and IX). Three of the six single-nucleotide non-synonymous substitutions (A66G, R68K, and V81F) occurred close to the N-terminal domain of protein VI, of which, two (A66G and R68K) were located within one of the two hexon binding domains. The other four non-synonymous substitutions were located in the 20kDa Protein, 55kDa Protein, Terminal Protein Precursor, and 14.9kDa Protein, respectively (Table 5). The 11 synonymous nucleotide substitutions were present in the 20kDa Protein, 55kDa Protein, DNA polymerase, 43kDa Protein, pV, pVI, Hexon, DNA Binding Protein, 100kDa Hexon-assembly Associated Protein, 11.7kDa Protein, and 20.1kDa Protein, respectively (Table 5). The capsid protein genes were highly conserved between the two strains isolated 5 years apart. There was only one single-nucleotide mutation (T to C) in the hexon gene, which resulted in a synonymous substitution (L to L). Both the penton base and fiber genes of the two strains were identical with each other.

FIGURE 4

Compared with the genome of strain QS-DLL, strain BJ01 had nucleotide insertions in three regions (TCCATATCCGTGT, AAAAAA, and A, respectively) and one deletion (AA) within L2 pX poly(A) region. All the insertions and deletions were located in non-coding regions (NCR). The insertions of “TCCATATCCGTGT” were located at the upstream of E1A poly(A) region, whereas the insertions of “AAAAAA” were located within IX poly(A) region, and the insertion of “A” were located within L2 pX poly(A) region (Table 5).

Discussion

Human adenoviruses were traditionally classified primarily according to their serological profile and hemagglutination properties (Hierholzer et al., 1991), and these methods are mainly based on the interaction between antibodies and the major capsid proteins. Presently, 51 HAdV serotypes have been recognized. However, serological methods seem not appropriate to characterize adenoviruses just by the epitopes. As an example, HAdV-D15, HAdV-D29, and HAdV-D56 were found to share the same virus neutralization (VN) epitopes (Stevens et al., 1967; Wigand, 1987; Robinson et al., 2011, 2013a), but all three viruses contain divergent genomes and biological properties (Robinson et al., 2011, 2013a). Currently, more and more traditional techniques have been replaced by genomic criteria. The genome data were proposed to be used to type HAdVs (Seto et al., 2011), i.e., HAdVs should be identified, typed, and characterized based on the complete genome sequences instead of serological methods (Zhang et al., 2006, 2012a; Robinson et al., 2009; Walsh et al., 2010; Zhou et al., 2012; Zhang Q. et al., 2016).

HAdV-B55 was first isolated from an ARD outbreak in a military training base in 1969, then it was named HAdV-B11a by partial characterization of its hexon and fiber genes (Hierholzer et al., 1974), as well as restriction enzyme patterns (Zhang Q. et al., 2016). Until 2006, along with several HAdV-B55-associated ARD outbreaks (Hierholzer et al., 1974; Chmielewicz et al., 2005; Yang et al., 2009; Kajon et al., 2010), it was ultimately revised as HAdV-B55 by Walsh et al. (2010) after bioinformatic analysis of the complete genome of strain QS-DLL. In 2011, 5 years after the emergence of the first HAdV-B55 strain QS-DLL in China (Walsh et al., 2010), HAdV-B55 re-emerged and caused a CAP outbreak in China (Gu et al., 2012; Zhang et al., 2012a; Cao et al., 2014). Since then, there have been more than 10 HAdV-B55-associated ARD outbreaks in China, for example, in Shanxi province, Chongqing, Beijing, Jiangsu province, Hebei province, Jinan, Tianjin, and Guangdong province (Yang et al., 2009; Walsh et al., 2010; Zhang et al., 2012a; Li et al., 2014; Chen et al., 2016; Gao et al., 2016; Tan et al., 2016; Yi et al., 2017), with severe and fatal cases (Zhu et al., 2009; Cao et al., 2014; Tan et al., 2016; Zhang S.Y. et al., 2016). Similar ARD outbreaks were also reported in other countries (Geng et al., 2013; Lafolie et al., 2016; Salama et al., 2016; Park et al., 2017).

Here presented is the detailed comparative genomic analysis of the re-emergent HAdV-B55 strain BJ01 with the first isolate QS-DLL. The genomic nucleotide identity between the two strains was 99.87%, along with 26 mutations, including 11 synonymous and 10 non-synonymous substitutions, 3 insertions, one deletion and one G to A mutation in L2 NCR (Table 5). Compared with other HAdV types, such as HAdV-B3, -C5, -B7, -B14, -D36, of which the genomes are usually conserved and stable (Mahadevan et al., 2010; Seto et al., 2010b; Nam et al., 2014; Alkhalaf et al., 2015; Yu et al., 2016; Zhang et al., 2017), HAdV-B55 showed similarly conserved genomes when looking at the two strains isolated five years apart (2006 vs 2011). However, the major population distribution of HAdV-B55 associated disease has been changed from juveniles to adults, whereas the disease severity remains similar. The major gene mutations between the two isolates are located in the L3 pVI protein, one of the structural proteins of HAdVs, in which six single-nucleotide non-synonymous substitutions and one synonymous substitution were identified (Figure 4). pVI and its processed form protein VI are special proteins because they have different functions at various stages of HAdV infections. For example, during adenoviral infection, protein VI aids the virion particle to escape from the endosome by inducing a pH-independent disruption of the membrane (Wiethoff et al., 2005). The N-terminal domain in pVI with predicted amphipathic α-helical structure is required for membrane lytic activity (Wiethoff et al., 2005). Additionally, it binds the hexon protein and facilitate nuclear import of hexon (Matthews and Russell, 1995). Two regions, between amino acid residues 48–74 and 233–239 of pVI, were required for the interaction with the hexon protein, where two non-synonymous substitutions were identified (Figure 4). Whether the non-synonymous substitutions in pVI can speed up the adenoviral maturation and further change the population distribution of HAdV-B55 infection remains unknown. Further wet-bench confirmation is necessary.

Phylogenetic analysis of genome sequences, as well as individual genes, showed a subclade composed of HAdV-B55 isolates (BJ01 included) with high bootstrap values (Figures 3A–D). HAdV-B55 was closely related to HAdV-B11 (Slobitski) and HAdV-B14 (de Wit), forming a clade (Figures 3A–D), which was consistent with the global pairwise alignment data (Figure 2A). Recombination analysis further revealed the mechanism that the two loops of the HAdV-B11 hexon were embedded in the whole genome of HAdV-B14 (Figures 2B–D), giving evidence that the hexon gene, particularly hypervariable regions, acts as an important hotspot for genome recombination. Although the penton base and shaft/knob boundary of the fiber gene were also known to be two recombination sites (Robinson et al., 2009; Liu et al., 2011), no recombination in these regions was identified in HAdV-B55 genomes. The neutralizing and protective antibodies in human are induced primarily by the hexon protein, partially by the fiber and penton base proteins (Gahery-Segard et al., 1998; Hong et al., 2003; Sumida et al., 2005). The lack of non-synonymous substitutions in the three capsid proteins of HAdV-B55 strains make the HAdV-B55 vaccine development easier when compared with other respiratory viruses, such as influenza viruses.

Genome recombination is an important strategy driving HAdV evolution (Robinson et al., 2009, 2011; Kajon et al., 2010; Walsh et al., 2010, 2011; Liu et al., 2011, 2012; Matsushima et al., 2012), and results in the changes in pathogenicity (Walsh et al., 2009; Singh et al., 2012; Zhou et al., 2012). For HAdVs, to some extent, recombination can be defined as an adaptive change resulting from environmental pressures. In the example of HAdV-B55, the initial HAdV-B14 infection may induce neutralizing antibodies against this type. Meanwhile the co-infection and co-replication of HAdV-B11 and HAdV-B14 in the intestine may lead to the recombination between HAdV-B11 and HAdV-B14 hexon genes, i.e., the resulting HAdV-B55 can escape from pre-existing immunity against HAdV-B14 from the host due to its acquisition of HAdV-B11 loops. The previous study has confirmed that the neutralizing and protective antibodies in human are induced primarily by the two hexon loops. The acquisition both loops of HAdV-B11 is beneficial to the survival and spreading of the recombinant HAdV-B55, because the previous absence of circulating HAdV-B11 in China (Han et al., 2013; Chen et al., 2016; Yu et al., 2016) led to a corresponding low level of herd immunity, which might be associated with the subsequent more and more ARD outbreaks caused by HAdV-B55 (Yang et al., 2009; Zhu et al., 2009; Walsh et al., 2010; Zhang et al., 2012a; Zhang S.Y. et al., 2016; Cao et al., 2014; Li et al., 2014; Chen et al., 2016; Gao et al., 2016; Tan et al., 2016; Yi et al., 2017). This is one of the immune escape capabilities of HAdVs. The detailed recombination mechanism of HAdV-B55 might be illustrated by in vitro co-infection and replication of the two types of adenoviruses, HAdV-B11 and -B14. In addition, we also infected A549, Hela, HEp-2, Hep-G2, and Vero cells with HAdV-B55, -B11, and -B14, respectively (data not shown). The replication efficiency of different types in the cells was compared. The three types of HAdVs replicated most efficiently in A549 cells. Unfortunately, no significant difference of replication efficiency among the three types of HAdVs was identified in all the tested cell clines. The prevalence and outbreak of HAdV-B55 might be not associated with the virus replication efficiency.

In summary, comparative genomic analysis of re-emergent HAdV-B55 pathogens revealed conserved capsid proteins and low non-synonymous substitutions per site of HAdV-B55, which provides the foundation for the development of effective vaccines against this pathogen. It also facilitates the further studies of the adenovirus epidemiology, molecular evolution, patterns of pathogen emergence and re-emergence (Zhang et al., 2017), and the predication of genome recombination between adenoviruses, as well as the development of public health strategies, including surveillance.

Statements

Author contributions

QZ, XW, and BC conceived and designed the study. ZC, YY, SJ, W-GL, W-WC, JZ, ML, ShZ, NC, JO, and SuZ performed the computational analyses and made the figures and tables. ZC, YY, and QZ wrote the manuscript. All authors reviewed the manuscript.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (31570155, 31370199, 81471961) and “Young Top-notch Talents” of the Guangdong Province Special Support Program (2014), as well as the Excellent Young Teacher Training Plan of Guangdong Province (Yq2013039) and Guangzhou Healthcare Collaborative Innovation Major Project (201803040004, 201803040007).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    Alissa AlkhalafM.Al QurashiY. M.GuiverM.CooperR. J. (2014). Genome sequences of three species d adenoviruses isolated from AIDS patients.Genome Announc.2:e01267-13.10.1128/genomeA.01267-13

  • 2

    AlkhalafM. A.GuiverM.CooperR. J. (2015). Genome stability of adenovirus types 3 and 7 during a simultaneous outbreak in Greater Manchester. UK.J. Med. Virol.87117124. 10.1002/jmv.23969

  • 3

    BrudnoM.DoC. B.CooperG. M.KimM. F.DavydovE.GreenE. D.et al (2003). LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA.Genome Res.13721731. 10.1101/gr.926603

  • 4

    CaoB.HuangG. H.PuZ. H.QuJ. X.YuX. M.ZhuZ.et al (2014). Emergence of community-acquired adenovirus type 55 as a cause of community-onset pneumonia.Chest1457986. 10.1378/chest.13-1186

  • 5

    ChenY.LiuF.WangC.ZhaoM.DengL.ZhongJ.et al (2016). Molecular identification and epidemiological features of human adenoviruses associated with acute respiratory infections in hospitalized children in Southern China, 2012-2013.PLoS One11:e0155412. 10.1371/journal.pone.0155412

  • 6

    ChmielewiczB.BenzlerJ.PauliG.KrauseG.BergmannF.SchweigerB. (2005). Respiratory disease caused by a species B2 adenovirus in a military camp in Turkey.J. Med. Virol.77232237. 10.1002/jmv.20441

  • 7

    ChuW.Pavan-LangstonD. (1979). Ocular surface manifestations of the major viruses.Int. Ophthalmol. Clin.19135167. 10.1097/00004397-197919020-00010

  • 8

    DehghanS.LiuE. B.SetoJ.TorresS. F.HudsonN. R.KajonA. E.et al (2012). Five genome sequences of subspecies B1 human adenoviruses associated with acute respiratory disease.J. Virol.86635636. 10.1128/JVI.06593-11

  • 9

    DhurandharN. V.KulkarniP.AjinkyaS. M.SherikarA. (1992). Effect of adenovirus infection on adiposity in chicken.Vet. Microbiol.31101107. 10.1016/0378-1135(92)90068-5

  • 10

    DingleJ. H.LangmuirA. D. (1968). Epidemiology of acute, respiratory disease in military recruits.Am. Rev. Respir. Dis.97(Suppl.), 165.

  • 11

    EchavarriaM. (2008). Adenoviruses in immunocompromised hosts.Clin. Microbiol. Rev.21704715. 10.1128/CMR.00052-07

  • 12

    Gahery-SegardH.FaraceF.GodfrinD.GastonJ.LengagneR.TurszT.et al (1998). Immune response to recombinant capsid proteins of adenovirus in humans: antifiber and anti-penton base antibodies have a synergistic effect on neutralizing activity.J. Virol.7223882397.

  • 13

    GaoH. W.WeiM. T.FanH. J.DingH.WeiW.LiuZ. Q.et al (2016). ). Dynamic changes in clinical characteristics during an outbreak of human adenovirus Serotype 55 in China.Disaster Med. Public Health Prep. [Epub ahead of print].

  • 14

    GengX.ZhangJ.YangG. (2013). Investigation of an adenovirus-induced respiratory disease outbreak.Adv. Infect. Dis.3257262. 10.4236/aid.2013.34039

  • 15

    GuL.LiuZ.LiX.QuJ.GuanW.LiuY.et al (2012). Severe community-acquired pneumonia caused by adenovirus type 11 in immunocompetent adults in Beijing.J. Clin. Virol.54295301. 10.1016/j.jcv.2012.04.018

  • 16

    HageE.LiebertU. G.BergsS.GanzenmuellerT.HeimA. (2015). Human Adenovirus type 70: a novel, multiple recombinant species D adenovirus isolated from diarrheal faeces of a haematopoietic stem cell transplantation recipient.J. Gen. Virol.9627342742. 10.1099/vir.0.000196

  • 17

    HakimF. A.TleyjehI. M. (2008). Severe adenovirus pneumonia in immunocompetent adults: a case report and review of the literature.Eur. J. Clin. Microbiol. Infect. Dis.27153158. 10.1007/s10096-007-0416-z

  • 18

    HanG.NiuH.ZhaoS.ZhuB.WangC.LiuY.et al (2013). Identification and typing of respiratory adenoviruses in Guangzhou, Southern China using a rapid and simple method.Virol. Sin.28103108. 10.1007/s12250-013-3308-7

  • 19

    HashimotoS.GonzalezG.HaradaS.OosakoH.HanaokaN.HinokumaR.et al (2018). Recombinant type Human mastadenovirus D85 associated with epidemic keratoconjunctivitis since 2015 in Japan.J. Med. Virol.90881889. 10.1002/jmv.25041

  • 20

    HierholzerJ. C.PumarolaA.Rodriguez-TorresA.BeltranM. (1974). Occurrence of respiratory illness due to an atypical strain of adenovirus type 11 during a large outbreak in Spanish military recruits.Am. J. Epidemiol.99434442. 10.1093/oxfordjournals.aje.a121632

  • 21

    HierholzerJ. C.StoneY. O.BrodersonJ. R. (1991). Antigenic relationships among the 47 human adenoviruses determined in reference horse antisera.Arch. Virol.121179197. 10.1007/BF01316753

  • 22

    HillemanM. R.WernerJ. H.GauldR. L. (1953). Influenza antibodies in the population of the USA; an epidemiological investigation.Bull. World Health Organ.8613631.

  • 23

    HongS. S.HabibN. A.FranquevilleL.JensenS.BoulangerP. A. (2003). Identification of adenovirus (ad) penton base neutralizing epitopes by use of sera from patients who had received conditionally replicative ad (addl1520) for treatment of liver tumors.J. Virol.771036610375. 10.1128/JVI.77.19.10366-10375.2003

  • 24

    IshikoH.ShimadaY.KonnoT.HayashiA.OhguchiT.TagawaY.et al (2008). Novel human adenovirus causing nosocomial epidemic keratoconjunctivitis.J. Clin. Microbiol.4620022008. 10.1128/JCM.01835-07

  • 25

    JonesM. S.IIHarrachB.GanacR. D.GozumM. M.Dela CruzW. P.RiedelB.et al (2007). New adenovirus species found in a patient presenting with gastroenteritis.J. Virol.8159785984. 10.1128/JVI.02650-06

  • 26

    KajonA. E.DicksonL. M.MetzgarD.HoungH. S.LeeV.TanB. H. (2010). Outbreak of febrile respiratory illness associated with adenovirus 11a infection in a Singapore military training cAMP.J. Clin. Microbiol.4814381441. 10.1128/JCM.01928-09

  • 27

    KojaoghlanianT.FlomenbergP.HorwitzM. S. (2003). The impact of adenovirus infection on the immunocompromised host.Rev. Med. Virol.13155171. 10.1002/rmv.386

  • 28

    LafolieJ.MirandA.SalmonaM.LautretteA.ArchimbaudC.BrebionA.et al (2016). Severe pneumonia associated with adenovirus type 55 infection, France, 2014.Emerg. Infect. Dis.2220122014. 10.3201/eid2211.160728

  • 29

    LeC. T.GrayG. C.PoddarS. K. (2001). A modified rapid method of nucleic acid isolation from suspension of matured virus: applied in restriction analysis of DNA from an adenovirus prototype strain and a patient isolate.J. Med. Microbiol.50571574. 10.1099/0022-1317-50-6-571

  • 30

    LiX.KongM.SuX.ZouM.GuoL.DongX.et al (2014). An outbreak of acute respiratory disease in China caused by human adenovirus type B55 in a physical training facility.Int. J. Infect. Dis.28117122. 10.1016/j.ijid.2014.06.019

  • 31

    LiuE. B.FerreyraL.FischerS. L.PavanJ. V.NatesS. V.HudsonN. R.et al (2011). Genetic analysis of a novel human adenovirus with a serologically unique hexon and a recombinant fiber gene.PLoS One6:e24491. 10.1371/journal.pone.0024491

  • 32

    LiuE. B.WadfordD. A.SetoJ.VuM.HudsonN. R.ThrasherL.et al (2012). Computational and serologic analysis of novel and known viruses in species human adenovirus D in which serology and genomics do not correlate.PLoS One7:e33212. 10.1371/journal.pone.0033212

  • 33

    LoleK. S.BollingerR. C.ParanjapeR. S.GadkariD.KulkarniS. S.NovakN. G.et al (1999). Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination.J. Virol.73152160.

  • 34

    LuQ. B.TongY. G.WoY.WangH. Y.LiuE. M.GrayG. C.et al (2014). Epidemiology of human adenovirus and molecular characterization of human adenovirus 55 in China, 2009-2012.Influenza Other Respir. Viruses8302308. 10.1111/irv.12232

  • 35

    MahadevanP.SetoJ.TibbettsC.SetoD. (2010). Natural variants of human adenovirus type 3 provide evidence for relative genome stability across time and geographic space.Virology397113118. 10.1016/j.virol.2009.10.052

  • 36

    MatsushimaY.ShimizuH.KanoA.NakajimaE.IshimaruY.DeyS. K.et al (2012). Novel human adenovirus strain, Bangladesh.Emerg. Infect. Dis.18846848. 10.3201/eid1805.111584

  • 37

    MatsushimaY.ShimizuH.KanoA.NakajimaE.IshimaruY.DeyS. K.et al (2013). Genome sequence of a novel virus of the species human adenovirus d associated with acute gastroenteritis.Genome Announc.1:e00068-12. 10.1128/genomeA.00068-12

  • 38

    MatsushimaY.ShimizuH.PhanT. G.UshijimaH. (2011). Genomic characterization of a novel human adenovirus type 31 recombinant in the hexon gene.J. Gen. Virol.9227702775. 10.1099/vir.0.034744-0

  • 39

    MatthewsD. A.RussellW. C. (1995). Adenovirus protein-protein interactions: molecular parameters governing the binding of protein VI to hexon and the activation of the adenovirus 23K protease.J. Gen. Virol. 76(Pt 8), 19591969. 10.1099/0022-1317-76-8-1959

  • 40

    NamJ. H.NaH. N.AtkinsonR. L.DhurandharN. V. (2014). Genomic stability of adipogenic human adenovirus 36.Int. J. Obes.38321324. 10.1038/ijo.2013.67

  • 41

    NeiM.GojoboriT. (1986). Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions.Mol. Biol. Evol.3418426.

  • 42

    ParkJ. Y.KimB. J.LeeE. J.ParkK. S.ParkH. S.JungS. S.et al (2017). Clinical Features and Courses of adenovirus pneumonia in healthy young adults during an outbreak among Korean military personnel.PLoS One12:e0170592. 10.1371/journal.pone.0170592

  • 43

    RobinsonC. M.RajaiyaJ.WalshM. P.SetoD.DyerD. W.JonesM. S.et al (2009). Computational analysis of human adenovirus type 22 provides evidence for recombination among species D human adenoviruses in the penton base gene.J. Virol.8389808985. 10.1128/JVI.00786-09

  • 44

    RobinsonC. M.SinghG.HenquellC.WalshM. P.Peigue-LafeuilleH.SetoD.et al (2011). Computational analysis and identification of an emergent human adenovirus pathogen implicated in a respiratory fatality.Virology409141147. 10.1016/j.virol.2010.10.020

  • 45

    RobinsonC. M.SinghG.LeeJ. Y.DehghanS.RajaiyaJ.LiuE. B.et al (2013a). Molecular evolution of human adenoviruses.Sci. Rep.3:1812. 10.1038/srep01812

  • 46

    RobinsonC. M.ZhouX.RajaiyaJ.YousufM. A.SinghG.DeserresJ. J.et al (2013b). Predicting the next eye pathogen: analysis of a novel adenovirus.mBio4:e00595-12. 10.1128/mBio.00595-12

  • 47

    RoweW. P.HuebnerR. J.GilmoreL. K.ParrottR. H.WardT. G. (1953). Isolation of a cytopathogenic agent from human adenoids undergoing spontaneous degeneration in tissue culture.Proc. Soc. Exp. Biol. Med.84570573. 10.3181/00379727-84-20714

  • 48

    SalamaM.AmitaiZ.AmirN.Gottesman-YekutieliT.SherbanyH.DroriY.et al (2016). Outbreak of adenovirus type 55 infection in Israel.J. Clin. Virol.783135. 10.1016/j.jcv.2016.03.002

  • 49

    SetoD.ChodoshJ.BristerJ. R.JonesM. S.Members of the Adenovirus Research Community (2011). Using the whole-genome sequence to characterize and name human adenoviruses.J. Virol.8557015702. 10.1128/JVI.00354-11

  • 50

    SetoJ.WalshM. P.MahadevanP.ZhangQ.SetoD. (2010a). Applying genomic and bioinformatic resources to human adenovirus genomes for use in vaccine development and for applications in vector development for gene delivery.Viruses2126. 10.3390/v2010001

  • 51

    SetoJ.WalshM. P.MetzgarD.SetoD. (2010b). Computational analysis of adenovirus serotype 5 (HAdV-C5) from an HAdV coinfection shows genome stability after 45 years of circulation.Virology404180186. 10.1016/j.virol.2010.05.010

  • 52

    SinghG.RobinsonC. M.DehghanS.JonesM. S.DyerD. W.SetoD.et al (2013). Homologous recombination in E3 genes of human adenovirus species D.J. Virol.871248112488. 10.1128/JVI.01927-13

  • 53

    SinghG.RobinsonC. M.DehghanS.SchmidtT.SetoD.JonesM. S.et al (2012). Overreliance on the hexon gene, leading to misclassification of human adenoviruses.J. Virol.8646934695. 10.1128/JVI.06969-11

  • 54

    StevensD. A.SchaefferM.FoxJ. P.BrandtC. D.RomanoM. (1967). Standardization and certification of reference antigens and antisera for 30 human adenovirus serotypes.Am. J. Epidemiol.86617633. 10.1093/oxfordjournals.aje.a120771

  • 55

    SumidaS. M.TruittD. M.LemckertA. A.VogelsR.CustersJ. H.AddoM. M.et al (2005). Neutralizing antibodies to adenovirus serotype 5 vaccine vectors are directed primarily against the adenovirus hexon protein.J. Immunol.17471797185. 10.4049/jimmunol.174.11.7179

  • 56

    TamuraK.DudleyJ.NeiM.KumarS. (2007). MEGA4: molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Mol. Biol. Evol.2415961599. 10.1093/molbev/msm092

  • 57

    TamuraK.PetersonD.PetersonN.StecherG.NeiM.KumarS. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.Mol. Biol. Evol.2827312739. 10.1093/molbev/msr121

  • 58

    TanD.ZhuH.FuY.TongF.YaoD.WallineJ.et al (2016). Severe community-acquired pneumonia caused by human adenovirus in immunocompetent adults: a multicenter case series.PLoS One11:e0151199. 10.1371/journal.pone.0151199

  • 59

    WadellG.HammarskjoldM. L.WinbergG.VarsanyiT. M.SundellG. (1980). Genetic variability of adenoviruses.Ann. N. Y. Acad. Sci.3541642. 10.1111/j.1749-6632.1980.tb27955.x

  • 60

    WalshM. P.ChintakuntlawarA.RobinsonC. M.MadischI.HarrachB.HudsonN. R.et al (2009). Evidence of molecular evolution driven by recombination events influencing tropism in a novel human adenovirus that causes epidemic keratoconjunctivitis.PLoS One4:e5635. 10.1371/journal.pone.0005635

  • 61

    WalshM. P.SetoJ.JonesM. S.ChodoshJ.XuW.SetoD. (2010). Computational analysis identifies human adenovirus type 55 as a re-emergent acute respiratory disease pathogen.J. Clin. Microbiol.48991993. 10.1128/JCM.01694-09

  • 62

    WalshM. P.SetoJ.LiuE. B.DehghanS.HudsonN. R.LukashevA. N.et al (2011). Computational analysis of two species C human adenoviruses provides evidence of a novel virus.J. Clin. Microbiol.4934823490. 10.1128/JCM.00156-11

  • 63

    WiethoffC. M.WodrichH.GeraceL.NemerowG. R. (2005). Adenovirus protein VI mediates membrane disruption following capsid disassembly.J. Virol.7919922000. 10.1128/JVI.79.4.1992-2000.2005

  • 64

    WigandR. (1987). Pitfalls in the identification of adenoviruses.J. Virol. Methods16161169. 10.1016/0166-0934(87)90001-2

  • 65

    WoodD. J. (1988). Adenovirus gastroenteritis.Br. Med. J. Clin. Res.296229230. 10.1136/bmj.296.6617.229-a

  • 66

    YangZ.ZhuZ.TangL.WangL.TanX.YuP.et al (2009). Genomic analyses of recombinant adenovirus type 11a in China.J. Clin. Microbiol.4730823090. 10.1128/JCM.00282-09

  • 67

    YiL.ZouL.LuJ.KangM.SongY.SuJ.et al (2017). A cluster of adenovirus type B55 infection in a neurosurgical inpatient department of a general hospital in Guangdong, China.Influenza Other Respir. Viruses11328336. 10.1111/irv.12457

  • 68

    YuZ.ZengZ.ZhangJ.PanY.ChenM.GuoY.et al (2016). Fatal community-acquired pneumonia in children caused by re-emergent human adenovirus 7d associated with higher severity of illness and fatality rate.Sci. Rep.6:37216. 10.1038/srep37216

  • 69

    ZhangS. Y.LuoY. P.HuangD. D.FanH.LuQ. B.WoY.et al (2016). Fatal pneumonia cases caused by human adenovirus 55 in immunocompetent adults.Infect. Dis.484047. 10.3109/23744235.2015.1055585

  • 70

    ZhangQ.DehghanS.SetoD. (2016). Pitfalls of restriction enzyme analysis in identifying, characterizing, typing, and naming viral pathogens in the era of whole genome data, as illustrated by HAdV type 55.Virologica Sinica31448453. 10.1007/s12250-016-3862-x

  • 71

    ZhangQ.JingS.ChengZ.YuZ.DehghanS.ShamsaddiniA.et al (2017). Comparative genomic analysis of two emergent human adenovirus type 14 respiratory pathogen isolates in China reveals similar yet divergent genomes.Emerg. Microbes Infect.6:e92. 10.1038/emi.2017.78

  • 72

    ZhangQ.SetoD.CaoB.ZhaoS.WanC. (2012a). Genome sequence of human adenovirus type 55, a re-emergent acute respiratory disease pathogen in China.J. Virol.861244112442. 10.1128/JVI.02225-12

  • 73

    ZhangQ.SetoD.ZhaoS.ZhuL.ZhaoW.WanC. (2012b). Genome sequence of the first human adenovirus type 14 isolated in China.J. Virol.8670197020. 10.1128/JVI.00814-12

  • 74

    ZhangQ.SuX.GongS.ZengQ.ZhuB.WuZ.et al (2006). Comparative genomic analysis of two strains of human adenovirus type 3 isolated from children with acute respiratory infection in southern China.J. Gen. Virol.8715311541. 10.1099/vir.0.81515-0

  • 75

    ZhaoS.WanC.KeC.SetoJ.DehghanS.ZouL.et al (2014). Re-emergent human adenovirus genome type 7d caused an acute respiratory disease outbreak in Southern China after a twenty-one year absence.Sci. Rep.4:7365. 10.1038/srep07365

  • 76

    ZhouX.RobinsonC. M.RajaiyaJ.DehghanS.SetoD.JonesM. S.et al (2012). Analysis of human adenovirus type 19 associated with epidemic keratoconjunctivitis and its reclassification as adenovirus type 64.Invest. Ophthalmol. Vis. Sci.5328042811. 10.1167/iovs.12-9656

  • 77

    ZhuZ.ZhangY.XuS.YuP.TianX.WangL.et al (2009). Outbreak of acute respiratory disease in China caused by B2 species of adenovirus type 11.J. Clin. Microbiol.47697703. 10.1128/JCM.01769-08

Summary

Keywords

human adenovirus type 55, recombination, comparative genomics, severe community-acquired pneumonia, China

Citation

Cheng Z, Yan Y, Jing S, Li W-G, Chen W-W, Zhang J, Li M, Zhao S, Cao N, Ou J, Zhao S, Wu X, Cao B and Zhang Q (2018) Comparative Genomic Analysis of Re-emergent Human Adenovirus Type 55 Pathogens Associated With Adult Severe Community-Acquired Pneumonia Reveals Conserved Genomes and Capsid Proteins. Front. Microbiol. 9:1180. doi: 10.3389/fmicb.2018.01180

Received

17 January 2018

Accepted

15 May 2018

Published

05 June 2018

Volume

9 - 2018

Edited by

Mei-Ru Chen, National Taiwan University, Taiwan

Reviewed by

Hongbin Song, Kirov Military Medical Academy, Russia; Rui Zhang, Xiamen University, China

Updates

Copyright

*Correspondence: Xianbo Wu, Bin Cao, Qiwei Zhang,

These authors have contributed equally to this work.

This article was submitted to Virology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics