Emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States.

The highly transmissible B.1.1.7 variant of SARS-CoV-2, first identified in the United Kingdom, has gained a foothold across the world. Using S gene target failure (SGTF) and SARS-CoV-2 genomic sequencing, we investigated the prevalence and dynamics of this variant in the United States (US), tracking it back to its early emergence. We found that, while the fraction of B.1.1.7 varied by state, the variant increased at a logistic rate with a roughly weekly doubling rate and an increased transmission of 40%-50%. We revealed several independent introductions of B.1.1.7 into the US as early as late November 2020, with community transmission spreading it to most states within months. We show that the US is on a similar trajectory as other countries where B.1.1.7 became dominant, requiring immediate and decisive action to minimize COVID-19 morbidity and mortality.

Evolution and Diversity of Bat and Rodent Paramyxoviruses from North America.

Paramyxoviruses are a diverse group of negative-sense, single-stranded RNA viruses of which several species cause significant mortality and morbidity. In recent years the collection of paramyxovirus sequences detected in wild mammals has substantially grown; however, little is known about paramyxovirus diversity in North American mammals. To better understand natural paramyxovirus diversity, host range, and host specificity, we sought to comprehensively characterize paramyxoviruses across a range of diverse cooccurring wild small mammals in southern Arizona. We used highly degenerate primers to screen fecal and urine samples and obtained a total of 55 paramyxovirus sequences from 12 rodent species and 6 bat species. We also performed Illumina transcriptome sequencing (RNA-seq) and de novo assembly on 14 of the positive samples to recover a total of 5 near-full-length viral genomes. We show there are at least two clades of rodent-borne paramyxoviruses in Arizona, while bat-associated paramyxoviruses formed a putative single clade. Using structural homology modeling of the viral attachment protein, we infer that three of the five novel viruses likely bind sialic acid in a manner similar to other respiroviruses, while the other two viruses from heteromyid rodents likely bind a novel host receptor. We find no evidence for cross-species transmission, even among closely related sympatric host species. Taken together, these data suggest paramyxoviruses are a common viral infection in some bat and rodent species present in North America and illuminate the evolution of these viruses. IMPORTANCE There are a number of viral lineages that are potential zoonotic threats to humans. One of these, paramyxoviruses have jumped into humans multiple times from wild and domestic animals. We conducted one of the largest viral surveys of wild mammals in the United States to better understand paramyxovirus diversity and evolution.

A near full-length HIV-1 genome from 1966 recovered from formalin-fixed paraffin-embedded tissue.

With very little direct biological data of HIV-1 from before the 1980s, far-reaching evolutionary and epidemiological inferences regarding the long prediscovery phase of this pandemic are based on extrapolations by phylodynamic models of HIV-1 genomic sequences gathered mostly over recent decades. Here, using a very sensitive multiplex RT-PCR assay, we screened 1,645 formalin-fixed paraffin-embedded tissue specimens collected for pathology diagnostics in Central Africa between 1958 and 1966. We report the near-complete viral genome in one HIV-1 positive specimen from Kinshasa, Democratic Republic of Congo (DRC), from 1966 ("DRC66")-a nonrecombinant sister lineage to subtype C that constitutes the oldest HIV-1 near full-length genome recovered to date. Root-to-tip plots showed the DRC66 sequence is not an outlier as would be expected if dating estimates from more recent genomes were systematically biased; and inclusion of the DRC66 sequence in tip-dated BEAST analyses did not significantly alter root and internal node age estimates based on post-1978 HIV-1 sequences. There was larger variation in divergence time estimates among datasets that were subsamples of the available HIV-1 genomes from 1978 to 2014, showing the inherent phylogenetic stochasticity across subsets of the real HIV-1 diversity. Our phylogenetic analyses date the origin of the pandemic lineage of HIV-1 to a time period around the turn of the 20th century (1881 to 1918). In conclusion, this unique archival HIV-1 sequence provides direct genomic insight into HIV-1 in 1960s DRC, and, as an ancient-DNA calibrator, it validates our understanding of HIV-1 evolutionary history.

Discovery of three cycloviruses in fecal samples from silver-haired bats (Lasionycteris noctivagans) in Arizona (USA).

Bats harbour a diverse array of viruses, some of which are zoonotic, and are one of the most speciose groups of mammals on earth. As part of an ongoing bat-associated viral diversity research project, we identified three cycloviruses (family Circoviridae) in fecal samples of silver-haired bats (Lasionycteris noctivagans) caught in Cave Creek Canyon of Arizona (USA). Two of the three identified genomes represent two new species in the genus Cyclovirus. Cycloviruses have been found in a wide range of environments and hosts; however, little is known about their biology. These new genomes of cycloviruses are the first from silver-haired bats, adding to the broader knowledge of cyclovirus diversity. With continuing studies, it is likely that additional viruses of the family Circoviridae will be identified in Arizona bat populations.

Increased HIV-1 vaccine efficacy against viruses with genetic signatures in Env V2.

The RV144 trial demonstrated 31% vaccine efficacy at preventing human immunodeficiency virus (HIV)-1 infection. Antibodies against the HIV-1 envelope variable loops 1 and 2 (Env V1 and V2) correlated inversely with infection risk. We proposed that vaccine-induced immune responses against V1/V2 would have a selective effect against, or sieve, HIV-1 breakthrough viruses. A total of 936 HIV-1 genome sequences from 44 vaccine and 66 placebo recipients were examined. We show that vaccine-induced immune responses were associated with two signatures in V2 at amino acid positions 169 and 181. Vaccine efficacy against viruses matching the vaccine at position 169 was 48% (confidence interval 18% to 66%; P = 0.0036), whereas vaccine efficacy against viruses mismatching the vaccine at position 181 was 78% (confidence interval 35% to 93%; P = 0.0028). Residue 169 is in a cationic glycosylated region recognized by broadly neutralizing and RV144-derived antibodies. The predicted distance between the two signature sites (21 ± 7 Å) and their match/mismatch dichotomy indicate that multiple factors may be involved in the protection observed in RV144. Genetic signatures of RV144 vaccination in V2 complement the finding of an association between high V1/V2-binding antibodies and reduced risk of HIV-1 acquisition, and provide evidence that vaccine-induced V2 responses plausibly had a role in the partial protection conferred by the RV144 regimen.

Fitness-Balanced Escape Determines Resolution of Dynamic Founder Virus Escape Processes in HIV-1 Infection.

UNLABELLED: To understand the interplay between host cytotoxic T-lymphocyte (CTL) responses and the mechanisms by which HIV-1 evades them, we studied viral evolutionary patterns associated with host CTL responses in six linked transmission pairs. HIV-1 sequences corresponding to full-length p17 and p24 gag were generated by 454 pyrosequencing for all pairs near the time of transmission, and seroconverting partners were followed for a median of 847 days postinfection. T-cell responses were screened by gamma interferon/interleukin-2 (IFN-γ/IL-2) FluoroSpot using autologous peptide sets reflecting any Gag variant present in at least 5% of sequence reads in the individual's viral population. While we found little evidence for the occurrence of CTL reversions, CTL escape processes were found to be highly dynamic, with multiple epitope variants emerging simultaneously. We found a correlation between epitope entropy and the number of epitope variants per response (r = 0.43; P = 0.05). In cases in which multiple escape mutations developed within a targeted epitope, a variant with no fitness cost became fixed in the viral population. When multiple mutations within an epitope achieved fitness-balanced escape, these escape mutants were each maintained in the viral population. Additional mutations found to confer escape but undetected in viral populations incurred high fitness costs, suggesting that functional constraints limit the available sites tolerable to escape mutations. These results further our understanding of the impact of CTL escape and reversion from the founder virus in HIV infection and contribute to the identification of immunogenic Gag regions most vulnerable to a targeted T-cell attack. IMPORTANCE: Rapid diversification of the viral population is a hallmark of HIV-1 infection, and understanding the selective forces driving the emergence of viral variants can provide critical insight into the interplay between host immune responses and viral evolution. We used deep sequencing to comprehensively follow viral evolution over time in six linked HIV transmission pairs. We then mapped T-cell responses to explore if mutations arose due to adaption to the host and found that escape processes were often highly dynamic, with multiple mutations arising within targeted epitopes. When we explored the impact of these mutations on replicative capacity, we found that dynamic escape processes only resolve with the selection of mutations that conferred escape with no fitness cost to the virus. These results provide further understanding of the complicated viral-host interactions that occur during early HIV-1 infection and may help inform the design of future vaccine immunogens.

Comprehensive sieve analysis of breakthrough HIV-1 sequences in the RV144 vaccine efficacy trial.

The RV144 clinical trial showed the partial efficacy of a vaccine regimen with an estimated vaccine efficacy (VE) of 31% for protecting low-risk Thai volunteers against acquisition of HIV-1. The impact of vaccine-induced immune responses can be investigated through sieve analysis of HIV-1 breakthrough infections (infected vaccine and placebo recipients). A V1/V2-targeted comparison of the genomes of HIV-1 breakthrough viruses identified two V2 amino acid sites that differed between the vaccine and placebo groups. Here we extended the V1/V2 analysis to the entire HIV-1 genome using an array of methods based on individual sites, k-mers and genes/proteins. We identified 56 amino acid sites or "signatures" and 119 k-mers that differed between the vaccine and placebo groups. Of those, 19 sites and 38 k-mers were located in the regions comprising the RV144 vaccine (Env-gp120, Gag, and Pro). The nine signature sites in Env-gp120 were significantly enriched for known antibody-associated sites (p = 0.0021). In particular, site 317 in the third variable loop (V3) overlapped with a hotspot of antibody recognition, and sites 369 and 424 were linked to CD4 binding site neutralization. The identified signature sites significantly covaried with other sites across the genome (mean = 32.1) more than did non-signature sites (mean = 0.9) (p < 0.0001), suggesting functional and/or structural relevance of the signature sites. Since signature sites were not preferentially restricted to the vaccine immunogens and because most of the associations were insignificant following correction for multiple testing, we predict that few of the genetic differences are strongly linked to the RV144 vaccine-induced immune pressure. In addition to presenting results of the first complete-genome analysis of the breakthrough infections in the RV144 trial, this work describes a set of statistical methods and tools applicable to analysis of breakthrough infection genomes in general vaccine efficacy trials for diverse pathogens.

Dendritic cells restore CD8+ T cell reactivity to autologous HIV-1.

UNLABELLED: Recall T cell responses to HIV-1 antigens are used as a surrogate for endogenous cellular immune responses generated during infection. Current methods of identifying antigen-specific T cell reactivity in HIV-1 infection use bulk peripheral blood mononuclear cells (PBMC) yet ignore professional antigen-presenting cells (APC) that could reveal otherwise hidden responses. In the present study, peptides representing autologous variants of major histocompatibility complex (MHC) class I-restricted epitopes from HIV-1 Gag and Env were used as antigens in gamma interferon (IFN-γ) enzyme-linked immunosorbent spot (ELISpot) and polyfunctional cytokine assays. Here we show that dendritic cells (DC) enhanced T cell reactivity at all stages of disease progression but specifically restored T cell reactivity after combination antiretroviral therapy (cART) to early infection levels. Type 1 cytokine secretion was also enhanced by DC and was most apparent late post-cART. We additionally show that DC reveal polyfunctional T cell responses after many years of treatment, when potential immunotherapies would be implemented. These data underscore the potential efficacy of DC immunotherapy that aims to awaken a dormant, autologous, HIV-1-specific CD8+ T cell response. IMPORTANCE: Assessment of endogenous HIV-1-specific T cell responses is critical for generating immunotherapies for subjects on cART. Current assays ignore the ability of dendritic cells to reveal these responses and may therefore underestimate the breadth and magnitude of T cell reactivity. As DC do not prime new responses in these assays, it can be assumed that the observed responses are not detected without appropriate stimulation. This is important because dogma states that HIV-1 mutates to evade host recognition and that CD8+ cytotoxic T lymphocyte (CTL) failure is due to the inability of T cells to recognize the autologous virus. The results presented here indicate that responses to autologous virus are generated during infection but may need additional stimulation to be effective. Detecting the breadth and magnitude of HIV-1-specific T cell reactivity generated in vivo is of the utmost importance for generating effective DC immunotherapies.

Analysis of HLA A*02 association with vaccine efficacy in the RV144 HIV-1 vaccine trial.

UNLABELLED: The RV144 HIV-1 vaccine trial demonstrated partial efficacy of 31% against HIV-1 infection. Studies into possible correlates of protection found that antibodies specific to the V1 and V2 (V1/V2) region of envelope correlated inversely with infection risk and that viruses isolated from trial participants contained genetic signatures of vaccine-induced pressure in the V1/V2 region. We explored the hypothesis that the genetic signatures in V1 and V2 could be partly attributed to selection by vaccine-primed T cells. We performed a T-cell-based sieve analysis of breakthrough viruses in the RV144 trial and found evidence of predicted HLA binding escape that was greater in vaccine versus placebo recipients. The predicted escape depended on class I HLA A*02- and A*11-restricted epitopes in the MN strain rgp120 vaccine immunogen. Though we hypothesized that this was indicative of postacquisition selection pressure, we also found that vaccine efficacy (VE) was greater in A*02-positive (A*02(+)) participants than in A*02(-) participants (VE = 54% versus 3%, P = 0.05). Vaccine efficacy against viruses with a lysine residue at site 169, important to antibody binding and implicated in vaccine-induced immune pressure, was also greater in A*02(+) participants (VE = 74% versus 15%, P = 0.02). Additionally, a reanalysis of vaccine-induced immune responses that focused on those that were shown to correlate with infection risk suggested that the humoral responses may have differed in A*02(+) participants. These exploratory and hypothesis-generating analyses indicate there may be an association between a class I HLA allele and vaccine efficacy, highlighting the importance of considering HLA alleles and host immune genetics in HIV vaccine trials. IMPORTANCE: The RV144 trial was the first to show efficacy against HIV-1 infection. Subsequently, much effort has been directed toward understanding the mechanisms of protection. Here, we conducted a T-cell-based sieve analysis, which compared the genetic sequences of viruses isolated from infected vaccine and placebo recipients. Though we hypothesized that the observed sieve effect indicated postacquisition T-cell selection, we also found that vaccine efficacy was greater for participants who expressed HLA A*02, an allele implicated in the sieve analysis. Though HLA alleles have been associated with disease progression and viral load in HIV-1 infection, these data are the first to suggest the association of a class I HLA allele and vaccine efficacy. While these statistical analyses do not provide mechanistic evidence of protection in RV144, they generate testable hypotheses for the HIV vaccine community and they highlight the importance of assessing the impact of host immune genetics in vaccine-induced immunity and protection. (This study has been registered at ClinicalTrials.gov under registration no. NCT00223080.).

HIV-1 conserved-element vaccines: relationship between sequence conservation and replicative capacity.

To overcome the problem of HIV-1 variability, candidate vaccine antigens have been designed to be composed of conserved elements of the HIV-1 proteome. Such candidate vaccines could be improved with a better understanding of both HIV-1 evolutionary constraints and the fitness cost of specific mutations. We evaluated the in vitro fitness cost of 23 mutations engineered in the HIV-1 subtype B Gag-p24 Center-of-Tree (COT) protein through fitness competition assays. While some mutations at conserved sites exacted a high fitness cost, as expected under the assumption that the most conserved residue confers the highest fitness, there was no overall strong relationship between sequence conservation and replicative capacity. By comparing sites that have evolved since the beginning of the epidemic to those that have remain unchanged, we found that sites that have evolved over time were more likely to correspond to HLA-associated sites and that their mutation had limited fitness costs. Our data showed no transcendent link between high conservation and high fitness cost, indicating that merely focusing on conserved segments of HIV-1 would not be sufficient for a successful vaccine strategy. Nonetheless, a subset of sites exacted a high fitness cost upon mutation--these sites have been under selective pressure to change since the beginning of the epidemic but have proved virtually nonmutable and could constitute preferred targets for vaccine design.

Indel and Carryforward Correction (ICC): a new analysis approach for processing 454 pyrosequencing data.

MOTIVATION: Pyrosequencing technology provides an important new approach to more extensively characterize diverse sequence populations and detect low frequency variants. However, the promise of this technology has been difficult to realize, as careful correction of sequencing errors is crucial to distinguish rare variants (∼1%) in an infected host with high sensitivity and specificity. RESULTS: We developed a new approach, referred to as Indel and Carryforward Correction (ICC), to cluster sequences without substitutions and locally correct only indel and carryforward sequencing errors within clusters to ensure that no rare variants are lost. ICC performs sequence clustering in the order of (i) homopolymer indel patterns only, (ii) indel patterns only and (iii) carryforward errors only, without the requirement of a distance cutoff value. Overall, ICC removed 93-95% of sequencing errors found in control datasets. On pyrosequencing data from a PCR fragment derived from 15 HIV-1 plasmid clones mixed at various frequencies as low as 0.1%, ICC achieved the highest sensitivity and similar specificity compared with other commonly used error correction and variant calling algorithms. AVAILABILITY AND IMPLEMENTATION: Source code is freely available for download at http://indra.mullins.microbiol.washington.edu/ICC. It is implemented in Perl and supported on Linux, Mac OS X and MS Windows.

Functional and antigenic landscape of the Nipah virus receptor binding protein.

Nipah virus recurrently spills over to humans, causing fatal infections. The viral receptor-binding protein (RBP or G) attaches to host receptors and is a major target of neutralizing antibodies. Here we use deep mutational scanning to measure how all amino-acid mutations to the RBP affect cell entry, receptor binding, and escape from neutralizing antibodies. We identify functionally constrained regions of the RBP, including sites involved in oligomerization, along with mutations that differentially modulate RBP binding to its two ephrin receptors. We map escape mutations for six anti-RBP antibodies, and find that few antigenic mutations are present in natural Nipah strains. Our findings offer insights into the potential for functional and antigenic evolution of the RBP that can inform the development of antibody therapies and vaccines.

Diverse DNA virus genomes identified in fecal samples of Mexican free-tailed bats (Tadarida brasiliensis) captured in Chiricahua Mountains of southeast Arizona (USA).

Bats (order Chiroptera) are some of the most abundant mammals on earth and their species ecology strongly influences zoonotic potential. While substantial research has been conducted on bat-associated viruses, particularly on those that can cause disease in humans and/or livestock, globally, limited research has focused on endemic bats in the USA. The southwest region of the US is of particular interest because of its high diversity of bat species. We identified 39 single-stranded DNA virus genomes in the feces of Mexican free-tailed bats (Tadarida brasiliensis) sampled in the Rucker Canyon (Chiricahua Mountains) of southeast Arizona (USA). Twenty-eight of these belong to the virus families Circoviridae (n = 6), Genomoviridae (n = 17), and Microviridae (n = 5). Eleven viruses cluster with other unclassified cressdnaviruses. Most of the viruses identified represent new species. Further research on identification of novel bat-associated cressdnaviruses and microviruses is needed to provide greater insights regarding their co-evolution and ecology relative to bats.

Using archived and biocollection samples towards deciphering the DNA virus diversity associated with rodent species in the families cricetidae and heteromyidae.

Rodentia is the most speciose order of mammals, and they are known to harbor a wide range of viruses. Although there has been significant research on zoonotic viruses in rodents, research on the diversity of other viruses has been limited, especially for rodents in the families Cricetidae and Heteromyidae. In fecal and liver samples of nine species of rodents, we identify 346 distinct circular DNA viral genomes. Of these, a large portion are circular, single-stranded DNA viruses in the families Anelloviridae (n = 3), Circoviridae (n = 5), Genomoviridae (n = 7), Microviridae (n = 297), Naryaviridae (n = 4), Vilyaviridae (n = 15) and in the phylum Cressdnaviricota (n = 13) that cannot be assigned established families. We also identified two large bacteriophages of 36 and 50 kb that are part of the class Caudoviricetes. Some of these viruses are clearly those that infect rodents, however, most of these likely infect various organisms associated with rodents, their environment or their diet.

Coevolutionary Analysis Implicates Toll-Like Receptor 9 in Papillomavirus Restriction.

Upon infection, DNA viruses can be sensed by pattern recognition receptors (PRRs), leading to the activation of type I and III interferons to block infection. Therefore, viruses must inhibit these signaling pathways, avoid being detected, or both. Papillomavirus virions are trafficked from early endosomes to the Golgi apparatus and wait for the onset of mitosis to complete nuclear entry. This unique subcellular trafficking strategy avoids detection by cytoplasmic PRRs, a property that may contribute to the establishment of infection. However, as the capsid uncoats within acidic endosomal compartments, the viral DNA may be exposed to detection by Toll-like receptor 9 (TLR9). In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we showed that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts. Finally, we demonstrated that the cancer-associated human papillomaviruses show a reduction in CpG dinucleotides within a TLR9 recognition complex. IMPORTANCE Viruses must avoid detection by the innate immune system. In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we demonstrated that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts.

Genomic epidemiology identifies emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States.

As of January of 2021, the highly transmissible B.1.1.7 variant of SARS-CoV-2, which was first identified in the United Kingdom (U.K.), has gained a strong foothold across the world. Because of the sudden and rapid rise of B.1.1.7, we investigated the prevalence and growth dynamics of this variant in the United States (U.S.), tracking it back to its early emergence and onward local transmission. We found that the RT-qPCR testing anomaly of S gene target failure (SGTF), first observed in the U.K., was a reliable proxy for B.1.1.7 detection. We sequenced 212 B.1.1.7 SARS-CoV-2 genomes collected from testing facilities in the U.S. from December 2020 to January 2021. We found that while the fraction of B.1.1.7 among SGTF samples varied by state, detection of the variant increased at a logistic rate similar to those observed elsewhere, with a doubling rate of a little over a week and an increased transmission rate of 35-45%. By performing time-aware Bayesian phylodynamic analyses, we revealed several independent introductions of B.1.1.7 into the U.S. as early as late November 2020, with onward community transmission enabling the variant to spread to at least 30 states as of January 2021. Our study shows that the U.S. is on a similar trajectory as other countries where B.1.1.7 rapidly became the dominant SARS-CoV-2 variant, requiring immediate and decisive action to minimize COVID-19 morbidity and mortality.

The emergence of SARS-CoV-2 in Europe and North America.

Accurate understanding of the global spread of emerging viruses is critical for public health responses and for anticipating and preventing future outbreaks. Here we elucidate when, where, and how the earliest sustained severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission networks became established in Europe and North America. Our results suggest that rapid early interventions successfully prevented early introductions of the virus from taking hold in Germany and the United States. Other, later introductions of the virus from China to both Italy and Washington state, United States, founded the earliest sustained European and North America transmission networks. Our analyses demonstrate the effectiveness of public health measures in preventing onward transmission and show that intensive testing and contact tracing could have prevented SARS-CoV-2 outbreaks from becoming established in these regions.

The emergence of SARS-CoV-2 in Europe and the US.

Accurate understanding of the global spread of emerging viruses is critically important for public health response and for anticipating and preventing future outbreaks. Here, we elucidate when, where and how the earliest sustained SARS-CoV-2 transmission networks became established in Europe and the United States (US). Our results refute prior findings erroneously linking cases in January 2020 with outbreaks that occurred weeks later. Instead, rapid interventions successfully prevented onward transmission of those early cases in Germany and Washington State. Other, later introductions of the virus from China to both Italy and Washington State founded the earliest sustained European and US transmission networks. Our analyses reveal an extended period of missed opportunity when intensive testing and contact tracing could have prevented SARS-CoV-2 from becoming established in the US and Europe.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity.

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2's unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

An Early Pandemic Analysis of SARS-CoV-2 Population Structure and Dynamics in Arizona.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.IMPORTANCE As the COVID-19 pandemic swept across the United States, there was great differential impact on local and regional communities. One of the earliest and hardest hit regions was in New York, while at the same time Arizona (for example) had low incidence. That situation has changed dramatically, with Arizona now having the highest rate of disease increase in the country. Understanding the roots of the pandemic during the initial months is essential as the pandemic continues and reaches new heights. Genomic analysis and phylogenetic modeling of SARS-COV-2 in Arizona can help to reconstruct population composition and predict the earliest undetected introductions. This foundational work represents the basis for future analysis and understanding as the pandemic continues.