Adaptive Long-Read Sequencing Reveals GGC Repeat Expansion in ZFHX3 Associated with Spinocerebellar Ataxia Type 4

BACKGROUND

Spinocerebellar ataxia type 4 (SCA4) is an autosomal dominant ataxia with invariable sensory neuropathy originally described in a family with Swedish ancestry residing in Utah more than 25 years ago. Despite tight linkage to the 16q22 region, the molecular diagnosis has since remained elusive.

OBJECTIVES

Inspired by pathogenic structural variation implicated in other 16q-ataxias with linkage to the same locus, we revisited the index SCA4 cases from the Utah family using novel technologies to investigate structural variation within the candidate region.

METHODS

We adopted a targeted long-read sequencing approach with adaptive sampling on the Oxford Nanopore Technologies (ONT) platform that enables the detection of segregating structural variants within a genomic region without a priori assumptions about any variant features.

RESULTS

Using this approach, we found a heterozygous (GGC)n repeat expansion in the last coding exon of the zinc finger homeobox 3 (ZFHX3) gene that segregates with disease, ranging between 48 and 57 GGC repeats in affected probands. This finding was replicated in a separate family with SCA4. Furthermore, the estimation of this GGC repeat size in short-read whole genome sequencing (WGS) data of 21,836 individuals recruited to the 100,000 Genomes Project in the UK and our in-house dataset of 11,258 exomes did not reveal any pathogenic repeats, indicating that the variant is ultrarare.

CONCLUSIONS

These findings support the utility of adaptive long-read sequencing as a powerful tool to decipher causative structural variation in unsolved cases of inherited neurological disease. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Radioimmunotherapy versus autologous hematopoietic stem cell transplantation in relapsed/refractory follicular lymphoma: a Fondazione Italiana Linfomi multicenter, randomized, phase III trial

BACKGROUND

Optimal consolidation for young patilents with relapsed/refractory (R/R) follicular lymphoma (FL) remains uncertain in the rituximab era, with an unclear benefit of autologous stem cell transplantation (ASCT). The multicenter, randomized, phase III FLAZ12 (NCT01827605) trial compared anti-CD20 radioimmunotherapy (RIT) with ASCT as consolidation after chemoimmunotherapy, both followed by rituximab maintenance.

PATIENTS AND METHODS

Patients (age 18-65 years) with R/R FL and without significant comorbidities were enrolled and treated with three courses of conventional, investigator-chosen chemoimmunotherapies. Those experiencing at least a partial response were randomized 1 : 1 to ASCT or RIT before CD34+ collection, and all received postconsolidation rituximab maintenance. Progression-free survival (PFS) was the primary endpoint. The target sample size was 210 (105/group).

RESULTS

Between August 2012 and September 2019, of 164 screened patients, 159 were enrolled [median age 57 (interquartile range 49-62) years, 55% male, 57% stage IV, 20% bulky disease]. The study was closed prematurely because of low accrual. Data were analyzed on 8 June 2023, on an intention-to-treat basis, with a 77-month median follow-up from enrollment. Of the 141 patients (89%), 70 were randomized to ASCT and 71 to RIT. The estimated 3-year PFS in both groups was 62% (hazard ratio 1.11, 95% confidence interval 0.69-1.80, P = 0.6662). The 3-year overall survival also was similar between the two groups. Rates of grade ≥3 hematological toxicity were 94% with ASCT versus 46% with RIT (P < 0.001), and grade ≥3 neutropenia occurred in 94% versus 41%, respectively (P < 0.001). Second cancers occurred in nine patients after ASCT and three after radioimmunotherapy (P = 0.189).

CONCLUSIONS

Even if prematurely discontinued, our study did not demonstrate the superiority of ASCT versus RIT. ASCT was more toxic and demanding for patients and health services. Both strategies yielded similar, favorable long-term outcomes, suggesting that consolidation programs milder than ASCT require further investigation in R/R FL.

Rare disease gene association discovery from burden analysis of the 100,000 Genomes Project data

To discover rare disease-gene associations, we developed a gene burden analytical framework and applied it to rare, protein-coding variants from whole genome sequencing of 35,008 cases with rare diseases and their family members recruited to the 100,000 Genomes Project (100KGP). Following in silico triaging of the results, 88 novel associations were identified including 38 with existing experimental evidence. We have published the confirmation of one of these associations, hereditary ataxia with UCHL1 , and independent confirmatory evidence has recently been published for four more. We highlight a further seven compelling associations: hypertrophic cardiomyopathy with DYSF and SLC4A3 where both genes show high/specific heart expression and existing associations to skeletal dystrophies or short QT syndrome respectively; monogenic diabetes with UNC13A with a known role in the regulation of β cells and a mouse model with impaired glucose tolerance; epilepsy with KCNQ1 where a mouse model shows seizures and the existing long QT syndrome association may be linked; early onset Parkinson's disease with RYR1 with existing links to tremor pathophysiology and a mouse model with neurological phenotypes; anterior segment ocular abnormalities associated with POMK showing expression in corneal cells and with a zebrafish model with developmental ocular abnormalities; and cystic kidney disease with COL4A3 showing high renal expression and prior evidence for a digenic or modifying role in renal disease. Confirmation of all 88 associations would lead to potential diagnoses in 456 molecularly undiagnosed cases within the 100KGP, as well as other rare disease patients worldwide, highlighting the clinical impact of a large-scale statistical approach to rare disease gene discovery.

A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 as a significant cause of intellectual disability

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites and underlie several congenital and late-onset disorders. Through a combination of DNA methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using PheWAS in 168,641 individuals from the UK Biobank, identifying 156 significant TRE:trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was linked with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared to controls. With a population prevalence that is at least 5-fold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a significant cause of neurodevelopmental delay.

Normal and pathogenic variation of RFC1 repeat expansions: implications for clinical diagnosis

Cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) is an autosomal recessive neurodegenerative disease, usually caused by biallelic AAGGG repeat expansions in RFC1. In this study, we leveraged whole genome sequencing data from nearly 10 000 individuals recruited within the Genomics England sequencing project to investigate the normal and pathogenic variation of the RFC1 repeat. We identified three novel repeat motifs, AGGGC (n = 6 from five families), AAGGC (n = 2 from one family) and AGAGG (n = 1), associated with CANVAS in the homozygous or compound heterozygous state with the common pathogenic AAGGG expansion. While AAAAG, AAAGGG and AAGAG expansions appear to be benign, we revealed a pathogenic role for large AAAGG repeat configuration expansions (n = 5). Long-read sequencing was used to characterize the entire repeat sequence, and six patients exhibited a pure AGGGC expansion, while the other patients presented complex motifs with AAGGG or AAAGG interruptions. All pathogenic motifs appeared to have arisen from a common haplotype and were predicted to form highly stable G quadruplexes, which have previously been demonstrated to affect gene transcription in other conditions. The assessment of these novel configurations is warranted in CANVAS patients with negative or inconclusive genetic testing. Particular attention should be paid to carriers of compound AAGGG/AAAGG expansions when the AAAGG motif is very large (>500 repeats) or the AAGGG motif is interrupted. Accurate sizing and full sequencing of the satellite repeat with long-read sequencing is recommended in clinically selected cases to enable accurate molecular diagnosis and counsel patients and their families.

Repeat expansions in NOP56 are a cause of spinocerebellar ataxia Type 36 in the British population

Spinocerebellar ataxias form a clinically and genetically heterogeneous group of neurodegenerative disorders characterized by progressive cerebellar ataxia. Their prevalence varies among populations and ethnicities. Spinocerebellar ataxia 36 is caused by a GGCCTG repeat expansion in the first intron of the NOP56 gene and is characterized by late-onset ataxia, sensorineural hearing loss and upper and lower motor neuron signs, including tongue fasciculations. Spinocerebellar ataxia 36 has been described mainly in East Asian and Western European patients and was thought to be absent in the British population. Leveraging novel bioinformatic tools to detect repeat expansions from whole-genome sequencing, we analyse the NOP56 repeat in 1257 British patients with hereditary ataxia and in 7506 unrelated controls. We identify pathogenic repeat expansions in five families (seven patients), representing the first cohort of White British descent patients with spinocerebellar ataxia 36. Employing in silico approaches using whole-genome sequencing data, we found an 87 kb shared haplotype in among the affected individuals from five families around the NOP56 repeat region, although this block was also shared between several controls, suggesting that the repeat arises on a permissive haplotype. Clinically, the patients presented with slowly progressive cerebellar ataxia with a low rate of hearing loss and variable rates of motor neuron impairment. Our findings show that the NOP56 expansion causes ataxia in the British population and that spinocerebellar ataxia 36 can be suspected in patients with a late-onset, slowly progressive ataxia, even without the findings of hearing loss and tongue fasciculation.

POPULATION FREQUENCY OF REPEAT EXPANSIONS INDICATES INCREASED DISEASE PREVALENCE ESTIMATES ACROSS DIFFERENT POPULATIONS

Repeat expansion disorders (REDs) are a devastating group of predominantly neurological diseases. Together they are common, affecting 1 in 3,000 people worldwide with population-specific differences. However, prevalence estimates of REDs are hampered by heterogeneous clinical presentation, variable geographic distributions, and technological limitations leading to under-ascertainment. Here, leveraging whole genome sequencing data from 82,176 individuals from different populations we found an overall carrier frequency of REDs of 1 in 340 individuals. Modelling disease prevalence using genetic data, age at onset and survival, we show that REDs are up to 3-fold more prevalent than currently reported figures. While some REDs are population-specific, e.g. Huntington's disease type 2, most REDs are represented in all broad genetic ancestries, including Africans and Asians, challenging the notion that some REDs are found only in European populations. These results have worldwide implications for local and global health communities in the diagnosis and management of REDs both at local and global levels.

Unexpected frequency of the pathogenic AR CAG repeat expansion in the general population

CAG repeat expansions in exon 1 of the AR gene on the X chromosome cause spinal and bulbar muscular atrophy, a male-specific progressive neuromuscular disorder associated with a variety of extra-neurological symptoms. The disease has a reported male prevalence of approximately 1:30 000 or less, but the AR repeat expansion frequency is unknown. We established a pipeline, which combines the use of the ExpansionHunter tool and visual validation, to detect AR CAG expansion on whole-genome sequencing data, benchmarked it to fragment PCR sizing, and applied it to 74 277 unrelated individuals from four large cohorts. Our pipeline showed sensitivity of 100% [95% confidence interval (CI) 90.8-100%], specificity of 99% (95% CI 94.2-99.7%), and a positive predictive value of 97.4% (95% CI 84.4-99.6%). We found the mutation frequency to be 1:3182 (95% CI 1:2309-1:4386, n = 117 734) X chromosomes-10 times more frequent than the reported disease prevalence. Modelling using the novel mutation frequency led to estimate disease prevalence of 1:6887 males, more than four times more frequent than the reported disease prevalence. This discrepancy is possibly due to underdiagnosis of this neuromuscular condition, reduced penetrance, and/or pleomorphic clinical manifestations.

Truncating Variants in RFC1 in Cerebellar Ataxia, Neuropathy, and Vestibular Areflexia Syndrome

BACKGROUND AND OBJECTIVE

Cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS) is an autosomal recessive neurodegenerative disease characterized by adult-onset and slowly progressive sensory neuropathy, cerebellar dysfunction, and vestibular impairment. In most cases, the disease is caused by biallelic (AAGGG)n repeat expansions in the second intron of the replication factor complex subunit 1 (RFC1). However, a small number of cases with typical CANVAS do not carry the common biallelic repeat expansion. The objective of this study was to expand the genotypic spectrum of CANVAS by identifying sequence variants in RFC1-coding region associated with this condition.

METHODS

Fifteen individuals diagnosed with CANVAS and carrying only 1 heterozygous (AAGGG)n expansion in RFC1 underwent whole-genome sequencing or whole-exome sequencing to test for the presence of a second variant in RFC1 or other unrelated gene. To assess the effect of truncating variants on RFC1 expression, we tested the level of RFC1 transcript and protein on patients' derived cell lines.

RESULTS

We identified 7 patients from 5 unrelated families with clinically defined CANVAS carrying a heterozygous (AAGGG)n expansion together with a second truncating variant in trans in RFC1, which included the following: c.1267C>T (p.Arg423Ter), c.1739_1740del (p.Lys580SerfsTer9), c.2191del (p.Gly731GlufsTer6), and c.2876del (p.Pro959GlnfsTer24). Patient fibroblasts containing the c.1267C>T (p.Arg423Ter) or c.2876del (p.Pro959GlnfsTer24) variants demonstrated nonsense-mediated mRNA decay and reduced RFC1 transcript and protein.

DISCUSSION

Our report expands the genotype spectrum of RFC1 disease. Full RFC1 sequencing is recommended in cases affected by typical CANVAS and carrying monoallelic (AAGGG)n expansions. In addition, it sheds further light on the pathogenesis of RFC1 CANVAS because it supports the existence of a loss-of-function mechanism underlying this complex neurodegenerative condition.

Functional genomics provide key insights to improve the diagnostic yield of hereditary ataxia

Improvements in functional genomic annotation have led to a critical mass of neurogenetic discoveries. This is exemplified in hereditary ataxia, a heterogeneous group of disorders characterised by incoordination from cerebellar dysfunction. Associated pathogenic variants in more than 300 genes have been described, leading to a detailed genetic classification partitioned by age-of-onset. Despite these advances, up to 75% of patients with ataxia remain molecularly undiagnosed even following whole genome sequencing, as exemplified in the 100,000 Genomes Project. This study aimed to understand whether we can improve our knowledge of the genetic architecture of hereditary ataxia by leveraging functional genomic annotations, and as a result, generate insights and strategies that raise the diagnostic yield. To achieve these aims, we used publicly-available multi-omics data to generate 294 genic features, capturing information relating to a gene's structure, genetic variation, tissue-specific, cell-type-specific and temporal expression, as well as protein products of a gene. We studied these features across genes typically causing childhood-onset, adult-onset or both types of disease first individually, then collectively. This led to the generation of testable hypotheses which we investigated using whole genome sequencing data from up to 2,182 individuals presenting with ataxia and 6,658 non-neurological probands recruited in the 100,000 Genomes Project. Using this approach, we demonstrated a high short tandem repeat (STR) density within childhood-onset genes suggesting that we may be missing pathogenic repeat expansions within this cohort. This was verified in both childhood- and adult-onset ataxia patients from the 100,000 Genomes Project who were unexpectedly found to have a trend for higher repeat sizes even at naturally-occurring STRs within known ataxia genes, implying a role for STRs in pathogenesis. Using unsupervised analysis, we found significant similarities in genomic annotation across the gene panels, which suggested adult- and childhood-onset patients should be screened using a common diagnostic gene set. We tested this within the 100,000 Genomes Project by assessing the burden of pathogenic variants among childhood-onset genes in adult-onset patients and vice versa. This demonstrated a significantly higher burden of rare, potentially pathogenic variants in conventional childhood-onset genes among individuals with adult-onset ataxia. Our analysis has implications for the current clinical practice in genetic testing for hereditary ataxia. We suggest that the diagnostic rate for hereditary ataxia could be increased by removing the age-of-onset partition, and through a modified screening for repeat expansions in naturally-occurring STRs within known ataxia-associated genes, in effect treating these regions as candidate pathogenic loci.

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 .

Genome sequencing reveals underdiagnosis of primary ciliary dyskinesia in bronchiectasis

BACKGROUND

Bronchiectasis can result from infectious, genetic, immunological and allergic causes. 60-80% of cases are idiopathic, but a well-recognised genetic cause is the motile ciliopathy, primary ciliary dyskinesia (PCD). Diagnosis of PCD has management implications including addressing comorbidities, implementing genetic and fertility counselling and future access to PCD-specific treatments. Diagnostic testing can be complex; however, PCD genetic testing is moving rapidly from research into clinical diagnostics and would confirm the cause of bronchiectasis.

METHODS

This observational study used genetic data from severe bronchiectasis patients recruited to the UK 100,000 Genomes Project and patients referred for gene panel testing within a tertiary respiratory hospital. Patients referred for genetic testing due to clinical suspicion of PCD were excluded from both analyses. Data were accessed from the British Thoracic Society audit, to investigate whether motile ciliopathies are underdiagnosed in people with bronchiectasis in the UK.

RESULTS

Pathogenic or likely pathogenic variants were identified in motile ciliopathy genes in 17 (12%) out of 142 individuals by whole-genome sequencing. Similarly, in a single centre with access to pathological diagnostic facilities, 5-10% of patients received a PCD diagnosis by gene panel, often linked to normal/inconclusive nasal nitric oxide and cilia functional test results. In 4898 audited patients with bronchiectasis, <2% were tested for PCD and <1% received genetic testing.

CONCLUSIONS

PCD is underdiagnosed as a cause of bronchiectasis. Increased uptake of genetic testing may help to identify bronchiectasis due to motile ciliopathies and ensure appropriate management.

Heterozygous UCHL1 loss-of-function variants cause a neurodegenerative disorder with spasticity, ataxia, neuropathy, and optic atrophy

PURPOSE

Biallelic variants in UCHL1 have been associated with a progressive early-onset neurodegenerative disorder, autosomal recessive spastic paraplegia type 79. In this study, we investigated heterozygous UCHL1 variants on the basis of results from cohort-based burden analyses.

METHODS

Gene-burden analyses were performed on exome and genome data of independent cohorts of patients with hereditary ataxia and spastic paraplegia from Germany and the United Kingdom in a total of 3169 patients and 33,141 controls. Clinical data of affected individuals and additional independent families were collected and evaluated. Patients' fibroblasts were used to perform mass spectrometry-based proteomics.

RESULTS

UCHL1 was prioritized in both independent cohorts as a candidate gene for an autosomal dominant disorder. We identified a total of 34 cases from 18 unrelated families, carrying 13 heterozygous loss-of-function variants (15 families) and an inframe insertion (3 families). Affected individuals mainly presented with spasticity (24/31), ataxia (28/31), neuropathy (11/21), and optic atrophy (9/17). The mass spectrometry-based proteomics showed approximately 50% reduction of UCHL1 expression in patients' fibroblasts.

CONCLUSION

Our bioinformatic analysis, in-depth clinical and genetic workup, and functional studies established haploinsufficiency of UCHL1 as a novel disease mechanism in spastic ataxia.

REViewer: haplotype-resolved visualization of read alignments in and around tandem repeats

BACKGROUND

Expansions of short tandem repeats are the cause of many neurogenetic disorders including familial amyotrophic lateral sclerosis, Huntington disease, and many others. Multiple methods have been recently developed that can identify repeat expansions in whole genome or exome sequencing data. Despite the widely recognized need for visual assessment of variant calls in clinical settings, current computational tools lack the ability to produce such visualizations for repeat expansions. Expanded repeats are difficult to visualize because they correspond to large insertions relative to the reference genome and involve many misaligning and ambiguously aligning reads.

RESULTS

We implemented REViewer, a computational method for visualization of sequencing data in genomic regions containing long repeat expansions and FlipBook, a companion image viewer designed for manual curation of large collections of REViewer images. To generate a read pileup, REViewer reconstructs local haplotype sequences and distributes reads to these haplotypes in a way that is most consistent with the fragment lengths and evenness of read coverage. To create appropriate training materials for onboarding new users, we performed a concordance study involving 12 scientists involved in short tandem repeat research. We used the results of this study to create a user guide that describes the basic principles of using REViewer as well as a guide to the typical features of read pileups that correspond to low confidence repeat genotype calls. Additionally, we demonstrated that REViewer can be used to annotate clinically relevant repeat interruptions by comparing visual assessment results of 44 FMR1 repeat alleles with the results of triplet repeat primed PCR. For 38 of these alleles, the results of visual assessment were consistent with triplet repeat primed PCR.

CONCLUSIONS

Read pileup plots generated by REViewer offer an intuitive way to visualize sequencing data in regions containing long repeat expansions. Laboratories can use REViewer and FlipBook to assess the quality of repeat genotype calls as well as to visually detect interruptions or other imperfections in the repeat sequence and the surrounding flanking regions. REViewer and FlipBook are available under open-source licenses at https://github.com/illumina/REViewer and https://github.com/broadinstitute/flipbook respectively.

022 Functional genomics and transcriptomics further characterise and potentially improve diagnostic yield of hereditary ataxias

zhongbo

chen@ucl.ac.uk|ABN Bursary

36

Whole genome sequencing for the diagnosis of neurological repeat expansion disorders in the UK: a retrospective diagnostic accuracy and prospective clinical validation study

BACKGROUND

Repeat expansion disorders affect about 1 in 3000 individuals and are clinically heterogeneous diseases caused by expansions of short tandem DNA repeats. Genetic testing is often locus-specific, resulting in underdiagnosis of people who have atypical clinical presentations, especially in paediatric patients without a previous positive family history. Whole genome sequencing is increasingly used as a first-line test for other rare genetic disorders, and we aimed to assess its performance in the diagnosis of patients with neurological repeat expansion disorders.

METHODS

We retrospectively assessed the diagnostic accuracy of whole genome sequencing to detect the most common repeat expansion loci associated with neurological outcomes (AR, ATN1, ATXN1, ATXN2, ATXN3, ATXN7, C9orf72, CACNA1A, DMPK, FMR1, FXN, HTT, and TBP) using samples obtained within the National Health Service in England from patients who were suspected of having neurological disorders; previous PCR test results were used as the reference standard. The clinical accuracy of whole genome sequencing to detect repeat expansions was prospectively examined in previously genetically tested and undiagnosed patients recruited in 2013-17 to the 100 000 Genomes Project in the UK, who were suspected of having a genetic neurological disorder (familial or early-onset forms of ataxia, neuropathy, spastic paraplegia, dementia, motor neuron disease, parkinsonian movement disorders, intellectual disability, or neuromuscular disorders). If a repeat expansion call was made using whole genome sequencing, PCR was used to confirm the result.

FINDINGS

The diagnostic accuracy of whole genome sequencing to detect repeat expansions was evaluated against 793 PCR tests previously performed within the NHS from 404 patients. Whole genome sequencing correctly classified 215 of 221 expanded alleles and 1316 of 1321 non-expanded alleles, showing 97·3% sensitivity (95% CI 94·2-99·0) and 99·6% specificity (99·1-99·9) across the 13 disease-associated loci when compared with PCR test results. In samples from 11 631 patients in the 100 000 Genomes Project, whole genome sequencing identified 81 repeat expansions, which were also tested by PCR: 68 were confirmed as repeat expansions in the full pathogenic range, 11 were non-pathogenic intermediate expansions or premutations, and two were non-expanded repeats (16% false discovery rate).

INTERPRETATION

In our study, whole genome sequencing for the detection of repeat expansions showed high sensitivity and specificity, and it led to identification of neurological repeat expansion disorders in previously undiagnosed patients. These findings support implementation of whole genome sequencing in clinical laboratories for diagnosis of patients who have a neurological presentation consistent with a repeat expansion disorder.

FUNDING

Medical Research Council, Department of Health and Social Care, National Health Service England, National Institute for Health Research, and Illumina.

Whole genome sequencing for the diagnosis of neurological repeat expansion disorders in the UK: a retrospective diagnostic accuracy and prospective clinical validation study

BACKGROUND

Repeat expansion disorders affect about 1 in 3000 individuals and are clinically heterogeneous diseases caused by expansions of short tandem DNA repeats. Genetic testing is often locus-specific, resulting in underdiagnosis of people who have atypical clinical presentations, especially in paediatric patients without a previous positive family history. Whole genome sequencing is increasingly used as a first-line test for other rare genetic disorders, and we aimed to assess its performance in the diagnosis of patients with neurological repeat expansion disorders.

METHODS

We retrospectively assessed the diagnostic accuracy of whole genome sequencing to detect the most common repeat expansion loci associated with neurological outcomes (AR, ATN1, ATXN1, ATXN2, ATXN3, ATXN7, C9orf72, CACNA1A, DMPK, FMR1, FXN, HTT, and TBP) using samples obtained within the National Health Service in England from patients who were suspected of having neurological disorders; previous PCR test results were used as the reference standard. The clinical accuracy of whole genome sequencing to detect repeat expansions was prospectively examined in previously genetically tested and undiagnosed patients recruited in 2013-17 to the 100 000 Genomes Project in the UK, who were suspected of having a genetic neurological disorder (familial or early-onset forms of ataxia, neuropathy, spastic paraplegia, dementia, motor neuron disease, parkinsonian movement disorders, intellectual disability, or neuromuscular disorders). If a repeat expansion call was made using whole genome sequencing, PCR was used to confirm the result.

FINDINGS

The diagnostic accuracy of whole genome sequencing to detect repeat expansions was evaluated against 793 PCR tests previously performed within the NHS from 404 patients. Whole genome sequencing correctly classified 215 of 221 expanded alleles and 1316 of 1321 non-expanded alleles, showing 97·3% sensitivity (95% CI 94·2-99·0) and 99·6% specificity (99·1-99·9) across the 13 disease-associated loci when compared with PCR test results. In samples from 11 631 patients in the 100 000 Genomes Project, whole genome sequencing identified 81 repeat expansions, which were also tested by PCR: 68 were confirmed as repeat expansions in the full pathogenic range, 11 were non-pathogenic intermediate expansions or premutations, and two were non-expanded repeats (16% false discovery rate).

INTERPRETATION

In our study, whole genome sequencing for the detection of repeat expansions showed high sensitivity and specificity, and it led to identification of neurological repeat expansion disorders in previously undiagnosed patients. These findings support implementation of whole genome sequencing in clinical laboratories for diagnosis of patients who have a neurological presentation consistent with a repeat expansion disorder.

FUNDING

Medical Research Council, Department of Health and Social Care, National Health Service England, National Institute for Health Research, and Illumina.

Structural-Equation-Modelling (SEM) to analyze climatic factor's role on COVID-19 spreading

Purpose

Climate seems to influence the COVID-19 spreading, but the results of the published studies are conflicting. Aim of this study was to perform a world-wide investigation to analyze the role of all the main climatic factors (CF), trying to identify the causes that led to the discrepancy of the results.

Methods & Materials

134,871 data (from 209 countries) were used for the analysis. These were extrapolated from an initial data-set of 1.200.000 data. To avoid biases present in most of the previously studies, a set of specific requirements was adopted: long observation period (16 weeks),

• the use of a relative time scale to synchronize the beginning of the outbreak among the countries,

• multiple data collection points (up to 4 cities/per country) to overcome the problem of climate variability within a country,

• the use of an appropriate technique to test the relationships among interdependent variables,

• the use of a lag-period to compensate the shift between the infection exposure and the diagnosis’ confirmation.

Data's analysis was performed with SEM, a flexible statistical technique for modeling causal chain of effects simultaneously. Using hypothesis-testing, this technique examines the relationships between observed variables and latent variables, in turn linked to observed variables, their indicators. With this statistical model it was possible to consider the integrated effects of all the CF on COVID-19 and, at the same time, to investigate the effects of population density (PD) too.

Results

The results of the analysis showed that both climate and population density significantly influence the spread of COVID-19 (p<0.001; p<0.01, respectively). Overall, climate outweighs population density (path coefficients: climate vs incidence=0.18, climate vs prevalence=0.11, PD vs incidence=0.04, PD vs prevalence=0.05). Among the climatic factors, irradiation plays the most relevant role, with a factor-loading of -0.77, followed by temperature (-0.56), humidity (0.52), precipitation (0.44), and pressure (0.073); for all p<0.001. Fit indices demonstrated a good fit of the model (determination-coefficient=0.826, Root-Mean-Square-Error-of-Approximation=0.088, Standardized-Root-Mean-Square-Residual=0.078).

Conclusion

This study demonstrates that CF significantly influence the spread of SARS-CoV-2. However, demographic factors, together with other determinants, can affect the transmission, overcoming the protective effect of climate, where favourable.

100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care - Preliminary Report

BACKGROUND

The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection.

METHODS

We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis.

RESULTS

Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives.

CONCLUSIONS

Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.).

Myoclonic status epilepticus and cerebellar hypoplasia associated with a novel variant in the GRIA3 gene

AMPA-type glutamate receptors (AMPARs) are postsynaptic ionotropic receptors which mediate fast excitatory currents. AMPARs have a heterotetrameric structure, variably composed by the four subunits GluA1-4 which are encoded by genes GRIA1-4. Increasing evidence support the role of pathogenic variants in GRIA1-4 genes as causative for syndromic intellectual disability (ID). We report an Italian pedigree where some male individuals share ID, seizures and facial dysmorphisms. The index subject was referred for severe ID, myoclonic seizures, cerebellar signs and short stature. Whole exome sequencing identified a novel variant in GRIA3, c.2360A > G, p.(Glu787Gly). The GRIA3 gene maps to chromosome Xq25 and the c.2360A > G variant was transmitted by his healthy mother. Subsequent analysis in the family showed a segregation pattern compatible with the causative role of this variant, further supported by preliminary functional insights. We provide a detailed description of the clinical evolution of the index subjects and stress the relevance of myoclonic seizures and cerebellar syndrome as cardinal features of his presentation.

Use of whole genome sequencing to determine genetic basis of suspected mitochondrial disorders: cohort study

OBJECTIVE

To determine whether whole genome sequencing can be used to define the molecular basis of suspected mitochondrial disease.

DESIGN

Cohort study.

SETTING

National Health Service, England, including secondary and tertiary care.

PARTICIPANTS

345 patients with suspected mitochondrial disorders recruited to the 100 000 Genomes Project in England between 2015 and 2018.

INTERVENTION

Short read whole genome sequencing was performed. Nuclear variants were prioritised on the basis of gene panels chosen according to phenotypes, ClinVar pathogenic/likely pathogenic variants, and the top 10 prioritised variants from Exomiser. Mitochondrial DNA variants were called using an in-house pipeline and compared with a list of pathogenic variants. Copy number variants and short tandem repeats for 13 neurological disorders were also analysed. American College of Medical Genetics guidelines were followed for classification of variants.

MAIN OUTCOME MEASURE

Definite or probable genetic diagnosis.

RESULTS

A definite or probable genetic diagnosis was identified in 98/319 (31%) families, with an additional 6 (2%) possible diagnoses. Fourteen of the diagnoses (4% of the 319 families) explained only part of the clinical features. A total of 95 different genes were implicated. Of 104 families given a diagnosis, 39 (38%) had a mitochondrial diagnosis and 65 (63%) had a non-mitochondrial diagnosis.

CONCLUSION

Whole genome sequencing is a useful diagnostic test in patients with suspected mitochondrial disorders, yielding a diagnosis in a further 31% after exclusion of common causes. Most diagnoses were non-mitochondrial disorders and included developmental disorders with intellectual disability, epileptic encephalopathies, other metabolic disorders, cardiomyopathies, and leukodystrophies. These would have been missed if a targeted approach was taken, and some have specific treatments.

Bi-allelic variants in SPATA5L1 lead to intellectual disability, spastic-dystonic cerebral palsy, epilepsy, and hearing loss

Spermatogenesis-associated 5 like 1 (SPATA5L1) represents an orphan gene encoding a protein of unknown function. We report 28 bi-allelic variants in SPATA5L1 associated with sensorineural hearing loss in 47 individuals from 28 (26 unrelated) families. In addition, 25/47 affected individuals (53%) presented with microcephaly, developmental delay/intellectual disability, cerebral palsy, and/or epilepsy. Modeling indicated damaging effect of variants on the protein, largely via destabilizing effects on protein domains. Brain imaging revealed diminished cerebral volume, thin corpus callosum, and periventricular leukomalacia, and quantitative volumetry demonstrated significantly diminished white matter volumes in several individuals. Immunofluorescent imaging in rat hippocampal neurons revealed localization of Spata5l1 in neuronal and glial cell nuclei and more prominent expression in neurons. In the rodent inner ear, Spata5l1 is expressed in the neurosensory hair cells and inner ear supporting cells. Transcriptomic analysis performed with fibroblasts from affected individuals was able to distinguish affected from controls by principal components. Analysis of differentially expressed genes and networks suggested a role for SPATA5L1 in cell surface adhesion receptor function, intracellular focal adhesions, and DNA replication and mitosis. Collectively, our results indicate that bi-allelic SPATA5L1 variants lead to a human disease characterized by sensorineural hearing loss (SNHL) with or without a nonprogressive mixed neurodevelopmental phenotype.

Production and persistence of specific antibodies in COVID-19 patients with hematologic malignancies: role of rituximab

The ability of patients with hematologic malignancies (HM) to develop an effective humoral immune response after COVID-19 is unknown. A prospective study was performed to monitor the immune response to SARS-CoV-2 of patients with follicular lymphoma (FL), diffuse large B-cell lymphoma (DLBCL), chronic lymphoproliferative disorders (CLD), multiple myeloma (MM), or myelodysplastic/myeloproliferative syndromes (MDS/MPN). Antibody (Ab) levels to the SARS-CoV-2 nucleocapsid (N) and spike (S) protein were measured at +1, +3, +6 months after nasal swabs became PCR-negative. Forty-five patients (9 FL, 8 DLBCL, 8 CLD, 10 MM, 10 MDS/MPS) and 18 controls were studied. Mean anti-N and anti-S-Ab levels were similar between HM patients and controls, and shared the same behavior, with anti-N Ab levels declining at +6 months and anti-S-Ab remaining stable. Seroconversion rates were lower in HM patients than in controls. In lymphoma patients mean Ab levels and seroconversion rates were lower than in other HM patients, primarily because all nine patients who had received rituximab within 6 months before COVID-19 failed to produce anti-N and anti-S-Ab. Only one patient requiring hematological treatment after COVID-19 lost seropositivity after 6 months. No reinfections were observed. These results may inform vaccination policies and clinical management of HM patients.

Biallelic variants in TSPOAP1, encoding the active-zone protein RIMBP1, cause autosomal recessive dystonia

Dystonia is a debilitating hyperkinetic movement disorder, which can be transmitted as a monogenic trait. Here, we describe homozygous frameshift, nonsense, and missense variants in TSPOAP1, which encodes the active-zone RIM-binding protein 1 (RIMBP1), as a genetic cause of autosomal recessive dystonia in 7 subjects from 3 unrelated families. Subjects carrying loss-of-function variants presented with juvenile-onset progressive generalized dystonia, associated with intellectual disability and cerebellar atrophy. Conversely, subjects carrying a pathogenic missense variant (p.Gly1808Ser) presented with isolated adult-onset focal dystonia. In mice, complete loss of RIMBP1, known to reduce neurotransmission, led to motor abnormalities reminiscent of dystonia, decreased Purkinje cell dendritic arborization, and reduced numbers of cerebellar synapses. In vitro analysis of the p.Gly1808Ser variant showed larger spike-evoked calcium transients and enhanced neurotransmission, suggesting that RIMBP1-linked dystonia can be caused by either reduced or enhanced rates of spike-evoked release in relevant neural networks. Our findings establish a direct link between dysfunction of the presynaptic active zone and dystonia and highlight the critical role played by well-balanced neurotransmission in motor control and disease pathogenesis.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.

Pathogenic Huntingtin Repeat Expansions in Patients with Frontotemporal Dementia and Amyotrophic Lateral Sclerosis

We examined the role of repeat expansions in the pathogenesis of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS) by analyzing whole-genome sequence data from 2,442 FTD/ALS patients, 2,599 Lewy body dementia (LBD) patients, and 3,158 neurologically healthy subjects. Pathogenic expansions (range, 40-64 CAG repeats) in the huntingtin (HTT) gene were found in three (0.12%) patients diagnosed with pure FTD/ALS syndromes but were not present in the LBD or healthy cohorts. We replicated our findings in an independent collection of 3,674 FTD/ALS patients. Postmortem evaluations of two patients revealed the classical TDP-43 pathology of FTD/ALS, as well as huntingtin-positive, ubiquitin-positive aggregates in the frontal cortex. The neostriatal atrophy that pathologically defines Huntington's disease was absent in both cases. Our findings reveal an etiological relationship between HTT repeat expansions and FTD/ALS syndromes and indicate that genetic screening of FTD/ALS patients for HTT repeat expansions should be considered.

Loss-of-Function Variants in HOPS Complex Genes VPS16 and VPS41 Cause Early Onset Dystonia Associated with Lysosomal Abnormalities

OBJECTIVES

The majority of people with suspected genetic dystonia remain undiagnosed after maximal investigation, implying that a number of causative genes have not yet been recognized. We aimed to investigate this paucity of diagnoses.

METHODS

We undertook weighted burden analysis of whole-exome sequencing (WES) data from 138 individuals with unresolved generalized dystonia of suspected genetic etiology, followed by additional case-finding from international databases, first for the gene implicated by the burden analysis (VPS16), and then for other functionally related genes. Electron microscopy was performed on patient-derived cells.

RESULTS

Analysis revealed a significant burden for VPS16 (Fisher's exact test p value, 6.9 × 10⁹ ). VPS16 encodes a subunit of the homotypic fusion and vacuole protein sorting (HOPS) complex, which plays a key role in autophagosome-lysosome fusion. A total of 18 individuals harboring heterozygous loss-of-function VPS16 variants, and one with a microdeletion, were identified. These individuals experienced early onset progressive dystonia with predominant cervical, bulbar, orofacial, and upper limb involvement. Some patients had a more complex phenotype with additional neuropsychiatric and/or developmental comorbidities. We also identified biallelic loss-of-function variants in VPS41, another HOPS-complex encoding gene, in an individual with infantile-onset generalized dystonia. Electron microscopy of patient-derived lymphocytes and fibroblasts from both patients with VPS16 and VPS41 showed vacuolar abnormalities suggestive of impaired lysosomal function.

INTERPRETATION

Our study strongly supports a role for HOPS complex dysfunction in the pathogenesis of dystonia, although variants in different subunits display different phenotypic and inheritance characteristics. ANN NEUROL 2020;88:867-877.

Neuronal intranuclear inclusion disease is genetically heterogeneous

Neuronal intranuclear inclusion disease (NIID) is a clinically heterogeneous neurodegenerative condition characterized by pathological intranuclear eosinophilic inclusions. A CGG repeat expansion in NOTCH2NLC was recently identified to be associated with NIID in patients of Japanese descent. We screened pathologically confirmed European NIID, cases of neurodegenerative disease with intranuclear inclusions and applied in silico-based screening using whole-genome sequencing data from 20 536 participants in the 100 000 Genomes Project. We identified a single European case harbouring the pathogenic repeat expansion with a distinct haplotype structure. Thus, we propose new diagnostic criteria as European NIID represents a distinct disease entity from East Asian cases.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions

SUMMARY

We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci.

AVAILABILITY AND IMPLEMENTATION

ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

De Novo SOX6 Variants Cause a Neurodevelopmental Syndrome Associated with ADHD, Craniosynostosis, and Osteochondromas

SOX6 belongs to a family of 20 SRY-related HMG-box-containing (SOX) genes that encode transcription factors controlling cell fate and differentiation in many developmental and adult processes. For SOX6, these processes include, but are not limited to, neurogenesis and skeletogenesis. Variants in half of the SOX genes have been shown to cause severe developmental and adult syndromes, referred to as SOXopathies. We here provide evidence that SOX6 variants also cause a SOXopathy. Using clinical and genetic data, we identify 19 individuals harboring various types of SOX6 alterations and exhibiting developmental delay and/or intellectual disability; the individuals are from 17 unrelated families. Additional, inconstant features include attention-deficit/hyperactivity disorder (ADHD), autism, mild facial dysmorphism, craniosynostosis, and multiple osteochondromas. All variants are heterozygous. Fourteen are de novo, one is inherited from a mosaic father, and four offspring from two families have a paternally inherited variant. Intragenic microdeletions, balanced structural rearrangements, frameshifts, and nonsense variants are predicted to inactivate the SOX6 variant allele. Four missense variants occur in residues and protein regions highly conserved evolutionarily. These variants are not detected in the gnomAD control cohort, and the amino acid substitutions are predicted to be damaging. Two of these variants are located in the HMG domain and abolish SOX6 transcriptional activity in vitro. No clear genotype-phenotype correlations are found. Taken together, these findings concur that SOX6 haploinsufficiency leads to a neurodevelopmental SOXopathy that often includes ADHD and abnormal skeletal and other features.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions

SUMMARY

We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci.

AVAILABILITY AND IMPLEMENTATION

ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.