REVIEW article

Front. Microbiol., 11 September 2018

Sec. Virology

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

Adenoviromics: Mining the Human Adenovirus Species D Genome

  • 1. Howe Laboratory, Massachusetts Eye and Ear, Harvard Medical School, Boston, MA, United States

  • 2. Molecular Virology Laboratory, Korea Zoonosis Research Institute, Jeonbuk National University, Jeonju, South Korea

  • 3. Department of Cell and Systems Biology, University of Toronto, Toronto, ON, Canada

  • 4. Department of Microbiology and Immunology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, United States

  • 5. Bioinformatics and Computational Biology Program, School of Systems Biology, George Mason University, Manassas, VI, United States

Abstract

Human adenovirus (HAdV) infections cause disease world-wide. Whole genome sequencing has now distinguished 90 distinct genotypes in 7 species (A-G). Over half of these 90 HAdVs fall within species D, with essentially all of the HAdV-D whole genome sequences generated in the last decade. Herein, we describe recent new findings made possible by mining of this expanded genome database, and propose future directions to elucidate new functional elements and new functions for previously known viral components.

Introduction

Human adenovirus (HAdV) infections represent a significant source of morbidity and mortality, world-wide and at all ages, through highly transmittable infections at mucosal sites, including the eye, and urinary, respiratory, and gastrointestinal tracts (Horwitz, ). HAdV causes fatal acute respiratory distress syndrome in healthy adults and is especially lethal in infants and the immune compromised (Bhanthumkosol, ; Ryu et al., ; Wallot et al., ; Engelmann et al., ; Tan et al., ; Zhang et al., ). No FDA-approved therapy for acute HAdV infection is available. At resolution of acute infection, persistence may develop within nasopharyngeal lymphoid tissue (Neumann et al., ; Garnett et al., , ; Zhang et al., ; Assadian et al., ), as yet uncharacterized cells in the gastrointestinal tract (Roy et al., ), and possibly the ocular surface (Kaye et al., ), permitting evolution of new HAdVs through homologous recombination between two or more HAdVs infecting the same cell(s) (Lee et al., , ; Echavarria et al., ; McCarthy et al., ; Seto et al., ).

HAdVs are divided phylogenetically into seven species (A-G), with a total of 90 recognized genotypes with whole genome sequences in GenBank, including the original 51 “serotypes”—determined by serum neutralization—which now all have been fully sequenced (Table 1) (Robinson et al., ) Human adenovirus species D (HAdV-D) is the largest and most rapidly growing among all HAdV species, and contains viruses associated with epidemic keratoconjunctivitis (EKC), a severe, hyperacute ocular surface infection (Butt and Chodosh, ). A collaboration funded by the American Recovery and Reinvestment Act of 2009 came to fruition with the complete whole genome sequencing and analysis of all previously unsequenced HAdV-D serotypes (Robinson et al., ), leading to a new understanding of adenovirus ontogeny (Jones et al., ; Robinson et al., , ,; Robinson et al., ,; Robinson et al., ,; Walsh et al., , ,; Arnold et al., ; Torres et al., ; Dehghan et al., , ,; Walsh et al., ; Liu et al., ; Seto et al., , ; Singh et al., , ; Zhou et al., )—including those HAdV-Ds associated with EKC (Robinson et al., , , ; Walsh et al., ; Zhou et al., )—and ultimately to a new typing system for HAdV based on genomics (Seto et al., ).

Table 1

TypeGenBank accession no.Genome lengthYear published
HAdV-C1AC_000017360012004
HAdV-C2AC_000007359372003
HAdV-B3AY599834353452006
HAdV-E4AY599837359642006
HAdV-C5AY601635359312006
HAdV-C6FJ349096357582011
HAdV-B7KP670856.2352392016
HAdV-D8AB448767349802009
HAdV-D9AJ854486350832008
HAdV-D10JN226746351052013
HAdV-B11AF532578347942003
HAdV-A12X73487341251979
HAdV-D13JN226747352092013
HAdV-B14JQ824845347672012
HAdV-D15KF268204351002013
HAdV-B16JN860680353842011
HAdV-D17HQ910407351392011
HAdV-A18GU191019341772010
HAdV-D19JQ326209351532011
HAdV-D20JN226749351812013
HAdV-B21AY601633353822006
HAdV-D22FJ619037351522009
HAdV-D23JN226750350502013
HAdV-D24JN226751351662013
HAdV-D25JN226752352482013
HAdV-D26EF153474351522007
HAdV-D27JN226753351542013
HAdV-D28FJ824826351302010
HAdV-D29JN226754352142013
HAdV-D30JN226755351782012
HAdV-A31AM749299337632005
HAdV-D32JN226756352482013
HAdV-D33JN226758351312013
HAdV-B34AY737797347752004
HAdV-B35AC_000019347942004
HAdV-D36GQ384080351522010
HAdV-D37AB448775352152009
HAdV-D38JN226759352212013
HAdV-D39JN226760351522013
HAdV-F40NC_001454342141993
HAdV-F41DQ315364.2341882007
HAdV-D42JN226761352312013
HAdV-D43JN226762350122013
HAdV-D44JN226763352142013
HAdV-D45JN226764351542013
HAdV-D46AY875648351782006
HAdV-D47JN226757351062013
HAdV-D48EF153473352062007
HAdV-D49DQ393829352152006
HAdV-B50AY737798353852007
HAdV-D51JN226765351142013

Species and type designations for the 51 human adenovirus (HAdV) serotypes.

Recent published work demonstrates how genome “mining,” in-depth analyses of the growing HAdV genome database, can bring about new realizations and add critical new information to prior ones. The trimeric fiber protein on adenoviruses mediates viral entry through interaction of the distal most “knob” structure on the fiber with host cell receptors. In a phylogenetic analysis of HAdV-D fiber genes, HAdV-D types associated with EKC were recently shown to form a unique clade (Ismail et al., ). By proteotyping, a new in silico methodology described in detail below, EKC virus-associated fiber knobs were uniquely shared, and signature amino acid positions distinguished EKC from non-EKC types. Remarkably, human corneal epithelial cell tropism could be predicted by the presence of a lysine or alanine at residue 240, and this amino acid residue in EKC viruses showed evidence for positive selection. These data added to the prior observation by Huang and coworkers that artificial mutation to a lysine at residue 240 in a non-EKC virus could confer infection of Chang cells, a conjunctiva derived continuous cell line (Huang et al., ). However, because Chang cells came later known to be contaminated by HeLa cells, the importance of residue 240 to ocular tropism was until this new observation, in some doubt.

Another recently published effort provided further evidence of the importance and potential for HAdV genome mining. Late adenoviral gene expression is initiated by the adenovirus major late promoter (Ramke et al., ), followed by splicing of mRNAs to the viral tripartite leader for translation (Chow et al., ; Akusjärvi and Pettersson, ; Chow and Broker, ; Logan and Shenk, ). The HAdV tripartite leader is a 200-nucleotide 5' noncoding region that circumvents the requirement for eukaryotic initiation factor 4F or cap binding protein complex (Ziff and Evans, ; Akusjärvi and Pettersson, ; Dolph et al., ; Zhang et al., ), and permits translation of HAdV mRNAs at late times in infection when cap-dependent translation is blocked due to shut down of host cellular cap-dependent mRNA translation. HAdV 5′ untranslated regions (5′UTRs) are critical for cap-independent initiation, and impact mRNA localization and stability. The HAdV tripartite leader (TPL), composed of three introns (TPL 1-3), drives translation of HAdV late mRNA. The annotation of 72 HAdV genotypes for the HAdV TPL and another previously described leader, the i-leader, let to identification of newly identified polycistronic mRNAs for RID-α and RID-β within the E3 transcription unit, and a potential new open reading frame (ORF) within the i-leader sequence, with termination of this potential protein in TPL3 (Ramke et al., ). In addition, the authors also identified a potential new leader sequence embedded within the E3 region, tentatively named the j-leader (Figure 1).

Figure 1

) with permission.

Structure and infection

The HAdV is non-enveloped, icosahedral in shape, and contains a double stranded DNA genome of ~36,000 base pairs (bp) with ~1 open reading frame (ORF) for every 1000 nucleotides. Viral DNA is associated with four (interior) core proteins including Mu, VII, V, and terminal protein. The histone-like protein (p) VII protects viral DNA from cellular DNA damage responses (Lischwe and Sung, ; Karen and Hearing, ; Avgousti et al., ). The outer protein coat (capsid) of the virus consists of 240 hexon capsomers and 12 penton capsomers, along with several minor capsid proteins. The latter include pVI, pIIIa, pVIII, and pIX and are important to capsid stability. Each penton capsomer contains a ring of five penton base proteins which bind and support the trimeric fiber protein with its distal fiber knob. During viral infection, the fiber knob binds to one of several host cell receptors (Nemerow, ; Goosney and Nemerow, ; Nemerow et al., ). The penton base protein contains two hypervariable loops. The interaction between fiber knob and a host cell receptor brings about secondary contact between the hypervariable loop 2 (HVL2) arginine-glycine-aspartic acid (RGD) motif in each penton base protein (five per penton base capsomer) with host cell integrins αvβ3, αvβ5, and αvβ1, that in turn induce endocytosis of the virus (Li et al., ,; Li et al., ). HAdV structural proteins can serve multiple functions. For example, the minor capsid structural protein VI (pVI) plays a critical role in at least three distinct aspects of the viral “life” cycle: endosomal escape during cell entry, nuclear assembly during viral replication, and stability of the intact, infectious virus outside the host (Wodrich et al., ; Wiethoff et al., ; Moyer et al., , ). These findings suggest that, as with pVI, other HAdV structural proteins may have multiple functions yet to be elucidated.

Genomics and evolution

The relatively large genome database for HAdV-D (over 50 unique viruses with available whole genome sequences) (Tables 1, 2) has permitted detailed analyses of genome relationships within this clinically important adenovirus species. HAdV-D genomes are highly conserved (>90%). However, whole genome analyses of HAdV-D have revealed specific loci of genetic hypervariability in the hexon, penton base, fiber, and E3 CR1α, β, and γ genes (Figure 2), dictating nonsynonymous amino acid changes in corresponding proteins (Figure 3). GC content confers genome stability and resistance to recombination (Gruss et al., ), and the genomes of HAdV-D have among the highest GC content among HAdV species (~56%). The hypervariable regions in HAdV-D were found to be sharply reduced in GC nucleotide content relative to the rest of the genome (Robinson et al., ). Mutations in HAdV are relatively infrequent, with genome stability now documented in some types across decades (Hofmayer et al., ; Mahadevan et al., ; Seto et al., ; Dehghan et al., ; Robinson et al., ; Alkhalaf et al., ). However, those regions of the genome shown to be hypervariable and relatively low in GC content are the very same also shown to undergo homologous recombination (Robinson et al., , ; Walsh et al., ; Zhou et al., ; Singh et al., ), driving the evolution of new genotypes.

Table 2

HAdV type#NameGenBank accession no.Genome lengthYear published
HAdV-G52P52H52F52/2003/USADQ923122.2342502007
HAdV-D53P37H22F8/2005/DEUFJ169625349092009
HAdV-D54P54H54F8/2000/JPNAB333801349202008
HAdV-B55P14H11F14/2006/CHNFJ643676347552010
HAdV-D56P56H15F9/2008/FRAHM770721350662011
HAdV-C57P1H57F6/2001/RUSHQ003817358182011
HAdV-D58P58H58F29/1996/ARGHQ883276352172011
HAdV-D59P64H25F56/2007/USAJF799911350722012
HAdV-D60P60H20F60/2009/CANHQ007053350502013
HAdV-A61P31H31F31/2004/JPNJF964962337762011
HAdV-D62P62H62F62/1993/GBRJN162671351272014
HAdV-D63P30H30F29/1959/USAJN935766351682012
HAdV-D64P22H19F37/1993/USAEF121005352312012
HAdV-D65P58H10F9/2004/BGDAP012285351722012
HAdV-B66P66H7F3/1987/ARGJN860676350802012
HAdV-D67P67H9F67/2005/BGDAP012302350752013
HAdV-B68P16H3F16/2004/ARGJN86067835538Unpublished
HAdV-D69P53H15F69/1955/SAUJN226748351242013
HAdV-D70P70H70F29/2014/DEUKP641339351862015
HAdV-D71P9H20F71/1987/DEUKF268207351922013
HAdV-D72P72H30F72/1985/DEUKF268335345532013
HAdV-D73P67H45F27/2015/DEUKY618676351902017
HAdV-D74P70H74F51/2015/DEUKY618677351552017
HAdV-D75P75H26F29/2015/DEUKY618678351042017
HAdV-B76P21H21F16/DEUKF633445355862013
HAdV-B77P35H34F7/1985/DEUKF268328346532013
HAdV-B78P11H11F7/2000/ARGKT97044034881Unpublished
HAdV-B79P11H34F11/2015/JPNLC177352347792017
HAdV-D80P19,23H28F22/2014/DEUTBA34909Unpublished
HAdV-D81P65H48F60/2012/JPNAB765926.1351982014
HAdV-D82P56H15F37/2011/JPNLC066535.135122Unpublished
HAdV-D83P83H9F15/2010/PARKX827426.1352072017
HAdV-D84P43H17F84/2011/PANMF416150352572017
HAdV-D85P37H19F8/2015/JPNLC314153352032018
HAdV-D86P9H25F25/1978/SWETBA35147Unpublished
HAdV-D87P9H15F25/1967/USAMF47684135159Unpublished
HAdV-D88P88H15F9/1963/USAMF47684235115Unpublished
HAdV-C89P89H2F2/2015/DEUTBA35998Unpublished
HAdV-D90P33H27F67/2017/BGDTBA34207*Unpublished

Species and molecular types of human adenovirus (HAdV) genotypes 52–90.

#

Name indicates molecular type (P, penton base; H, hexon; F, fiber)/year of isolation/country of isolation.

TBA: GenBank accession number, to be assigned.

*

Metagenomics project, missing the inverted terminal repeat sequences.

Figure 2

) with permission.

Figure 3

) with permission.

Adenoviruses recombine specifically during viral replication (Williams et al., ; Meinschad and Winnacker, ; Munz et al., ), and do so by both homologous and heterologous mechanisms (Young et al., ; Epstein and Young, ; Crawford-Miksza and Schnurr, ). However, the evidence for homologous recombination as the major mechanism driving HAdV-D evolution is unassailable (Robinson et al., ; Singh et al., ). Specifically, recombination occurs in the two penton base hypervariable regions (these code for two hypervariable loops (HVLs) on the penton base protein, separated from one another by ~125 conserved amino acids), seven hexon hypervariable regions (these are closely adjacent in the hexon gene and determine two adjacent HVLs on the hexon protein), fiber (fiber gene and protein are entirely hypervariable), and E3 CR1α, β, and γ (each also entirely hypervariable). For homologous recombination between two HAdVs to occur, at least two virus types with high nucleotide sequence homology at corresponding locations in both genomes must co-infect the same cell, and viral DNA replication should be ongoing. Co-infection by two or more HAdVs has been well documented (Lee et al., ; Echavarria et al., ; Vora et al., ; McCarthy et al., ; Halstead et al., ; Seto et al., ), as has the presence of two HAdV types in archived clinical samples (Singh et al., ).

“Proteotyping” is a novel approach to the study of genome evolution (Obenauer et al., ), and has been applied to characterize recombination among HAdV-D (Robinson et al., ; Singh et al., ). In this method, maximum likelihood trees are used to align amino acid sequences of hypervariable, frequently recombined proteins. Each amino acid is assigned a unique, arbitrary color. Consensus residues are colored white, and gaps in the alignment are colored black. A threshold of <10% sequence divergence is used to distinguish unique proteotypes. An example of proteotyping is shown (Figure 4), comparing an amino acid alignment from E3 14.7K, a highly conserved gene with one distinct proteotype, with the hypervariable E3 CR1α (Singh et al., ), with six distinct proteotypes observed among 38 HAdV-Ds. E3 14.7K is therefore not hypervariable and not recombinant. E3 CR1α is hypervariable and recombinant.

Figure 4

) with permission.

Another way to interpret the analyses for those proteins like E3 CR1α, with more than one proteotype is that those proteotypes containing more than one HAdV type have previously recombined in nature, while those proteotypes with only one HAdV type are those that have not (yet) been shown to recombine in nature. HAdV-D37 and 29 fall within different hexon proteotypes (Figure 5). HAdV-D37 shares a hexon proteotype with HAdV-D13 and 30 (Robinson et al., ), while HAdV-D29 shares a hexon proteotype with HAdV-D15, 56, and 69 (Singh et al., ). These two hexon proteotypes therefore have undergone prior homologous recombination. In contrast, the hexon proteins of HAdV-D10 and 28 are each in a proteotype with only one member; hexon recombination for these two viruses has therefore not yet been documented in nature. In sum, these data show by independent means that homologous recombination within HAdV-D is common, and confirm previously recognized patterns of homologous recombination among HAdV-D (Robinson et al., , , ,; Walsh et al., , ; Singh et al., , ; Zhou et al., ; Gonzalez et al., ).

Figure 5

). For HAdV-D29 and –D15, the recombination event extended through the E3 CR1β ORF gene and then ended.

The local sequence and/or structure of DNA in regions flanking recombinogenic sites is significant for directing cellular recombination machinery to those regions. In bacteria, a signal for recombination between homologous DNA is the crossover hotspot instigator, or Chi nucleotide sequence. This was first discovered in bacteriophage lambda, then in bacterial DNA, and later shown to mediate recombination between them (Stahl, ). The Chi sequence in E. coli (ChiEC) is 5′-GCTGGTGG-3′ (Smith et al., ; Bianco and Kowalczykowski, ), and its presence induces the conversion of the RecBCD enzyme from a helicase to an exonuclease, producing ssDNA that can invade homologous dsDNA during recombination (Taylor et al., ). The RecA enzyme of E. coli is loaded onto unwound ssDNA by RecBCD and promotes ssDNA exchange/recombination with homologous dsDNA (Cox, ; Smith, ). RecA has significant homology to eukaryotic Rad51 and its paralogs (Suwaki et al., ), enzymes that repair dsDNA breaks in human cells, facilitate homologous recombination, and during adenovirus infection, bind to the E2 DNA binding protein (Tookman et al., ). In our study of the region just 5′ to HVL2 on the penton base gene, a recombination hot-spot for HAdV-D (Robinson et al., ), we found Chi-like sequences (ChiAD), e.g., 5′-ACTTCTGA-3′ in the proteotype containing HAdV-D64, and 5′-TCTCCTGA-3′ in the proteotype including HAdV-D37 (Lee et al., ). The putative ChiAD sequences we identified in HAdV-D were found within the GC-rich component of GC/AT transition zones that precede and include HVL2, and were conserved within each proteotype. In vitro, E. coli lysates containing RecA protein increased recombination of two HAdV-D genotypes with the same penton base HVL2 proteotype. RecA was shown by ChIP to bind specifically to ChiAD nucleotide sequence in the same regions, and also colocalize with adenovirus DNA within infected cell nuclei. These data suggest that Chi-like nucleotide sequences adjacent to the junction of conserved and hypervariable gene segments in HAdV-D may be an important signal for homologous recombination, and provide evidence in support of the idea that local bacterial flora might enhance natural recombination through Chi-like nucleotide sequences at HAdV-D recombination hotspots.

Another explanation for homologous recombination between HAdV, not exclusive of a role for ChiAD, is the potential for GC-low (AT-rich) single stranded DNA (ssDNA) to form hairpin loops (Nagy and Bujarski, ; Ohshima et al., ), a physical nonlinearity that would facilitate binding of ssDNA of one HAdV-D type to a homologous segment of ss or dsDNA from a physically adjacent but different HAdV-D type during co-infection of the same cell. Hairpin loops and other alterations in the physical configuration of ssDNA during DNA replication might also contribute to polymerase jumping (Jennings et al., ; Spaan et al., ; Pääbo et al., ; Viswanathan et al., ), in which physical constraints to polymerization lead to translocation of the DNA polymerase to an adjacent DNA from a different virus, resulting in a recombined DNA. Polymerase jumping has been shown to occur during HAdV DNA replication (King et al., ; de Jong et al., ), although it has not been suggested previously as a mechanism for HAdV-D evolution. Analysis of 38 HAdV-D whole genome sequences identified instances of 15 nucleotide-long GC-rich sequence adjacent to 15 nucleotide-long AT-rich sequence (sometimes with a 15 or 30 nucleotide-long GC-moderate sequence intervening), located just 5′ and 3′ to frequently recombined gene segments, and which were shown by in silico analysis of their corresponding ssDNA to form hairpin loops (Robinson et al., ). Taken together, these data suggest covariant effects of nucleotide sequence and ssDNA secondary structures on homologous recombination between two HAdV-Ds.

Transcriptome

Regions of the HAdV-D genome currently thought to be “noncoding,” may contain functional elements. Because viruses exist on the nano-scale, viral genomes are by necessity constrained by size, and “junk” nucleotide sequences represent an extravagance. The National Human Genome Research Institute project to identify functional elements in the human genome (Encyclopedia of DNA Elements, or ENCODE) identified functionality in much of the human genome previously without known utility (Consortium et al., ; Qu and Fang, ; Kellis et al., ) The double-stranded DNA genomes of HAdV also contain regions with no known function. Transcriptional profiling of host gene expression has been studied after HAdV infection (Dorer et al., ) However, although viral transcriptomes have been reported for several viruses, most notably dengue, varicella zoster, and Epstein-Barr viruses (Ortmann et al., ; Ertl et al., ; Nagel et al., , ; Arvey et al., ; Sujayanont et al., ), a de novo HAdV transcriptome has not been reported. Wu and coworkers used deep RNA sequencing to confirm known bat AdV transcripts (Wu et al., ), but did not investigate “noncoding” regions. In silico ORF prediction in HAdV can be difficult due to splice variants and inconsistencies in banked gene annotations (Davison et al., ), but in a prior annotation of HAdV-D37, ~60 new additional ORFs were predicted using in combination, the NCBI ORF finder, TIGR annotation engine, and GeneMark Heuristic model (Robinson et al., ) Putative genes were found within the large regions of noncoding DNA on the complementary strand opposite to established HAdV genes (Figure 6), in smaller regions on the coding strand within established transcription units but between confirmed genes, and overlapping or completely within established genes. Work is in progress in our laboratories to identify putative new genomic elements in HAdV by high-throughput sequencing of the viral transcriptome of HAdV-D37.

Figure 6

) with permission.

Transcription factor binding sites

HAdV uses host TFs nuclear factor I and III (NF-I and NF-III) as part of the viral DNA replication complex (Pruijn et al., ; Mul et al., ; Hatfield and Hearing, , ) Simian AdVs typically lack the NF-I binding site, while human viruses express it. It was previously reported that HAdV-E4, originally isolated in 1952, is a product of recombination between HAdV-B16 and the simian AdV, SAdV-E26. Clinical strains of HAdV-E4 isolated recently contain a NF-I binding site in the inverted terminal repeat (Houng et al., ; Dehghan et al., ,) that is absent in the original isolate (Purkayastha et al., ), suggesting that NF-I binding may be important to viral fitness in humans. To further elucidate mechanisms of viral gene expression, we are exploring novel TF binding sites on HAdV-D DNA, using ENCODE validated methodologies (Gerstein et al., ; Landt et al., ).

Protein interactomes

The major HAdV capsid structural proteins—hexon, penton base, and fiber—interact directly with extracellular mediators of host immunity (Gahéry-Ségard et al., ; Molinier-Frenkel et al., ; Schoggins and Falck-Pedersen, ; Tamanini et al., ; Kalyuzhniy et al., ; Chintakuntlawar et al., ; Bradshaw et al., ; Flatt et al., ). The hexon, penton base, and fiber proteins also exhibit distinct amino acid signatures, characterizing discrete proteotypes (Robinson et al., ). Gene products from the E3 transcription unit of HAdV function in viral immune evasion (Horwitz, ; Lichtenstein et al., ; Windheim et al., ). In HAdV-D, the open reading frames for three of eight E3 genes—CR1α, CR1β, and CR1γ–are uniquely hypervariable compared to the other ORFs within the E3 transcription unit, and also segregate into discrete proteotypes (Singh et al., ). Highly conserved genes, such as DNA binding protein, DNA polymerase, and E3 14.7K, show no such variability (Robinson et al., ; Singh et al., ).

While it may be assumed that hypervariablity in major capsid and E3 proteins is driven through evolutionary selection by the extracellular interactome, amino acid differences in a hypervariable protein can also lead to differences in that protein's intracellular interactome, the set of intrinsic host cell proteins which network with the viral protein, as was recently confirmed for E3 CR1 genes across HAdV species (Martinez-Martin et al., ). Viral capsid structural proteins are critical to virion stability. For the nonenveloped HAdV, fiber and penton base proteins on the external surface of the capsid serve as ligands for attachment to the host cell (Huang et al., ) and initiate viral entry (Wickham et al., ), respectively. HAdVs are typically internalized via endosomes. Endosomal acidification leads to structural instability of the capsid and endosomal release into the cytosol. HAdV capsid is then transported by microtubules to the nuclear membrane. Viral DNA then enters the nucleus through nuclear pores, leaving almost all the viral structural proteins in the cytosol (Henaff et al., ). Viral capsid proteins within the cell are eventually targeted for ubiquitination (Ko et al., ; Marvin and Wiethoff, ; Horan et al., ; Li et al., ) and degraded (Greber et al., ), but there are many opportunities for interaction with intracellular host cellular proteins during entry, trafficking, translation, assembly, and egress. Penton base HVL2, with its RGD motif, is critical to viral internalization through the interaction with host cell integrins,; (Wickham et al., ) but function of penton base HVL1 is unknown, and might be revealed though knowledge of its protein interactome. The closely adjacent hexon HVL1 and 2 form the epsilon epitope that determines serum neutralization, and interactions between the hexon protein and serum coagulation factor X confers liver tropism to HAdV-C5 (Sumarheni et al., ). However, nothing is known about potential hexon interactions with intracellular proteins during infection.

The E3 transcription unit of HAdV codes for proteins that mediate immune evasion by the virus (Horwitz, ). Although E3 is labeled as an early transcription region, its transcripts are expressed both early and late during viral infection (Chow et al., ; Chow and Broker, ; Bhat and Wold, ), and there is evidence for at least one E3 protein that late transcripts are translated (Robinson et al., ). E3 gene products are not required for viral replication in cultured cells (Morin et al., ), but inhibit cellular and cytokine mediated host immune responses to infection (Horwitz, ; Lichtenstein et al., ; Windheim et al., ). Almost all of what is known about the function of specific E3 proteins derives from studies on HAdV-C. For example, HAdV-C2 E3 CR1α directs another E3 protein (19K) to the endoplasmic reticulum of cytotoxic T cells (Wilson-Rawls et al., ), where 19K binds to and retains MHC class I proteins (Jefferies and Burgert, ), preventing presentation of viral peptides within MHC class I at the cell surface (Burgert and Kvist, , ; Andersson et al., ; Burgert et al., ; Cox et al., ). CR1α, RIDα, and RIDβ proteins cooperate to evade TNFα-related apoptosis through TRAIL (Elsing and Burgert, ; Tollefson et al., ; Benedict et al., ; Lichtenstein et al., ). CR1β (Wold et al., ), also called the adenovirus death protein (Tollefson et al., ), is required for cell lysis (Tollefson et al., ) and viral spread (Doronin et al., ). The ORF size of each E3 gene varies across HAdV species (Figure 7) (Robinson et al., ). Similarly, immune evasion functions of E3 gene products may not be the same across HAdV species, or function similarly in all cell types (Routes and Cook, ). Windheim and coworkers recently showed that the CR1β protein of the eye pathogen HAdV-D64 suppresses natural killer cell function (Windheim et al., ). The E3 CR1 genes are uniquely hypervariable within HAdV-D, and as predicted, overlapping but distinguishable intracellular interactomes across proteotypes were recently reported by Martinez-Martin and colleagues, who used protein microarrays to identify novel CR1β binding partners (Martinez-Martin et al., ).

Figure 7

) with permission.

Conclusions

HAdV was critical to the dual discoveries of viral oncogenesis and RNA splicing (Berget et al., ; Chow et al., ; Whyte et al., ). HAdV is also a significant agent of disease for which there is no approved treatment. Recent mining of HAdV genomes has been highly productive, and there is ample evidence to suggest that further whole genome analyses will elucidate new and fundamental mechanisms in HAdV biology. In the last decade, of 27 newly identified HAdVs, 19 were HAdV-Ds, suggesting the continuing evolution of new pathogens from species D. Analyses of 38 fully sequenced HAdV-D whole genomes identified homologous recombination of specific regions within the hexon, penton base, fiber, and E3 CR1 genes as the major mechanism behind HAdV-D evolution, a new finding (Robinson et al., ; Singh et al., ). Stereotypical reductions in GC content at the junction of conserved and hypervariable regions, along with Chi-like sequence motifs (also a new finding), appear likely to augment the intrinsic tendency of HAdV to undergo homologous recombination in vivo (Lee et al., ).

Recently, the whole genome sequences of 85 HAdVs from archives and current collections were determined, including both historical and circulating strains, respectively (Ismail et al., ). Of these, 3 novel recombinants within HAdV-B and 15 within HAdV-D were identified. Only two of the 15 HAdV-Ds were found to contain novel genes (penton base and fiber); these were subsequently typed as HAdV-D71 and 72. Isolates of HAdV-D53 and HAdV-D58, two novel genotypes recently recognized, were also identified, adding confidence in their clinical importance. Fully genotyped HAdVs now number 90, with more awaiting type numbers, and the scientific community has a 10-fold larger database of unique HAdV genomes than available only 15 years ago. Published and validated ENCODE methodologies can now be applied, and comparisons made across disparate HAdV genomes. We suggest that the HAdV genome contains previously uncharacterized functional elements, and that every HAdV protein has pleiotropic interactions. Current technologies should afford a wave of new and important discoveries that may lead to needed therapies against adenoviral diseases.

Statements

Author contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

This work was funded by National Institutes of Health grants EY013124, EY021558, and EY014104, a Senior Scientific Investigator Award grant (to JC) from Research to Prevent Blindness, Inc., New York, NY, The Falk Foundation, and the Massachusetts Lions Eye Research Fund.

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

Summary

Keywords

adenovirus, genome, evolution, transcription factor, interactome

Citation

Ismail AM, Lee JS, Lee JY, Singh G, Dyer DW, Seto D, Chodosh J and Rajaiya J (2018) Adenoviromics: Mining the Human Adenovirus Species D Genome. Front. Microbiol. 9:2178. doi: 10.3389/fmicb.2018.02178

Received

20 July 2018

Accepted

24 August 2018

Published

11 September 2018

Volume

9 - 2018

Edited by

Qiwei Zhang, Southern Medical University, China

Reviewed by

Jun Hang, Walter Reed Army Institute of Research, United States; Gabriel Gonzalez, Hokkaido University, Japan

Updates

Copyright

*Correspondence: Jaya Rajaiya

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