H19 influenza A virus exhibits species-specific MHC class II receptor usage.

Avian influenza A virus (IAV) surveillance in Northern California, USA, revealed unique IAV hemagglutinin (HA) genome sequences in cloacal swabs from lesser scaups. We found two closely related HA sequences in the same duck species in 2010 and 2013. Phylogenetic analyses suggest that both sequences belong to the recently discovered H19 subtype, which thus far has remained uncharacterized. We demonstrate that H19 does not bind the canonical IAV receptor sialic acid (Sia). Instead, H19 binds to the major histocompatibility complex class II (MHC class II), which facilitates viral entry. Unlike the broad MHC class II specificity of H17 and H18 from bat IAV, H19 exhibits a species-specific MHC class II usage that suggests a limited host range and zoonotic potential. Using cell lines overexpressing MHC class II, we rescued recombinant H19 IAV. We solved the H19 crystal structure and identified residues within the putative Sia receptor binding site (RBS) that impede Sia-dependent entry.

Lack of selective resistance of influenza A virus in presence of host-targeted antiviral, UV-4B.

Development of antiviral drug resistance is a continuous concern for viruses with high mutation rates such as influenza. The use of antiviral drugs targeting host proteins required for viral replication is less likely to result in the selection of resistant viruses than treating with direct-acting antivirals. The iminosugar UV-4B is a host-targeted glucomimetic that inhibits endoplasmic reticulum α-glucosidase I and II enzymes resulting in improper glycosylation and misfolding of viral glycoproteins. UV-4B has broad-spectrum antiviral activity against diverse viruses including dengue and influenza. To examine the ability of influenza virus to develop resistance against UV-4B, mouse-adapted influenza virus was passaged in mice in the presence or absence of UV-4B and virus isolated from lungs was used to infect the next cohort of mice, for five successive passages. Deep sequencing was performed to identify changes in the viral genome during passaging in the presence or absence of UV-4B. Relatively few minor variants were identified within each virus and the ratio of nonsynonymous to synonymous (dN/dS) substitutions of minor variants confirmed no apparent positive selection following sustained exposure to UV-4B. Three substitutions (one synonymous in PB2, one nonsynonymous in M and PA each) were specifically enriched (>3%) in UV-4B-treated groups at passage five. Recombinant viruses containing each individual or combinations of these nonsynonymous mutations remained sensitive to UV-4B treatment in mice. Overall, these data provide evidence that there is a high genetic barrier to the generation and selection of escape mutants following exposure to host-targeted iminosugar antivirals.

Author Correction: Genomic analysis of a large set of currently-and historically-important human adenovirus pathogens.

The original publication of this article [1] was missing the below author that made contributions to the research and the published article.

Avian Influenza Viruses in Wild Birds: Virus Evolution in a Multihost Ecosystem.

Wild ducks and gulls are the major reservoirs for avian influenza A viruses (AIVs). The mechanisms that drive AIV evolution are complex at sites where various duck and gull species from multiple flyways breed, winter, or stage. The Republic of Georgia is located at the intersection of three migratory flyways: the Central Asian flyway, the East Africa/West Asia flyway, and the Black Sea/Mediterranean flyway. For six complete study years (2010 to 2016), we collected AIV samples from various duck and gull species that breed, migrate, and overwinter in Georgia. We found a substantial subtype diversity of viruses that varied in prevalence from year to year. Low-pathogenic AIV (LPAIV) subtypes included H1N1, H2N3, H2N5, H2N7, H3N8, H4N2, H6N2, H7N3, H7N7, H9N1, H9N3, H10N4, H10N7, H11N1, H13N2, H13N6, H13N8, and H16N3, and two highly pathogenic AIVs (HPAIVs) belonging to clade 2.3.4.4, H5N5 and H5N8, were found. Whole-genome phylogenetic trees showed significant host species lineage restriction for nearly all gene segments and significant differences in observed reassortment rates, as defined by quantification of phylogenetic incongruence, and in nucleotide sequence diversity for LPAIVs among different host species. Hemagglutinin clade 2.3.4.4 H5N8 viruses, which circulated in Eurasia during 2014 and 2015, did not reassort, but analysis after their subsequent dissemination during 2016 and 2017 revealed reassortment in all gene segments except NP and NS. Some virus lineages appeared to be unrelated to AIVs in wild bird populations in other regions, with maintenance of local AIVs in Georgia, whereas other lineages showed considerable genetic interrelationships with viruses circulating in other parts of Eurasia and Africa, despite relative undersampling in the area.IMPORTANCE Waterbirds (e.g., gulls and ducks) are natural reservoirs of avian influenza viruses (AIVs) and have been shown to mediate the dispersal of AIVs at intercontinental scales during seasonal migration. The segmented genome of influenza viruses enables viral RNA from different lineages to mix or reassort when two viruses infect the same host. Such reassortant viruses have been identified in most major human influenza pandemics and several poultry outbreaks. Despite their importance, we have only recently begun to understand AIV evolution and reassortment in their natural host reservoirs. This comprehensive study illustrates AIV evolutionary dynamics within a multihost ecosystem at a stopover site where three major migratory flyways intersect. Our analysis of this ecosystem over a 6-year period provides a snapshot of how these viruses are linked to global AIV populations. Understanding the evolution of AIVs in the natural host is imperative to mitigating both the risk of incursion into domestic poultry and the potential risk to mammalian hosts, including humans.

A draft genome sequence for the Ixodes scapularis cell line, ISE6.

Background: The tick cell line ISE6, derived from Ixodes scapularis, is commonly used for amplification and detection of arboviruses in environmental or clinical samples. Methods: To assist with sequence-based assays, we sequenced the ISE6 genome with single-molecule, long-read technology. Results: The draft assembly appears near complete based on gene content analysis, though it appears to lack some instances of repeats in this highly repetitive genome. The assembly appears to have separated the haplotypes at many loci. DNA short read pairs, used for validation only, mapped to the cell line assembly at a higher rate than they mapped to the Ixodes scapularis reference genome sequence. Conclusions: The assembly could be useful for filtering host genome sequence from sequence data obtained from cells infected with pathogens.

Large-Scale Complete-Genome Sequencing and Phylodynamic Analysis of Eastern Equine Encephalitis Virus Reveals Source-Sink Transmission Dynamics in the United States.

Eastern equine encephalitis virus (EEEV) has a high case-fatality rate in horses and humans, and Florida has been hypothesized to be the source of EEEV epidemics for the northeastern United States. To test this hypothesis, we sequenced complete genomes of 433 EEEV strains collected within the United States from 1934 to 2014. Phylogenetic analysis suggested EEEV evolves relatively slowly and that transmission is enzootic in Florida, characterized by higher genetic diversity and long-term local persistence. In contrast, EEEV strains in New York and Massachusetts were characterized by lower genetic diversity, multiple introductions, and shorter local persistence. Our phylogeographic analysis supported a source-sink model in which Florida is the major source of EEEV compared to the other localities sampled. In sum, this study revealed the complex epidemiological dynamics of EEEV in different geographic regions in the United States and provided general insights into the evolution and transmission of other avian mosquito-borne viruses in this region.IMPORTANCE Eastern equine encephalitis virus (EEEV) infections are severe in horses and humans on the east coast of the United States with a >90% mortality rate in horses, an ∼33% mortality rate in humans, and significant brain damage in most human survivors. However, little is known about the evolutionary characteristics of EEEV due to the lack of genome sequences. By generating large collection of publicly available complete genome sequences, this study comprehensively determined the evolution of the virus, described the epidemiological dynamics of EEEV in different states in the United States, and identified Florida as one of the major sources. These results may have important implications for the control and prevention of other mosquito-borne viruses in the Americas.

Biological Sexing of a 4000-Year-Old Egyptian Mummy Head to Assess the Potential of Nuclear DNA Recovery from the Most Damaged and Limited Forensic Specimens.

High throughput sequencing (HTS) has been used for a number of years in the field of paleogenomics to facilitate the recovery of small DNA fragments from ancient specimens. Recently, these techniques have also been applied in forensics, where they have been used for the recovery of mitochondrial DNA sequences from samples where traditional PCR-based assays fail because of the very short length of endogenous DNA molecules. Here, we describe the biological sexing of a ~4000-year-old Egyptian mummy using shotgun sequencing and two established methods of biological sex determination (RX and RY), by way of mitochondrial genome analysis as a means of sequence data authentication. This particular case of historical interest increases the potential utility of HTS techniques for forensic purposes by demonstrating that data from the more discriminatory nuclear genome can be recovered from the most damaged specimens, even in cases where mitochondrial DNA cannot be recovered with current PCR-based forensic technologies. Although additional work remains to be done before nuclear DNA recovered via these methods can be used routinely in operational casework for individual identification purposes, these results indicate substantial promise for the retrieval of probative individually identifying DNA data from the most limited and degraded forensic specimens.

Poor Immunogenicity, Not Vaccine Strain Egg Adaptation, May Explain the Low H3N2 Influenza Vaccine Effectiveness in 2012-2013.

Background: Influenza vaccination aims to prevent infection by influenza virus and reduce associated morbidity and mortality; however, vaccine effectiveness (VE) can be modest, especially for subtype A(H3N2). Low VE has been attributed to mismatches between the vaccine and circulating influenza strains and to the vaccine's elicitation of protective immunity in only a subset of the population. The low H3N2 VE in the 2012-2013 season was attributed to egg-adaptive mutations that created antigenic mismatch between the actual vaccine strain (IVR-165) and both the intended vaccine strain (A/Victoria/361/2011) and the predominant circulating strains (clades 3C.2 and 3C.3). Methods: We investigated the basis of low VE in 2012-2013 by determining whether vaccinated and unvaccinated individuals were infected by different viral strains and by assessing the serologic responses to IVR-165, A/Victoria/361/2011, and 3C.2 and 3C.3 strains in an adult cohort before and after vaccination. Results: We found no significant genetic differences between the strains that infected vaccinated and unvaccinated individuals. Vaccination increased titers to A/Victoria/361/2011 and 3C.2 and 3C.3 representative strains as much as to IVR-165. These results are consistent with the hypothesis that vaccination boosted cross-reactive immune responses instead of specific responses against unique vaccine epitopes. Only approximately one-third of the cohort achieved a ≥4-fold increase in titer. Conclusions: In contrast to analyses based on ferret studies, low H3N2 VE in 2012-2013 in adults does not appear to be due to egg adaptation of the vaccine strain. Instead, low VE might have been caused by low vaccine immunogenicity in a subset of the population.

Genomic analysis of a large set of currently-and historically-important human adenovirus pathogens.

Human adenoviruses (HAdVs) are uniquely important "model organisms" as they have been used to elucidate fundamental biological processes, are recognized as complex pathogens, and are used as remedies for human health. As pathogens, HAdVs may effect asymptomatic or mild and severe symptomatic disease upon their infection of respiratory, ocular, gastrointestinal, and genitourinary systems. High-resolution genomic data have enhanced the understanding of HAdV epidemiology, with recombination recognized as an important and major pathway in the molecular evolution and genesis of emergent HAdV pathogens. To support this view and to actualize an algorithm for identifying, characterizing, and typing novel HAdVs, we determined the DNA sequence of 95 isolates from archives containing historically important pathogens and collections housing currently circulating strains to be sequenced. Of the 85 samples that were completely sequenced, 18 novel recombinants within species HAdV-B and D were identified. Two HAdV-D genomes were found to contain novel penton base and fiber genes with significant divergence from known molecular types. In this data set, we found additional isolates of HAdV-D53 and HAdV-D58, two novel genotypes recognized recently using genomics. This supports the thesis that novel HAdV genotypes are not limited to "one-time" appearances of the prototype but are of importance in HAdV epidemiology. These data underscore the significance of lateral genomic transfer in HAdV evolution and reinforce the potential public health impact of novel genotypes of HAdVs emerging in the population.

Genetic stability of foot-and-mouth disease virus during long-term infections in natural hosts.

Foot-and-mouth disease (FMD) is a severe infection caused by a picornavirus that affects livestock and wildlife. Persistence in ruminants is a well-documented feature of Foot-and-mouth disease virus (FMDV) pathogenesis and a major concern for disease control. Persistently infected animals harbor virus for extended periods, providing a unique opportunity to study within-host virus evolution. This study investigated the genetic dynamics of FMDV during persistent infections of naturally infected Asian buffalo. Using next-generation sequencing (NGS) we obtained 21 near complete FMDV genome sequences from 12 sub-clinically infected buffalo over a period of one year. Four animals yielded only one virus isolate and one yielded two isolates of different serotype suggesting a serial infection. Seven persistently infected animals yielded more than one virus of the same serotype showing a long-term intra-host viral genetic divergence at the consensus level of less than 2.5%. Quasi-species analysis showed few nucleotide variants and non-synonymous substitutions of progeny virus despite intra-host persistence of up to 152 days. Phylogenetic analyses of serotype Asia-1 VP1 sequences clustered all viruses from persistent animals with Group VII viruses circulating in Pakistan in 2011, but distinct from those circulating on 2008-2009. Furthermore, signature amino acid (aa) substitutions were found in the antigenically relevant VP1 of persistent viruses compared with viruses from 2008-2009. Intra-host purifying selective pressure was observed, with few codons in structural proteins undergoing positive selection. However, FMD persistent viruses did not show a clear pattern of antigenic selection. Our findings provide insight into the evolutionary dynamics of FMDV populations within naturally occurring subclinical and persistent infections that may have implications to vaccination strategies in the region.

Analysis of the Aedes albopictus C6/36 genome provides insight into cell line utility for viral propagation.

Background: The 50-year-old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome. Results: The C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells were mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of nonhost reads. Conclusions: The C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease.

A host subtraction database for virus discovery in human cell line sequencing data.

The human cell lines HepG2, HuH-7, and Jurkat are commonly used for amplification of the RNA viruses present in environmental samples. To assist with assays by RNAseq, we sequenced these cell lines and developed a subtraction database that contains sequences expected in sequence data from uninfected cells. RNAseq data from cell lines infected with Sendai virus were analyzed to test host subtraction. The process of mapping RNAseq reads to our subtraction database vastly reduced the number non-viral reads in the dataset to allow for efficient secondary analyses.

Multiplex Reverse Transcription-PCR for Simultaneous Surveillance of Influenza A and B Viruses.

Influenza A and B viruses are the causative agents of annual influenza epidemics that can be severe, and influenza A viruses intermittently cause pandemics. Sequence information from influenza virus genomes is instrumental in determining mechanisms underpinning antigenic evolution and antiviral resistance. However, due to sequence diversity and the dynamics of influenza virus evolution, rapid and high-throughput sequencing of influenza viruses remains a challenge. We developed a single-reaction influenza A/B virus (FluA/B) multiplex reverse transcription-PCR (RT-PCR) method that amplifies the most critical genomic segments (hemagglutinin [HA], neuraminidase [NA], and matrix [M]) of seasonal influenza A and B viruses for next-generation sequencing, regardless of viral type, subtype, or lineage. Herein, we demonstrate that the strategy is highly sensitive and robust. The strategy was validated on thousands of seasonal influenza A and B virus-positive specimens using multiple next-generation sequencing platforms.

Cloning, Assembly, and Modification of the Primary Human Cytomegalovirus Isolate Toledo by Yeast-Based Transformation-Associated Recombination.

Genetic engineering of cytomegalovirus (CMV) currently relies on generating a bacterial artificial chromosome (BAC) by introducing a bacterial origin of replication into the viral genome using in vivo recombination in virally infected tissue culture cells. However, this process is inefficient, results in adaptive mutations, and involves deletion of viral genes to avoid oversized genomes when inserting the BAC cassette. Moreover, BAC technology does not permit the simultaneous manipulation of multiple genome loci and cannot be used to construct synthetic genomes. To overcome these limitations, we adapted synthetic biology tools to clone CMV genomes in Saccharomyces cerevisiae. Using an early passage of the human CMV isolate Toledo, we first applied transformation-associated recombination (TAR) to clone 16 overlapping fragments covering the entire Toledo genome in Saccharomyces cerevisiae. Then, we assembled these fragments by TAR in a stepwise process until the entire genome was reconstituted in yeast. Since next-generation sequence analysis revealed that the low-passage-number isolate represented a mixture of parental and fibroblast-adapted genomes, we selectively modified individual DNA fragments of fibroblast-adapted Toledo (Toledo-F) and again used TAR assembly to recreate parental Toledo (Toledo-P). Linear, full-length HCMV genomes were transfected into human fibroblasts to recover virus. Unlike Toledo-F, Toledo-P displayed characteristics of primary isolates, including broad cellular tropism in vitro and the ability to establish latency and reactivation in humanized mice. Our novel strategy thus enables de novo cloning of CMV genomes, more-efficient genome-wide engineering, and the generation of viral genomes that are partially or completely derived from synthetic DNA. IMPORTANCE The genomes of large DNA viruses, such as human cytomegalovirus (HCMV), are difficult to manipulate using current genetic tools, and at this time, it is not possible to obtain, molecular clones of CMV without extensive tissue culture. To overcome these limitations, we used synthetic biology tools to capture genomic fragments from viral DNA and assemble full-length genomes in yeast. Using an early passage of the HCMV isolate Toledo containing a mixture of wild-type and tissue culture-adapted virus. we directly cloned the majority sequence and recreated the minority sequence by simultaneous modification of multiple genomic regions. Thus, our novel approach provides a paradigm to not only efficiently engineer HCMV and other large DNA viruses on a genome-wide scale but also facilitates the cloning and genetic manipulation of primary isolates and provides a pathway to generating entirely synthetic genomes.

Evolution and spread of Venezuelan equine encephalitis complex alphavirus in the Americas.

Venezuelan equine encephalitis (VEE) complex alphaviruses are important re-emerging arboviruses that cause life-threatening disease in equids during epizootics as well as spillover human infections. We conducted a comprehensive analysis of VEE complex alphaviruses by sequencing the genomes of 94 strains and performing phylogenetic analyses of 130 isolates using complete open reading frames for the nonstructural and structural polyproteins. Our analyses confirmed purifying selection as a major mechanism influencing the evolution of these viruses as well as a confounding factor in molecular clock dating of ancestors. Times to most recent common ancestors (tMRCAs) could be robustly estimated only for the more recently diverged subtypes; the tMRCA of the ID/IAB/IC/II and IE clades of VEE virus (VEEV) were estimated at ca. 149-973 years ago. Evolution of the IE subtype has been characterized by a significant evolutionary shift from the rest of the VEEV complex, with an increase in structural protein substitutions that are unique to this group, possibly reflecting adaptation to its unique enzootic mosquito vector Culex (Melanoconion) taeniopus. Our inferred tree topologies suggest that VEEV is maintained primarily in situ, with only occasional spread to neighboring countries, probably reflecting the limited mobility of rodent hosts and mosquito vectors.

Identification of Dezidougou Virus in a DAK AR 41524 Zika Virus Stock.

We report here the complete genome of a Dezidougou virus (DEZV) isolated from a passaged culture of the Zika virus strain DAK AR 41524. The consensus DEZV sequence we recovered shows 99% nucleotide similarity using BLASTN to a previously reported DEZV (accession no. JQ675604.1). The current sequence has additional repeat regions as well as a deleted repeat region, which we confirmed by Sanger sequencing, that were not present in the originally published sequence, JQ675604.1.

Initial genome sequencing of the sugarcane CP 96-1252 complex hybrid.

The CP 96-1252 cultivar of sugarcane is a complex hybrid of commercial importance. DNA was extracted from lab-grown leaf tissue and sequenced. The raw Illumina DNA sequencing results provide 101 Gbp of genome sequence reads. The dataset is available from https://www.ncbi.nlm.nih.gov/bioproject/PRJNA345486/.

The Ebola virus VP35 protein binds viral immunostimulatory and host RNAs identified through deep sequencing.

Ebola virus and Marburg virus are members of the Filovirdae family and causative agents of hemorrhagic fever with high fatality rates in humans. Filovirus virulence is partially attributed to the VP35 protein, a well-characterized inhibitor of the RIG-I-like receptor pathway that triggers the antiviral interferon (IFN) response. Prior work demonstrates the ability of VP35 to block potent RIG-I activators, such as Sendai virus (SeV), and this IFN-antagonist activity is directly correlated with its ability to bind RNA. Several structural studies demonstrate that VP35 binds short synthetic dsRNAs; yet, there are no data that identify viral immunostimulatory RNAs (isRNA) or host RNAs bound to VP35 in cells. Utilizing a SeV infection model, we demonstrate that both viral isRNA and host RNAs are bound to Ebola and Marburg VP35s in cells. By deep sequencing the purified VP35-bound RNA, we identified the SeV copy-back defective interfering (DI) RNA, previously identified as a robust RIG-I activator, as the isRNA bound by multiple filovirus VP35 proteins, including the VP35 protein from the West African outbreak strain (Makona EBOV). Moreover, RNAs isolated from a VP35 RNA-binding mutant were not immunostimulatory and did not include the SeV DI RNA. Strikingly, an analysis of host RNAs bound by wild-type, but not mutant, VP35 revealed that select host RNAs are preferentially bound by VP35 in cell culture. Taken together, these data support a model in which VP35 sequesters isRNA in virus-infected cells to avert RIG-I like receptor (RLR) activation.

The effective rate of influenza reassortment is limited during human infection.

We characterise the evolutionary dynamics of influenza infection described by viral sequence data collected from two challenge studies conducted in human hosts. Viral sequence data were collected at regular intervals from infected hosts. Changes in the sequence data observed across time show that the within-host evolution of the virus was driven by the reversion of variants acquired during previous passaging of the virus. Treatment of some patients with oseltamivir on the first day of infection did not lead to the emergence of drug resistance variants in patients. Using an evolutionary model, we inferred the effective rate of reassortment between viral segments, measuring the extent to which randomly chosen viruses within the host exchange genetic material. We find strong evidence that the rate of effective reassortment is low, such that genetic associations between polymorphic loci in different segments are preserved during the course of an infection in a manner not compatible with epistasis. Combining our evidence with that of previous studies we suggest that spatial heterogeneity in the viral population may reduce the extent to which reassortment is observed. Our results do not contradict previous findings of high rates of viral reassortment in vitro and in small animal studies, but indicate that in human hosts the effective rate of reassortment may be substantially more limited.

A Universal Next-Generation Sequencing Protocol To Generate Noninfectious Barcoded cDNA Libraries from High-Containment RNA Viruses.

Several biosafety level 3 and/or 4 (BSL-3/4) pathogens are high-consequence, single-stranded RNA viruses, and their genomes, when introduced into permissive cells, are infectious. Moreover, many of these viruses are select agents (SAs), and their genomes are also considered SAs. For this reason, cDNAs and/or their derivatives must be tested to ensure the absence of infectious virus and/or viral RNA before transfer out of the BSL-3/4 and/or SA laboratory. This tremendously limits the capacity to conduct viral genomic research, particularly the application of next-generation sequencing (NGS). Here, we present a sequence-independent method to rapidly amplify viral genomic RNA while simultaneously abolishing both viral and genomic RNA infectivity across multiple single-stranded positive-sense RNA (ssRNA+) virus families. The process generates barcoded DNA amplicons that range in length from 300 to 1,000 bp, which cannot be used to rescue a virus and are stable to transport at room temperature. Our barcoding approach allows for up to 288 barcoded samples to be pooled into a single library and run across various NGS platforms without potential reconstitution of the viral genome. Our data demonstrate that this approach provides full-length genomic sequence information not only from high-titer virion preparations but it can also recover specific viral sequence from samples with limited starting material in the background of cellular RNA, and it can be used to identify pathogens from unknown samples. In summary, we describe a rapid, universal standard operating procedure that generates high-quality NGS libraries free of infectious virus and infectious viral RNA. IMPORTANCE This report establishes and validates a standard operating procedure (SOP) for select agents (SAs) and other biosafety level 3 and/or 4 (BSL-3/4) RNA viruses to rapidly generate noninfectious, barcoded cDNA amenable for next-generation sequencing (NGS). This eliminates the burden of testing all processed samples derived from high-consequence pathogens prior to transfer from high-containment laboratories to lower-containment facilities for sequencing. Our established protocol can be scaled up for high-throughput sequencing of hundreds of samples simultaneously, which can dramatically reduce the cost and effort required for NGS library construction. NGS data from this SOP can provide complete genome coverage from viral stocks and can also detect virus-specific reads from limited starting material. Our data suggest that the procedure can be implemented and easily validated by institutional biosafety committees across research laboratories.