Skip to main content

Genome sequences of human cytomegalovirus strain TB40/E variants propagated in fibroblasts and epithelial cells

Abstract

The advent of whole genome sequencing has revealed that common laboratory strains of human cytomegalovirus (HCMV) have major genetic deficiencies resulting from serial passage in fibroblasts. In particular, tropism for epithelial and endothelial cells is lost due to mutations disrupting genes UL128, UL130, or UL131A, which encode subunits of a virion-associated pentameric complex (PC) important for viral entry into these cells but not for entry into fibroblasts. The endothelial cell-adapted strain TB40/E has a relatively intact genome and has emerged as a laboratory strain that closely resembles wild-type virus. However, several heterogeneous TB40/E stocks and cloned variants exist that display a range of sequence and tropism properties. Here, we report the use of PacBio sequencing to elucidate the genetic changes that occurred, both at the consensus level and within subpopulations, upon passaging a TB40/E stock on ARPE-19 epithelial cells. The long-read data also facilitated examination of the linkage between mutations. Consistent with inefficient ARPE-19 cell entry, at least 83% of viral genomes present before adaptation contained changes impacting PC subunits. In contrast, and consistent with the importance of the PC for entry into endothelial and epithelial cells, genomes after adaptation lacked these or additional mutations impacting PC subunits. The sequence data also revealed six single noncoding substitutions in the inverted repeat regions, single nonsynonymous substitutions in genes UL26, UL69, US28, and UL122, and a frameshift truncating gene UL141. Among the changes affecting protein-coding regions, only the one in UL122 was strongly selected. This change, resulting in a D390H substitution in the encoded protein IE2, has been previously implicated in rendering another viral protein, UL84, essential for viral replication in fibroblasts. This finding suggests that IE2, and perhaps its interactions with UL84, have important functions unique to HCMV replication in epithelial cells.

Passage of human cytomegalovirus (HCMV; species Human betaherpesvirus 5) in cell culture results in mutations ranging from small substitutions, insertions, or deletions to large-scale deletions, duplications, and rearrangements. Some mutations appear to be stochastic, whereas others consistently disrupt specific genes and confer growth advantages in certain cell types or growth conditions [1,2,3]. Of particular importance are mutations that invariably arise during passage in fibroblasts and alter one or more of three contiguous genes, UL128, UL130, and UL131A [1]. These mutations prevent assembly of a pentameric complex (PC) consisting of the UL128, UL130, and UL131A proteins complexed with glycoproteins H and L on the virion surface that is important for entry into epithelial and endothelial cells but dispensable for entry into fibroblasts [4,5,6,7,8].

Many commonly used laboratory strains have large deletions as well as various additional mutations that disrupt PC expression, thereby rendering the virus non-epitheliotropic and non-endotheliotropic [9]. Consequently, research has shifted toward the use of genetically more authentic strains such as the endothelial cell-adapted strain TB40/E [10]. However, the designation “TB40/E” has been applied indiscriminately to heterogeneous stocks propagated from the original TB40/E stock [10], as well as to a variety of viruses derived from bacterial artificial chromosome (BAC) clones generated from such stocks [11, 12]. The sequences and tropisms of these viruses can differ significantly from each other and from the original TB40/E strain. For example, clone Lisa, a virus that was plaque-purified from a TB40/E stock [13], has a 1 bp insertion in UL128 causing a frameshift after codon 69 [13], whereas the widely utilized BAC clone TB40-BAC4 carries a single nucleotide substitution in the intron between UL128 exons 2 and 3 that reduces splicing efficiency, lowers levels of the encoded protein, and reduces infection efficiency in epithelial cells [3]. In contrast, a more recently constructed BAC clone, TB40-KL7-SE, has no obvious mutations impacting PC expression and is both endotheliotropic and epitheliotropic [12].

To begin addressing the diversity of TB40/E stocks and the impact of propagation using different cell types, a TB40/E stock amplified twice in primary human foreskin fibroblasts (HFF; TB40/EF, stock 31,519) and entering retinal pigmented epithelial cells (ARPE-19; ATCC® CRL-2302) with poor efficiency [14] was passaged five times in ARPE-19 cells, generating a stock (TB40/EE, stock SE2) capable of infecting HFF and ARPE-19 cells with similar efficiencies [14]. Despite efficient entry, the amounts of cell-free virus generated by ARPE-19 cells infected with TB40/EE were consistently 100- to 1,000-fold lower than those produced by HFF cells, revealing the existence of epithelial cell-specific post-entry restrictions. TB40/EE also exhibited an increased propensity to form multinucleated syncytia in ARPE-19 cell populations, suggesting an enhancement in the ability to induce cell–cell fusion during infection [14].

In the current study, we used long-read PacBio sequencing to examine in detail the genetic changes potentially associated with adaptation to epithelial cells. DNA was isolated from TB40/EF and TB40/EE cell-free virions as described previously [15]. HiFi SMRTbell library construction and sequencing were performed at the Genomics Core at Virginia Commonwealth University as described previously [16]. The data were processed using tools from the PacBio SMRT-Link command-line package (https://www.pacb.com/support/software-downloads/) with default settings. Two-modal polymerase reads for TB40/EF (25,124) or TB40/EE (8,364) were indexed using pbindex, XML files of the subread counts were produced using dataset, and 16,920 HiFi reads were generated for TB40/EF and 6,985 HiFi reads for TB40/EE using CCS. The approximately three-fold difference in TB40/EF compared to TB40/EE reads is consistent with a three-fold difference in DNA concentration, which in turn reflects reduced levels of cell-free virus released from TB40/EE-infected ARPE-19 cells [14]. Final HCMV genome assemblies were made by reference-guided de novo assembly using LoReTTA v0.1 (https:/bioinformatics.cvr.ac.uk/software/; [16] with default settings and the sequence of strain TB40/E clone Lisa (13; GenBank accession no. KF297339.1) as the reference. HiFi reads were then mapped to the respective final assemblies using minimap v2.17-r941 [17], and the read alignments were visualized using the Integrative Genomics Viewer [18]. The consensus genome sequences were deposited in GenBank under accession numbers MW439038 (TB40/EF) and MW439039 (TB40/EE), and had median coverage depths of 202 and 43 reads/nucleotide, respectively. Differences between these sequences were identified, and these and other major heterogeneities noted during examination of the read alignments were quantified by counting their occurrence in the reads.

The HCMV genome (236 kbp) consists of unique long (UL) and unique short (US) regions, each of which is flanked by inverted repeats in the arrangement ab-UL-bʹaʹcʹ-US-ca (the primes denote the inverted repeats of a, b, and c). Comparison of the TB40/EF and TB40/EE consensus genome sequences identified 12 single nucleotide substitutions and one single nucleotide deletion within noncoding regions in the a, b, or c inverted repeats (Table 1). As seven of these were replicated in the inverted repeats (one each in a/aʹ/a, b/bʹ, and c/cʹ), these loci represent only six unique differences. Although these mutations are in noncoding regions, the significant levels of enrichment suggest that they may provide a selective advantage in ARPE-19 cells. However, the inverted repeats have been reported to be particularly prone to mutation during passage of HCMV, with changes being generally replicated in all copies presumably as the result of recombination [1].

Table 1 Nucleotide differences between TB40/EF and TB40/EE within non-coding regions

Comparison of the TB40/EF and TB40/EE consensus genome sequences and examination of the read alignments also identified changes in coding regions (Table 2). Given the low efficiency of epithelial cell entry observed for TB40/EF, mutations disrupting UL128, UL130, or UL131A (encoding the PC subunit proteins UL128, UL130, and UL131A, respectively) were anticipated. Indeed, targeted sequencing of TB40/EF had previously identified a suppressor substitution converting the UL128 stop codon (TGA) to TTA (encoding leucine), thereby extending UL128 by 19 residues [14]. Although the consensus genome sequence of TB40/EF did not reflect this mutation, examination of read frequencies revealed the existence of a subpopulation in which 30% of TB40/EF genomes contain this suppressor mutation. Further examination identified two additional subpopulations: one containing a 2 bp deletion causing a frameshift in UL128 and resulting in truncation of UL128 (44% of genomes), and one containing a single nucleotide substitution in UL130 resulting in a C207S substitution in UL130 (9% of genomes) (Table 2). There was no evidence for subpopulations with mutations in UL131A. The long length of PacBio reads connected not only the two loci in UL128, which are separated by107 bp, but also the UL130 locus 814 bp beyond; among the connecting reads, those containing one mutation did not contain the others. Thus, consistent with the low epithelial cell entry efficiency of the TB40/EF stock [14], these findings suggest that cumulatively 83% of TB40/EF genomes contain mutations potentially impacting PC assembly or function. In contrast, and consistent with efficient epithelial cell entry [14], mutations impacting PC subunits were absent from TB40/EE.

Table 2 Nucleotide differences between TB40/EF and TB40/EE within coding regions

Although it has not been demonstrated that the UL128 suppressor substitution disrupts the function of the PC, indirect evidence indicates that the UL130 C207S substitution is likely to have a negative effect. The crystal structure of the PC indicates the presence of a disulfide bond between C207 and C172 [19], and the converse mutation, C172W, has been reported to occur during serial fibroblast passage of HCMV strain IgKG-H2 in conjunction with loss of epithelial cell tropism [16]. Moreover, in HCMV strain Towne, a frameshift in UL130 after codon 203 replaces 11 C-terminal residues (including C207) with 26 novel residues, resulting in rapid degradation of the mutant protein and loss of endothelial cell tropism [20]. These findings suggest that the C172-C207 disulfide bond is critical for the function or stability of UL130, and for its essential role in PC formation and epithelial cell entry. Curiously, the three mutations impacting PC subunit genes in TB40/EF are different from the mutations in clone Lisa or TB40-BAC4, and were not detected in reads from TB40/EE. Thus, within the available TB40/E lineages, at least five distinct mutations targeting PC subunit genes have been identified thus far.

Five other changes impacting coding regions were also identified (Table 2). These included single nonsynonymous substitutions in genes UL26 (resulting in E98K in UL26), UL69 (H492Q in UL69), UL122 (D390H in IE2), and US28 (C320W in US28), and a two nucleotide insertion in gene UL141 introducing a frameshift truncating UL141. Examination of read frequencies revealed that most of these changes were enriched to a marginal or modest level: UL26 (from 53 to 100%), UL69 (from 44 to 52%), UL141 (from 44 to 69%) and US28 (from 88 to 100%). Thus, although these changes may be associated with improved replication in ARPE-19 cells, they may have been the consequence of stochastic effects. Moreover, both variants of UL69 and UL141 have been reported previously in consensus sequences of strain TB40/E, namely a partial TB40/E sequence (GenBank accession number AY446866.1), clone Lisa (KF297339.1), and the BAC clones TB40-BAC4 and TB40-KL7-SE (EF999921.1 and MF871618.1, respectively) [9, 11,12,13]. In UL69, clone Lisa and TB40-BAC4 encode Q492, whereas TB40-KL7-SE encodes H492. In UL141, the frameshift is absent from clone Lisa but present in the partial TB40/E sequence, TB40-BAC4, and TB40-KL7-SE. Thus, it appears that parental TB40/E stocks contained two variants of both genes, with capture of one or the other allele in the genomes of clone Lisa, TB40-BAC4, and TB40-KL7-SE resulting from cloning. In contrast, the prevalence of the D390H substitution in IE2 increased markedly from 5 to 87% (Table 2), suggesting strong selective pressure favoring this allele during ARPE-19 adaptation. The D390 allele is unique to strain TB40/E and is present in all currently reported TB40/E-derived sequences, whereas the H390 allele is conserved among all other HCMV strains for which sequences are publicly available. This observation is all the more interesting given that gene UL84 has been identified as being essential for replication in vitro in fibroblasts in the presence of the IE2 H390 allele, but non-essential in the presence of the D390 allele [21].

In summary, whole genome sequencing identified variants impacting IE2 and PC subunits UL128 and UL130 as being potentially selected during adaptation of HCMV strain TB40/E for growth in epithelial cells. Enrichment of viral genomes lacking disruptive mutations in UL128 and UL130 is consistent with the detected improvement in efficiency of epithelial cell entry [14], and, as the PC has been associated with increased cell–cell fusion [22,23,24,25], may also explain the reported increase in syncytium formation [14]. It is not known how the D390H polymorphism in IE2 determines the requirement for UL84 during fibroblast replication, or whether this phenomenon also extends to epithelial cells, but selection of genomes encoding the H390 allele suggests that IE2 and, perhaps, its interplay with UL84 provide important functions that are unique to HCMV replication in epithelial cells. Construction and phenotypic characterization of viral mutants containing these genetic changes in isolation are in progress to further elucidate the role of UL84 in the context of IE2 H390 or IE2 D390.

Availability of data and materials

Genome sequences are available from GenBank under accession numbers MW439038 (TB40/EF) and MW439039 (TB40/EE).

Abbreviations

HCMV:

Human cytomegalovirus

PC:

Pentameric complex

BAC:

Bacterial artificial chromosome

HFF:

Human foreskin fibroblast

UL :

Unique long

US :

Unique short

References

  1. 1.

    Dargan DJ, Douglas E, Cunningham C, Jamieson F, Stanton RJ, Baluchova K, McSharry BP, Tomasec P, Emery VC, Percivalle E, et al. Sequential mutations associated with adaptation of human cytomegalovirus to growth in cell culture. J Gen Virol. 2010;91:1535–46.

    CAS  Article  Google Scholar 

  2. 2.

    Murrell I, Wilkie GS, Davison AJ, Statkute E, Fielding CA, Tomasec P, Wilkinson GW, Stanton RJ. Genetic stability of bacterial artificial chromosome-derived human cytomegalovirus during culture in vitro. J Virol. 2016;90:3929–43.

    CAS  Article  Google Scholar 

  3. 3.

    Murrell I, Tomasec P, Wilkie GS, Dargan DJ, Davison AJ, Stanton RJ. Impact of sequence variation in the UL128 locus on production of human cytomegalovirus in fibroblast and epithelial cells. J Virol. 2013;87:10489–500.

    CAS  Article  Google Scholar 

  4. 4.

    Hahn G, Revello MG, Patrone M, Percivalle E, Campanini G, Sarasini A, Wagner M, Gallina A, Milanesi G, Koszinowski U, et al. Human cytomegalovirus UL131-128 genes are indispensable for virus growth in endothelial cells and virus transfer to leukocytes. J Virol. 2004;78:10023–33.

    CAS  Article  Google Scholar 

  5. 5.

    Wang D, Shenk T. Human cytomegalovirus virion protein complex required for epithelial and endothelial cell tropism. Proc Natl Acad Sci U S A. 2005;102:18153–8.

    CAS  Article  Google Scholar 

  6. 6.

    Adler B, Scrivano L, Ruzcics Z, Rupp B, Sinzger C, Koszinowski U. Role of human cytomegalovirus UL131A in cell type-specific virus entry and release. J Gen Virol. 2006;87:2451–60.

    CAS  Article  Google Scholar 

  7. 7.

    Ryckman BJ, Rainish BL, Chase MC, Borton JA, Nelson JA, Jarvis MA, Johnson DC. Characterization of the human cytomegalovirus gH/gL/UL128-131 complex that mediates entry into epithelial and endothelial cells. J Virol. 2008;82:60–70.

    CAS  Article  Google Scholar 

  8. 8.

    Freed DC, Tang Q, Tang A, Li F, He X, Huang Z, Meng W, Xia L, Finnefrock AC, Durr E, et al. Pentameric complex of viral glycoprotein H is the primary target for potent neutralization by a human cytomegalovirus vaccine. Proc Natl Acad Sci USA. 2013;110:E4997-5005.

    CAS  Article  Google Scholar 

  9. 9.

    Dolan A, Cunningham C, Hector RD, Hassan-Walker AF, Lee L, Addison C, Dargan DJ, McGeoch DJ, Gatherer D, Emery VC, et al. Genetic content of wild-type human cytomegalovirus. J Gen Virol. 2004;85:1301–12.

    CAS  Article  Google Scholar 

  10. 10.

    Sinzger C, Schmidt K, Knapp J, Kahl M, Beck R, Waldman J, Hebart H, Einsele H, Jahn G. Modification of human cytomegalovirus tropism through propagation in vitro is associated with changes in the viral genome. J Gen Virol. 1999;80:2867–77.

    CAS  Article  Google Scholar 

  11. 11.

    Sinzger C, Hahn G, Digel M, Katona R, Sampaio KL, Messerle M, Hengel H, Koszinowski U, Brune W, Adler B. Cloning and sequencing of a highly productive, endotheliotropic virus strain derived from human cytomegalovirus TB40/E. J Gen Virol. 2008;89:359–68.

    CAS  Article  Google Scholar 

  12. 12.

    Sampaio KL, Weyell A, Subramanian N, Wu Z, Sinzger C. A TB40/E-derived human cytomegalovirus genome with an intact US-gene region and a self-excisable BAC cassette for immunological research. Biotechniques. 2017;63:205–14.

    CAS  Article  Google Scholar 

  13. 13.

    Tomasec P, Wang EC, Davison AJ, Vojtesek B, Armstrong M, Griffin C, McSharry BP, Morris RJ, Llewellyn-Lacey S, Rickards C, et al. Downregulation of natural killer cell-activating ligand CD155 by human cytomegalovirus UL141. Nat Immunol. 2005;6:181–8.

    CAS  Article  Google Scholar 

  14. 14.

    Vo M, Aguiar A, McVoy MA, Hertel L. Cytomegalovirus Strain TB40/E Restrictions and Adaptations to Growth in ARPE-19 Epithelial Cells. Microorganisms 2020, 8.

  15. 15.

    Ourahmane A, Cui X, He L, Catron M, Dittmer DP, Al Qaffasaa A, Schleiss MR, Hertel L, McVoy MA. Inclusion of Antibodies to Cell Culture Media Preserves the Integrity of Genes Encoding RL13 and the Pentameric Complex Components During Fibroblast Passage of Human Cytomegalovirus. Viruses 2019, 11.

  16. 16.

    Qaffas AA, Nichols J, Davison AJ, Ourahmane A, Hertel L, McVoy MA, Camiolo S. LoReTTA, a user-friendly tool for assembling viral genomes from PacBio sequence data. Virus Evolution 2021.

  17. 17.

    Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34:3094–100.

    CAS  Article  Google Scholar 

  18. 18.

    Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.

    CAS  Article  Google Scholar 

  19. 19.

    Chandramouli S, Malito E, Nguyen T, Luisi K, Donnarumma D, Xing Y, Norais N, Yu D, Carfi A. Structural basis for potent antibody-mediated neutralization of human cytomegalovirus. Sci Immunol 2017, 2.

  20. 20.

    Patrone M, Secchi M, Fiorina L, Ierardi M, Milanesi G, Gallina A. Human cytomegalovirus UL130 protein promotes endothelial cell infection through a producer cell modification of the virion. J Virol. 2005;79:8361–73.

    CAS  Article  Google Scholar 

  21. 21.

    Spector DJ. UL84-independent replication of human cytomegalovirus strains conferred by a single codon change in UL122. Virology. 2015;476:345–54.

    CAS  Article  Google Scholar 

  22. 22.

    Cui X, Freed DC, Wang D, Qiu P, Li F, Fu TM, Kauvar LM, McVoy MA: Impact of Antibodies and Strain Polymorphisms on Cytomegalovirus Entry and Spread in Fibroblasts and Epithelial Cells. J Virol 2017, 91.

  23. 23.

    Chiuppesi F, Wussow F, Johnson E, Bian C, Zhuo M, Rajakumar A, Barry PA, Britt WJ, Chakraborty R, Diamond DJ. Vaccine-derived neutralizing antibodies to the human cytomegalovirus gH/gL pentamer potently block primary cytotrophoblast infection. J Virol. 2015;89:11884–98.

    CAS  Article  Google Scholar 

  24. 24.

    Ciferri C, Chandramouli S, Donnarumma D, Nikitin PA, Cianfrocco MA, Gerrein R, Feire AL, Barnett SW, Lilja AE, Rappuoli R, et al. Structural and biochemical studies of HCMV gH/gL/gO and Pentamer reveal mutually exclusive cell entry complexes. Proc Natl Acad Sci USA. 2015;112:1767–72.

    CAS  Article  Google Scholar 

  25. 25.

    Gerna G, Percivalle E, Perez L, Lanzavecchia A, Lilleri D. Monoclonal antibodies to different components of the human cytomegalovirus (HCMV) pentamer gH/gL/pUL128L and trimer gH/gL/gO as well as antibodies elicited during primary HCMV infection prevent epithelial cell syncytium formation. J Virol. 2016;90:6216–23.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

The authors thank Christian Sinzger of Ulm University Medical Center, Ulm, Germany for providing a seed stock of TB40/E, and Gregory Buck and Vladimir Lee of the Genomics Core at Virginia Commonwealth University for conducting the PacBio sequencing.

Funding

Financial support was provided by a PacBio SMRT Grant (MAM), National Institutes of Health Grant R01AI128912 (MAM and LH), and Wellcome Grant 204870/Z/16/Z (AJD). Funding bodies had no role in study design, data collection, analysis, interpretation, or writing of the manuscript.

Author information

Affiliations

Authors

Contributions

LH and MAM conceived the study, LH, MAM, MS, and AJD supervised the work, MV, AA, and AO prepared and provided the materials, AAQ and SC assembled, annotated, and analyzed the sequence data, LH and MAM wrote the draft manuscript, and all authors contributed to and approved the final manuscript.

Corresponding authors

Correspondence to Michael A. McVoy or Laura Hertel.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Al Qaffas, A., Camiolo, S., Vo, M. et al. Genome sequences of human cytomegalovirus strain TB40/E variants propagated in fibroblasts and epithelial cells. Virol J 18, 112 (2021). https://doi.org/10.1186/s12985-021-01583-3

Download citation

Keywords

  • Human cytomegalovirus
  • Epithelial cell adaptation
  • PacBio
  • Whole genome sequence
  • TB40/E
  • UL122