Clinical Whole Genome Sequencing Research Articles

Resolving unsolved whole-genome sequencing data in paediatric neurological disorders: a cohort study

ObjectiveTo resolve unsolved whole-genome sequencing (WGS) data in individuals with paediatric neurological disorders.DesignA cohort study method using updated bioinformatic tools, new analysis targets, clinical information and literature databases was employed...

Whole-genome sequencing in prenatally detected congenital malformations: prospective cohort study in clinical setting.

To investigate the diagnostic yield of trio whole-genome sequencing (WGS) in fetuses with various congenital malformations referred to a tertiary center for prenatal diagnosis. In this prospective study, 50 pregnancies with different congenital malformations, negative for trisomies and causative copy-number variants, were analyzed further with fetal-parental trio WGS analysis. Parents were eligible for inclusion if they accepted further investigation following the detection of isolated or multiple malformations on prenatal ultrasound. Cases with isolated increased nuchal translucency, gamete donation or multiple pregnancy were excluded. WGS with the Illumina Inc. 30× polymerase-chain-reaction-free short-read sequencing included analysis of single-nucleotide variants, insertions and deletions, structural variants, short tandem repeats and copy-number identification of SMN1 and SMN2 genes. A molecular diagnosis was achieved in 13/50 (26%) cases. Causative sequence variants were identified in 12 genes: FGFR3 (n = 2), ACTA1 (n = 1), CDH2 (n = 1), COL1A2 (n = 1), DHCR7 (n = 1), EYA1 (n = 1), FBXO11 (n = 1), FRAS1 (n = 1), L1CAM (n = 1), OFD1 (n = 1), PDHA1 (n = 1) and SOX9 (n = 1). The phenotypes of the cases were divided into different groups, with the following diagnostic yields: skeletal malformation (4/9 (44%)), multisystem malformation (3/7 (43%)), central nervous system malformation (5/15 (33%)) and thoracic malformation (1/10 (10%)). Additionally, two cases carried variants that were considered potentially clinically relevant, even though they were assessed as variants of uncertain significance, according to the guidelines provided by the American College of Medical Genetics and Genomics. Overall, we identified a causative or potentially clinically relevant variant in 15/50 (30%) cases. We demonstrate a diagnostic yield of 26% with clinical WGS in prenatally detected congenital malformations. This study emphasizes the benefits that WGS can bring to the diagnosis of fetal structural anomalies. It is important to note that causative chromosomal aberrations were excluded from our cohort before WGS. As chromosomal aberrations are a well-known cause of prenatally detected congenital malformations, future studies using WGS as a primary diagnostic test, including assessment of chromosomal aberrations, may show that the detection rate exceeds the diagnostic yield of this study. WGS can add clinically relevant information, explaining the underlying cause of the fetal anomaly, which will provide information concerning the specific prognosis of the condition, as well as estimate the risk of recurrence. A genetic diagnosis can also provide more reproductive choice for future pregnancies. © 2024 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

Impact of deep phenotyping: high diagnostic yield in a diverse pediatric population of 172 patients through clinical whole-genome sequencing at a single center.

Background: Pediatric patients with undiagnosed conditions, particularly those suspected of having Mendelian genetic disorders, pose a significant challenge in healthcare. This study investigates the diagnostic yield of whole-genome sequencing (WGS) in a pediatric cohort with diverse phenotypes, particularly focusing on the role of clinical expertise in interpreting WGS results. Methods: A retrospective cohort study was conducted at Acibadem University's Maslak Hospital in Istanbul, Turkey, involving pediatric patients (0-18years) who underwent diagnostic WGS testing. Clinical assessments, family histories, and previous laboratory and imaging studies were analyzed. Variants were classified and interpreted in conjunction with clinical findings. Results: The cohort comprised 172 pediatric patients, aged 0-5years (62.8%). International patients (28.5%) were from 20 different countries. WGS was used as a first-tier approach in 61.6% of patients. The diagnostic yield of WGS reached 61.0%, enhanced by reclassification of variants of uncertain significance (VUS) through reverse phenotyping by an experienced clinical geneticist. Consanguinity was 18.6% of the overall cohort. Dual diagnoses were carried out for 8.5% of solved patients. Discussion: Our study particularly advocates for the selection of WGS as a first-tier testing approach in infants and children with rare diseases, who were under 5years of age, thereby potentially shortening the duration of the diagnostic odyssey. The results also emphasize the critical role of a single clinical geneticist's expertise in deep phenotyping and reverse phenotyping, which contributed significantly to the high diagnostic yield.

Abstract A034: Development of a novel class of genomic biomarkers through mathematical modeling of cancer genome contiguity

Abstract Introduction: Chromosomal instability (CIN) is a hallmark of cancer aggressiveness. However, in lung adenocarcinoma (LUAD) the most common methods for measuring genomic instability, such as the weighted Genome Instability Index (wGII), are not reliable predictors of patient outcomes. The wGII measures the portion of the genome altered by copy number variation (CNV). The wGII quantifies aneuploidy but ignores other properties of CNVs such as their length distribution. We hypothesized that quantification of genome contiguity may overcome the limitations of wGII and help identify novel prognostic biomarkers. Method: We develop a novel measure of CIN based on mathematical modeling of genome breakage patterns. Utilizing segmented CNV profiles we infer two properties of the genome breakage process – its intensity and dispersion. Intensity is proportional to the number of genome breakpoints and is estimated by the total number of CNV segments. Dispersion quantifies spatial dependence of breakpoints over a wide range of genomic scales. It is computed by Ripley’s K-function which compares the distribution of breakpoints with a random Poisson point process. Jointly, intensity and dispersion, which we term the weighted Genome Contiguity Index (wGCI), estimate the contiguity of the genome. We contrast the wGCI with a heuristic score based on quantiles of segment lengths. Results: We sought to evaluate the association between the wGCI and patients’ survival in a proteogenomic cohort of 215 primary LUAD tumors profiled using whole-genome sequencing (WGS). In support of our hypothesis, we observed that patients had significantly better outcomes if their tumor genomes had longer segments compared to those with a high degree of fragmentation (P=0.014). Utilizing the wGCI, we found that tumor genomes with dispersed breakage (and low contiguity) were associated with worst outcomes (P=0.01). To identify mechanisms associated with dispersed breakage patterns we compared genomic features of the two groups. We found that LUAD tumors with highly fragmented genomes have a higher rate of TERT amplifications (P=0.005) and increased TERT RNA expression (P=2e-6). These tumors are also enriched for TP53 mutations (P=2e-16) and MYC amplifications (P=0.01). Given the cost and practical challenges of clinical WGS we leveraged mass-spectrometry proteomics to identify protein biomarkers. We found IGF2BP3 to be significantly upregulated (P=0.0006) in highly fragmented tumors and to be associated with worse outcomes (P=0.011). We use immunohistochemistry (IHC) to develop and validate IGF2BP3 as a prognostic biomarker in independent lung cancer cohorts. Conclusion: We introduced a novel quantitative measure of CIN that uniquely captures the impact of chromosome breakage on LUAD patients' clinical outcomes. Further, we identified upregulation of IGF2BP3 as a proxy to identify fragmented genomes using IHC. Our work underscores the value of mathematical modeling of genome breakage patterns to develop complex genomic biomarkers. Citation Format: Noshad Hosseini, Rahul Mannan, Alexey Nesvizhskii, Saravana Dhanasekaran, Marcin Cieslik. Development of a novel class of genomic biomarkers through mathematical modeling of cancer genome contiguity [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr A034.

Computational pharmacogenotype extraction from clinical next-generation sequencing.

Next-generation sequencing (NGS), including whole genome sequencing (WGS) and whole exome sequencing (WES), is increasingly being used for clinic care. While NGS data have the potential to be repurposed to support clinical pharmacogenomics (PGx), current computational approaches have not been widely validated using clinical data. In this study, we assessed the accuracy of the Aldy computational method to extract PGx genotypes from WGS and WES data for 14 and 13 major pharmacogenes, respectively. Germline DNA was isolated from whole blood samples collected for 264 patients seen at our institutional molecular solid tumor board. DNA was used for panel-based genotyping within our institutional Clinical Laboratory Improvement Amendments- (CLIA-) certified PGx laboratory. DNA was also sent to other CLIA-certified commercial laboratories for clinical WGS or WES. Aldy v3.3 and v4.4 were used to extract PGx genotypes from these NGS data, and results were compared to the panel-based genotyping reference standard that contained 45 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, TPMT, and VKORC1. Mean WGS read depth was >30x for all variant regions except for G6PD (average read depth was 29 reads), and mean WES read depth was >30x for all variant regions. For 94 patients with WGS, Aldy v3.3 diplotype calls were concordant with those from the genotyping reference standard in 99.5% of cases when excluding diplotypes with additional major star alleles not tested by targeted genotyping, ambiguous phasing, and CYP2D6 hybrid alleles. Aldy v3.3 identified 15 additional clinically actionable star alleles not covered by genotyping within CYP2B6, CYP2C19, DPYD, SLCO1B1, and NUDT15. Within the WGS cohort, Aldy v4.4 diplotype calls were concordant with those from genotyping in 99.7% of cases. When excluding patients with CYP2D6 copy number variation, all Aldy v4.4 diplotype calls except for one CYP3A4 diplotype call were concordant with genotyping for 161 patients in the WES cohort. Aldy v3.3 and v4.4 called diplotypes for major pharmacogenes from clinical WES and WGS data with >99% accuracy. These findings support the use of Aldy to repurpose clinical NGS data to inform clinical PGx.

The Need for Speed and Energy Efficiency in Genome Analysis

The Need for Speed and Energy Efficiency in Genome Analysis

Increased diagnostic yield from negative whole genome-slice panels using automated reanalysis.

We evaluated the diagnostic yield using genome-slice panel reanalysis in the clinical setting using an automated phenotype/gene ranking system. We analyzed whole genome sequencing (WGS) data produced from clinically ordered panels built as bioinformatic slices for 16 clinically diverse, undiagnosed cases referred to the Pediatric Mendelian Genomics Research Center, an NHGRI-funded GREGoR Consortium site. Genome-wide reanalysis was performed using Moon™, a machine-learning-based tool for variant prioritization. In five out of 16 cases, we discovered a potentially clinically significant variant. In four of these cases, the variant was found in a gene not included in the original panel due to phenotypic expansion of a disorder or incomplete initial phenotyping of the patient. In the fifth case, the gene containing the variant was included in the original panel, but being a complex structural rearrangement with intronic breakpoints outside the clinically analyzed regions, it was not initially identified. Automated genome-wide reanalysis of clinical WGS data generated during targeted panels testing yielded a 25% increase in diagnostic findings and a possibly clinically relevant finding in one additional case, underscoring the added value of analyses versus those routinely performed in the clinical setting.

An improved method for determining frequency of multiple variants of SARS-CoV-2 in wastewater using qPCR assays

Wastewater-based surveillance has become an effective tool around the globe for indirect monitoring of COVID-19 in communities. Variants of Concern (VOCs) have been detected in wastewater by use of reverse transcription polymerase chain reaction (RT-PCR) or whole genome sequencing (WGS). Rapid, reliable RT-PCR assays continue to be needed to determine the relative frequencies of VOCs and sub-lineages in wastewater-based surveillance programs. The presence of multiple mutations in a single region of the N-gene allowed for the design of a single amplicon, multiple probe assay, that can distinguish among several VOCs in wastewater RNA extracts. This approach which multiplexes probes designed to target mutations associated with specific VOC's along with an intra-amplicon universal probe (non-mutated region) was validated in singleplex and multiplex. The prevalence of each mutation (i.e. VOC) is estimated by comparing the abundance of the targeted mutation with a non-mutated and highly conserved region within the same amplicon. This is advantageous for the accurate and rapid estimation of variant frequencies in wastewater. The N200 assay was applied to monitor frequencies of VOCs in wastewater extracts from several communities in Ontario, Canada in near real time from November 28, 2021 to January 4, 2022. This includes the period of the rapid replacement of the Delta variant with the introduction of the Omicron variant in these Ontario communities in early December 2021. The frequency estimates using this assay were highly reflective of clinical WGS estimates for the same communities. This style of qPCR assay, which simultaneously measures signal from a non-mutated comparator probe and multiple mutation-specific probes contained within a single qPCR amplicon, can be applied to future assay development for rapid and accurate estimations of variant frequencies.

Test development, optimization and validation of a WGS pipeline for genetic disorders

BackgroundWith advances in massive parallel sequencing (MPS) technology, whole-genome sequencing (WGS) has gradually evolved into the first-tier diagnostic test for genetic disorders. However, deployment practice and pipeline testing for clinical WGS are lacking.MethodsIn this study, we introduced a whole WGS pipeline for genetic disorders, which included the entire process from obtaining a sample to clinical reporting. All samples that underwent WGS were constructed using polymerase chain reaction (PCR)-free library preparation protocols and sequenced on the MGISEQ-2000 platform. Bioinformatics pipelines were developed for the simultaneous detection of various types of variants, including single nucleotide variants (SNVs), insertions and deletions (indels), copy number variants (CNVs) and balanced rearrangements, mitochondrial (MT) variants, and other complex variants such as repeat expansion, pseudogenes and absence of heterozygosity (AOH). A semiautomatic pipeline was developed for the interpretation of potential SNVs and CNVs. Forty-five samples (including 14 positive commercially available samples, 23 laboratory-held positive cell lines and 8 clinical cases) with known variants were used to validate the whole pipeline.ResultsIn this study, a whole WGS pipeline for genetic disorders was developed and optimized. Forty-five samples with known variants (6 with SNVs and Indels, 3 with MT variants, 5 with aneuploidies, 1 with triploidy, 23 with CNVs, 5 with balanced rearrangements, 2 with repeat expansions, 1 with AOHs, and 1 with exon 7–8 deletion of SMN1 gene) validated the effectiveness of our pipeline.ConclusionsThis study has been piloted in test development, optimization, and validation of the WGS pipeline for genetic disorders. A set of best practices were recommended using our pipeline, along with a dataset of positive samples for benchmarking.

Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method.

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

Scalable, high quality, whole genome sequencing from archived, newborn, dried blood spots

Universal newborn screening (NBS) is a highly successful public health intervention. Archived dried bloodspots (DBS) collected for NBS represent a rich resource for population genomic studies. To fully harness this resource in such studies, DBS must yield high-quality genomic DNA (gDNA) for whole genome sequencing (WGS). In this pilot study, we hypothesized that gDNA of sufficient quality and quantity for WGS could be extracted from archived DBS up to 20 years old without PCR (Polymerase Chain Reaction) amplification. We describe simple methods for gDNA extraction and WGS library preparation from several types of DBS. We tested these methods in DBS from 25 individuals who had previously undergone diagnostic, clinical WGS and 29 randomly selected DBS cards collected for NBS from the California State Biobank. While gDNA from DBS had significantly less yield than from EDTA blood from the same individuals, it was of sufficient quality and quantity for WGS without PCR. All samples DBS yielded WGS that met quality control metrics for high-confidence variant calling. Twenty-eight variants of various types that had been reported clinically in 19 samples were recapitulated in WGS from DBS. There were no significant effects of age or paper type on WGS quality. Archived DBS appear to be a suitable sample type for WGS in population genomic studies.

105. Clinical whole-genome sequencing identifies NSD3 as the correct fusion partner of NUP98 in a patient with acute myeloid

<h3>Background</h3> Only rare cases of acute myeloid leukemia (AML) have been shown to harbor a t(8;11)(p11.2;p15.4). This translocation is believed to involve the fusion of NSD3 or FGFR1 with NUP98; however, apart from targeted mRNA quantitative PCR analysis, no molecular approaches have been utilized to define the chimeric fusions present in these rare cases. <h3>Case Presentation</h3> Here we present the case of a 51-year-old female with AML with myelodysplastic-related morphologic changes, 13q deletion and t(8;11), where initial fluorescence in situ hybridization (FISH) assays were consistent with the presence of NUP98 and FGFR1 rearrangements, and suggestive of NUP98/FGFR1 fusion. Using a streamlined clinical whole-genome sequencing approach, we resolved the breakpoints of this translocation to intron 4 of NSD3 and intron 12 of NUP98, indicating NUP98/NSD3 rearrangement as the underlying aberration. Furthermore, our approach identified small variants in WT1 and STAG2, as well as an interstitial deletion on the short arm of chromosome 12, which were cryptic in G-banded chromosomes. <h3>Conclusions</h3> NUP98 fusions in acute leukemia are predictive of poor prognosis. The associated fusion partner and the presence of co-occurring mutations, such as WT1, further refine this prognosis with potential clinical implications. Using a clinical whole-genome sequencing analysis, we resolved t(8;11) breakpoints to NSD3 and NUP98, ruling out the involvement of FGFR1 suggested by FISH while also identifying multiple chromosomal and sequence level aberrations.

Epidemiological investigation of Pseudomonas aeruginosa isolates including multidrug-resistant serogroup O12 isolates, by use of a rapid and simplified multiple-locus variable-number of tandem repeats analysis and whole genome sequencing

Clustered cases of Pseudomonas aeruginosa infection in immunocompromised patients' wards require rapid characterization of a potential epidemic to guide investigations and identify the potential source of contamination. To design and evaluate a rapid and simple typing method for P.aeruginosa in comparison to whole genome sequencing (WGS). A simplified polymerase chain reaction based on multiple-locus variable-number of tandem repeats analysis (MLVA) was designed and used to investigate cases of P.aeruginosa infection and colonization in a paediatric haematology department. The method was compared to WGS by using the Illumina method. On the 17 isolates recovered from 15 children (eight from blood cultures, three from urinary tract infections, one from sputum and five stool isolates), MLVA distinguished 10 different profiles, and seven isolates from six children shared the same profile. Analysis by WGS revealed that these seven isolates belonged to sequence type ST111 and serotype O12, allowing at least three different genotypes to be distinguished among them. Five environmental strains had three MLVA profiles; one was shared with a clinical isolate but WGS excluded any relationship. The simplified and inexpensive MLVA method enabled the exclusion, in less than 5h, of most of the unrelated isolates and thus to focus investigations on a small number of cases, whereas WGS, taking several days of work, drew definitive conclusions concerning the outbreak and the genetic relationships of the ST111 isolates circulating in the department. We conclude that sequential use of both methods is the optimal strategy to investigate clustered cases of P.aeruginosa infections.

Diagnostic yield using whole-genome sequencing and in-silico panel of 281 genes associated with non-immune hydrops fetalis in clinical setting.

To investigate the diagnostic yield of clinical whole-genome sequencing (WGS) in prenatally diagnosed non-immune hydrops fetalis (NIHF). This was a retrospective study of 23 fetuses with prenatally diagnosed NIHF, negative for trisomies and copy-number variants, referred for analysis by WGS with an in-silico panel of 281 genes associated with hydrops fetalis. Due to identification of a high proportion of causative variants in the HRAS gene in the main cohort, Sanger sequencing of HRAS was performed in a replication cohort, consisting of 24 additional fetuses with NIHF that were negative for trisomies and copy-number variants and had not undergone WGS. Of the 23 fetuses in the main cohort, a molecular diagnosis was achieved in 12 (52.2%). Pathogenic or likely pathogenic variants were identified in seven genes: HRAS (n = 5), RIT1 (n = 2), FOXP3 (n = 1), GLB1 (n = 1), MAP2K1 (n = 1), PTPN11 (n = 1) and RASA1 (n = 1). The inheritance pattern of the 12 causative variants was autosomal dominant in 10 cases (HRAS, MAP2K1, PTPN11, RASA1, RIT1), autosomal recessive in one (GLB1) and X-linked recessive in one (FOXP3). Of the 24 fetuses in the replication cohort, a pathogenic variant in HRAS was identified in one, resulting in an overall frequency of causative HRAS variants of 12.8% (6/47) in our two cohorts. We demonstrate a diagnostic yield of 52% with clinical WGS in NIHF using an in-silico panel of 281 genes. However, the high diagnostic yield may be attributed to the small sample size and possible over-representation of severe phenotypes in the included fetuses. Bearing in mind that chromosomal abnormalities were excluded in our cohorts, a detection rate of up to 75% is possible in prenatally diagnosed NIHF when WGS analysis includes calling of chromosomal aberrations. © 2022 The Authors. Ultrasound in Obstetrics &amp; Gynecology published by John Wiley &amp; Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

Personalised selection of experimental treatment in patients with advanced solid cancer is feasible using whole-genome sequencing

BackgroundBiomarker-guided therapy in an experimental setting has been suggested to improve patient outcomes. However, trial-specific pre-screening tests are time and tissue consuming and complicate the personalised treatment of patients eligible for early-phase clinical trials. In this study the feasibility of whole-genome sequencing (WGS) as a one-test-for-all for guided inclusion in early-phase trials was investigated.MethodsPhase I Molecular Tumor Board (MTB) at the Erasmus MC Cancer Institute reviewed patients with advanced cancer without standard-of-care treatment (SOC) options for a ‘fresh-frozen’ (FF) tumour biopsy for WGS based on clinical-pathological features. Clinical grade WGS was performed by Hartwig Medical Foundation. MTB matched the patient with a trial, if available.ResultsFrom September 2019–March 2021, 31 patients with highly diverse tumour types underwent a tumour biopsy for WGS. The median turnaround time (TAT) was 15 days [10–42 days]. At least one actionable event was found in 84% of the patients (26/31). One-third of the patients (11/31) received matched experimental treatment.ConclusionsWGS on fresh FF biopsies is a feasible tool for the selection of personalised experimental therapy in patients with advanced cancer without SOC options. WGS is now possible in an acceptable TAT and thus could fulfil the role of a universal genomic pre-screening test.

Performance characterization of PCR-free whole genome sequencing for clinical diagnosis.

To evaluate the performance of polymerase chain reaction (PCR)-free whole genome sequencing (WGS) for clinical diagnosis, and thereby revealing how experimental parameters affect variant detection.Five NA12878 samples were sequenced using MGISEQ-2000. NA12878 samples underwent WGS with differing deoxyribonucleic acid (DNA) input and library preparation protocol (PCR-based vs PCR-free protocols for library preparation). The depth of coverage and genotype quality of each sample were compared. The performance of each sample was measured for sensitivity, coverage of depth and breadth of coverage of disease-related genes, and copy number variants. We also developed a systematic WGS pipeline (PCR-free) for the analysis of 11 clinical cases.In general, NA12878-2 (PCR-free WGS) showed better depth of coverage and genotype quality distribution than NA12878-1 (PCR-based WGS). With a mean depth of ∼40×, the sensitivity of homozygous and heterozygous single nucleotide polymorphisms (SNPs) of NA12878-2 showed higher sensitivity (>99.77% and >99.82%) than NA12878-1, and positive predictive value exceeded 99.98% and 99.07%. The sensitivity and positive predictive value of homozygous and heterozygous indels for NA12878-2 (PCR-free WGS) showed great improvement than NA128878-1. The breadths of coverage for disease-related genes and copy number variants are slightly better for samples with PCR-free library preparation protocol than the sample with PCR-based library preparation protocol. DNA input also influences the performance of variant detection in samples with PCR-free WGS. All the 19 previously confirmed variants in 11 clinical cases were successfully detected by our WGS pipeline (PCR free).Different experimental parameters may affect variant detection for clinical WGS. Clinical scientists should know the range of sensitivity of variants for different methods of WGS, which would be useful when interpreting and delivering clinical reports.

The benefit of diagnostic whole genome sequencing in schizophrenia and other psychotic disorders.

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families.

Hyper-IgE Syndrome due to an Elusive Novel Intronic Homozygous Variant in DOCK8.

Rare, biallelic loss-of-function mutations in DOCK8 result in a combined immune deficiency characterized by severe and recurrent cutaneous infections, eczema, allergies, and susceptibility to malignancy, as well as impaired humoral and cellular immunity and hyper-IgE. The advent of next-generation sequencing technologies has enabled the rapid molecular diagnosis of rare monogenic diseases, including inborn errors of immunity. These advances have resulted in the implementation of gene-guided treatments, such as hematopoietic stem cell transplant for DOCK8 deficiency. However, putative disease-causing variants revealed by next-generation sequencing need rigorous validation to demonstrate pathogenicity. Here, we report the eventual diagnosis of DOCK8 deficiency in a consanguineous family due to a novel homozygous intronic deletion variant that caused aberrant exon splicing and subsequent loss of expression of DOCK8 protein. Remarkably, the causative variant was not initially detected by clinical whole-genome sequencing but was subsequently identified and validated by combining advanced genomic analysis, RNA-seq, and flow cytometry. This case highlights the need to adopt multipronged confirmatory approaches to definitively solve complex genetic cases that result from variants outside protein-coding exons and conventional splice sites.

Effect of Whole-Genome Sequencing on the Clinical Management of Acutely Ill Infants With Suspected Genetic Disease

Whole-genome sequencing (WGS) shows promise as a first-line genetic test for acutely ill infants, but widespread adoption and implementation requires evidence of an effect on clinical management. To determine the effect of WGS on clinical management in a racially and ethnically diverse and geographically distributed population of acutely ill infants in the US. This randomized, time-delayed clinical trial enrolled participants from September 11, 2017, to April 30, 2019, with an observation period extending to July 2, 2019. The study was conducted at 5 US academic medical centers and affiliated children's hospitals. Participants included infants aged between 0 and 120 days who were admitted to an intensive care unit with a suspected genetic disease. Data were analyzed from January 14 to August 20, 2020. Patients were randomized to receive clinical WGS results 15 days (early) or 60 days (delayed) after enrollment, with the observation period extending to 90 days. Usual care was continued throughout the study. The main outcome was the difference in the proportion of infants in the early and delayed groups who received a change of management (COM) 60 days after enrollment. Additional outcome measures included WGS diagnostic efficacy, within-group COM at 90 days, length of hospital stay, and mortality. A total of 354 infants were randomized to the early (n = 176) or delayed (n = 178) arms. The mean participant age was 15 days (IQR, 7-32 days); 201 participants (56.8%) were boys; 19 (5.4%) were Asian; 47 (13.3%) were Black; 250 (70.6%) were White; and 38 (10.7%) were of other race. At 60 days, twice as many infants in the early group vs the delayed group received a COM (34 of 161 [21.1%; 95% CI, 15.1%-28.2%] vs 17 of 165 [10.3%; 95% CI, 6.1%-16.0%]; P = .009; odds ratio, 2.3; 95% CI, 1.22-4.32) and a molecular diagnosis (55 of 176 [31.0%; 95% CI, 24.5%-38.7%] vs 27 of 178 [15.0%; 95% CI, 10.2%-21.3%]; P < .001). At 90 days, the delayed group showed a doubling of COM (to 45 of 161 [28.0%; 95% CI, 21.2%-35.6%]) and diagnostic efficacy (to 56 of 178 [31.0%; 95% CI, 24.7%-38.8%]). The most frequent COMs across the observation window were subspecialty referrals (39 of 354; 11%), surgery or other invasive procedures (17 of 354; 4%), condition-specific medications (9 of 354; 2%), or other supportive alterations in medication (12 of 354; 3%). No differences in length of stay or survival were observed. In this randomized clinical trial, for acutely ill infants in an intensive care unit, introduction of WGS was associated with a significant increase in focused clinical management compared with usual care. Access to first-line WGS may reduce health care disparities by enabling diagnostic equity. These data support WGS adoption and implementation in this population. ClinicalTrials.gov Identifier: NCT03290469.

Low-Level Rifampin Resistance and rpoB Mutations in Mycobacterium tuberculosis: an Analysis of Whole-Genome Sequencing and Drug Susceptibility Test Data in New York.

Rapid and reliable detection of rifampin (RIF) resistance is critical for the diagnosis and treatment of drug-resistant and multidrug-resistant (MDR) tuberculosis. Discordant RIF phenotype/genotype susceptibility results remain a challenge due to the presence of rpoB mutations that do not confer high levels of RIF resistance, as have been exhibited in strains with mutations such as Ser450Leu. These strains, termed low-level RIF resistant, exhibit elevated RIF MICs compared to fully susceptible strains but remain phenotypically susceptible by mycobacterial growth indicator tube (MGIT) testing and have been associated with poor patient outcomes. Here, we assess RIF resistance prediction by whole-genome sequencing (WGS) among a set of 1,779 prospectively tested strains by both prevalence of rpoB gene mutation and phenotype as part of routine clinical testing during a 2.5-year period. During this time, 139 strains were found to have nonsynonymous rpoB mutations, 53 of which were associated with RIF resistance, including both low-level and high-level resistance. Resistance to RIF (1.0 μg/ml in MGIT) was identified in 43 (81.1%) isolates. The remaining 10 (18.9%) strains were susceptible by MGIT but were confirmed to be low-level RIF resistant by MIC testing. Full rpoB gene sequencing overcame the limitations of critical concentration phenotyping, probe-based genotyping, and partial gene sequencing methods. Universal clinical WGS with concurrent phenotypic testing provided a more complete understanding of the prevalence and type of rpoB mutations and their association with RIF resistance in New York.