Evaluating the performance of a clinical genome sequencing program for diagnosis of rare genetic disease, seen through the lens of craniosynostosis

A Comparison of Structural Variant Calling from Short-Read and Nanopore-Based Whole-Genome Sequencing Using Optical Genome Mapping as a Benchmark

The identification of structural variants (SVs) in genomic data represents an ongoing challenge because of difficulties in reliable SV calling leading to reduced sensitivity and specificity. We prepared high-quality DNA from 9 parent-child trios, who had previously undergone short-read whole-genome sequencing (Illumina platform) as part of the Genomics England 100,000 Genomes Project. We reanalysed the genomes using both Bionano optical genome mapping (OGM; 8 probands and one trio) and Nanopore long-read sequencing (Oxford Nanopore Technologies [ONT] platform; all samples). To establish a "truth" dataset, we asked whether rare proband SV calls (n = 234) made by the Bionano Access (version 1.6.1)/Solve software (version 3.6.1_11162020) could be verified by individual visualisation using the Integrative Genomics Viewer with either or both of the Illumina and ONT raw sequence. Of these, 222 calls were verified, indicating that Bionano OGM calls have high precision (positive predictive value 95%). We then asked what proportion of the 222 true Bionano SVs had been identified by SV callers in the other two datasets. In the Illumina dataset, sensitivity varied according to variant type, being high for deletions (115/134; 86%) but poor for insertions (13/58; 22%). In the ONT dataset, sensitivity was generally poor using the original Sniffles variant caller (48% overall) but improved substantially with use of Sniffles2 (36/40; 90% and 17/23; 74% for deletions and insertions, respectively). In summary, we show that the precision of OGM is very high. In addition, when applying the Sniffles2 caller, the sensitivity of SV calling using ONT long-read sequence data outperforms Illumina sequencing for most SV types.

The phenotype of MEGF8-related Carpenter syndrome (CRPT2) is refined through the identification of eight new patients

Carpenter syndrome (CRPTS) is a rare autosomal recessive condition caused by biallelic variants in genes that encode negative regulators of hedgehog signalling (RAB23 [CRPT1] or, more rarely, MEGF8 [CRPT2]), and is characterised by craniosynostosis, polysyndactyly, and other congenital abnormalities. We describe a further six families comprising eight individuals with MEGF8-associated CRPT2, increasing the total number of reported cases to fifteen, and refine the phenotype of CRPT2 compared to CRPT1. The core features of craniosynostosis, polysyndactyly and (in males) cryptorchidism are almost universal in both CRPT1 and CRPT2. However, laterality defects are present in nearly half of those with MEGF8-associated CRPT2, but are rare in RAB23-associated CRPT1. Craniosynostosis in CRPT2 commonly involves a single midline suture in comparison to the multi-suture craniosynostosis characteristic of CRPT1. No patient to date has carried two MEGF8 gene alterations that are both predicted to lead to complete loss-of-function, suggesting that a variable degree of residual MEGF8 activity may be essential for viability and potentially contributing to variable phenotypic severity. These data refine the phenotypic spectrum of CRPT2 in comparison to CRPT1 and more than double the number of likely pathogenic MEGF8 variants in this rare disorder.

The impact of inversions across 33,924 families with rare disease from a national genome sequencing project

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Increasing the diagnostic yield of childhood glaucoma cases recruited into the 100,000 Genomes Project

Childhood glaucoma (CG) encompasses a heterogeneous group of genetic eye disorders that is responsible for approximately 5% of childhood blindness worldwide. Understanding the molecular aetiology is key to improving diagnosis, prognosis and unlocking the potential for optimising clinical management. In this study, we investigated 86 CG cases from 78 unrelated families of diverse ethnic backgrounds, recruited into the Genomics England 100,000 Genomes Project (GE100KGP) rare disease cohort, to improve the genetic diagnostic yield. Using the Genomics England/Genomic Medicine Centres (GE/GMC) diagnostic pipeline, 13 unrelated families were solved (13/78, 17%). Further interrogation using an expanded gene panel yielded a molecular diagnosis in 7 more unrelated families (7/78, 9%). This analysis effectively raises the total number of solved CG families in the GE100KGP to 26% (20/78 families). Twenty-five percent (5/20) of the solved families had primary congenital glaucoma (PCG), while 75% (15/20) had secondary CG; 53% of this group had non-acquired ocular anomalies (including iris hypoplasia, megalocornea, ectopia pupillae, retinal dystrophy, and refractive errors) and 47% had non-acquired systemic diseases such as cardiac abnormalities, hearing impairment, and developmental delay. CYP1B1 was the most frequently implicated gene, accounting for 55% (11/20) of the solved families. We identified two novel likely pathogenic variants in the TEK gene, in addition to one novel pathogenic copy number variant (CNV) in FOXC1. Variants that passed undetected in the GE100KGP diagnostic pipeline were likely due to limitations of the tiering process, the use of smaller gene panels during analysis, and the prioritisation of coding SNVs and indels over larger structural variants, CNVs, and non-coding variants.

The value of genome-wide analysis in craniosynostosis

Background: This study assessed the diagnostic yield of high-throughput sequencing methods in a cohort of craniosynostosis (CS) patients not presenting causal variants identified through previous targeted analysis. Methods: Whole-genome or whole-exome sequencing (WGS/WES) was performed in a cohort of 59 patients (from 57 families) assessed by retrospective phenotyping as having syndromic or nonsyndromic CS. Results: A syndromic form was identified in 51% of the unrelated cases. A genetic cause was identified in 38% of syndromic cases, with novel variants detected in FGFR2 (a rare Alu insertion), TWIST1, TCF12, KIAA0586, HDAC9, FOXP1, and NSD2. Additionally, we report two patients with rare recurrent variants in KAT6A and YY1 as well as two patients with structural genomic aberrations: one with a 22q13 duplication and one with a complex rearrangement involving chromosome 2 (2p25 duplication including SOX11 and deletion of 2q22). Moreover, we identified potentially relevant variants in 87% of the remaining families with no previously detected causal variants, including novel variants in ADAMTSL4, ASH1L, ATRX, C2CD3, CHD5, ERF, H4C5, IFT122, IFT140, KDM6B, KMT2D, LTBP1, MAP3K7, NOTCH2, NSD1, SOS1, SPRY1, POLR2A, PRRX1, RECQL4, TAB2, TAOK1, TET3, TGFBR1, TCF20, and ZBTB20. Conclusion: These results confirm WGS/WES as a powerful diagnostic tool capable of either targeted in silico or broad genomic analysis depending on phenotypic presentation (e.g., classical or unusual forms of syndromic CS).

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND

Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome.

METHODS

We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants.

RESULTS

Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving.

CONCLUSIONS

Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Genetic diagnostic yield in an 11-year cohort of craniosynostosis patients

Craniosynostosis may present in isolation, 'non-syndromic', or with additional congenital anomalies/neurodevelopmental disorders, 'syndromic'. Clinical focus shifted from confirming classical syndromic cases to offering genetic testing to all craniosynostosis patients. This retrospective study assesses diagnostic yield of molecular testing by investigating prevalences of chromosomal and monogenic (likely) pathogenic variants in an 11-year cohort of 1020 craniosynostosis patients. 502 children underwent genetic testing. Pathogenic variants were identified in 174 patients (35%). Diagnostic yield was significantly higher in syndromic craniosynostosis (62%) than in non-syndromic craniosynostosis (6%). Before whole exome sequencing (WES) emerged, single-gene testing was performed using Sanger sequencing or multiplex ligation-dependent probe amplification (MLPA). Diagnostic yield was 11% and was highest for EFNB1, FGFR2, FGFR3, and IL11RA. Diagnostic yield for copy number variant analysis using microarray was 8%. From 2015 onwards, the WES craniosynostosis panel was implemented, with a yield of 10%. In unsolved, mainly syndromic, cases suspected of a genetic cause, additional WES panels (multiple congenital anomalies (MCA)/intellectual disability (ID)) or open exome analysis were performed with an 18% diagnostic yield. To conclude, microarray and the WES craniosynostosis panel are key to identifying pathogenic variants. in craniosynostosis patients. Given the advances in genetic diagnostics, we should look beyond the scope of the WES craniosynostosis panel and consider extensive genetic diagnostics (e.g. open exome sequencing, whole genome sequencing, RNA sequencing and episignature analysis) if no diagnosis is obtained through microarray and/or WES craniosynostosis panel. If parents are uncomfortable with more extensive diagnostics, MCA or ID panels may be considered.

Craniosynostosis, inner ear, and renal anomalies in a child with complete loss of SPRY1 (sprouty homolog 1) function

INTRODUCTION

SPRY1 encodes protein sprouty homolog 1 (Spry-1), a negative regulator of receptor tyrosine kinase signalling. Null mutant mice display kidney/urinary tract abnormalities and altered size of the skull; complete loss-of-function of Spry-1 in humans has not been reported.

METHODS

Analysis of whole-genome sequencing data from individuals with craniosynostosis enrolled in the 100,000 Genomes Project identified a likely pathogenic variant within SPRY1. Reverse-transcriptase PCR and western blot analysis were used to investigate the effect of the variant on SPRY1 mRNA and protein, in lymphoblastoid cell lines from the patient and both parents.

RESULTS

A nonsense variant in SPRY1, encoding p.(Leu27*), was confirmed to be heterozygous in the unaffected parents and homozygous in the child. The child's phenotype, which included sagittal craniosynostosis, subcutaneous cystic lesions overlying the lambdoid sutures, hearing loss associated with bilateral cochlear and vestibular dysplasia and a unilateral renal cyst, overlapped the features reported in Spry1-/- null mice. Functional studies supported escape from nonsense-mediated decay, but western blot analysis demonstrated complete absence of full-length protein in the affected child and a marked reduction in both parents.

CONCLUSION

This is the first report of complete loss of Spry-1 function in humans, associated with abnormalities of the cranial sutures, inner ear, and kidneys.

Review of Recurrently Mutated Genes in Craniosynostosis Supports Expansion of Diagnostic Gene Panels

Craniosynostosis, the premature fusion of the cranial sutures, affects ~1 in 2000 children. Although many patients with a genetically determined cause harbor a variant in one of just seven genes or have a chromosomal abnormality, over 60 genes are known to be recurrently mutated, thus comprising a long tail of rarer diagnoses. Genome sequencing for the diagnosis of rare diseases is increasingly used in clinical settings, but analysis of the data is labor intensive and involves a trade-off between achieving high sensitivity or high precision. PanelApp, a crowd-sourced disease-focused set of gene panels, was designed to enable prioritization of variants in known disease genes for a given pathology, allowing enhanced identification of true-positives. For heterogeneous disorders like craniosynostosis, these panels must be regularly updated to ensure that diagnoses are not being missed. We provide a systematic review of genetic literature on craniosynostosis over the last 5 years, including additional results from resequencing a 42-gene panel in 617 affected individuals. We identify 16 genes (representing a 25% uplift) that should be added to the list of bona fide craniosynostosis disease genes and discuss the insights that these new genes provide into pathophysiological mechanisms of craniosynostosis.

Case report: A third variant in the 5' UTR of TWIST1 creates a novel upstream translation initiation site in a child with Saethre-Chotzen syndrome

Introduction: Saethre-Chotzen syndrome, a craniosynostosis syndrome characterized by the premature closure of the coronal sutures, dysmorphic facial features and limb anomalies, is caused by haploinsufficiency of TWIST1. Although the majority of variants localize in the coding region of the gene, two variants in the 5' UTR have been recently reported to generate novel upstream initiation codons. Methods: Skeletal dysplasia Next-generation sequencing (NGS) panel was used for genetic analysis in a patient with bicoronal synostosis, facial dysmorphisms and limb anomalies. The variant pathogenicity was assessed by a luciferase reporter promoter assay. Results: Here, we describe the identification of a third ATG-creating de novo variant, c.-18C>T, in the 5' UTR of TWIST1 in the patient with a clinical diagnosis of Saethre-Chotzen syndrome. It was predicted to create an out-of-frame new upstream translation initiation codon resulting in a 40 amino acid larger functionally inactive protein. We performed luciferase reporter promoter assays to demonstrate that the variant does indeed reduce translation from the main open reading frame. Conclusion: This is the third variant identified in this region and confirms the introduction of upstream ATGs in the 5' UTR of TWIST1 as a pathogenic mechanism in Saethre-Chotzen syndrome. This case report shows the necessity for performing functional characterization of variants of unknown significance within national health services.

Conclusion of diagnostic odysseys due to inversions disrupting GLI3 and FBN1

Many genetic testing methodologies are biased towards picking up structural variants (SVs) that alter copy number. Copy-neutral rearrangements such as inversions are therefore likely to suffer from underascertainment. In this study, manual review prompted by a virtual multidisciplinary team meeting and subsequent bioinformatic prioritisation of data from the 100K Genomes Project was performed across 43 genes linked to well-characterised skeletal disorders. Ten individuals from three independent families were found to harbour diagnostic inversions. In two families, inverted segments of 1.2/14.8 Mb unequivocally disrupted GLI3 and segregated with skeletal features consistent with Greig cephalopolysyndactyly syndrome. For one family, phenotypic blending was due to the opposing breakpoint lying ~45 kb from HOXA13 In the third family, long suspected to have Marfan syndrome, a 2.0 Mb inversion disrupting FBN1 was identified. These findings resolved lengthy diagnostic odysseys of 9-20 years and highlight the importance of direct interaction between clinicians and data-analysts. These exemplars of a rare mutational class inform future SV prioritisation strategies within the NHS Genomic Medicine Service and similar genome sequencing initiatives. In over 30 years since these two disease-gene associations were identified, large inversions have yet to be described and so our results extend the mutational spectra linked to these conditions.

Specialist multidisciplinary input maximises rare disease diagnoses from whole genome sequencing

Diagnostic whole genome sequencing (WGS) is increasingly used in rare diseases. However, standard, semi-automated WGS analysis may overlook diagnoses in complex disorders. Here, we show that specialist multidisciplinary analysis of WGS, following an initial 'no primary findings' (NPF) report, improves diagnostic rates and alters management. We undertook WGS in 102 adults with diagnostically challenging primary mitochondrial disease phenotypes. NPF cases were reviewed by a genomic medicine team, thus enabling bespoke informatic approaches, co-ordinated phenotypic validation, and functional work. We enhanced the diagnostic rate from 16.7% to 31.4%, with management implications for all new diagnoses, and detected strong candidate disease-causing variants in a further 3.9% of patients. This approach presents a standardised model of care that supports mainstream clinicians and enhances diagnostic equity for complex disorders, thereby facilitating access to the potential benefits of genomic healthcare. This research was made possible through access to the data and findings generated by the 100,000 Genomes Project: http://www.genomicsengland.co.uk .

An additional whole-exome sequencing study in 102 panel-undiagnosed patients: A retrospective study in a Chinese craniosynostosis cohort

Craniosynostosis (CRS) is a disease with prematurely fused cranial sutures. In the last decade, the whole-exome sequencing (WES) was widely used in Caucasian populations. The WES largely contributed in genetic diagnosis and exploration on new genetic mechanisms of CRS. In this study, we enrolled 264 CRS patients in China. After a 17-gene-panel sequencing designed in the previous study, 139 patients were identified with pathogenic/likely pathogenic (P/LP) variants according to the ACMG guideline as positive genetic diagnosis. WES was then performed on 102 patients with negative genetic diagnosis by panel. Ten P/LP variants were additionally identified in ten patients, increasing the genetic diagnostic yield by 3.8% (10/264). The novel variants in ANKH, H1-4, EIF5A, SOX6, and ARID1B expanded the mutation spectra of CRS. Then we designed a compatible research pipeline (RP) for further exploration. The RP could detect all seven P/LP SNVs and InDels identified above, in addition to 15 candidate variants found in 13 patients with worthy of further study. In sum, the 17-gene panel and WES identified positive genetic diagnosis for 56.4% patients (149/264) in 16 genes. At last, in our estimation, the genetic testing strategy of “Panel-first” saves 24.3% of the cost compared with “WES only”, suggesting the “Panel-first” is an economical strategy.

A systematic analysis of splicing variants identifies new diagnoses in the 100,000 Genomes Project

Background

Genomic variants which disrupt splicing are a major cause of rare genetic diseases. However, variants which lie outside of the canonical splice sites are difficult to interpret clinically. Improving the clinical interpretation of non-canonical splicing variants offers a major opportunity to uplift diagnostic yields from whole genome sequencing data.

Methods

Here, we examine the landscape of splicing variants in whole-genome sequencing data from 38,688 individuals in the 100,000 Genomes Project and assess the contribution of non-canonical splicing variants to rare genetic diseases. We use a variant-level constraint metric (the mutability-adjusted proportion of singletons) to identify constrained functional variant classes near exon–intron junctions and at putative splicing branchpoints. To identify new diagnoses for individuals with unsolved rare diseases in the 100,000 Genomes Project, we identified individuals with de novo single-nucleotide variants near exon–intron boundaries and at putative splicing branchpoints in known disease genes. We identified candidate diagnostic variants through manual phenotype matching and confirmed new molecular diagnoses through clinical variant interpretation and functional RNA studies.

Results

We show that near-splice positions and splicing branchpoints are highly constrained by purifying selection and harbour potentially damaging non-coding variants which are amenable to systematic analysis in sequencing data. From 258 de novo splicing variants in known rare disease genes, we identify 35 new likely diagnoses in probands with an unsolved rare disease. To date, we have confirmed a new diagnosis for six individuals, including four in whom RNA studies were performed.

Conclusions

Overall, we demonstrate the clinical value of examining non-canonical splicing variants in individuals with unsolved rare diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01087-x.

Navigating an early career in genomics and data science, written from the perspective of a current PhD student

Rebecca Tooze, a PhD student in the Oxford Clinical Genetics’ group, discusses the importance and application for data science in her field. Using bioinformatic approaches, she analyzes whole-genome sequencing data from patients with craniosynostosis. In this paper, she comments on her current research and her opportunity as an editorial intern with Patterns.