SwabExpress: An End-to-End Protocol for Extraction-Free COVID-19 Testing.

BACKGROUND: The urgent need for massively scaled clinical testing for SARS-CoV-2, along with global shortages of critical reagents and supplies, has necessitated development of streamlined laboratory testing protocols. Conventional nucleic acid testing for SARS-CoV-2 involves collection of a clinical specimen with a nasopharyngeal swab in transport medium, nucleic acid extraction, and quantitative reverse-transcription PCR (RT-qPCR). As testing has scaled across the world, the global supply chain has buckled, rendering testing reagents and materials scarce. To address shortages, we developed SwabExpress, an end-to-end protocol developed to employ mass produced anterior nares swabs and bypass the requirement for transport media and nucleic acid extraction. METHODS: We evaluated anterior nares swabs, transported dry and eluted in low-TE buffer as a direct-to-RT-qPCR alternative to extraction-dependent viral transport media. We validated our protocol of using heat treatment for viral inactivation and added a proteinase K digestion step to reduce amplification interference. We tested this protocol across archived and prospectively collected swab specimens to fine-tune test performance. RESULTS: After optimization, SwabExpress has a low limit of detection at 2-4 molecules/µL, 100% sensitivity, and 99.4% specificity when compared side by side with a traditional RT-qPCR protocol employing extraction. On real-world specimens, SwabExpress outperforms an automated extraction system while simultaneously reducing cost and hands-on time. CONCLUSION: SwabExpress is a simplified workflow that facilitates scaled testing for COVID-19 without sacrificing test performance. It may serve as a template for the simplification of PCR-based clinical laboratory tests, particularly in times of critical shortages during pandemics.

Comparison of Symptoms and RNA Levels in Children and Adults With SARS-CoV-2 Infection in the Community Setting.

Importance: The association between COVID-19 symptoms and SARS-CoV-2 viral levels in children living in the community is not well understood. Objective: To characterize symptoms of pediatric COVID-19 in the community and analyze the association between symptoms and SARS-CoV-2 RNA levels, as approximated by cycle threshold (Ct) values, in children and adults. Design, Setting, and Participants: This cross-sectional study used a respiratory virus surveillance platform in persons of all ages to detect community COVID-19 cases from March 23 to November 9, 2020. A population-based convenience sample of children younger than 18 years and adults in King County, Washington, who enrolled online for home self-collection of upper respiratory samples for SARS-CoV-2 testing were included. Exposures: Detection of SARS-CoV-2 RNA by reverse transcription-polymerase chain reaction (RT-PCR) from participant-collected samples. Main Outcomes and Measures: RT-PCR-confirmed SARS-CoV-2 infection, with Ct values stratified by age and symptoms. Results: Among 555 SARS-CoV-2-positive participants (mean [SD] age, 33.7 [20.1] years; 320 were female [57.7%]), 47 of 123 children (38.2%) were asymptomatic compared with 31 of 432 adults (7.2%). When symptomatic, fewer symptoms were reported in children compared with adults (mean [SD], 1.6 [2.0] vs 4.5 [3.1]). Symptomatic individuals had lower Ct values (which corresponded to higher viral RNA levels) than asymptomatic individuals (adjusted estimate for children, -3.0; 95% CI, -5.5 to -0.6; P = .02; adjusted estimate for adults, -2.9; 95% CI, -5.2 to -0.6; P = .01). The difference in mean Ct values was neither statistically significant between symptomatic children and symptomatic adults (adjusted estimate, -0.7; 95% CI, -2.2 to 0.9; P = .41) nor between asymptomatic children and asymptomatic adults (adjusted estimate, -0.6; 95% CI, -4.0 to 2.8; P = .74). Conclusions and Relevance: In this community-based cross-sectional study, SARS-CoV-2 RNA levels, as determined by Ct values, were significantly higher in symptomatic individuals than in asymptomatic individuals and no significant age-related differences were found. Further research is needed to understand the role of SARS-CoV-2 RNA levels and viral transmission.

Viral genomes reveal patterns of the SARS-CoV-2 outbreak in Washington State.

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has gravely affected societies around the world. Outbreaks in different parts of the globe have been shaped by repeated introductions of new viral lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State (USA) to characterize how the spread of SARS-CoV-2 in Washington State in early 2020 was shaped by differences in timing of mitigation strategies across counties and by repeated introductions of viral lineages into the state. In addition, we show that the increase in frequency of a potentially more transmissible viral variant (614G) over time can potentially be explained by regional mobility differences and multiple introductions of 614G but not the other variant (614D) into the state. At an individual level, we observed evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we did not find any evidence that the 614G variant affects clinical severity or patient outcomes. Overall, this suggests that with regard to D614G, the behavior of individuals has been more important in shaping the course of the pandemic in Washington State than this variant of the virus.

Effects of weather-related social distancing on city-scale transmission of respiratory viruses: a retrospective cohort study.

BACKGROUND: Unusually high snowfall in western Washington State in February 2019 led to widespread school and workplace closures. We assessed the impact of social distancing caused by this extreme weather event on the transmission of respiratory viruses. METHODS: Residual specimens from patients evaluated for acute respiratory illness at hospitals in the Seattle metropolitan area were screened for a panel of respiratory viruses. Transmission models were fit to each virus to estimate the magnitude reduction in transmission due to weather-related disruptions. Changes in contact rates and care-seeking were informed by data on local traffic volumes and hospital visits. RESULTS: Disruption in contact patterns reduced effective contact rates during the intervention period by 16 to 95%, and cumulative disease incidence through the remainder of the season by 3 to 9%. Incidence reductions were greatest for viruses that were peaking when the disruption occurred and least for viruses in an early epidemic phase. CONCLUSION: High-intensity, short-duration social distancing measures may substantially reduce total incidence in a respiratory virus epidemic if implemented near the epidemic peak. For SARS-CoV-2, this suggests that, even when SARS-CoV-2 spread is out of control, implementing short-term disruptions can prevent COVID-19 deaths.

Evaluating Specimen Quality and Results from a Community-Wide, Home-Based Respiratory Surveillance Study.

While influenza and other respiratory pathogens cause significant morbidity and mortality, the community-based burden of these infections remains incompletely understood. The development of novel methods to detect respiratory infections is essential for mitigating epidemics and developing pandemic-preparedness infrastructure. From October 2019 to March 2020, we conducted a home-based cross-sectional study in the greater Seattle, WA, area, utilizing electronic consent and data collection instruments. Participants received nasal swab collection kits via rapid delivery within 24 hours of self-reporting respiratory symptoms. Samples were returned to the laboratory and were screened for 26 respiratory pathogens and a housekeeping gene. Participant data were recorded via online survey at the time of sample collection and 1 week later. Of the 4,572 consented participants, 4,359 (95.3%) received a home swab kit and 3,648 (83.7%) returned a nasal specimen for respiratory pathogen screening. The 3,638 testable samples had a mean RNase P relative cycle threshold (Crt ) value of 19.0 (SD, 3.4), and 1,232 (33.9%) samples had positive results for one or more pathogens, including 645 (17.7%) influenza-positive specimens. Among the testable samples, the median time between shipment of the home swab kit and completion of laboratory testing was 8.0 days (interquartile range [IQR], 7.0 to 14.0). A single adverse event occurred and did not cause long-term effects or require medical attention. Home-based surveillance using online participant enrollment and specimen self-collection is a safe and feasible method for community-level monitoring of influenza and other respiratory pathogens, which can readily be adapted for use during pandemics.

SwabExpress: An end-to-end protocol for extraction-free COVID-19 testing.

BACKGROUND: The urgent need for massively scaled clinical testing for SARS-CoV-2, along with global shortages of critical reagents and supplies, has necessitated development of streamlined laboratory testing protocols. Conventional nucleic acid testing for SARS-CoV-2 involves collection of a clinical specimen with a nasopharyngeal swab in transport medium, nucleic acid extraction, and quantitative reverse transcription PCR (RT-qPCR) (1). As testing has scaled across the world, the global supply chain has buckled, rendering testing reagents and materials scarce (2). To address shortages, we developed SwabExpress, an end-to-end protocol developed to employ mass produced anterior nares swabs and bypass the requirement for transport media and nucleic acid extraction. METHODS: We evaluated anterior nares swabs, transported dry and eluted in low-TE buffer as a direct-to-RT-qPCR alternative to extraction-dependent viral transport media. We validated our protocol of using heat treatment for viral activation and added a proteinase K digestion step to reduce amplification interference. We tested this protocol across archived and prospectively collected swab specimens to fine-tune test performance. RESULTS: After optimization, SwabExpress has a low limit of detection at 2-4 molecules/uL, 100% sensitivity, and 99.4% specificity when compared side-by-side with a traditional RT-qPCR protocol employing extraction. On real-world specimens, SwabExpress outperforms an automated extraction system while simultaneously reducing cost and hands-on time. CONCLUSION: SwabExpress is a simplified workflow that facilitates scaled testing for COVID-19 without sacrificing test performance. It may serve as a template for the simplification of PCR-based clinical laboratory tests, particularly in times of critical shortages during pandemics.

Comparable specimen collection from both ends of at-home mid-turbinate swabs.

Unsupervised upper respiratory specimen collection is a key factor in the ability to massively scale SARS-CoV-2 testing. But there is concern that unsupervised specimen collection may produce inferior samples. Across two studies that included unsupervised at-home mid-turbinate specimen collection, ~1% of participants used the wrong end of the swab. We found that molecular detection of respiratory pathogens and a human biomarker were comparable between specimens collected from the handle of the swab and those collected correctly. Older participants were more likely to use the swab backwards. Our results suggest that errors made during home-collection of nasal specimens do not preclude molecular detection of pathogens and specialized swabs may be an unnecessary luxury during a pandemic.

Viral genomes reveal patterns of the SARS-CoV-2 outbreak in Washington State.

The rapid spread of SARS-CoV-2 has gravely impacted societies around the world. Outbreaks in different parts of the globe are shaped by repeated introductions of new lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State to characterize how the spread of SARS-CoV-2 in Washington State (USA) was shaped by differences in timing of mitigation strategies across counties, as well as by repeated introductions of viral lineages into the state. Additionally, we show that the increase in frequency of a potentially more transmissible viral variant (614G) over time can potentially be explained by regional mobility differences and multiple introductions of 614G, but not the other variant (614D) into the state. At an individual level, we see evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we do not find any evidence that the 614G variant impacts clinical severity or patient outcomes. Overall, this suggests that at least to date, the behavior of individuals has been more important in shaping the course of the pandemic than changes in the virus.

Cryptic transmission of SARS-CoV-2 in Washington state.

After its emergence in Wuhan, China, in late November or early December 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus rapidly spread globally. Genome sequencing of SARS-CoV-2 allows the reconstruction of its transmission history, although this is contingent on sampling. We analyzed 453 SARS-CoV-2 genomes collected between 20 February and 15 March 2020 from infected patients in Washington state in the United States. We find that most SARS-CoV-2 infections sampled during this time derive from a single introduction in late January or early February 2020, which subsequently spread locally before active community surveillance was implemented.

Cryptic transmission of SARS-CoV-2 in Washington State.

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

Evidence for Limited Early Spread of COVID-19 Within the United States, January-February 2020.

From January 21 through February 23, 2020, public health agencies detected 14 U.S. cases of coronavirus disease 2019 (COVID-19), all related to travel from China (1,2). The first nontravel-related U.S. case was confirmed on February 26 in a California resident who had become ill on February 13 (3). Two days later, on February 28, a second nontravel-related case was confirmed in the state of Washington (4,5). Examination of four lines of evidence provides insight into the timing of introduction and early transmission of SARS-CoV-2, the virus that causes COVID-19, into the United States before the detection of these two cases. First, syndromic surveillance based on emergency department records from counties affected early by the pandemic did not show an increase in visits for COVID-19-like illness before February 28. Second, retrospective SARS-CoV-2 testing of approximately 11,000 respiratory specimens from several U.S. locations beginning January 1 identified no positive results before February 20. Third, analysis of viral RNA sequences from early cases suggested that a single lineage of virus imported directly or indirectly from China began circulating in the United States between January 18 and February 9, followed by several SARS-CoV-2 importations from Europe. Finally, the occurrence of three cases, one in a California resident who died on February 6, a second in another resident of the same county who died February 17, and a third in an unidentified passenger or crew member aboard a Pacific cruise ship that left San Francisco on February 11, confirms cryptic circulation of the virus by early February. These data indicate that sustained, community transmission had begun before detection of the first two nontravel-related U.S. cases, likely resulting from the importation of a single lineage of virus from China in late January or early February, followed by several importations from Europe. The widespread emergence of COVID-19 throughout the United States after February highlights the importance of robust public health systems to respond rapidly to emerging infectious threats.

Annotation of pseudogenic gene segments by massively parallel sequencing of rearranged lymphocyte receptor loci.

BACKGROUND: The adaptive immune system generates a remarkable range of antigen-specific T-cell receptors (TCRs), allowing the recognition of a diverse set of antigens. Most of this diversity is encoded in the complementarity determining region 3 (CDR3) of the ß chain of the αß TCR, which is generated by somatic recombination of noncontiguous variable (V), diversity (D), and joining (J) gene segments. Deletion and non-templated insertion of nucleotides at the D-J and V-DJ junctions further increases diversity. Many of these gene segments are annotated as non-functional owing to defects in their primary sequence, the absence of motifs necessary for rearrangement, or chromosomal locations outside the TCR locus. METHODS: We sought to utilize a novel method, based on high-throughput sequencing of rearranged TCR genes in a large cohort of individuals, to evaluate the use of functional and non-functional alleles. We amplified and sequenced genomic DNA from the peripheral blood of 587 healthy volunteers using a multiplexed polymerase chain reaction assay that targets the variable region of the rearranged TCRß locus, and we determined the presence and the proportion of productive rearrangements for each TCRß V gene segment in each individual. We then used this information to annotate the functional status of TCRß V gene segments in this cohort. RESULTS: For most TCRß V gene segments, our method agrees with previously reported functional annotations. However, we identified novel non-functional alleles for several gene segments, some of which were used exclusively in our cohort to the detriment of reported functional alleles. We also saw that some gene segments reported to have both functional and non-functional alleles consistently behaved in our cohort as either functional or non-functional, suggesting that some reported alleles were not present in the population studied. CONCLUSIONS: In this proof-of-principle study, we used high-throughput sequencing of the TCRß locus of a large cohort of healthy volunteers to evaluate the use of functional and non-functional alleles of individual TCRß V gene segments. With some modifications, our method has the potential to be extended to gene segments in the α, γ, and δ TCR loci, as well as the genes encoding for B-cell receptor chains.

High-throughput pairing of T cell receptor α and ß sequences.

The T cell receptor (TCR) protein is a heterodimer composed of an α chain and a ß chain. TCR genes undergo somatic DNA rearrangements to generate the diversity of T cell binding specificities needed for effective immunity. Recently, high-throughput immunosequencing methods have been developed to profile the TCR α (TCRA) and TCR ß (TCRB) repertoires. However, these methods cannot determine which TCRA and TCRB chains combine to form a specific TCR, which is essential for many functional and therapeutic applications. We describe and validate a method called pairSEQ, which can leverage the diversity of TCR sequences to accurately pair hundreds of thousands of TCRA and TCRB sequences in a single experiment. Our TCR pairing method uses standard laboratory consumables and equipment without the need for single-cell technologies. We show that pairSEQ can be applied to T cells from both blood and solid tissues, such as tumors.

Exome Sequencing of Phenotypic Extremes Identifies CAV2 and TMC6 as Interacting Modifiers of Chronic Pseudomonas aeruginosa Infection in Cystic Fibrosis.

Discovery of rare or low frequency variants in exome or genome data that are associated with complex traits often will require use of very large sample sizes to achieve adequate statistical power. For a fixed sample size, sequencing of individuals sampled from the tails of a phenotype distribution (i.e., extreme phenotypes design) maximizes power and this approach was recently validated empirically with the discovery of variants in DCTN4 that influence the natural history of P. aeruginosa airway infection in persons with cystic fibrosis (CF; MIM219700). The increasing availability of large exome/genome sequence datasets that serve as proxies for population-based controls affords the opportunity to test an alternative, potentially more powerful and generalizable strategy, in which the frequency of rare variants in a single extreme phenotypic group is compared to a control group (i.e., extreme phenotype vs. control population design). As proof-of-principle, we applied this approach to search for variants associated with risk for age-of-onset of chronic P. aeruginosa airway infection among individuals with CF and identified variants in CAV2 and TMC6 that were significantly associated with group status. These results were validated using a large, prospective, longitudinal CF cohort and confirmed a significant association of a variant in CAV2 with increased age-of-onset of P. aeruginosa airway infection (hazard ratio = 0.48, 95% CI=[0.32, 0.88]) and variants in TMC6 with diminished age-of-onset of P. aeruginosa airway infection (HR = 5.4, 95% CI=[2.2, 13.5]) A strong interaction between CAV2 and TMC6 variants was observed (HR=12.1, 95% CI=[3.8, 39]) for children with the deleterious TMC6 variant and without the CAV2 protective variant. Neither gene showed a significant association using an extreme phenotypes design, and conditions for which the power of an extreme phenotype vs. control population design was greater than that for the extreme phenotypes design were explored.

Measures of exposure impact genetic association studies: an example in vitamin K levels and VKORC1.

Studies assessing the impact of gene-environment interactions on common human diseases and traits have been relatively few for many reasons. One often acknowledged reason is that it is difficult to accurately measure the environment or exposure. Indeed, most large-scale epidemiologic studies use questionnaires to assess and measure past and current exposure levels. While questionnaires may be cost-effective, the data may or may not accurately represent the exposure compared with more direct measurements (e.g., self-reported current smoking status versus direct measurement for cotinine levels). Much like phenotyping, the choice in how an exposure is measured may impact downstream tests of genetic association and gene-environment interaction studies. As a case study, we performed tests of association between five common VKORC1 SNPs and two different measurements of vitamin K levels, dietary (n=5,725) and serum (n=348), in the Third National Health and Nutrition Examination Studies (NHANES III). We did not replicate previously reported associations between VKORC1 and vitamin K levels using either measure. Furthermore, the suggestive associations and estimated genetic effect sizes identified in this study differed depending on the vitamin K measurement. This case study of VKORC1 and vitamin K levels serves as a cautionary example of the downstream consequences that the type of exposure measurement choices will have on genetic association and possibly gene-environment studies.

Rare variation facilitates inferences of fine-scale population structure in humans.

Understanding the genetic structure of human populations has important implications for the design and interpretation of disease mapping studies and reconstructing human evolutionary history. To date, inferences of human population structure have primarily been made with common variants. However, recent large-scale resequencing studies have shown an abundance of rare variation in humans, which may be particularly useful for making inferences of fine-scale population structure. To this end, we used an information theory framework and extensive coalescent simulations to rigorously quantify the informativeness of rare and common variation to detect signatures of fine-scale population structure. We show that rare variation affords unique insights into patterns of recent population structure. Furthermore, to empirically assess our theoretical findings, we analyzed high-coverage exome sequences in 6,515 European and African American individuals. As predicted, rare variants are more informative than common polymorphisms in revealing a distinct cluster of European-American individuals, and subsequent analyses demonstrate that these individuals are likely of Ashkenazi Jewish ancestry. Our results provide new insights into the population structure using rare variation, which will be an important factor to account for in rare variant association studies.

Quantifying rare, deleterious variation in 12 human cytochrome P450 drug-metabolism genes in a large-scale exome dataset.

The study of genetic influences on drug response and efficacy ('pharmacogenetics') has existed for over 50 years. Yet, we still lack a complete picture of how genetic variation, both common and rare, affects each individual's responses to medications. Exome sequencing is a promising alternative method for pharmacogenetic discovery as it provides information on both common and rare variation in large numbers of individuals. Using exome data from 2203 AA and 4300 Caucasian individuals through the NHLBI Exome Sequencing Project, we conducted a survey of coding variation within 12 Cytochrome P450 (CYP) genes that are collectively responsible for catalyzing nearly 75% of all known Phase I drug oxidation reactions. In addition to identifying many polymorphisms with known pharmacogenetic effects, we discovered over 730 novel nonsynonymous alleles across the 12 CYP genes of interest. These alleles include many with diverse functional effects such as premature stop codons, aberrant splicesites and mutations at conserved active site residues. Our analysis considering both novel, predicted functional alleles as well as known, actionable CYP alleles reveals that rare, deleterious variation contributes markedly to the overall burden of pharmacogenetic alleles within the populations considered, and that the contribution of rare variation to this burden is over three times greater in AA individuals as compared with Caucasians. While most of these impactful alleles are individually rare, 7.6-11.7% of individuals interrogated in the study carry at least one newly described potentially deleterious alleles in a major drug-metabolizing CYP.