Genomic surveillance reveals dynamic shifts in the connectivity of COVID-19 epidemics.

The maturation of genomic surveillance in the past decade has enabled tracking of the emergence and spread of epidemics at an unprecedented level. During the COVID-19 pandemic, for example, genomic data revealed that local epidemics varied considerably in the frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage importation and persistence, likely due to a combination of COVID-19 restrictions and changing connectivity. Here, we show that local COVID-19 epidemics are driven by regional transmission, including across international boundaries, but can become increasingly connected to distant locations following the relaxation of public health interventions. By integrating genomic, mobility, and epidemiological data, we find abundant transmission occurring between both adjacent and distant locations, supported by dynamic mobility patterns. We find that changing connectivity significantly influences local COVID-19 incidence. Our findings demonstrate a complex meaning of "local" when investigating connected epidemics and emphasize the importance of collaborative interventions for pandemic prevention and mitigation.

The Suitability of RNA from Positive SARS-CoV-2 Rapid Antigen Tests for Whole Virus Genome Sequencing and Variant Identification to Maintain Genomic Surveillance.

The COVID-19 pandemic has transformed laboratory management, with a surge in demand for diagnostic tests prompting the adoption of new diagnostic assays and the spread of variant surveillance tools. Rapid antigen tests (RATs) were initially used only for screening and later as suitable infection assessment tools. This study explores the feasibility of sequencing the SARS-CoV-2 genome from the residue of the nasopharyngeal swab extraction buffers of rapid antigen tests (RATs) to identify different COVID-19 lineages and sub-lineages. METHODS: Viral RNA was extracted from the residue of the nasopharyngeal swab extraction buffers of RATs and, after a confirmation of positivity through a reaction of RT-PCR, viral genome sequencing was performed. RESULTS: Overall, the quality of the sequences obtained from the RNA extracted from the residue of the nasopharyngeal swab extraction buffers of RATs was adequate and allowed us to identify the SARS-CoV-2 variants' circulation and distribution in a period when the use of molecular swabs had been drastically reduced. CONCLUSIONS: This study demonstrates the potential for genomic surveillance by sequencing SARS-CoV-2 from the residue of the nasopharyngeal swab extraction buffers of RATs, highlighting alternative possibilities for tracking variants.

A laboratory framework for ongoing optimization of amplification-based genomic surveillance programs.

IMPORTANCE: This study provides a laboratory framework to ensure ongoing relevance and performance of amplification-based whole genome sequencing to strengthen public health surveillance during extended outbreaks or pandemics. The framework integrates regular reviews of the performance of a genomic surveillance system and highlights the importance of ongoing monitoring and the identification and implementation of improvements to whole genome sequencing methods to enhance public health responses to pathogen outbreaks.

Genotypic testing improves detection of antiviral resistance in human herpes simplex virus.

BACKGROUND: Antiviral resistance in human herpes simplex viruses (HSV) remains a significant clinical challenge in immunocompromised populations. Although molecular tests have largely replaced viral culture for HSV diagnosis and molecular antiviral resistance testing is available for many viruses, HSV resistance testing continues to rely on phenotypic, viral culture-based methods, requiring weeks for results. Consequently, treatment of suspected HSV resistance remains largely empiric. METHODS: We used HSV whole genome sequencing and a database of previously characterized HSV acyclovir and foscarnet resistance mutations to evaluate the performance of genotypic antiviral resistance testing among 19 control strains compared to in-house plaque reduction assay (PRA) and 25 clinical isolates sent for reference lab PRA antiviral resistance testing. RESULTS: Among control strains, 23/29 (79.3%) results were concordant, 5 (17.2%) were indeterminate, and 1 (3.4%) was discordant. Indeterminate results were caused by variants of uncertain significance (VUS), including mutations without published phenotypes and mutations with contradictory results. Among clinical isolates, 14/40 (35%) results were concordant, 17 (42.5%) were indeterminate, and 9 (22.5%) were discordant. All discordant results were in reportedly phenotypically-susceptible HSV-1 strains yet possessed resistance mutations. Three contained resistant subpopulations. 6/8 (75%) discordant phenotypes were concordant with resistant genotypes upon repeat PRA. CONCLUSIONS: These data support the combination of genotypic and phenotypic testing to diagnose HSV resistance more accurately and likely more rapidly than phenotypic testing alone. Genotypic context of resistance mutations and the ability of viral strains to form plaques in culture may affect phenotypic resistance results, highlighting the limitations of PRA alone as a gold standard method.

Inspecting the interaction between human immunodeficiency virus and the immune system through genetic turnover.

Chronic infections of the human immunodeficiency virus (HIV) create a very complex coevolutionary process, where the virus tries to escape the continuously adapting host immune system. Quantitative details of this process are largely unknown and could help in disease treatment and vaccine development. Here we study a longitudinal dataset of ten HIV-infected people, where both the B-cell receptors and the virus are deeply sequenced. We focus on simple measures of turnover, which quantify how much the composition of the viral strains and the immune repertoire change between time points. At the single-patient level, the viral-host turnover rates do not show any statistically significant correlation, however, they correlate if one increases the amount of statistics by aggregating the information across patients. We identify an anti-correlation: large changes in the viral pool composition come with small changes in the B-cell receptor repertoire. This result seems to contradict the naïve expectation that when the virus mutates quickly, the immune repertoire needs to change to keep up. However, a simple model of antagonistically evolving populations can explain this signal. If it is sampled at intervals comparable with the sweep time, one population has had time to sweep while the second cannot start a counter-sweep, leading to the observed anti-correlation. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA.

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Detection and sequencing of Zika virus in normocephalic newborns with congenital Zika infection.

Fourteen asymptomatic normocephalic newborns with confirmed congenital Zika infection were investigated. All newborns presented Zika virus (ZIKV) positivity on reverse transcriptase polymerase chain reaction. Following ZIKV-specific NS5 gene fragment sequencing in one child, phylogenetic analysis revealed that this isolate belonged to the Asian genotype, and clustered closely with other sequences previously isolated in north-east and northern regions of Brazil.

Distinguishing Features Common to Dual Fatal Herpes Simplex Virus Infections That Occur in Both a Pregnant Woman and Her Newborn Infant.

Deaths from herpes simplex virus type 1 (HSV-1) and herpes simplex virus type 2 (HSV-2) are rare. A major exception is perinatally acquired HSV-1 or HSV-2 infection where the neonatal death rate is substantial. Fatal HSV infection also occurs occasionally in pregnant women. The goal of this review is to enumerate the reports that describe dual deaths of both a pregnant woman and her newborn from a herpesvirus infection. A total of 15 reports were found in the medical literature, of which five described pregnant women with HSV encephalitis and 10 described women with disseminated HSV infection. When the virus was typed, most cases of dual mother/newborn deaths were caused by HSV-2. Of interest, in two situations caused by HSV-1, the pregnant woman probably acquired her primary HSV-1 infection from one of her children and not by sexual transmission. Complete genomic sequencing was performed on one set of HSV-1 isolates collected from mother (blood) and newborn (blood and skin). The mother's strain and the newborn's skin strain were 98.9% identical. When the newborn's two strains were compared, they were 97.4% identical. Only one mother was tested by the HerpeSelect IgG antibody kit. During the nine days of her undiagnosed disseminated infection preceding her death, her serology was negative. In summary, although dual mother/newborn deaths from HSV infection are rare, they continue to be reported as recently as 2017.

Combined single-cell transcriptional, translational, and genomic profiling reveals HIV-1 reservoir diversity.

Although understanding the diversity of HIV-1 reservoirs is key to achieving a cure, their study at the single-cell level in primary samples remains challenging. We combine flow cytometric multiplexed fluorescent in situ RNA hybridization for different viral genes with HIV-1 p24 protein detection, cell phenotyping, and downstream near-full-length single-cell vDNA sequencing. Stimulation-induced viral RNA-positive (vRNA+) cells from viremic and antiretroviral-therapy (ART)-suppressed individuals differ in their ability to produce p24. In participants on ART, latency-reversing agents (LRAs) induce a wide variety of viral gene transcription and translation patterns with LRA class-specific differences in reactivation potency. Reactivated proviruses, including in p24+ cells, are mostly defective. Although LRAs efficiently induce transcription in all memory cell subsets, we observe induction of translation mostly in effector memory cells, rather than in the long-lived central memory pool. We identify HIV-1 clones with diverse transcriptional and translational patterns between individual cells, and this finding suggests that cell-intrinsic factors influence reservoir persistence and heterogeneity.

Emergence of an early SARS-CoV-2 epidemic in the United States.

The emergence of the COVID-19 epidemic in the United States (U.S.) went largely undetected due to inadequate testing. New Orleans experienced one of the earliest and fastest accelerating outbreaks, coinciding with Mardi Gras. To gain insight into the emergence of SARS-CoV-2 in the U.S. and how large-scale events accelerate transmission, we sequenced SARS-CoV-2 genomes during the first wave of the COVID-19 epidemic in Louisiana. We show that SARS-CoV-2 in Louisiana had limited diversity compared to other U.S. states and that one introduction of SARS-CoV-2 led to almost all of the early transmission in Louisiana. By analyzing mobility and genomic data, we show that SARS-CoV-2 was already present in New Orleans before Mardi Gras, and the festival dramatically accelerated transmission. Our study provides an understanding of how superspreading during large-scale events played a key role during the early outbreak in the U.S. and can greatly accelerate epidemics.

Field and Molecular Epidemiology: How Viral Sequencing Changed Transmission Inferences in the First Portuguese SARS-CoV-2 Infection Cluster.

Field epidemiology and viral sequencing provide a comprehensive characterization of transmission chains and allow a better identification of superspreading events. However, very few examples have been presented to date during the COVID-19 pandemic. We studied the first COVID-19 cluster detected in Portugal (59 individuals involved amongst extended family and work environments), following the return of four related individuals from work trips to Italy. The first patient to introduce the virus would be misidentified following the traditional field inquiry alone, as shown by the viral sequencing in isolates from 23 individuals. The results also pointed out family, and not work environment, as the primary mode of transmission.

High Prevalence of SARS-CoV-2 Genetic Variation and D614G Mutation in Pediatric Patients With COVID-19.

BACKGROUND: The full spectrum of the disease phenotype and viral genotype of coronavirus disease 2019 (COVID-19) have yet to be thoroughly explored in children. Here, we analyze the relationships between viral genetic variants and clinical characteristics in children. METHODS: Whole-genome sequencing was performed on respiratory specimens collected for all SARS-CoV-2-positive children (n = 141) between March 13 and June 16, 2020. Viral genetic variations across the SARS-CoV-2 genome were identified and investigated to evaluate genomic correlates of disease severity. RESULTS: Higher viral load was detected in symptomatic patients (P = .0007) and in children <5 years old (P = .0004). Genomic analysis revealed a mean pairwise difference of 10.8 single nucleotide variants (SNVs), and the majority (55.4%) of SNVs led to an amino acid change in the viral proteins. The D614G mutation in the spike protein was present in 99.3% of the isolates. The calculated viral mutational rate of 22.2 substitutions/year contrasts the 13.5 substitutions/year observed in California isolates without the D614G mutation. Phylogenetic clade 20C was associated with severe cases of COVID-19 (odds ratio, 6.95; P = .0467). Epidemiological investigation revealed major representation of 3 of 5 major Nextstrain clades (20A, 20B, and 20C) consistent with multiple introductions of SARS-CoV-2 in Southern California. CONCLUSIONS: Genomic evaluation demonstrated greater than expected genetic diversity, presence of the D614G mutation, increased mutation rate, and evidence of multiple introductions of SARS-CoV-2 into Southern California. Our findings suggest a possible association of phylogenetic clade 20C with severe disease, but small sample size precludes a definitive conclusion. Our study warrants larger and multi-institutional genomic evaluation and has implications for infection control practices.

Viral Sequencing to Investigate Sources of SARS-CoV-2 Infection in US Healthcare Personnel.

BACKGROUND: Healthcare personnel (HCP) are at increased risk of infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We posit that current infection control guidelines generally protect HCP from SARS-CoV-2 infection in a healthcare setting. METHODS: In this retrospective case series, we used viral genomics to investigate the likely source of SARS-CoV-2 infection in HCP at a major academic medical institution in the Upper Midwest of the United States between 25 March and 27 December 2020. We obtained limited epidemiological data through informal interviews and review of the electronic health record and combined this information with healthcare-associated viral sequences and viral sequences collected in the broader community to infer the most likely source of infection in HCP. RESULTS: We investigated SARS-CoV-2 infection clusters involving 95 HCP and 137 possible patient contact sequences. The majority of HCP infections could not be linked to a patient or coworker (55 of 95 [57.9%]) and were genetically similar to viruses circulating concurrently in the community. We found that 10.5% of HCP infections (10 of 95) could be traced to a coworker. Strikingly, only 4.2% (4 of 95) could be traced to a patient source. CONCLUSIONS: Infections among HCP add further strain to the healthcare system and put patients, HCP, and communities at risk. We found no evidence for healthcare-associated transmission in the majority of HCP infections evaluated. Although we cannot rule out the possibility of cryptic healthcare-associated transmission, it appears that HCP most commonly become infected with SARS-CoV-2 via community exposure. This emphasizes the ongoing importance of mask wearing, physical distancing, robust testing programs, and rapid distribution of vaccines.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings.

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Molecular Detection of SARS-CoV-2 Infection in FFPE Samples and Histopathologic Findings in Fatal SARS-CoV-2 Cases.

OBJECTIVES: To report methods and findings of 2 autopsies with molecular evaluation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) positive individuals. METHODS: Postmortem examination was completed following Centers for Disease Control and Prevention public guidelines. Numerous formalin-fixed paraffin-embedded (FFPE) tissue types from each case were surveyed for SARS-CoV-2 RNA by quantitative reverse transcription polymerase chain reaction (qRT-PCR). SARS-CoV-2 viral genome was sequenced by next-generation sequencing (NGS) from FFPE lung tissue blocks. RESULTS: Postmortem examinations revealed diffuse alveolar damage, while no viral-associated hepatic, cardiac, or renal damage was observed. Viral RNA was detected in lungs, bronchi, lymph nodes, and spleen in both cases using qRT-PCR method. RNA sequencing using NGS in case 1 revealed mutations most consistent with Western European Clade A2a with ORF1a L3606F mutation. CONCLUSIONS: SARS-CoV-2 testing and viral sequencing can be performed from FFPE tissue. Detection and sequencing of SARS-CoV-2 in combination with morphological findings from postmortem tissue examination can aid in gaining a better understanding of the virus's pathophysiologic effects on human health.

Next-generation sequencing of dsRNA is greatly improved by treatment with the inexpensive denaturing reagent DMSO.

dsRNA is the genetic material of important viruses and a key component of RNA interference-based immunity in eukaryotes. Previous studies have noted difficulties in determining the sequence of dsRNA molecules that have affected studies of immune function and estimates of viral diversity in nature. DMSO has been used to denature dsRNA prior to the reverse-transcription stage to improve reverse transcriptase PCR and Sanger sequencing. We systematically tested the utility of DMSO to improve the sequencing yield of a dsRNA virus (Φ6) in a short-read next-generation sequencing platform. DMSO treatment improved sequencing read recovery by over two orders of magnitude, even when RNA and cDNA concentrations were below the limit of detection. We also tested the effects of DMSO on a mock eukaryotic viral community and found that dsRNA virus reads increased with DMSO treatment. Furthermore, we provide evidence that DMSO treatment does not adversely affect recovery of reads from a ssRNA viral genome (influenza A/California/07/2009). We suggest that up to 50 % DMSO treatment be used prior to cDNA synthesis when samples of interest are composed of or may contain dsRNA.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar.

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

FATAL RANAVIRUS INFECTION IN A GROUP OF ZOO-HOUSED MELLER'S CHAMELEONS (TRIOCEROS MELLERI).

A group of five juvenile Meller's chameleons (Trioceros melleri) experienced 100% mortality over a period of 1 mo due to ranavirus infection. The index case was found dead without premonitory signs. The three subsequent cases presented with nonspecific clinical signs (lethargy, decreased appetite, ocular discharge) and were ultimately euthanatized. The final case died after initially presenting with skin lesions. Postmortem examination revealed thin body condition in all five animals and mild coelomic effusion and petechiae affecting the tongue and kidneys of one animal. Microscopically, all animals had multifocal necrosis of the spleen, liver, and kidney; four of five animals had necrosis of the nasal cavity; and two of five had necrosis of adrenal tissue, bone marrow, and skin. Numerous basophilic intracytoplasmic inclusions were present in the liver of all animals and nasal mucosa of three of the five animals. Consensus polymerase chain reaction for herpesvirus and adenovirus were negative, whereas ranavirus quantitative polymerase chain reaction was positive. Virus isolation followed by whole genome sequencing and Bayesian phylogenetic analysis classified the isolates as a strain of frog virus 3 (FV3) most closely related to an FV3 isolate responsible for a previous outbreak in the zoo's eastern box turtle (Terrapene carolina carolina) group. This case series documents the first known occurrence of ranavirus-associated disease in chameleons and demonstrates the potential for interspecies transmission between chelonian and squamate reptiles.

Diagnosis of Fatal Human Case of St. Louis Encephalitis Virus Infection by Metagenomic Sequencing, California, 2016.

We used unbiased metagenomic next-generation sequencing to diagnose a fatal case of meningoencephalitis caused by St. Louis encephalitis virus in a patient from California in September 2016. This case is associated with the recent 2015-2016 reemergence of this virus in the southwestern United States.

Rapid Metagenomic Next-Generation Sequencing during an Investigation of Hospital-Acquired Human Parainfluenza Virus 3 Infections.

Metagenomic next-generation sequencing (mNGS) is increasingly used for the unbiased detection of viruses, bacteria, fungi, and eukaryotic parasites in clinical samples. Whole-genome sequencing (WGS) of clinical bacterial isolates has been shown to inform hospital infection prevention practices, but this technology has not been utilized during potential respiratory virus outbreaks. Here, we report on the use of mNGS to inform the real-time infection prevention response to a cluster of hospital-acquired human parainfluenza 3 virus (HPIV3) infections at a children's hospital. Samples from 3 patients with hospital-acquired HPIV3 identified over a 12-day period on a general medical unit and 10 temporally associated samples from patients with community-acquired HPIV3 were analyzed. Our sample-to-sequencer time was <24 h, while our sample-to-answer turnaround time was <60 h with a hands-on time of approximately 6 h. Eight (2 cases and 6 controls) of 13 samples had sufficient sequencing coverage to yield the whole genome for HPIV3, while 10 (2 cases and 8 controls) of 13 samples gave partial genomes and all 13 samples had >1 read for HPIV3. Phylogenetic clustering revealed the presence of identical HPIV3 genomic sequence in the two of the cases with hospital-acquired infection, consistent with the concern for recent transmission within the medical unit. Adequate sequence coverage was not recovered for the third case. This work demonstrates the promise of mNGS for providing rapid information for infection prevention in addition to microbial detection.