Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

A framework for the clinical implementation of optical genome mapping in hematologic malignancies.

Optical Genome Mapping (OGM) is rapidly emerging as an exciting cytogenomic technology both for research and clinical purposes. In the last 2 years alone, multiple studies have demonstrated that OGM not only matches the diagnostic scope of conventional standard of care cytogenomic clinical testing but it also adds significant new information in certain cases. Since OGM consolidates the diagnostic benefits of multiple costly and laborious tests (e.g., karyotyping, fluorescence in situ hybridization, and chromosomal microarrays) in a single cost-effective assay, many clinical laboratories have started to consider utilizing OGM. In 2021, an international working group of early adopters of OGM who are experienced with routine clinical cytogenomic testing in patients with hematological neoplasms formed a consortium (International Consortium for OGM in Hematologic Malignancies, henceforth "the Consortium") to create a consensus framework for implementation of OGM in a clinical setting. The focus of the Consortium is to provide guidance for laboratories implementing OGM in three specific areas: validation, quality control and analysis and interpretation of variants. Since OGM is a complex technology with many variables, we felt that by consolidating our collective experience, we could provide a practical and useful tool for uniform implementation of OGM in hematologic malignancies with the ultimate goal of achieving globally accepted standards.

A complex structural variant near SOX3 causes X-linked split-hand/foot malformation.

Split-hand/foot malformation (SHFM) is a congenital limb defect most typically presenting with median clefts in hands and/or feet, that can occur in a syndromic context as well as in isolated form. SHFM is caused by failure to maintain normal apical ectodermal ridge function during limb development. Although several genes and contiguous gene syndromes are implicated in the monogenic etiology of isolated SHFM, the disorder remains genetically unexplained for many families and associated genetic loci. We describe a family with isolated X-linked SHFM, for which the causative variant could be detected after a diagnostic journey of 20 years. We combined well-established approaches including microarray-based copy number variant analysis and fluorescence in situ hybridization coupled with optical genome mapping and whole genome sequencing. This strategy identified a complex structural variant (SV) comprising a 165-kb gain of 15q26.3 material ([GRCh37/hg19] chr15:99795320-99960362dup) inserted in inverted position at the site of a 38-kb deletion on Xq27.1 ([GRCh37/hg19] chrX:139481061-139518989del). In silico analysis suggested that the SV disrupts the regulatory framework on the X chromosome and may lead to SOX3 misexpression. We hypothesize that SOX3 dysregulation in the developing limb disturbed the fine balance between morphogens required for maintaining AER function, resulting in SHFM in this family.

Biallelic variants in CRIPT cause a Rothmund-Thomson-like syndrome with increased cellular senescence.

PURPOSE: Rothmund-Thomson syndrome (RTS) is characterized by poikiloderma, sparse hair, small stature, skeletal defects, cancer, and cataracts, resembling features of premature aging. RECQL4 and ANAPC1 are the 2 known disease genes associated with RTS in >70% of cases. We describe RTS-like features in 5 individuals with biallelic variants in CRIPT (OMIM 615789). METHODS: Two newly identified and 4 published individuals with CRIPT variants were systematically compared with those with RTS using clinical data, computational analysis of photographs, histologic analysis of skin, and cellular studies on fibroblasts. RESULTS: All CRIPT individuals fulfilled the diagnostic criteria for RTS and additionally had neurodevelopmental delay and seizures. Using computational gestalt analysis, CRIPT individuals showed greatest facial similarity with individuals with RTS. Skin biopsies revealed a high expression of senescence markers (p53/p16/p21) and the senescence-associated ß-galactosidase activity was elevated in CRIPT-deficient fibroblasts. RECQL4- and CRIPT-deficient fibroblasts showed an unremarkable mitotic progression and unremarkable number of mitotic errors and no or only mild sensitivity to genotoxic stress by ionizing radiation, mitomycin C, hydroxyurea, etoposide, and potassium bromate. CONCLUSION: CRIPT causes an RTS-like syndrome associated with neurodevelopmental delay and epilepsy. At the cellular level, RECQL4- and CRIPT-deficient cells display increased senescence, suggesting shared molecular mechanisms leading to the clinical phenotypes.

Precision oncology using organoids of a secretory carcinoma of the salivary gland treated with TRK-inhibitors.

The use of anticancer drugs targeting specific molecular tumor characteristics is rapidly increasing in clinical practice, but selecting patients to benefit from these remains a challenge. It has been suggested that organoid cultures would be ideally suited to test drug responses in vitro. Here we describe and characterize in depth a case of ETV6-NTRK3 gene fusion-positive secretory carcinoma of the salivary glands and corresponding organoid cultures that responded and subsequently acquired resistance to TRK targeting therapy with larotrectinib. This case-culture-characterization illustrates the advances made in precision oncology, but also exposes important caveats in using organoids to predict treatment response.

Long-Read Sequencing Identifies Novel Pathogenic Intronic Variants in Gitelman Syndrome.

BACKGROUND: Gitelman syndrome is a salt-losing tubulopathy characterized by hypokalemic alkalosis and hypomagnesemia. It is caused by homozygous recessive or compound heterozygous pathogenic variants in SLC12A3 , which encodes the Na + -Cl - cotransporter (NCC). In up to 10% of patients with Gitelman syndrome, current genetic techniques detect only one specific pathogenic variant. This study aimed to identify a second pathogenic variant in introns, splice sites, or promoters to increase the diagnostic yield. METHODS: Long-read sequencing of SLC12A3 was performed in 67 DNA samples from individuals with suspected Gitelman syndrome in whom a single likely pathogenic or pathogenic variant was previously detected. In addition, we sequenced DNA samples from 28 individuals with one variant of uncertain significance or no candidate variant. Midigene splice assays assessed the pathogenicity of novel intronic variants. RESULTS: A second likely pathogenic/pathogenic variant was identified in 45 (67%) patients. Those with two likely pathogenic/pathogenic variants had a more severe electrolyte phenotype than other patients. Of the 45 patients, 16 had intronic variants outside of canonic splice sites (nine variants, mostly deep intronic, six novel), whereas 29 patients had an exonic variant or canonic splice site variant. Midigene splice assays of the previously known c.1670-191C>T variant and intronic candidate variants demonstrated aberrant splicing patterns. CONCLUSION: Intronic pathogenic variants explain an important part of the missing heritability in Gitelman syndrome. Long-read sequencing should be considered in diagnostic workflows for Gitelman syndrome.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes.

PURPOSE: Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs. METHODS: Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria. RESULTS: In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses. CONCLUSION: Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

Lessons learned from rapid exome sequencing for 575 critically ill patients across the broad spectrum of rare disease.

Introduction: Rapid exome sequencing (rES) has become the first-choice genetic test for critically ill patients, mostly neonates, young infants, or fetuses in prenatal care, in time-sensitive situations and when it is expected that the genetic test result may guide clinical decision making. The implementation of rES has revolutionized medicine by enabling timely identification of genetic causes for various rare diseases. The utilization of rES has increasingly been recognized as an essential diagnostic tool for the identification of complex and undiagnosed genetic disorders. Methods: We conducted a retrospective evaluation of our experiences with rES performed on 575 critically ill patients from various age groups (prenatal to adulthood), over a four-year period (2016-2019). These patients presented with a wide spectrum of rare diseases, including but not limited to neurological disorders, severe combined immune deficiency, and cancer. Results: During the study period, there was a significant increase in rES referrals, with a rise from a total of two referrals in Q1-2016 to 10 referrals per week in Q4-2019. The median turnaround time (TAT) decreased from 17 to 11 days in the period 2016-2019, with an overall median TAT of 11 days (IQR 8-15 days). The overall diagnostic yield for this cohort was 30.4%, and did not significantly differ between the different age groups (e.g. adults 22.2% vs children 31.0%; p-value 0.35). However, variability in yield was observed between clinical entities: craniofacial anomalies yielded 58.3%, while for three clinical entities (severe combined immune deficiency, aneurysm, and hypogonadotropic hypogonadism) no diagnoses were obtained. Discussion: Importantly, whereas clinical significance is often only attributed to a conclusive diagnosis, we also observed impact on clinical decision-making for individuals in whom no genetic diagnosis was established. Hence, our experience shows that rES has an important role for patients of all ages and across the broad spectrum of rare diseases to impact clinical outcomes.

Diagnostic analysis of the highly complex OPN1LW/OPN1MW gene cluster using long-read sequencing and MLPA.

Pathogenic variants in the OPN1LW/OPN1MW gene cluster are causal for a range of mild to severe visual impairments with color deficiencies. The widely utilized short-read next-generation sequencing (NGS) is inappropriate for the analysis of the OPN1LW/OPN1MW gene cluster and many patients with pathogenic variants stay underdiagnosed. A diagnostic genetic assay was developed for the OPN1LW/OPN1MW gene cluster, consisting of copy number analysis via multiplex ligation-dependent probe amplification and sequence analysis via long-read circular consensus sequencing. Performance was determined on 50 clinical samples referred for genetic confirmation of the clinical diagnosis (n = 43) or carrier status analysis (n = 7). A broad range of pathogenic haplotypes were detected, including deletions, hybrid genes, single variants and combinations of variants. The developed genetic assay for the OPN1LW/OPN1MW gene cluster is a diagnostic test that can detect both structural and nucleotide variants with a straightforward analysis, improving diagnostic care of patients with visual impairment.

Scrutinizing pathogenicity of the USH2A c.2276 G > T; p.(Cys759Phe) variant.

The USH2A variant c.2276 G > T (p.(Cys759Phe)) has been described by many authors as a frequent cause of autosomal recessive retinitis pigmentosa (arRP). However, this is in contrast with the description of two asymptomatic individuals homozygous for this variant. We therefore assessed pathogenicity of the USH2A c.2276 G > T variant using extensive genetic and functional analyses. Whole genome sequencing and optical genome mapping were performed for three arRP cases homozygous for USH2A c.2276 G > T to exclude alternative genetic causes. A minigene splice assay was designed to investigate the effect of c.2276 G > T on pre-mRNA splicing, in presence or absence of the nearby c.2256 T > C variant. Moreover, an ush2ap.(Cys771Phe) zebrafish knock-in model mimicking human p.(Cys759Phe) was generated and characterized using functional and immunohistochemical analyses. Besides the homozygous c.2276 G > T USH2A variant, no alternative genetic causes were identified. Evaluation of the ush2ap.(Cys771Phe) zebrafish model revealed strongly reduced levels of usherin expression at the photoreceptor periciliary membrane, increased levels of rhodopsin localization in the photoreceptor cell body and decreased electroretinogram (ERG) b-wave amplitudes compared to wildtype controls. In conclusion, we confirmed pathogenicity of USH2A c.2276 G > T (p.(Cys759Phe)). Consequently, cases homozygous for c.2276 G > T can now receive a definite genetic diagnosis and can be considered eligible for receiving future QR-421a-mediated exon 13 skipping therapy.

Optical genome mapping for structural variation analysis in hematologic malignancies.

Optical genome mapping (OGM) is a technology that is rapidly being adopted in clinical genetics laboratories for its ability to detect structural variation (e.g., translocations, inversions, deletions, duplications, etc.) and replace several concurrent standard of care techniques (karyotype, fluorescence in situ hybridization, and chromosomal microarray). OGM can dramatically simplify lab workflow by reducing multiple tests (conventional karyotype, fluorescence in situ hybridization [FISH], and chromosomal microarray) into one test. The superior ability to detect structural variation across the genome removes the need for reflex FISH studies, which can dramatically reduce cost and turnaround time per sample. Parallel studies of OGM versus standard of care testing have demonstrated it can detect and resolve more abnormalities than karyotyping or FISH. However, like many molecular tests that normalize copy number it can have difficulty with non-diploid karyotypes. This Test of the Month review will summarize how the technique works, review the strengths and weaknesses of OGM compared to standard of care techniques and illustrate how the technique is likely to change front line testing in many hematologic malignancies-including summarizing the clinical utility in acute myeloid leukemia, myelodysplastic syndromes, and B cell acute lymphoblastic leukemia.

RNA Sequencing Elucidates Drug-Specific Mechanisms of Antibiotic Tolerance and Resistance in Mycobacterium abscessus.

Mycobacterium abscessus is an opportunistic pathogen notorious for its resistance to most classes of antibiotics and low cure rates. M. abscessus carries an array of mostly unexplored defense mechanisms. A deeper understanding of antibiotic resistance and tolerance mechanisms is pivotal in development of targeted therapeutic regimens. We provide the first description of all major transcriptional mechanisms of tolerance to all antibiotics recommended in current guidelines, using RNA sequencing-guided experiments. M. abscessus ATCC 19977 bacteria were subjected to subinhibitory concentrations of clarithromycin (CLR), amikacin (AMK), tigecycline (TIG), cefoxitin (FOX), and clofazimine (CFZ) for 4 and 24 h, followed by RNA sequencing. To confirm key mechanisms of tolerance suggested by transcriptomic responses, we performed time-kill kinetic analysis using bacteria after preexposure to CLR, AMK, or TIG for 24 h and constructed isogenic knockout and knockdown strains. To assess strain specificity, pan-genome analysis of 35 strains from all three subspecies was performed. Mycobacterium abscessus shows both drug-specific and common transcriptomic responses to antibiotic exposure. Ribosome-targeting antibiotics CLR, AMK, and TIG elicit a common response characterized by upregulation of ribosome structural genes, the WhiB7 regulon and transferases, accompanied by downregulation of respiration through NuoA-N. Exposure to any of these drugs decreases susceptibility to ribosome-targeting drugs from multiple classes. The cytochrome bd-type quinol oxidase contributes to CFZ tolerance in M. abscessus, and the sigma factor sigH but not antisigma factor MAB_3542c is involved in TIG resistance. The observed transcriptomic responses are not strain-specific, as all genes involved in tolerance, except erm(41), are found in all included strains.

Exploring the missing heritability in subjects with hearing loss, enlarged vestibular aqueducts, and a single or no pathogenic SLC26A4 variant.

Pathogenic variants in SLC26A4 have been associated with autosomal recessive hearing loss (arHL) and a unilateral or bilateral enlarged vestibular aqueduct (EVA). SLC26A4 is the second most frequently mutated gene in arHL. Despite the strong genotype-phenotype correlation, a significant part of cases remains genetically unresolved. In this study, we investigated a cohort of 28 Dutch index cases diagnosed with HL in combination with an EVA but without (M0) or with a single (M1) pathogenic variant in SLC26A4. To explore the missing heritability, we first determined the presence of the previously described EVA-associated haplotype (Caucasian EVA (CEVA)), characterized by 12 single nucleotide variants located upstream of SLC26A4. We found this haplotype and a delimited V1-CEVA haplotype to be significantly enriched in our M1 patient cohort (10/16 cases). The CEVA haplotype was also present in two M0 cases (2/12). Short- and long-read whole genome sequencing and optical genome mapping could not prioritize any of the variants present within the CEVA haplotype as the likely pathogenic defect. Short-read whole-genome sequencing of the six M1 cases without this haplotype and the two M0/CEVA cases only revealed previously overlooked or misinterpreted splice-altering SLC26A4 variants in two cases, who are now genetically explained. No deep-intronic or structural variants were identified in any of the M1 subjects. With this study, we have provided important insights that will pave the way for elucidating the missing heritability in M0 and M1 SLC26A4 cases. For pinpointing the pathogenic effect of the CEVA haplotype, additional analyses are required addressing defect(s) at the RNA, protein, or epigenetic level.

Characterization of a genomic region 8 kb downstream of GFI1B associated with myeloproliferative neoplasms.

A genomic locus 8 kb downstream of the transcription factor GFI1B (Growth Factor Independence 1B) predisposes to clonal hematopoiesis and myeloproliferative neoplasms. One of the most significantly associated polymorphisms in this region is rs621940-G. GFI1B auto-represses GFI1B, and altered GFI1B expression contributes to myeloid neoplasms. We studied whether rs621940-G affects GFI1B expression and growth of immature cells. GFI1B ChIP-seq showed clear binding to the rs621940 locus. Preferential binding of various hematopoietic transcription factors to either the rs621940-C or -G allele was observed, but GFI1B showed no preference. In gene reporter assays the rs621940 region inhibited GFI1B promoter activity with the G-allele having less suppressive effects compared to the C-allele. However, CRISPR-Cas9 mediated deletion of the locus in K562 cells did not alter GFI1B expression nor auto-repression. In healthy peripheral blood mononuclear cells GFI1B expression did not differ consistently between the rs621940 alleles. Long range and targeted deep sequencing did not detect consistent effects of rs621940-G on allelic GFI1B expression either. Finally, we observed that myeloid colony formation was not significantly affected by either rs621940 allele in 193 healthy donors. Together, these findings show no evidence that rs621940 or its locus affect GFI1B expression, auto-repression or growth of immature myeloid cells.

Optical genome mapping enables constitutional chromosomal aberration detection.

Chromosomal aberrations including structural variations (SVs) are a major cause of human genetic diseases. Their detection in clinical routine still relies on standard cytogenetics. Drawbacks of these tests are a very low resolution (karyotyping) and the inability to detect balanced SVs or indicate the genomic localization and orientation of duplicated segments or insertions (copy number variant [CNV] microarrays). Here, we investigated the ability of optical genome mapping (OGM) to detect known constitutional chromosomal aberrations. Ultra-high-molecular-weight DNA was isolated from 85 blood or cultured cells and processed via OGM. A de novo genome assembly was performed followed by structural variant and CNV calling and annotation, and results were compared to known aberrations from standard-of-care tests (karyotype, FISH, and/or CNV microarray). In total, we analyzed 99 chromosomal aberrations, including seven aneuploidies, 19 deletions, 20 duplications, 34 translocations, six inversions, two insertions, six isochromosomes, one ring chromosome, and four complex rearrangements. Several of these variants encompass complex regions of the human genome involved in repeat-mediated microdeletion/microduplication syndromes. High-resolution OGM reached 100% concordance compared to standard assays for all aberrations with non-centromeric breakpoints. This proof-of-principle study demonstrates the ability of OGM to detect nearly all types of chromosomal aberrations. We also suggest suited filtering strategies to prioritize clinically relevant aberrations and discuss future improvements. These results highlight the potential for OGM to provide a cost-effective and easy-to-use alternative that would allow comprehensive detection of chromosomal aberrations and structural variants, which could give rise to an era of "next-generation cytogenetics."

Next-generation cytogenetics: Comprehensive assessment of 52 hematological malignancy genomes by optical genome mapping.

Somatic structural variants (SVs) are important drivers of cancer development and progression. In a diagnostic set-up, especially for hematological malignancies, the comprehensive analysis of all SVs in a given sample still requires a combination of cytogenetic techniques, including karyotyping, FISH, and CNV microarrays. We hypothesize that the combination of these classical approaches could be replaced by optical genome mapping (OGM). Samples from 52 individuals with a clinical diagnosis of a hematological malignancy, divided into simple (<5 aberrations, n = 36) and complex (≥5 aberrations, n = 16) cases, were processed for OGM, reaching on average: 283-fold genome coverage. OGM called a total of 918 high-confidence SVs per sample, of which, on average, 13 were rare and >100 kb. In addition, on average, 73 CNVs were called per sample, of which six were >5 Mb. For the 36 simple cases, all clinically reported aberrations were detected, including deletions, insertions, inversions, aneuploidies, and translocations. For the 16 complex cases, results were largely concordant between standard-of-care and OGM, but OGM often revealed higher complexity than previously recognized. Detailed technical comparison with standard-of-care tests showed high analytical validity of OGM, resulting in a sensitivity of 100% and a positive predictive value of >80%. Importantly, OGM resulted in a more complete assessment than any previous single test and most likely reported the most accurate underlying genomic architecture (e.g., for complex translocations, chromoanagenesis, and marker chromosomes). In conclusion, the excellent concordance of OGM with diagnostic standard assays demonstrates its potential to replace classical cytogenetic tests as well as to rapidly map novel leukemia drivers.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors.

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat-rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long-range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE-VNTR-Alu, subfamily E (SVA-E) retrotransposon element, which was present in a mosaic state in the mother. This SVA-E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long-read sequencing to identify genomic variations in high-risk cancer-predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.