Bedside Back to Bench: Building Bridges between Basic and Clinical Genomic Research.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

Calibration of variant effect predictors on genome-wide data masks heterogeneous performance across genes.

In silico variant effect predictions are available for nearly all missense variants but played a minimal role in clinical variant classification because they were deemed to provide only supporting evidence. Recently, the ClinGen Sequence Variant Interpretation (SVI) Working Group updated recommendations for variant effect prediction use. By analyzing control pathogenic and benign variants across all genes, they were able to compute evidence strength for predictor score intervals with some intervals generating moderate, strong, or even very strong evidence. However, this genome-wide approach could obscure heterogeneous predictor performance in different genes. We quantified the gene-by-gene performance of two top predictors, REVEL and BayesDel, by analyzing control variants in each predictor score interval in 3,668 disease-relevant genes. Approximately 10% of intervals had sufficient control variants for analysis, and ∼70% of these intervals exceeded the maximum number of incorrect predictions implied by the SVI recommendations. These trending discordant intervals arose owing to the divergence of the gene-specific distribution of predictions from the genome-wide distribution, suggesting that gene-specific calibration is needed in many cases. Approximately 22% of ClinVar missense variants of uncertain significance in genes we analyzed (REVEL = 100,629, BayesDel = 71,928) had predictions in trending discordant intervals. Thus, genome-wide calibrations could result in many variants receiving inappropriate evidence strength. To facilitate a review of the SVI's calibrations, we developed a web application enabling visualization of gene-specific predictions and trending concordant and discordant intervals.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions.

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

Local-scale phylodynamics reveal differential community impact of SARS-CoV-2 in a metropolitan US county.

SARS-CoV-2 transmission is largely driven by heterogeneous dynamics at a local scale, leaving local health departments to design interventions with limited information. We analyzed SARS-CoV-2 genomes sampled between February 2020 and March 2022 jointly with epidemiological and cell phone mobility data to investigate fine scale spatiotemporal SARS-CoV-2 transmission dynamics in King County, Washington, a diverse, metropolitan US county. We applied an approximate structured coalescent approach to model transmission within and between North King County and South King County alongside the rate of outside introductions into the county. Our phylodynamic analyses reveal that following stay-at-home orders, the epidemic trajectories of North and South King County began to diverge. We find that South King County consistently had more reported and estimated cases, COVID-19 hospitalizations, and longer persistence of local viral transmission when compared to North King County, where viral importations from outside drove a larger proportion of new cases. Using mobility and demographic data, we also find that South King County experienced a more modest and less sustained reduction in mobility following stay-at-home orders than North King County, while also bearing more socioeconomic inequities that might contribute to a disproportionate burden of SARS-CoV-2 transmission. Overall, our findings suggest a role for local-scale phylodynamics in understanding the heterogeneous transmission landscape.

DNA repair function scores for 2172 variants in the BRCA1 amino-terminus.

Single nucleotide variants are the most frequent type of sequence changes detected in the genome and these are frequently variants of uncertain significance (VUS). VUS are changes in DNA for which disease risk association is unknown. Thus, methods that classify the functional impact of a VUS can be used as evidence for variant interpretation. In the case of the breast and ovarian cancer specific tumor suppressor protein, BRCA1, pathogenic missense variants frequently score as loss of function in an assay for homology-directed repair (HDR) of DNA double-strand breaks. We previously published functional results using a multiplexed assay for 1056 amino acid substitutions residues 2-192 in the amino terminus of BRCA1. In this study, we have re-assessed the data from this multiplexed assay using an improved analysis pipeline. These new analysis methods yield functional scores for more variants in the first 192 amino acids of BRCA1, plus we report new results for BRCA1 amino acid residues 193-302. We now present the functional classification of 2172 BRCA1 variants in BRCA1 residues 2-302 using the multiplexed HDR assay. Comparison of the functional determinations of the missense variants with clinically known benign or pathogenic variants indicated 93% sensitivity and 100% specificity for this assay. The results from BRCA1 variants tested in this assay are a resource for clinical geneticists for evidence to evaluate VUS in BRCA1.

The functional impact of BRCA1 BRCT domain variants using multiplexed DNA double-strand break repair assays.

Pathogenic variants in BRCA1 are associated with a greatly increased risk of hereditary breast and ovarian cancer (HBOC). With the increased availability and affordability of genetic testing, many individuals have been identified with BRCA1 variants of uncertain significance (VUSs), which are individually detected in the population too infrequently to ascertain a clinical risk. Functional assays can be used to experimentally assess the effects of these variants. In this study, we used multiplexed DNA repair assays of variants in the BRCA1 carboxyl terminus to functionally characterize 2,271 variants for homology-directed repair function (HDR) and 1,427 variants for cisplatin resistance (CR). We found a high level of consistent results (Pearson's r = 0.74) in the two multiplexed functional assays with non-functional variants located within regions of the BRCA1 protein necessary for its tumor suppression activity. In addition, functional categorizations of variants tested in the multiplex HDR and CR assays correlated with known clinical significance and with other functional assays for BRCA1 (Pearson's r = 0.53 to 0.71). The results of the multiplex HDR and CR assays are useful resources for characterizing large numbers of BRCA1 VUSs.

Comparative Diagnostic Utility of SARS-CoV-2 Rapid Antigen and Molecular Testing in a Community Setting.

BACKGROUND: SARS-CoV-2 antigen-detection rapid diagnostic tests (Ag-RDTs) have become widely utilized but longitudinal characterization of their community-based performance remains incompletely understood. METHODS: This prospective longitudinal study at a large public university in Seattle, WA utilized remote enrollment, online surveys, and self-collected nasal swab specimens to evaluate Ag-RDT performance against real-time reverse transcription polymerase chain reaction (rRT-PCR) in the context of SARS-CoV-2 Omicron. Ag-RDT sensitivity and specificity within 1 day of rRT-PCR were evaluated by symptom status throughout the illness episode and Orf1b cycle threshold (Ct). RESULTS: From February to December 2022, 5757 participants reported 17 572 Ag-RDT results and completed 12 674 rRT-PCR tests, of which 995 (7.9%) were rRT-PCR positive. Overall sensitivity and specificity were 53.0% (95% confidence interval [CI], 49.6%-56.4%) and 98.8% (95% CI, 98.5%-99.0%), respectively. Sensitivity was comparatively higher for Ag-RDTs used 1 day after rRT-PCR (69.0%), 4-7 days after symptom onset (70.1%), and Orf1b Ct ≤20 (82.7%). Serial Ag-RDT sensitivity increased with repeat testing ≥2 (68.5%) and ≥4 (75.8%) days after an initial Ag-RDT-negative result. CONCLUSIONS: Ag-RDT performance varied by clinical characteristics and temporal testing patterns. Our findings support recommendations for serial testing following an initial Ag-RDT-negative result, especially among recently symptomatic persons or those at high risk for SARS-CoV-2 infection.

Assigning credit where it is due: an information content score to capture the clinical value of multiplexed assays of variant effect.

BACKGROUND: A variant can be pathogenic or benign with relation to a human disease. Current classification categories from benign to pathogenic reflect a probabilistic summary of the current understanding. A primary metric of clinical utility for multiplexed assays of variant effect (MAVE) is the number of variants that can be reclassified from uncertain significance (VUS). However, a gap in this measure of utility is that it underrepresents the information gained from MAVEs. The aim of this study was to develop an improved quantification metric for MAVE utility. We propose adopting an information content approach that includes data that does not reclassify variants will better reflect true information gain. We adopted an information content approach to evaluate the information gain, in bits, for MAVEs of BRCA1, PTEN, and TP53. Here, one bit represents the amount of information required to completely classify a single variant starting from no information. RESULTS: BRCA1 MAVEs produced a total of 831.2 bits of information, 6.58% of the total missense information in BRCA1 and a 22-fold increase over the information that only contributed to VUS reclassification. PTEN MAVEs produced 2059.6 bits of information which represents 32.8% of the total missense information in PTEN and an 85-fold increase over the information that contributed to VUS reclassification. TP53 MAVEs produced 277.8 bits of information which represents 6.22% of the total missense information in TP53 and a 3.5-fold increase over the information that contributed to VUS reclassification. CONCLUSIONS: An information content approach will more accurately portray information gained through MAVE mapping efforts than by counting the number of variants reclassified. This information content approach may also help define the impact of guideline changes that modify the information definitions used to classify groups of variants.

Clinical and Genomic Epidemiology of Coxsackievirus A21 and Enterovirus D68 in Homeless Shelters, King County, Washington, USA, 2019-2021.

Congregate homeless shelters are disproportionately affected by infectious disease outbreaks. We describe enterovirus epidemiology across 23 adult and family shelters in King County, Washington, USA, during October 2019-May 2021, by using repeated cross-sectional respiratory illness and environmental surveillance and viral genome sequencing. Among 3,281 participants >3 months of age, we identified coxsackievirus A21 (CVA21) in 39 adult residents (3.0% [95% CI 1.9%-4.8%] detection) across 7 shelters during October 2019-February 2020. We identified enterovirus D68 (EV-D68) in 5 adult residents in 2 shelters during October-November 2019. Of 812 environmental samples, 1 was EV-D68-positive and 5 were CVA21-positive. Other enteroviruses detected among residents, but not in environmental samples, included coxsackievirus A6/A4 in 3 children. No enteroviruses were detected during April 2020-May 2021. Phylogenetically clustered CVA21 and EV-D68 cases occurred in some shelters. Some shelters also hosted multiple CVA21 lineages.

Closing the gap: Systematic integration of multiplexed functional data resolves variants of uncertain significance in BRCA1, TP53, and PTEN.

Clinical interpretation of missense variants is challenging because the majority identified by genetic testing are rare and their functional effects are unknown. Consequently, most variants are of uncertain significance and cannot be used for clinical diagnosis or management. Although not much can be done to ameliorate variant rarity, multiplexed assays of variant effect (MAVEs), where thousands of single-nucleotide variant effects are simultaneously measured experimentally, provide functional evidence that can help resolve variants of unknown significance (VUSs). However, a rigorous assessment of the clinical value of multiplexed functional data for variant interpretation is lacking. Thus, we systematically combined previously published BRCA1, TP53, and PTEN multiplexed functional data with phenotype and family history data for 324 VUSs identified by a single diagnostic testing laboratory. We curated 49,281 variant functional scores from MAVEs for these three genes and integrated four different TP53 multiplexed functional datasets into a single functional prediction for each variant by using machine learning. We then determined the strength of evidence provided by each multiplexed functional dataset and reevaluated 324 VUSs. Multiplexed functional data were effective in driving variant reclassification when combined with clinical data, eliminating 49% of VUSs for BRCA1, 69% for TP53, and 15% for PTEN. Thus, multiplexed functional data, which are being generated for numerous genes, are poised to have a major impact on clinical variant interpretation.

Tiny swabs: nasal swabs integrated into tube caps facilitate large-scale self-collected SARS-CoV-2 testing.

The COVID-19 pandemic spurred the development of innovative solutions for specimen collection and molecular detection for large-scale community testing. Among these developments is the RHINOstic nasal swab, a plastic anterior nares swab built into the cap of a standard matrix tube that facilitates automated processing of up to 96 specimens at a time. In a study of unsupervised self-collection utilizing these swabs, we demonstrate comparable analytic performance and shipping stability compared to traditional anterior nares swabs, as well as significant improvements in laboratory processing efficiency. The use of these swabs may allow laboratories to accommodate large numbers of sample collections during periods of high testing demand. Automation-friendly nasal swabs are an important tool for high-throughput processing of samples that may be adopted in response to future respiratory viral pandemics.

Remote surveillance and detection of SARS-CoV-2 transmission among household members in King County, Washington.

BACKGROUND: Early during the COVID-19 pandemic, it was important to better understand transmission dynamics of SARS-CoV-2, the virus that causes COVID-19. Household contacts of infected individuals are particularly at risk for infection, but delays in contact tracing, delays in testing contacts, and isolation and quarantine posed challenges to accurately capturing secondary household cases. METHODS: In this study, 346 households in the Seattle region were provided with respiratory specimen collection kits and remotely monitored using web-based surveys for respiratory illness symptoms weekly between October 1, 2020, and June 20, 2021. Symptomatic participants collected respiratory specimens at symptom onset and mailed specimens to the central laboratory in Seattle. Specimens were tested for SARS-CoV-2 using RT-PCR with whole genome sequencing attempted when positive. SARS-CoV-2-infected individuals were notified, and their household contacts submitted specimens every 2 days for 14 days. RESULTS: In total, 1371 participants collected 2029 specimens that were tested; 16 individuals (1.2%) within 6 households tested positive for SARS-CoV-2 during the study period. Full genome sequences were generated from 11 individuals within 4 households. Very little genetic variation was found among SARS-CoV-2 viruses sequenced from different individuals in the same household, supporting transmission within the household. CONCLUSIONS: This study indicates web-based surveillance of respiratory symptoms, combined with rapid and longitudinal specimen collection and remote contact tracing, provides a viable strategy to monitor households and detect household transmission of SARS-CoV-2. TRIAL REGISTRATION IDENTIFIER: NCT04141930, Date of registration 28/10/2019.

Accurate classification of BRCA1 variants with saturation genome editing.

Variants of uncertain significance fundamentally limit the clinical utility of genetic information. The challenge they pose is epitomized by BRCA1, a tumour suppressor gene in which germline loss-of-function variants predispose women to breast and ovarian cancer. Although BRCA1 has been sequenced in millions of women, the risk associated with most newly observed variants cannot be definitively assigned. Here we use saturation genome editing to assay 96.5% of all possible single-nucleotide variants (SNVs) in 13 exons that encode functionally critical domains of BRCA1. Functional effects for nearly 4,000 SNVs are bimodally distributed and almost perfectly concordant with established assessments of pathogenicity. Over 400 non-functional missense SNVs are identified, as well as around 300 SNVs that disrupt expression. We predict that these results will be immediately useful for the clinical interpretation of BRCA1 variants, and that this approach can be extended to overcome the challenge of variants of uncertain significance in additional clinically actionable genes.

CRaTER enrichment for on-target gene editing enables generation of variant libraries in hiPSCs.

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CRISPRa On-Target Editing Retrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of variants in MYH7, a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different variants. We differentiated these hiPSCs to cardiomyocytes and show MHC-ß fusion proteins can localize as expected. Additionally, single-cell contractility analyses revealed cardiomyocytes with a pathogenic, hypertrophic cardiomyopathy-associated MYH7 variant exhibit salient HCM physiology relative to isogenic controls. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of functional transgenic cell lines at unprecedented scale.

Multiplexed Functional Assessment of Genetic Variants in CARD11.

Genetic testing has increased the number of variants identified in disease genes, but the diagnostic utility is limited by lack of understanding variant function. CARD11 encodes an adaptor protein that expresses dominant-negative and gain-of-function variants associated with distinct immunodeficiencies. Here, we used a "cloning-free" saturation genome editing approach in a diploid cell line to simultaneously score 2,542 variants for decreased or increased function in the region of CARD11 associated with immunodeficiency. We also described an exon-skipping mechanism for CARD11 dominant-negative activity. The classification of reported clinical variants was sensitive (94.6%) and specific (88.9%), which rendered the data immediately useful for interpretation of seven coding and splicing variants implicated in immunodeficiency found in our clinic. This approach is generalizable for variant interpretation in many other clinically actionable genes, in any relevant cell type.

Phylogenomics of Severe Acute Respiratory Syndrome Coronavirus 2 in Emergency Shelters for People Experiencing Homelessness.

BACKGROUND: Residents and staff of emergency shelters for people experiencing homelessness (PEH) are at high risk of infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The importance of shelter-related transmission of SARS-CoV-2 in this population remains unclear. It is also unknown whether there is significant spread of shelter-related viruses into surrounding communities. METHODS: We analyzed genome sequence data for 28 SARS-CoV-2-positive specimens collected from 8 shelters in King County, Washington between March and October, 2020. RESULTS: We identified at least 12 separate SARS-CoV-2 introduction events into these 8 shelters and estimated that 57% (16 of 28) of the examined cases of SARS-CoV-2 infection were the result of intrashelter transmission. However, we identified just a few SARS-CoV-2 specimens from Washington that were possible descendants of shelter viruses. CONCLUSIONS: Our data suggest that SARS-CoV-2 spread in shelters is common, but we did not observe evidence of widespread transmission of shelter-related viruses into the general population.

The Clinical and Genomic Epidemiology of Rhinovirus in Homeless Shelters-King County, Washington.

BACKGROUND: Rhinovirus (RV) is a common cause of respiratory illness in all people, including those experiencing homelessness. RV epidemiology in homeless shelters is unknown. METHODS: We analyzed data from a cross-sectional homeless shelter study in King County, Washington, October 2019-May 2021. Shelter residents or guardians aged ≥3 months reporting acute respiratory illness completed questionnaires and submitted nasal swabs. After 1 April 2020, enrollment expanded to residents and staff regardless of symptoms. Samples were tested by multiplex RT-PCR for respiratory viruses. A subset of RV-positive samples was sequenced. RESULTS: There were 1066 RV-positive samples with RV present every month of the study period. RV was the most common virus before and during the coronavirus disease 2019 (COVID-19) pandemic (43% and 77% of virus-positive samples, respectively). Participants from family shelters had the highest prevalence of RV. Among 131 sequenced samples, 33 RV serotypes were identified with each serotype detected for ≤4 months. CONCLUSIONS: RV infections persisted through community mitigation measures and were most prevalent in shelters housing families. Sequencing showed a diversity of circulating RV serotypes, each detected over short periods of time. Community-based surveillance in congregate settings is important to characterize respiratory viral infections during and after the COVID-19 pandemic. CLINICAL TRIALS REGISTRATION: NCT04141917.

Associations Between Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants and Risk of Coronavirus Disease 2019 (COVID-19) Hospitalization Among Confirmed Cases in Washington State: A Retrospective Cohort Study.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic is dominated by variant viruses; the resulting impact on disease severity remains unclear. Using a retrospective cohort study, we assessed the hospitalization risk following infection with 7 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants. METHODS: Our study includes individuals with positive SARS-CoV-2 reverse transcription polymerase chain reaction (RT-PCR) in the Washington Disease Reporting System with available viral genome data, from 1 December 2020 to 14 January 2022. The analysis was restricted to cases with specimens collected through sentinel surveillance. Using a Cox proportional hazards model with mixed effects, we estimated hazard ratios (HR) for hospitalization risk following infection with a variant, adjusting for age, sex, calendar week, and vaccination. RESULTS: In total, 58 848 cases were sequenced through sentinel surveillance, of which 1705 (2.9%) were hospitalized due to COVID-19. Higher hospitalization risk was found for infections with Gamma (HR 3.20, 95% confidence interval [CI] 2.40-4.26), Beta (HR 2.85, 95% CI 1.56-5.23), Delta (HR 2.28 95% CI 1.56-3.34), or Alpha (HR 1.64, 95% CI 1.29-2.07) compared to infections with ancestral lineages; Omicron (HR 0.92, 95% CI .56-1.52) showed no significant difference in risk. Following Alpha, Gamma, or Delta infection, unvaccinated patients show higher hospitalization risk, while vaccinated patients show no significant difference in risk, both compared to unvaccinated, ancestral lineage cases. Hospitalization risk following Omicron infection is lower with vaccination. CONCLUSIONS: Infection with Alpha, Gamma, or Delta results in a higher hospitalization risk, with vaccination attenuating that risk. Our findings support hospital preparedness, vaccination, and genomic surveillance.

Human Parainfluenza Virus in Homeless Shelters before and during the COVID-19 Pandemic, Washington, USA.

To determine the epidemiology of human parainfluenza virus in homeless shelters during the COVID-19 pandemic, we analyzed data and sequences from respiratory specimens collected in 23 shelters in Washington, USA, during 2019-2021. Two clusters in children were genetically similar by shelter of origin. Shelter-specific interventions are needed to reduce these infections.

Simultaneous monitoring of HIV viral load and screening of SARS-CoV-2 employing a low-cost RT-qPCR test workflow.

The COVID-19 pandemic interrupted routine care for individuals living with HIV, putting them at risk of virologic failure and HIV-associated illness. Often this population is at high risk for exposure to SARS-CoV-2 infection, and once infected, for severe disease. Therefore, close monitoring of HIV plasma viral load (VL) and screening for SARS-CoV-2 infection are needed. We developed a non-proprietary method to isolate RNA from plasma, nasal secretions (NS), or both. The extracted RNA is then submitted to RT-qPCR to estimate the VL and classify HIV/SARS-CoV-2 status (i.e., HIV virologic failure or suppressed; SARS-CoV-2 as positive, presumptive positive, negative, or indeterminate). In contrived samples, the in-house RNA extraction workflow achieved a detection limit of 200-copies per mL for HIV RNA in plasma and 100-copies per mL for SARS-CoV-2 RNA in NS. Similar detection limits were observed for HIV and SARS-CoV-2 in pooled plasma/NS contrived samples. When comparing in-house with standard extraction methods, we found high agreement (>0.91) between input and measured RNA copies for HIV LTR in contrived plasma; SARS-CoV-2 N1/N2 in contrived NS; and LTR, N1, and N2 in pooled plasma/NS samples. We further evaluated this workflow on 133 clinical specimens: 40 plasma specimens (30 HIV-positive), 67 NS specimens (31 SARS-CoV-2-positive), and 26 combined plasma/NS specimens (26 HIV-positive with 10 SARS-CoV-2-positive), and compared the results obtained using the in-house RNA extraction to those using a commercial kit (standard extraction method). The in-house extraction and standard extraction of clinical specimens were positively correlated: plasma HIV VL (R2 of 0.81) and NS SARS-CoV-2 VL (R2 of 0.95 and 0.99 for N1 and N2 genes, respectively); and pooled plasma/NS HIV VL (R2 of 0.71) and SARS-CoV-2 VL (R2 of 1 both for N1 and N2 genes). Our low-cost molecular test workflow ($1.85 per pooled sample extraction) for HIV RNA and SARS-CoV-2 RNA could serve as an alternative to current standard assays ($12 per pooled sample extraction) for laboratories in low-resource settings.