Original Research ARTICLE
Whole-Genome Approach to Assessing Human Cytomegalovirus Dynamics in Transplant Patients Undergoing Antiviral Therapy
- 1MRC-University of Glasgow Centre for Virus Research, Glasgow, United Kingdom
- 2Westmead Institute for Medical Research, Sydney, NSW, Australia
- 3Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, Australia
- 4Blood and Bone Marrow Transplant Unit, Westmead Hospital, Sydney, NSW, Australia
- 5Sydney Cellular Therapies Laboratory, Westmead Hospital, Sydney, NSW, Australia
- 6Institute for Medical Virology, University Hospital Tübingen, Tübingen, Germany
- 7Institute of Virology, Hannover Medical School, Hanover, Germany
- 8Department of Nephrology, Hannover Medical School, Hanover, Germany
- 9Public Health Agency of Lower Saxony, Hanover, Germany
- 10German Center for Infection Research, Hanover, Germany
- 11West of Scotland Specialist Virology Centre, Glasgow Royal Infirmary, Glasgow, United Kingdom
- 12Discipline of Infectious Diseases and Immunology, Faculty of Medicine and Health, Charles Perkins Centre, The University of Sydney, Sydney, NSW, Australia
Human cytomegalovirus (HCMV) is the most frequent cause of opportunistic viral infection following transplantation. Viral factors of potential clinical importance include the selection of mutants resistant to antiviral drugs and the occurrence of infections involving multiple HCMV strains. These factors are typically addressed by analyzing relevant HCMV genes by PCR and Sanger sequencing, which involves independent assays of limited sensitivity. To assess the dynamics of viral populations with high sensitivity, we applied high-throughput sequencing coupled with HCMV-adapted target enrichment to samples collected longitudinally from 11 transplant recipients (solid organ, n = 9, and allogeneic hematopoietic stem cell, n = 2). Only the latter presented multiple-strain infections. Four cases presented resistance mutations (n = 6), two (A594V and L595S) at high (100%) and four (V715M, V781I, A809V, and T838A) at low (<25%) frequency. One allogeneic hematopoietic stem cell transplant recipient presented up to four resistance mutations, each at low frequency. The use of high-throughput sequencing to monitor mutations and strain composition in people at risk of HCMV disease is of potential value in helping clinicians implement the most appropriate therapy.
Despite continuous advances in diagnostics and therapy, human cytomegalovirus (HCMV) remains the most frequent opportunistic viral infection following transplantation, contributing significantly to patient morbidity, and mortality (Boeckh and Ljungman, 2009). The risk of developing HCMV disease after transplantation depends on a number of host factors, including the HCMV serological status of the donor and the recipient, the type of transplant and the level of immunosuppression. In addition, viral factors are associated with poorer outcomes, including the selection of mutants resistant to antiviral drugs and the occurrence of infections involving multiple HCMV strains (Wu et al., 2010; Lisboa et al., 2012).
Viral factors are usually addressed by analyzing individual HCMV genes by polymerase chain reaction (PCR) and Sanger sequencing. Resistance mutations can be screened in the genes responsible for antiviral drug activity (typically UL54 encoding the viral DNA polymerase, and UL97 encoding a phosphotransferase), whereas multiple-strain infections may be detected by genotyping hypervariable genes (commonly UL73 encoding glycoprotein N [gN], UL74 [gO], and UL55 [gB]). This strategy involves an independent assay for each gene analyzed, and its sensitivity is generally limited to the detection of subpopulations exceeding 20% of the total viral population (Schuurman et al., 1999; Sahoo et al., 2013).
Monitoring the complexity of viral populations and the dynamics of resistance mutations is critical to understanding the evolutionary mechanisms operating on the virus during antiviral selection. This is also clinically relevant because the presence of cells infected with multiple strains provides an opportunity for the virus to recombine (Rasmussen et al., 2003), which might lead to an increase of viral strain pathogenicity. Moreover, the presence of drug-resistant populations at low frequencies (<5%) during the early stages of antiviral treatment may signal eventual treatment failure when these populations come to predominate at later stages (Chou et al., 2014; Houldcroft et al., 2016). Indeed, continuation with the antiviral drug to which resistance has developed can lead to the development of additional resistance mutations (Lurain and Chou, 2010). In contrast, therapy modification can accelerate viral clearance (Castagnola et al., 2004). Thus, rapid, comprehensive and sensitive molecular diagnostic approaches would be beneficial for assisting clinicians in improving the management of HCMV infections.
The application of commercial target enrichment methods coupled with the use of high-throughput sequencing (HTS) has dramatically surpassed the sensitivity of PCR methods for characterizing HCMV in clinical material (Hage et al., 2017; Cudini et al., 2019; Suárez et al., 2019a,b). This approach has the great advantage of allowing resistant populations and strain numbers to be assessed across the whole viral genome in a single assay. Studies employing analysis of individual HCMV genes by PCR and HTS to detect antiviral resistance mutations in HCMV infections have been reported (Sahoo et al., 2013; Chou et al., 2014; Guermouche et al., 2020), but none has applied a whole-genome approach. Moreover, sophisticated bioinformatic tools have recently been developed to facilitate accurate interrogation of HTS data for the presence of multiple HCMV strains (Suárez et al., 2019a,b). Here, we used HTS and HCMV-adapted bioinformatics resources to assess the dynamics of HCMV viral populations in 11 transplant recipients undergoing antiviral therapy.
Materials and Methods
Plasma (n = 42) and whole blood (n = 20) samples from confirmed HCMV-infected transplant recipients (nine solid organ transplants and two allogeneic hematopoietic stem cell transplants) were collected longitudinally during episodes of viremia. They were collected at times following transplantation or after HCMV reactivation ranging from 16 to 336 days, and had viral loads ranging from 6.94 × 102 to 8.13 × 106 HCMV international units per ml (IU/ml) of plasma or whole blood. All transplant patients underwent antiviral therapy with valganciclovir, ganciclovir, foscarnet or cidofovir alone or in combination with T cell therapy (Table 1 and Figure 1). Samples were collected with the approval of the University of Sydney Human Research Ethics Committee (reference 2014/440), the National Health Service research Scotland Greater Glasgow and Clyde Biorepository (reference 405) and the Institutional review boards of Hannover Medical School (reference 2527-2014).
Figure 1. Information on the cohort (n = 11): clinical data, level of viremia, and presence of multiple-strain infections and resistance mutations. Data on type of transplant and HCMV serological status (positive, +; and negative, –) of the donor (D) and recipient (R) is given on the left side of each panel (NA, not available). Viral loads (traces) are expressed in log10 [IU/ml] (left Y-axis). Dashed traces indicate lack of HTS data for the corresponding time points. Chronological data (X-axis) are represented in days after transplant, or in days after reactivation (*) if the time between transplant and sample collection was >3 years (patients GLA-2, GLA-3, and GLA-4). Sample collection points lacking corresponding HTS data are shown in square brackets. Sanger sequencing results and application of virus-specific T cell (VST) infusions are indicated at the corresponding time points under the X-axis. Application and duration of antiviral therapy are represented by colored blocks (light green, ganciclovir; dark green, valganciclovir; purple, cidofovir; and orange, foscarnet). Single-strain and multiple-strain infections are represented by single and double traces, respectively. Red traces indicate the presence of fixed resistance mutations (frequency 100%; A594V in patient HAN-4 and L595S in patient HAN-5) in the UL97 gene, and a blue trace indicates the presence of resistance mutations in the UL54 gene (patient SYD-2). In patients SYD-1 and SYD-2, the right Y-axis shows the frequency of low-level resistance mutations (frequency >2%), which are represented by colored triangles embedded in the panels. The resistance mutations detected in each of these patients are indicated on the right side of the panels.
DNA Extraction, Genomic Library Preparation, and DNA Sequencing
Total DNA was extracted from 200 μl of clinical material using the QIAamp MinElute virus spin kit for plasma samples and the QIAamp DNA blood mini kit for whole blood samples (QIAGEN, Crawley, UK). An aliquot of 50 μl of extracted DNA was sheared acoustically using an LE220 sonicator (Covaris, Woburn, MA, USA), aiming to achieve an average fragment size of 500 bp. The fragmented DNA was prepared for sequencing using a library preparation kit (KAPA Biosystems, London, UK) and following the SureSelectXT version 1.7 target enrichment system (Agilent Technologies, Santa Clara, CA, USA) as described previously (Hage et al., 2017; Suárez et al., 2019b). Libraries were indexed using ultrapure (TruGrade) oligonucleotides (Integrated DNA Technologies, Leuven, Belgium), and loaded onto a NextSeq (Illumina, San Diego, CA, USA) DNA sequencer to generate 2 × 150 bp paired-ended reads.
Genome Assembly and Data Analysis
Sequencing data quality was enhanced by removing low-quality regions from the reads using Trim Galore (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) (length = 21, quality = 10, and stringency = 3). The remaining paired-ended reads were screened for the presence of multiple HCMV strains using a method based on detecting signatures unique to the individual genotypes of hypervariable genes (Suárez et al., 2019b). Briefly, multiple strains were assigned if at least two genotypic signatures were found for at least two of the 12 hypervariable genes tested, and if these represented ≥2% of the total number of reads having any genotypic signature of these genes. In the datasets reflecting the presence of a single strain (i.e., presenting single genotypic signatures in at least ten hypervariable genes), the paired-ended reads were assembled de novo using SPAdes 3.5.0 (Bankevich et al., 2012) after removing reads mapping to the Genome Reference Consortium Human Reference 38 sequence (http://genome.ucsc.edu/) using Bowtie2 v.2.2.6 (Langmead and Salzberg, 2012). The contigs were then assembled into genomes using Scaffold_builder (Silva et al., 2013) and the publicly available sequence of HCMV reference strain Merlin (GenBank accession no. AY446894.2). Subsequently, the reads were mapped to the respective assembled genome using Bowtie2, and the alignments were inspected visually using Tablet (Milne et al., 2013). Drug resistance mutations within the antiviral target genes (UL54 and UL97) were screened at two frequency levels: (1) at high frequency by interrogating the UL54 and UL97 gene consensus sequences of the assembled genomes using the online resource Mutation Resistance Analyzer (Chevillotte et al., 2010), which employs HCMV strain TB40/E as reference, and (2) at low frequency by interrogating high-quality datasets (criteria are defined in Results) mapped to corresponding assemblies using LoFreq (Wilm et al., 2012). The latter analysis was conducted using a conservative multi-step filtering approach aimed at minimizing the detection of false positive variants. Briefly, read quality was enhanced using Trim Galore with default parameters, trimmed reads were deduplicated (i.e., duplicated reads were removed) using FastUniq (Xu et al., 2012), and read quality was further improved using PRINSEQ (Schmieder and Edwards, 2011) (minimum quality mean = 25, trimming quality left = 30, trimming quality window = 5, trimming quality step = 1, minimum length = 80, and trimming 3' Ns = 20). Imperfect poly-G tracts located at the 3' ends of reads were trimmed when a segment of 40 nucleotides contained >30 G residues. Filtered sequencing reads were mapped to their respective assemblies using Bowtie2 in end-to-end alignment mode. Reads were further deduplicated on the basis of alignment coordinates using Picard (http://broadinstitute.github.io/picard/). Variants with a quality score of ≥30 were called from the deduplicated alignments.
Datasets were purged of human reads and deposited in the European Nucleotide Archive (ENA; project no. PRJEB36759), and the annotated consensus genome sequences were deposited in GenBank (accession nos. MT044476-MT044485).
Quality Assessment of Sequencing Libraries and Enumeration of Strains
A total of 62 sequencing libraries from 11 transplant patients (2–10 samples per patient) were generated, with input HCMV values ranging from 39 to 1.6 × 106 IU per library. Sequencing yielded 7.6 × 106 to 24.5 × 106 trimmed reads per sample, with 0–87% mapping to the HCMV reference strain Merlin genome (Table S1, rows 14 and 15). When reads were mapped to the strain Merlin genome, the average sequencing depth ranged from 5 to 12,016 reads/nucleotide (nt), covering between 35 and 100% of the Merlin reference genome (Table S1, rows 16 and 17).
To examine library fragment diversity, the average sequencing depth of reads mapped to the strain Merlin genome was calculated after removing duplicates. The average sequencing depth of unique reads ranged from 0.001 (low diversity) to 4,458 (high diversity) reads/nt (Table S1, row 18), and generally reflected the number of HCMV genome copies in the input material (Figure 2).
Figure 2. Scatter plot showing the effect of the number of genome copies used to make a sequencing library on the diversity of the sequencing data obtained. The number of genome copies (IU) is represented on the Y-axis, and the average sequencing depth (reads/nt) of unique HCMV reads is represented on the X-axis. Unique reads were identified as described in Methods. A significant (p < 0.01) positive correlation (Spearman) was noted.
Using the previously published method for distinguishing viral strains by monitoring the genotypes of hypervariable genes (Suárez et al., 2019b), multiple HCMV strains were detected in 15 datasets derived from three patients (SYD-1, SYD-2, and GLA-3). Multiple strains were detected in 4/6 datasets derived from patient SYD-1, 10/10 datasets derived from patient SYD-2, and 1/9 datasets derived from patient GLA-3 (Table 2).
Eight complete and two incomplete HCMV genome sequences were assembled. Complete genome sizes ranged from 235,306 to 235,917 bp. Three genomes presented the complete set of HCMV genes. The remaining seven genomes presented mutations leading to premature termination of at least one gene (most frequently RL5A), a proportion that is in line with previous reports (Sijmons et al., 2015; Suárez et al., 2019b). One genome presented four mutated genes, including UL111A (encoding viral interleukin-10), one genome presented three mutated genes, one genome presented two mutated genes, and the remaining four genomes presented a single mutated gene (Table 3).
Minor Variants and Resistance Mutations
HTS is highly PCR-dependent, which might lead to the incorporation of artifactual variability into sequence datasets. Consequently, only high-quality datasets, that is those with (1) an average sequencing depth of all reads of ≥500 reads/nt, (2) an average sequencing depth of unique reads of ≥10 reads/nt, and (3) ≥95% coverage of the strain Merlin genome, were considered for screening minor variants. In addition, only variants occurring at a frequency of >2% were scored (Xu et al., 2017). Of the 62 sequencing libraries generated, 46 met all quality thresholds. One dataset (SYD-1_T03) met criteria (2) and (3) but exhibited slightly lower average coverage depth (486). With this caveat, this dataset was included in the analysis. The total number of variants across the genome ranged from 0 to 2,151 (Table 4), with low diversity libraries (Figure 3) and multiple-strain infections (Figure 4) showing the highest values.
Table 4. Low frequency variants (>2%) in whole genome sequences and in antiviral target genes UL54 and UL97.
Figure 3. Scatter plot showing the effect of sequencing library diversity and presence of multiple strains on the number of low frequency variants detected. Unique HCMV reads were identified as described in Methods, and the average sequencing depth (reads/nt) of these reads is represented on the Y-axis. The numbers of low frequency (>2%) variants were determined from deduplicated datasets as described in Methods, and are represented on the X-axis. A significant (p < 0.01) negative correlation (Spearman) was noted.
Figure 4. Box-and-whisker plot created using ggplot2 (https://ggplot2.tidyverse.org) showing the number of low frequency (>2%) variants detected by LoFreq in single-strain infections (n = 33 datasets) and multiple-strain infections (n = 14 datasets). Datasets were classified as single-strain or multiple-strain infections using a genotype-based method described previously (Suárez et al., 2019b). Each box encompasses the first to third quartiles (Q1–Q3) and shows the median as a thick line. For each box, the horizontal line at the end of the upper dashed whisker marks the upper extreme (defined as the smaller of Q3+1.5 [Q3–Q1] and the highest single value), and the horizontal line at the end of the lower dashed whisker marks indicates the lower extreme (the greater of Q1–1.5 [Q3–Q1] and the lowest single value). A significant difference was noted between the number of low frequency variants in libraries with multiple strains and those in libraries with single strains (Mann-Whitney test, p < 0.01).
Mutations in the UL97 gene that are known to confer resistance to ganciclovir and valganciclovir were fixed (i.e., present in 100% of the population) in the assembled genomes and at all time points from two solid organ transplant patients, HAN-4 (A594V) and HAN-5 (L595S) (Table 5). The remaining seven solid organ transplant patients did not present resistance mutations in either the UL97 or the UL54 gene. The assembled genomes from both allogeneic hematopoietic stem cell transplant patients presented mutations at low levels (<25%) in the UL54 gene that are known to confer resistance to foscarnet and ganciclovir (patient SYD-1: A809V; and patient SYD-2: V715M, V781I, A809V, and T838A) (Table 6). However, only three of these mutations were observed in samples collected from at least two different time points from patient SYD-2 (V715M from day 166 to 219; V781I at days 177, 205, and 219; and A809V at days 198 and 205), whereas the rest were observed in samples taken at a single time point (Table 6 and Figure 1).
Table 5. High frequency non-synonymous variants in genes UL54 and UL97 identified in assembled genomes.
The emergence of resistance mutations is one of the major risks of antiviral therapy after transplantation (Emery and Griffiths, 2000). In addition, the presence of multiple-strain infections has been linked to poorer outcomes in transplant recipients (Coaquette et al., 2004; Manuel et al., 2009; Lisboa et al., 2012). Therefore, recognizing the presence of resistance mutations and multiple HCMV strains in patients, even at low levels, would provide potentially important information to clinicians for implementing the most appropriate antiviral therapy. Whole-genome HTS makes it possible to screen for the presence of minor variants (including resistance mutations) and the number of viral strains in a single assay with an unprecedented level of sensitivity. However, a comprehensive assessment of the data is required, as sound interpretation is dependent on understanding and monitoring a range of biological and technological factors (Suárez et al., 2019b).
Here, the dynamics of HCMV genome variation were examined in 11 transplant recipients, including solid organ and allogeneic hematopoietic stem cell transplant recipients. As demonstrated in recent studies (Cudini et al., 2019; Suárez et al., 2019b), the presence of multiple strains in clinical samples (i.e., multiple-strain infections) results in gross overestimation of HCMV genome variability within an individual host, as the variants detected almost entirely reflect genetic differences between the strains present rather than intrahost evolution of any single strain. This conclusion was supported in the present study by the much lower level of genome variability observed in confirmed single-strain infections. In addition, sequencing libraries with low diversity, which were usually derived from samples with low viral loads, also presented a high number of variants, possible due to the prolific incorporation of PCR errors. These findings highlight the importance of ensuring sound interpretation of data, particularly as HTS becomes more widely applied to genotypic resistance testing. It is notable that far more progress in applying HTS in the clinical setting has been achieved in the HIV field, and yet numerous quality control issues remain (Weber et al., 2019). In this study, we opted for a conservative approach in assessing the data, by only analyzing datasets meeting specific quality criteria, and scoring variants at a frequency of >2%. In this regard, further studies focused on evaluating accurate quality thresholds, assessing detection limits, and exploring additional sources of misinterpretation are needed in order to allow HTS technology to be applied confidently to HCMV in the clinical setting.
As observed in this and previous studies (Hage et al., 2017; Cudini et al., 2019; Suárez et al., 2019a,b), the HCMV genome is highly stable in patients over time and during different reactivation episodes (patients HAN-1, HAN-2, and HAN-3), as indicated by the low number of variants detected in longitudinal samples from single-strain infections (Figure 4). This finding is consistent with the expression of a proof-reading DNA polymerase by HCMV (Nishiyama et al., 1983), and contrasts with previous studies characterizing HCMV as a fast-evolving virus that readily generates variants in an individual host at a frequency approaching that of RNA viruses such as HIV (Renzette et al., 2013). These discrepancies, which may have important clinical implications, may be explained by differences in the sequencing technologies used, the analytical approaches taken, and the clinical conditions analyzed.
The occurrence of infections involving multiple HCMV strains may disrupt genome stability by providing opportunities for recombination (Rasmussen et al., 2003; Sijmons et al., 2015; Suárez et al., 2019b). Previous studies have demonstrated that multiple-strain infections are common in allogeneic hematopoietic stem cell and solid organ transplant recipients and are associated with a worse outcome, including a higher prevalence of HCMV disease and a higher rate of graft rejection. Both allogeneic hematopoietic stem cell transplant recipients analyzed in the present study presented multiple-strain (n = 2) infections that may have resulted from the positive serological status of both donor and recipient in each case or from the higher level of immunosuppression that such patients typically undergo in comparison to solid organ transplant recipients. However, in order to derive an accurate incidence of multiple-strain infections in transplantation, further studies are warranted that would increase the sample size and represent the various combinations of the serological status of the donor and recipient in both solid organ and allogeneic hematopoietic stem cell transplant cases.
Antiviral drugs may also promote disruption of genome stability by stimulating the selection of variants in genes that render the virus susceptible to those drugs. Approximately 40% of the cases analyzed in the present study developed resistance mutations at various levels in UL54 and UL97. The HCMV genomes sequenced from two solid organ transplant cases (patients HAN-4 and HAN-5) contained fixed mutations in UL97 (A594V and L595S, respectively) that confer resistance to ganciclovir and valganciclovir. This finding is in line with a previous report that the category of solid organ transplants involving a seropositive donor and a seronegative recipient (D+/R–), to which these two cases belonged, account for the highest incidence of resistance (Lurain and Chou, 2010). These mutations were also detected in these patients by Sanger sequencing (Figure 1), but only at specific time points following a recommended procedure (i.e., when the viral load had persisted or increased after 2 weeks of antiviral therapy). However, these mutations stayed fixed in the viral population, even after a change of therapy (patient HAN-5). No resistance mutations were detected in the remaining seven solid organ transplant patients, even though these included three additional D+/R- cases. The incidence of resistance mutations obtained in the present study may not necessarily reflect the general situation, as the number of cases analyzed in this study was limited and included some cases (patients HAN-4 and HAN-5) with previously recognized antiviral resistance mutations detected by Sanger sequencing. In addition, our findings might also have been influenced by the generation of a number of low quality sequencing libraries, such as for patients HAN-2 and HAN-3, or by a shorter time frame of sample collection, such as for the GLA patients (maximum time frame of 31 days between samples). The latter factor is particularly relevant because detection of resistance mutations, except in the immunodeficient pediatric population (Wolf et al., 1998; Houldcroft et al., 2016), is unusual during the first 6 weeks of antiviral therapy (Springer et al., 2005; Lurain and Chou, 2010).
In the present study, both allogeneic hematopoietic stem cell transplant recipients presented resistance mutations (V715M, V781I, A809V, and T838A) at a low level (<25%) in UL54. Mutation V715M, which confers resistance to foscarnet (Baldanti et al., 1996), was detected in patient SYD-2 by Sanger sequencing at day 114 post-transplant after a period of 41 days receiving this antiviral drug. However, this method failed to detect the same mutation 52 days later, as it was present at a low level (6%) (Figure 1 and Table 6). Indeed, this mutation was never cleared despite the change of therapy to ganciclovir or cidofovir. There are several possible explanations for the resilience of this mutation, including the administration of a suboptimal drug dosage, the presence of compensatory mutations in a different region of the HCMV genome, and host factors associated with the level of immunosuppression. Also, in the same patient, mutation V781I, which confers resistance to foscarnet (Cihlar et al., 1998) and ganciclovir (Chou, 2011), was first detected by HTS at 205 days after transplantation at a frequency of <5%, soon after initiation of cidofovir. This mutation was also detected by Sanger sequencing 25 days later (a time point at which no HTS data were available), after completion of cidofovir treatment and initiation of ganciclovir administration. In contrast, mutation A809V, which confers resistance to ganciclovir and foscarnet (Chou et al., 2007), was detected by HTS at 198 days after transplantation and remained for 7 days until being effectively cleared with the initiation of cidofovir. The presence of multiple resistance mutations simultaneously in the same patient, even at low frequency, is of particular interest, as these mutations are unlikely to be detected by traditional non-HTS methods (i.e., Sanger sequencing). Multiple resistance mutations occurring at sub-fixation levels (i.e., present in <100% of the viral population) can lead to a drug resistant phenotype (Chou et al., 2014), a phenomenon that has been observed in fatally immunocompromised pediatric patients (Houldcroft et al., 2016) and solid organ transplant patients (Guermouche et al., 2020). In such circumstances, novel HCMV-specific cell therapies may offer an alternative intervention for patients with clinically resistant disease, especially those found to have clinically resistant HCMV associated with known resistance mutations (Withers et al., 2017). These findings highlight the importance of monitoring resistance mutations with high sensitivity during the course of an antiviral regime, as those present at low level may become prevalent at later stages due to a change to a regime that selects them.
A number of other non-synonymous mutations in UL54 and UL97 were also detected in the cohort. Their capacity to contribute to antiviral resistance is not known and is worthy of further investigation. Moreover, we focused on resistance mutations in UL54 and UL97 as these are the target genes of the antivirals used in the cases analyzed. The application of alternative drugs (Griffiths, 2019), targeting different loci, such as maribavir (UL27) and letermovir (the terminase complex consisting of UL56 and UL89), the potential release of new antiviral drugs acting at other essential HCMV genes, and the possible development of compensatory mutations in other genomic regions, indicates the usefulness of an assay capable of detecting mutations throughout the whole genome, such as the one applied in this study.
Improved treatment of HCMV-infected transplant patients would benefit from more effective monitoring of host and viral determinants. In this regard, HTS provides an opportunity to carry out comprehensive analyses of the evolution of resistance mutants and other minor variants in the HCMV populations present in transplant recipients. The stability of the viral genome in such patients may be perturbed by the application of antiviral drugs or the presence of multiple strains. Monitoring the development of resistance mutations, even at a low frequency, is potentially of use in assisting clinicians to implement the most appropriate antiviral therapy at a time when it is most likely to be effective. However, it is necessary to have a firm understanding of the data, in terms of its quality and the ways in which it is analyzed, in order to ensure that the conclusions are sound and meaningful.
Data Availability Statement
The datasets generated for this study were deposited in the European Nucleotide Archive (ENA; project no. PRJEB36759), and the annotated consensus genome sequences were deposited in GenBank (Accession Nos.: MT044476-MT044485).
The studies involving human samples were collected with the approval of: The University of Sydney human research ethics committee (reference 2014/440), National Health Service research Scotland Greater Glasgow and Clyde biorepository (reference 405), and The Institutional review boards of Hannover Medical School (reference 2527-2014). Written informed consent for participation was obtained for samples collected in the University of Sydney. The samples analysed from Glasgow and Hanover were derived from an anonymised archive of samples. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
NS and AJD contributed to the conception and design of the study. EB, KL, TG, TS, RG, and DG recruited subjects. NS, EB, KL, TG, SA, BW, SL-H, WG, AD, AH, and BS collected and organized the database. NS wrote the first draft of the manuscript. SC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
This work was supported by funding from Tenovus Scotland (grant S17-22 to NS), the Medical Research Council (grant number MC_UU_12014/3 to AJD and KL), the Wellcome Trust (grant number 204870/Z/16/Z to AJD), the Deutsche Forschungsgemeinschaft (Collaborative Research Center 900 to TFS) and the German Center for Infection Research, Thematic Area Infections of the Immunocompromised Host (to TS and TG). BS and EB received Sydney Medical School BioMed-Connect grant funding from the University of Sydney.
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.00267/full#supplementary-material
Table S1. Information on samples, sequencing libraries, genomes assembled and genotyping results.
Baldanti, F., Underwood, M. R., Stanat, S. C., Biron, K. K., Chou, S., Sarasini, A., et al. (1996). Single amino acid changes in the DNA polymerase confer foscarnet resistance and slow-growth phenotype, while mutations in the UL97-encoded phosphotransferase confer ganciclovir resistance in three double-resistant human cytomegalovirus strains recovered from patients with AIDS. J. Virol. 70, 1390–1395. doi: 10.1128/JVI.70.3.1390-1395.1996
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Castagnola, E., Cappelli, B., Erba, D., Rabagliati, A., Lanino, E., and Dini, G. (2004). Cytomegalovirus infection after bone marrow transplantation in children. Hum. Immunol. 65, 416–422. doi: 10.1016/j.humimm.2004.02.013
Chevillotte, M., von Einem, J., Meier, B. M., Lin, F. M., Kestler, H. A., and Mertens, T. (2010). A new tool linking human cytomegalovirus drug resistance mutations to resistance phenotypes. Antiviral Res. 85, 318–327. doi: 10.1016/j.antiviral.2009.10.004
Chou, S., Ercolani, R. J., Sahoo, M. K., Lefterova, M. I., Strasfeld, L. M., and Pinsky, B. A. (2014). Improved detection of emerging drug-resistant mutant cytomegalovirus subpopulations by deep sequencing. Antimicrob. Agents Chemother. 58, 4697–4702. doi: 10.1128/AAC.03214-14
Chou, S., Marousek, G. I., Van Wechel, L. C., Li, S., and Weinberg, A. (2007). Growth and drug resistance phenotypes resulting from cytomegalovirus DNA polymerase region III mutations observed in clinical specimens. Antimicrob. Agents Chemother. 51, 4160–4162. doi: 10.1128/AAC.00736-07
Cihlar, T., Fuller, M. D., and Cherrington, J. M. (1998). Characterization of drug resistance-associated mutations in the human cytomegalovirus DNA polymerase gene by using recombinant mutant viruses generated from overlapping DNA fragments. J. Virol. 72, 5927–5936. doi: 10.1128/JVI.72.7.5927-5936.1998
Coaquette, A., Bourgeois, A., Dirand, C., Varin, A., Chen, W., and Herbein, G. (2004). Mixed cytomegalovirus glycoprotein B genotypes in immunocompromised patients. Clin. Infect. Dis. 39, 155–161. doi: 10.1086/421496
Cudini, J., Roy, S., Houldcroft, C. J., Bryant, J. M., Depledge, D., Tutill, H., et al. (2019). Human cytomegalovirus haplotype reconstruction reveals high diversity due to superinfection and evidence of within-host recombination. Proc. Natl. Acad. Sci. U.S.A. 116, 5693–5698. doi: 10.1073/pnas.1818130116
Guermouche, H., Burrel, S., Mercier-Darty, M., Kofman, T., Rogier, O., Pawlotsky, J. M., et al. (2020). Characterization of the dynamics of human cytomegalovirus resistance to antiviral drugs by ultra-deep sequencing. Antiviral Res. 173:104647. doi: 10.1016/j.antiviral.2019.104647
Hage, E., Wilkie, G. S., Linnenweber-Held, S., Dhingra, A., Suárez, N. M., Schmidt, J. J., et al. (2017). Characterization of human cytomegalovirus genome diversity in immunocompromised hosts by whole-genome sequencing directly from clinical specimens. J. Infect. Dis. 215, 1673–1683. doi: 10.1093/infdis/jix157
Houldcroft, C. J., Bryant, J. M., Depledge, D., Margetts, B. K., Simmonds, J., Nicolaou, S., et al. (2016). Detection of low frequency multi-drug resistance and novel putative maribavir resistance in immunocompromised pediatric patients with cytomegalovirus. Front. Microbiol. 7:1317. doi: 10.3389/fmicb.2016.01317
Lisboa, L. F., Tong, Y., Kumar, D., Pang, X. L., Asberg, A., Hartmann, A., et al. (2012). Analysis and clinical correlation of genetic variation in cytomegalovirus. Transpl. Infect Dis. 14, 132–140. doi: 10.1111/j.1399-3062.2011.00685.x
Manuel, O., Asberg, A., Pang, X., Rollag, H., Emery, C., Preiksaitis, J. K., et al. (2009). Impact of genetic polymorphisms in cytomegalovirus glycoprotein B on outcomes in solid-organ transplant recipients with cytomegalovirus disease. Clin. Infect. Dis. 49, 1160–1166. doi: 10.1086/605633
Milne, I., Stephen, G., Bayer, M., Cock, J., Pritchard, L., Cardle, L., et al. (2013). Using tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 14, 193–202. doi: 10.1093/bib/bbs012
Nishiyama, Y., Maeno, K., and Yoshida, S. (1983). Characterization of human cytomegalovirus-induced DNA polymerase and the associated 3'-to-5', exonuclease. Virology 124, 221–231. doi: 10.1016/0042-6822(83)90339-2
Renzette, N., Gibson, L., Bhattacharjee, B., Fisher, D., Schleiss, M. R., Jensen, J. D., et al. (2013). Rapid intrahost evolution of human cytomegalovirus is shaped by demography and positive selection. PLoS Genet. 9:e1003735. doi: 10.1371/journal.pgen.1003735
Sahoo, M. K., Lefterova, M. I., Yamamoto, F., Waggoner, J. J., Chou, S., Holmes, S., et al. (2013). Detection of cytomegalovirus drug resistance mutations by next-generation sequencing. J. Clin. Microbiol. 51, 3700–3710. doi: 10.1128/JCM.01605-13
Schuurman, R., Demeter, L., Reichelderfer, P., Tijnagel, J., de Groot, T., and Boucher, C. (1999). Worldwide evaluation of DNA sequencing approaches for identification of drug resistance mutations in the human immunodeficiency virus type 1 reverse transcriptase. J. Clin. Microbiol. 37, 2291–2296. doi: 10.1128/JCM.37.7.2291-2296.1999
Scott, G. M., Weinberg, A., Rawlinson, W. D., and Chou, S. (2007). Multidrug resistance conferred by novel DNA polymerase mutations in human cytomegalovirus isolates. Antimicrob. Agents Chemother. 51, 89–94. doi: 10.1128/AAC.00633-06
Sijmons, S., Thys, K., Mbong Ngwese, M., Van Damme, E., Dvorak, J., Van Loock, M., et al. (2015). High-throughput analysis of human cytomegalovirus genome diversity highlights the widespread occurrence of gene-disrupting mutations and pervasive recombination. J. Virol. 89, 7673–7695. doi: 10.1128/JVI.00578-15
Silva, G. G., Dutilh, B. E., Matthews, T. D., Elkins, K., Schmieder, R., Dinsdale, E. A., et al. (2013). Combining de novo and reference-guided assembly with scaffold_builder. Source Code Biol. Med. 8:23. doi: 10.1186/1751-0473-8-23
Springer, K. L., Chou, S., Li, S., Giller, R. H., Quinones, R., Shira, J. E., et al. (2005). How evolution of mutations conferring drug resistance affects viral dynamics and clinical outcomes of cytomegalovirus-infected hematopoietic cell transplant recipients. J. Clin. Microbiol. 43, 208–213. doi: 10.1128/JCM.43.1.208-213.2005
Suárez, N. M., Musonda, K. G., Escriva, E., Njenga, M., Agbueze, A., Camiolo, S., et al. (2019a). Multiple-strain infections of human cytomegalovirus with high genomic diversity are common in breast milk from human immunodeficiency virus-infected women in Zambia. J. Infect. Dis. 220, 792–801. doi: 10.1093/infdis/jiz209
Suárez, N. M., Wilkie, G. S., Hage, E., Camiolo, S., Holton, M., Hughes, J., et al. (2019b). Human cytomegalovirus genomes sequenced directly from clinical material. variation, multiple-strain infection, recombination, and gene loss. J. Infect. Dis. 220, 781–791. doi: 10.1093/infdis/jiz208
Weber, J., Volkova, I., Sahoo, M. K., Tzou, L., Shafer, R. W., and Pinsky, B. A. (2019). Prospective evaluation of the vela diagnostics next-generation sequencing platform for HIV-1 genotypic resistance testing. J. Mol. Diagn. 21, 961–970. doi: 10.1016/j.jmoldx.2019.06.003
Wilm, A., Aw, B. D., Yeo, G. H., Ong, S. H., Wong, C. H., Khor, C. C., et al. (2012). LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res. 40, 11189–11201. doi: 10.1093/nar/gks918
Withers, B., Blyth, E., Clancy, L. E., Yong, A., Fraser, C., Burgess, J., et al. (2017). Long-term control of recurrent or refractory viral infections after allogeneic HSCT with third-party virus-specific T cells. Blood Adv. 1, 2193–2205. doi: 10.1182/bloodadvances.2017010223
Wolf, D. G., Yaniv, I., Honigman, A., Kassis, I., Schonfeld, T., and Ashkenazi, S. (1998). Early emergence of ganciclovir-resistant human cytomegalovirus strains in children with primary combined immunodeficiency. J. Infect. Dis. 178, 535–538. doi: 10.1086/517468
Wu, X., Wang, Y., Xu, Y., Wu, D., Sun, A., Zhu, Z., et al. (2010). Cytomegalovirus glycoprotein B genotype in hematopoietic stem cell transplant patients from China. Biol. Blood Marrow Transplant. 16, 647–652. doi: 10.1016/j.bbmt.2009.12.001
Xu, C., Nezami Ranjbar, M. R., Wu, Z., DiCarlo, J., and Wang, Y. (2017). Detecting very low allele fraction variants using targeted DNA sequencing and a novel molecular barcode-aware variant caller. BMC Genomics. 18:5. doi: 10.1186/s12864-016-3425-4
Keywords: human cytomegalovirus, transplantation, genome sequence, target enrichment, multiple-strain infection, antiviral therapy, resistance mutation
Citation: Suárez NM, Blyth E, Li K, Ganzenmueller T, Camiolo S, Avdic S, Withers B, Linnenweber-Held S, Gwinner W, Dhingra A, Heim A, Schulz TF, Gunson R, Gottlieb D, Slobedman B and Davison AJ (2020) Whole-Genome Approach to Assessing Human Cytomegalovirus Dynamics in Transplant Patients Undergoing Antiviral Therapy. Front. Cell. Infect. Microbiol. 10:267. doi: 10.3389/fcimb.2020.00267
Received: 14 February 2020; Accepted: 05 May 2020;
Published: 15 June 2020.
Edited by:Michael Nevels, University of St. Andrews, United Kingdom
Reviewed by:David Navarro, Hospital Clínic Universitari de València, Spain
Pilar Perez-Romero, Carlos III Health Institute (ISCIII), Spain
Copyright © 2020 Suárez, Blyth, Li, Ganzenmueller, Camiolo, Avdic, Withers, Linnenweber-Held, Gwinner, Dhingra, Heim, Schulz, Gunson, Gottlieb, Slobedman and Davison. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.