Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program.

BACKGROUND: The Centers for Disease Control and Prevention contracted with laboratories to sequence the SARS-CoV-2 genome from positive samples across the United States to enable public health officials to investigate the impact of variants on disease severity as well as the effectiveness of vaccines and treatment. Herein we present the initial results correlating RT-PCR quality control metrics with sample collection and sequencing methods from full SARS-CoV-2 viral genomic sequencing of 24,441 positive patient samples between April and June 2021. METHODS: RT-PCR confirmed (N Gene Ct value < 30) positive patient samples, with nucleic acid extracted from saliva, nasopharyngeal and oropharyngeal swabs were selected for viral whole genome SARS-CoV-2 sequencing. Sequencing was performed using Illumina COVIDSeq™ protocol on either the NextSeq550 or NovaSeq6000 systems. Informatic variant calling, and lineage analysis were performed using DRAGEN COVID Lineage applications on Illumina's Basespace cloud analytical system. All sequence data and variant calls were uploaded to NCBI and GISAID. RESULTS: An association was observed between higher sequencing coverage, quality, and samples with a lower Ct value, with < 27 being optimal, across both sequencing platforms and sample collection methods. Both nasopharyngeal swabs and saliva samples were found to be optimal samples of choice for SARS-CoV-2 surveillance sequencing studies, both in terms of strain identification and sequencing depth of coverage, with NovaSeq 6000 providing higher coverage than the NextSeq 550. The most frequent variants identified were the B.1.617.2 Delta (India) and P.1 Gamma (Brazil) variants in the samples sequenced between April 2021 and June 2021. At the time of submission, the most common variant > 99% of positives sequenced was Omicron. CONCLUSION: These initial analyses highlight the importance of sequencing platform, sample collection methods, and RT-PCR Ct values in guiding surveillance efforts. These surveillance studies evaluating genetic changes of SARS-CoV-2 have been identified as critical by the CDC that can affect many aspects of public health including transmission, disease severity, diagnostics, therapeutics, and vaccines.

Rapid Whole Genome Sequencing Has Clinical Utility in Children in the PICU.

OBJECTIVES: Genetic disorders are a leading contributor to mortality in the neonatal ICU and PICU in the United States. Although individually rare, there are over 6,200 single-gene diseases, which may preclude a genetic diagnosis prior to ICU admission. Rapid whole genome sequencing is an emerging method of diagnosing genetic conditions in time to affect ICU management of neonates; however, its clinical utility has yet to be adequately demonstrated in critically ill children. This study evaluates next-generation sequencing in pediatric critical care. DESIGN: Retrospective cohort study. SETTING: Single-center PICU in a tertiary children's hospital. PATIENTS: Children 4 months to 18 years admitted to the PICU who were nominated between July 2016 and May 2018. INTERVENTIONS: Rapid whole genome sequencing with targeted phenotype-driven analysis was performed on patients and their parents, when parental samples were available. MEASUREMENTS AND MAIN RESULTS: A molecular diagnosis was made by rapid whole genome sequencing in 17 of 38 children (45%). In four of the 17 patients (24%), the genetic diagnoses led to a change in management while in the PICU, including genome-informed changes in pharmacotherapy and transition to palliative care. Nine of the 17 diagnosed children (53%) had no dysmorphic features or developmental delay. Eighty-two percent of diagnoses affected the clinical management of the patient and/or family after PICU discharge, including avoidance of biopsy, administration of factor replacement, and surveillance for disorder-related sequelae. CONCLUSIONS: This study demonstrates a retrospective evaluation for undiagnosed genetic disease in the PICU and clinical utility of rapid whole genome sequencing in a portion of critically ill children. Further studies are needed to identify PICU patients who will benefit from rapid whole genome sequencing early in PICU admission when the underlying etiology is unclear.

MDS-associated somatic mutations and clonal hematopoiesis are common in idiopathic cytopenias of undetermined significance.

Establishing a diagnosis in patients suspected of having a myelodysplastic syndrome (MDS) can be challenging and could be informed by the identification of somatic mutations. We performed a prospective study to examine the frequency and types of mutations encountered in 144 patients with unexplained cytopenias. Based on bone marrow findings, 17% were diagnosed with MDS, 15% with idiopathic cytopenias of undetermined significance (ICUS) and some evidence of dysplasia, and 69% with ICUS and no dysplasia. Bone marrow DNA was sequenced for mutations in 22 frequently mutated myeloid malignancy genes. Somatic mutations were identified in 71% of MDS patients, 62% of patients with ICUS and some dysplasia, and 20% of ICUS patients and no dysplasia. In total, 35% of ICUS patients carried a somatic mutation or chromosomal abnormality indicative of clonal hematopoiesis. We validated these results in a cohort of 91 lower-risk MDS and 249 ICUS cases identified over a 6-month interval. Mutations were found in 79% of those with MDS, in 45% of those with ICUS with dysplasia, and in 17% of those with ICUS without dysplasia. The spectrum of mutated genes was similar with the exception of SF3B1 which was rarely mutated in patients without dysplasia. Variant allele fractions were comparable between clonal ICUS (CCUS) and MDS as were mean age and blood counts. We demonstrate that CCUS is a more frequent diagnosis than MDS in cytopenic patients. Clinical and mutational features are similar in these groups and may have diagnostic utility once outcomes in CCUS patients are better understood.

ATM-dependent MiR-335 targets CtIP and modulates the DNA damage response.

ATM plays a critical role in cellular responses to DNA double-strand breaks (DSBs). We describe a new ATM-mediated DSB-induced DNA damage response pathway involving microRNA (miRNA): irradiation (IR)-induced DSBs activate ATM, which leads to the downregulation of miR-335, a miRNA that targets CtIP, which is an important trigger of DNA end resection in homologous recombination repair (HRR). We demonstrate that CREB is responsible for a large portion of miR-335 expression by binding to the promoter region of miR-335. CREB binding is greatly reduced after IR, corroborating with previous studies that IR-activated ATM phosphorylates CREB to reduce its transcription activity. Overexpression of miR-335 in HeLa cells resulted in reduced CtIP levels and post-IR colony survival and BRCA1 foci formation. Further, in two patient-derived lymphoblastoid cell lines with decreased post-IR colony survival, a "radiosensitive" phenotype, we demonstrated elevated miR-335 expression, reduced CtIP levels, and reduced BRCA1 foci formation. Colony survival, BRCA1 foci, and CtIP levels were partially rescued by miRNA antisense AMO-miR-335 treatment. Taken together, these findings strongly suggest that an ATM-dependent CREB-miR-335-CtIP axis influences the selection of HRR for repair of certain DSB lesions.

Defective DNA double-strand break repair in pediatric systemic lupus erythematosus.

OBJECTIVE: Previous reports of cells from patients with systemic lupus erythematosus (SLE) note that repair of single-strand breaks is delayed, and these lesions may be converted to double-strand breaks (DSBs) at DNA replication forks. We undertook this study to assess the integrity of DSB recognition, signaling, and repair mechanisms in B lymphoblastoid cell lines derived from patients with pediatric SLE. METHODS: Nine assays were used to interrogate DSB repair and recognition in lymphoblastoid cell lines from patients with pediatric SLE, including the neutral comet assay (NCA), colony survival assay (CSA), irradiation-induced foci formation for γ-H2AX and 53BP1 proteins, kinetics of phosphorylation of structural maintenance of chromosomes protein 1 (SMC1), postirradiation bromodeoxyuridine incorporation to evaluate S phase checkpoint integrity, monoubiquitination of Fanconi protein D2, ATM protein expression, and non-homologous DNA end joining protein expression and function. RESULTS: Three of the 9 assays revealed abnormal patterns of response to irradiation-induced DNA damage. The NCA and CSA yielded aberrant results in the majority of SLE lymphoblastoid cell lines. Abnormal prolongation of SMC1 phosphorylation was also noted in 2 of 16 SLE lymphoblastoid cell lines. CONCLUSION: Our data suggest that DSB repair is defective in some lymphoblastoid cell lines from pediatric patients with SLE, especially when assessed by both NCA and CSA. Since these studies are nonspecific, further studies of DNA repair and kinetics are indicated to further delineate the underlying pathogenesis of SLE and possibly identify therapeutic targets.

Common copy number variations in fifty radiosensitive cell lines.

Hypersensitivity to radiation exposure is a major challenge to radiotherapy in the treatment of cancer patients. Copy number variations (CNVs) are believed to identify genomic regions of functional significance for radiosensitivity (RS) but have yet to be systematically investigated. We used Affymetrix 6.0 SNP arrays to survey common CNVs in a cohort of 50 radiosensitive lymphoblastoid cell lines (RS-LCLs) derived from patients with undiagnosed diseases. A total of 317 CNVs that were present in at least 10% of the studied cell lines were identified. Three hundred and eight CNVs overlapped with polymorphic CNVs, 13 of which were significantly enriched in the RS-LCLs compared to the reference. The remaining 9 CNVs were novel. The majority of these enriched and novel CNVs were chromosomal gains. The dominance of the chromosomal gains over losses is inconsistent with the traditional concept of molecular basis of RS and suggests more complex genetic mechanisms for RS.

An automated 13.5 hour system for scalable diagnosis and acute management guidance for genetic diseases.

While many genetic diseases have effective treatments, they frequently progress rapidly to severe morbidity or mortality if those treatments are not implemented immediately. Since front-line physicians frequently lack familiarity with these diseases, timely molecular diagnosis may not improve outcomes. Herein we describe Genome-to-Treatment, an automated, virtual system for genetic disease diagnosis and acute management guidance. Diagnosis is achieved in 13.5 h by expedited whole genome sequencing, with superior analytic performance for structural and copy number variants. An expert panel adjudicated the indications, contraindications, efficacy, and evidence-of-efficacy of 9911 drug, device, dietary, and surgical interventions for 563 severe, childhood, genetic diseases. The 421 (75%) diseases and 1527 (15%) effective interventions retained are integrated with 13 genetic disease information resources and appended to diagnostic reports ( https://gtrx.radygenomiclab.com ). This system provided correct diagnoses in four retrospectively and two prospectively tested infants. The Genome-to-Treatment system facilitates optimal outcomes in children with rapidly progressive genetic diseases.

Loss of FOCAD, operating via the SKI messenger RNA surveillance pathway, causes a pediatric syndrome with liver cirrhosis.

Cirrhosis is usually a late-onset and life-threatening disease characterized by fibrotic scarring and inflammation that disrupts liver architecture and function. While it is typically the result of alcoholism or hepatitis viral infection in adults, its etiology in infants is much less understood. In this study, we report 14 children from ten unrelated families presenting with a syndromic form of pediatric liver cirrhosis. By genome/exome sequencing, we found recessive variants in FOCAD segregating with the disease. Zebrafish lacking focad phenocopied the human disease, revealing a signature of altered messenger RNA (mRNA) degradation processes in the liver. Using patient's primary cells and CRISPR-Cas9-mediated inactivation in human hepatic cell lines, we found that FOCAD deficiency compromises the SKI mRNA surveillance pathway by reducing the levels of the RNA helicase SKIC2 and its cofactor SKIC3. FOCAD knockout hepatocytes exhibited lowered albumin expression and signs of persistent injury accompanied by CCL2 overproduction. Our results reveal the importance of FOCAD in maintaining liver homeostasis and disclose a possible therapeutic intervention point via inhibition of the CCL2/CCR2 signaling axis.

Ending a diagnostic odyssey: Moving from exome to genome to identify cockayne syndrome.

BACKGROUND: Cockayne syndrome (CS) is a rare autosomal recessive disorder characterized by growth failure and multisystemic degeneration. Excision repair cross-complementation group 6 (ERCC6 OMIM: *609413) is the gene most frequently mutated in CS. METHODS: A child with pre and postnatal growth failure and progressive neurologic deterioration with multisystem involvement, and with nondiagnostic whole-exome sequencing, was screened for causal variants with whole-genome sequencing (WGS). RESULTS: WGS identified biallelic ERCC6 variants, including a previously unreported intronic variant. Pathogenicity of these variants was established by demonstrating reduced levels of ERCC6 mRNA and protein expression, normal unscheduled DNA synthesis, and impaired recovery of RNA synthesis in patient fibroblasts following UV-irradiation. CONCLUSION: The study confirms the pathogenicity of a previously undescribed upstream intronic variant, highlighting the power of genome sequencing to identify noncoding variants. In addition, this report provides evidence for the utility of a combination approach of genome sequencing plus functional studies to provide diagnosis in a child for whom a lengthy diagnostic odyssey, including exome sequencing, was previously unrevealing.

Rapid whole genome sequencing impacts care and resource utilization in infants with congenital heart disease.

Congenital heart disease (CHD) is the most common congenital anomaly and a major cause of infant morbidity and mortality. While morbidity and mortality are highest in infants with underlying genetic conditions, molecular diagnoses are ascertained in only ~20% of cases using widely adopted genetic tests. Furthermore, cost of care for children and adults with CHD has increased dramatically. Rapid whole genome sequencing (rWGS) of newborns in intensive care units with suspected genetic diseases has been associated with increased rate of diagnosis and a net reduction in cost of care. In this study, we explored whether the clinical utility of rWGS extends to critically ill infants with structural CHD through a retrospective review of rWGS study data obtained from inpatient infants < 1 year with structural CHD at a regional children's hospital. rWGS diagnosed genetic disease in 46% of the enrolled infants. Moreover, genetic disease was identified five times more frequently with rWGS than microarray ± gene panel testing in 21 of these infants (rWGS diagnosed 43% versus 10% with microarray ± gene panels, p = 0.02). Molecular diagnoses ranged from syndromes affecting multiple organ systems to disorders limited to the cardiovascular system. The average daily hospital spending was lower in the time period post blood collection for rWGS compared to prior (p = 0.003) and further decreased after rWGS results (p = 0.000). The cost was not prohibitive to rWGS implementation in the care of this cohort of infants. rWGS provided timely actionable information that impacted care and there was evidence of decreased hospital spending around rWGS implementation.

Functional Precision Medicine Identifies New Therapeutic Candidates for Medulloblastoma.

Medulloblastoma is among the most common malignant brain tumors in children. Recent studies have identified at least four subgroups of the disease that differ in terms of molecular characteristics and patient outcomes. Despite this heterogeneity, most patients with medulloblastoma receive similar therapies, including surgery, radiation, and intensive chemotherapy. Although these treatments prolong survival, many patients still die from the disease and survivors suffer severe long-term side effects from therapy. We hypothesize that each patient with medulloblastoma is sensitive to different therapies and that tailoring therapy based on the molecular and cellular characteristics of patients' tumors will improve outcomes. To test this, we assembled a panel of orthotopic patient-derived xenografts (PDX) and subjected them to DNA sequencing, gene expression profiling, and high-throughput drug screening. Analysis of DNA sequencing revealed that most medulloblastomas do not have actionable mutations that point to effective therapies. In contrast, gene expression and drug response data provided valuable information about potential therapies for every tumor. For example, drug screening demonstrated that actinomycin D, which is used for treatment of sarcoma but rarely for medulloblastoma, was active against PDXs representing Group 3 medulloblastoma, the most aggressive form of the disease. Functional analysis of tumor cells was successfully used in a clinical setting to identify more treatment options than sequencing alone. These studies suggest that it should be possible to move away from a one-size-fits-all approach and begin to treat each patient with therapies that are effective against their specific tumor. SIGNIFICANCE: These findings show that high-throughput drug screening identifies therapies for medulloblastoma that cannot be predicted by genomic or transcriptomic analysis.

Rapid flow cytometry-based structural maintenance of chromosomes 1 (SMC1) phosphorylation assay for identification of ataxia-telangiectasia homozygotes and heterozygotes.

BACKGROUND: No rapid reliable method exists for identifying ataxia-telangiectasia (A-T) homozygotes or heterozygotes. Heterozygotes are at an increased risk of cancer and are more sensitive to the effects of ionizing radiation (IR) than the general population. We report a rapid flow cytometry (FC)-based ataxia-telangiectasia mutated (ATM) kinase assay that measures ATM- dependent phosphorylation of structural maintenance of chromosomes 1 (SMC1) following DNA damage (FC-pSMC1 assay). METHODS: After optimizing conditions with lymphoblastoid cell lines (LCLs), we studied peripheral blood mononuclear cells (PBMCs) isolated from 16 healthy donors (unknowns), 10 obligate A-T heterozygotes, and 6 unrelated A-T patients. One hour after DNA damage (by either IR or bleomycin), the cells were fixed and incubated with a primary antibody to SMC1pSer966. We analyzed the stained cells by FC to determine the difference in geometric mean fluorescence intensity (DeltaGMFI) of untreated and treated cells; this difference was expressed as a percentage of daily experimental controls. RESULTS: The FC-pSMC1 assay reliably distinguished ATM heterozygotes and homozygotes from controls. Average DeltaGMFI percentages (SD) of daily controls were, for unknowns, 106.1 (37.6); for A-T heterozygotes, 37.0 (18.7); and for A-T homozygotes; -8.73 (16.2). Values for heterozygotes and homozygotes were significantly different from those of controls (P < 0.0001). CONCLUSIONS: The FC-pSMC1 assay shortens the turnaround time for diagnosing A-T homozygotes from approximately 3 months to approximately 3 h. It also identifies A-T heterozygotes and can be used for prenatal counseling or for screening individuals in large study cohorts for potential ATM heterozygosity, which can then be confirmed by sequencing.

The case for early use of rapid whole-genome sequencing in management of critically ill infants: late diagnosis of Coffin-Siris syndrome in an infant with left congenital diaphragmatic hernia, congenital heart disease, and recurrent infections.

Congenital diaphragmatic hernia (CDH) results from incomplete formation of the diaphragm leading to herniation of abdominal organs into the thoracic cavity. CDH is associated with pulmonary hypoplasia, congenital heart disease, and pulmonary hypertension. Genetically, it is associated with aneuploidies, chromosomal copy-number variants, and single gene mutations. CDH is the most expensive noncardiac congenital defect. Management frequently requires implementation of extracorporeal membrane oxygenation (ECMO), which increases management expenditures 2.4-3.5-fold. The cost of management of CDH has been estimated to exceed $250 million per year. Despite in-hospital survival of 80%-90%, current management is imperfect, as a great proportion of surviving children have long-term functional deficits. We report the case of a premature infant prenatally diagnosed with CDH and congenital heart disease, who had a protracted and complicated course in the intensive care unit with multiple surgical interventions, including postcardiac surgery ECMO, gastrostomy tube placement with Nissen fundoplication, tracheostomy for respiratory failure, recurrent infections, and developmental delay. Rapid whole-genome sequencing (rWGS) identified a de novo, likely pathogenic, c.3096_ 3100delCAAAG (p.Lys1033Argfs*32) variant in ARID1B, providing a diagnosis of Coffin-Siris syndrome. Her parents elected palliative care and she died later that day.

Rapid whole-genome sequencing identifies a novel GABRA1 variant associated with West syndrome.

A 9-mo-old infant was admitted with infantile spasms that improved on administration of topiramate and steroids. He also had developmental delay, esotropia, and hypsarrhythmia on interictal electroencephalogram (EEG), and normal brain magnetic resonance imaging (MRI). West syndrome is the triad of infantile spasms, interictal hypsarrhythmia, and mental retardation. Rapid trio whole-genome sequencing (WGS) revealed a novel, likely pathogenic, de novo variant in the gene encoding γ-aminobutyric acid (GABA) type A receptor, α1 polypeptide (GABRA1 c.789G>A, p.Met263Ile) in the proband. GABRA1 mutations have been associated with early infantile epileptic encephalopathy type 19 (EIEE19). We suggest that GABRA1 p.Met263Ile is associated with a distinct West syndrome phenotype.

Post-irradiation phosphorylation of structural maintenance chromosome 1 (SMC1) is independent of the Fanconi protein pathway.

PURPOSE: To confirm the sensitivity of cells from patients with Fanconi anemia (FA) to ionizing radiation, and to determine whether the phosphorylation of structural maintenance chromosome 1 (SMC1) was associated with radiosensitivity, as it is in other DNA repair disorders. METHODS AND MATERIALS: Using lymphoblastoid cell lines from FA patients before and after exposure to ionizing radiation, the colony survival assay, radioresistant DNA synthesis, and SMC1 phosphorylation were measured. FA lymphoblastoid cell lines that had been transfected with the wild-type FANC gene were used as controls. RESULTS: Cells from FA patients of six complementation groups were radiosensitive. Despite this, SMC1 phosphorylation was normal in each case; radioresistant DNA synthesis, a measure of S phase checkpoint integrity, was defective in FANCD2 lymphoblastoid cell lines and was corrected in FANCD2 + D2 cells. CONCLUSIONS: The data indicate that the FANC pathway proteins play a major role in the cellular responses to ionizing radiation, but not in SMC1 phosphorylation or in the S phase checkpoint of FANCD2-deficient cells. Thus, SMC1 activation is not a common denominator of radiosensitivity, as has been suggested by radiation responses of cells from ataxia-telangiectasia, Nijmegen breakage syndrome, or Mre11 deficiency patients.