Rare variant associations with plasma protein levels in the UK Biobank.

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

Rare variant contribution to human disease in 281,104 UK Biobank exomes.

Genome-wide association studies have uncovered thousands of common variants associated with human disease, but the contribution of rare variants to common disease remains relatively unexplored. The UK Biobank contains detailed phenotypic data linked to medical records for approximately 500,000 participants, offering an unprecedented opportunity to evaluate the effect of rare variation on a broad collection of traits1,2. Here we study the relationships between rare protein-coding variants and 17,361 binary and 1,419 quantitative phenotypes using exome sequencing data from 269,171 UK Biobank participants of European ancestry. Gene-based collapsing analyses revealed 1,703 statistically significant gene-phenotype associations for binary traits, with a median odds ratio of 12.4. Furthermore, 83% of these associations were undetectable via single-variant association tests, emphasizing the power of gene-based collapsing analysis in the setting of high allelic heterogeneity. Gene-phenotype associations were also significantly enriched for loss-of-function-mediated traits and approved drug targets. Finally, we performed ancestry-specific and pan-ancestry collapsing analyses using exome sequencing data from 11,933 UK Biobank participants of African, East Asian or South Asian ancestry. Our results highlight a significant contribution of rare variants to common disease. Summary statistics are publicly available through an interactive portal ( http://azphewas.com/ ).

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank.

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders.

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

Whole-genome sequencing of patients with rare diseases in a national health system.

Most patients with rare diseases do not receive a molecular diagnosis and the aetiological variants and causative genes for more than half such disorders remain to be discovered1. Here we used whole-genome sequencing (WGS) in a national health system to streamline diagnosis and to discover unknown aetiological variants in the coding and non-coding regions of the genome. We generated WGS data for 13,037 participants, of whom 9,802 had a rare disease, and provided a genetic diagnosis to 1,138 of the 7,065 extensively phenotyped participants. We identified 95 Mendelian associations between genes and rare diseases, of which 11 have been discovered since 2015 and at least 79 are confirmed to be aetiological. By generating WGS data of UK Biobank participants2, we found that rare alleles can explain the presence of some individuals in the tails of a quantitative trait for red blood cells. Finally, we identified four novel non-coding variants that cause disease through the disruption of transcription of ARPC1B, GATA1, LRBA and MPL. Our study demonstrates a synergy by using WGS for diagnosis and aetiological discovery in routine healthcare.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants.

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Uncovering variable neoplasms between ATM protein-truncating and common missense variants using 394 694 UK Biobank exomes.

As an essential regulator of DNA damage, ataxia-telangiectasia mutated (ATM) gene has been widely studied in oncology. However, the independent effects of ATM missense variants and protein-truncating variants (PTVs) on neoplasms have not been heavily studied. Whole-exome sequencing data and the clinical health records of 394,694 UK Biobank European participants were used in this analysis. We mined genetic associations from gene-level and variant-level phenome-wide association studies, and conducted a variant-level conditional association study to test whether the effects of ATM missense variants on neoplasms were independent of ATM PTV carrier status. The gene-level PTV collapsing analysis was consistent with established ATM PTV literature showing that the aggregated impact of 286 ATM PTVs significantly (p < 2 × 10-9 ) associated with 31 malignant neoplasm phenotypes. Of 773 distinct protein-coding variants in ATM, three individual missense variants significantly (p < 2 × 10-9 ) associated with nine phenotypes. Remarkably, although the nine phenotypes were tumor-related, none overlapped the established ATM PTV-linked malignancies. A subsequent conditional analysis identified that the missense signals were acting independently of the known clinically relevant ATM PTVs.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia.

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Cav2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca2+ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Cav2.2 in normal human neurodevelopment.

MED27, SLC6A7, and MPPE1 Variants in a Complex Neurodevelopmental Disorder with Severe Dystonia.

BACKGROUND: Despite advances in next generation sequencing technologies, the identification of variants of uncertain significance (VUS) can often hinder definitive diagnosis in patients with complex neurodevelopmental disorders. OBJECTIVE: The objective of this study was to identify and characterize the underlying cause of disease in a family with two children with severe developmental delay associated with generalized dystonia and episodic status dystonicus, chorea, epilepsy, and cataracts. METHODS: Candidate genes identified by autozygosity mapping and whole-exome sequencing were characterized using cellular and vertebrate model systems. RESULTS: Homozygous variants were found in three candidate genes: MED27, SLC6A7, and MPPE1. Although the patients had features of MED27-related disorder, the SLC6A7 and MPPE1 variants were functionally investigated. SLC6A7 variant in vitro overexpression caused decreased proline transport as a result of reduced cell-surface expression, and zebrafish knockdown of slc6a7 exhibited developmental delay and fragile motor neuron morphology that could not be rescued by L-proline transporter-G396S RNA. Lastly, patient fibroblasts displayed reduced cell-surface expression of glycophosphatidylinositol-anchored proteins linked to MPPE1 dysfunction. CONCLUSIONS: We report a family harboring a homozygous MED27 variant with additional loss-of-function SLC6A7 and MPPE1 gene variants, which potentially contribute to a blended phenotype caused by multilocus pathogenic variants. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

De Novo Truncating Mutations in WASF1 Cause Intellectual Disability with Seizures.

Next-generation sequencing has been invaluable in the elucidation of the genetic etiology of many subtypes of intellectual disability in recent years. Here, using exome sequencing and whole-genome sequencing, we identified three de novo truncating mutations in WAS protein family member 1 (WASF1) in five unrelated individuals with moderate to profound intellectual disability with autistic features and seizures. WASF1, also known as WAVE1, is part of the WAVE complex and acts as a mediator between Rac-GTPase and actin to induce actin polymerization. The three mutations connected by Matchmaker Exchange were c.1516C>T (p.Arg506Ter), which occurs in three unrelated individuals, c.1558C>T (p.Gln520Ter), and c.1482delinsGCCAGG (p.Ile494MetfsTer23). All three variants are predicted to partially or fully disrupt the C-terminal actin-binding WCA domain. Functional studies using fibroblast cells from two affected individuals with the c.1516C>T mutation showed a truncated WASF1 and a defect in actin remodeling. This study provides evidence that de novo heterozygous mutations in WASF1 cause a rare form of intellectual disability.

Comprehensive Cancer-Predisposition Gene Testing in an Adult Multiple Primary Tumor Series Shows a Broad Range of Deleterious Variants and Atypical Tumor Phenotypes.

Multiple primary tumors (MPTs) affect a substantial proportion of cancer survivors and can result from various causes, including inherited predisposition. Currently, germline genetic testing of MPT-affected individuals for variants in cancer-predisposition genes (CPGs) is mostly targeted by tumor type. We ascertained pre-assessed MPT individuals (with at least two primary tumors by age 60 years or at least three by 70 years) from genetics centers and performed whole-genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment and molecular investigations, pathogenic variants in moderate- and high-risk CPGs were detected in 67/440 (15.2%) probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies, including low-frequency structural variants (6/440 [1.4%] probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant), the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice-site CPG variants and at least one discordant tumor type was significantly higher than in a control population (χ2 = 43.642; p ≤ 0.0001). 2/67 (3%) probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Together with variant detection rates from a previous series of similarly ascertained MPT-affected individuals, the present results suggest that first-line comprehensive CPG analysis in an MPT cohort referred to clinical genetics services would detect a deleterious variant in about a third of individuals.

Ceramide synthase TLCD3B is a novel gene associated with human recessive retinal dystrophy.

PURPOSE: Previous studies suggest that ceramide is a proapoptotic lipid as high levels of ceramides can lead to apoptosis of neuronal cells, including photoreceptors. However, no pathogenic variant in ceramide synthases has been identified in human patients and knockout of various ceramide synthases in mice has not led to photoreceptor degeneration. METHODS: Exome sequencing was used to identify candidate disease genes in patients with vision loss as confirmed by standard evaluation methods, including electroretinography (ERG) and optical coherence tomography. The vision loss phenotype in mice was evaluated by ERG and histological analyses. RESULTS: Here we have identified four patients with cone-rod dystrophy or maculopathy from three families carrying pathogenic variants in TLCD3B. Consistent with the phenotype observed in patients, the Tlcd3bKO/KO mice exhibited a significant reduction of the cone photoreceptor light responses, thinning of the outer nuclear layer, and loss of cone photoreceptors across the retina. CONCLUSION: Our results provide a link between loss-of-function variants in a ceramide synthase gene and human retinal dystrophy. Establishment of the Tlcd3b knockout murine model, an in vivo photoreceptor cell degeneration model due to loss of a ceramide synthase, will provide a unique opportunity in probing the role of ceramide in survival and function of photoreceptor cells.

Large-Scale Whole-Genome Sequencing Reveals the Genetic Architecture of Primary Membranoproliferative GN and C3 Glomerulopathy.

BACKGROUND: Primary membranoproliferative GN, including complement 3 (C3) glomerulopathy, is a rare, untreatable kidney disease characterized by glomerular complement deposition. Complement gene mutations can cause familial C3 glomerulopathy, and studies have reported rare variants in complement genes in nonfamilial primary membranoproliferative GN. METHODS: We analyzed whole-genome sequence data from 165 primary membranoproliferative GN cases and 10,250 individuals without the condition (controls) as part of the National Institutes of Health Research BioResource-Rare Diseases Study. We examined copy number, rare, and common variants. RESULTS: Our analysis included 146 primary membranoproliferative GN cases and 6442 controls who were unrelated and of European ancestry. We observed no significant enrichment of rare variants in candidate genes (genes encoding components of the complement alternative pathway and other genes associated with the related disease atypical hemolytic uremic syndrome; 6.8% in cases versus 5.9% in controls) or exome-wide. However, a significant common variant locus was identified at 6p21.32 (rs35406322) (P=3.29×10-8; odds ratio [OR], 1.93; 95% confidence interval [95% CI], 1.53 to 2.44), overlapping the HLA locus. Imputation of HLA types mapped this signal to a haplotype incorporating DQA1*05:01, DQB1*02:01, and DRB1*03:01 (P=1.21×10-8; OR, 2.19; 95% CI, 1.66 to 2.89). This finding was replicated by analysis of HLA serotypes in 338 individuals with membranoproliferative GN and 15,614 individuals with nonimmune renal failure. CONCLUSIONS: We found that HLA type, but not rare complement gene variation, is associated with primary membranoproliferative GN. These findings challenge the paradigm of complement gene mutations typically causing primary membranoproliferative GN and implicate an underlying autoimmune mechanism in most cases.

Biallelic Mutation of ARHGEF18, Involved in the Determination of Epithelial Apicobasal Polarity, Causes Adult-Onset Retinal Degeneration.

Mutations in more than 250 genes are implicated in inherited retinal dystrophy; the encoded proteins are involved in a broad spectrum of pathways. The presence of unsolved families after highly parallel sequencing strategies suggests that further genes remain to be identified. Whole-exome and -genome sequencing studies employed here in large cohorts of affected individuals revealed biallelic mutations in ARHGEF18 in three such individuals. ARHGEF18 encodes ARHGEF18, a guanine nucleotide exchange factor that activates RHOA, a small GTPase protein that is a key component of tight junctions and adherens junctions. This biological pathway is known to be important for retinal development and function, as mutation of CRB1, encoding another component, causes retinal dystrophy. The retinal structure in individuals with ARHGEF18 mutations resembled that seen in subjects with CRB1 mutations. Five mutations were found on six alleles in the three individuals: c.808A>G (p.Thr270Ala), c.1617+5G>A (p.Asp540Glyfs∗63), c.1996C>T (p.Arg666∗), c.2632G>T (p.Glu878∗), and c.2738_2761del (p.Arg913_Glu920del). Functional tests suggest that each disease genotype might retain some ARHGEF18 activity, such that the phenotype described here is not the consequence of nullizygosity. In particular, the p.Thr270Ala missense variant affects a highly conserved residue in the DBL homology domain, which is required for the interaction and activation of RHOA. Previously, knock-out of Arhgef18 in the medaka fish has been shown to cause larval lethality which is preceded by retinal defects that resemble those seen in zebrafish Crumbs complex knock-outs. The findings described here emphasize the peculiar sensitivity of the retina to perturbations of this pathway, which is highlighted as a target for potential therapeutic strategies.

Comprehensive Rare Variant Analysis via Whole-Genome Sequencing to Determine the Molecular Pathology of Inherited Retinal Disease.

Inherited retinal disease is a common cause of visual impairment and represents a highly heterogeneous group of conditions. Here, we present findings from a cohort of 722 individuals with inherited retinal disease, who have had whole-genome sequencing (n = 605), whole-exome sequencing (n = 72), or both (n = 45) performed, as part of the NIHR-BioResource Rare Diseases research study. We identified pathogenic variants (single-nucleotide variants, indels, or structural variants) for 404/722 (56%) individuals. Whole-genome sequencing gives unprecedented power to detect three categories of pathogenic variants in particular: structural variants, variants in GC-rich regions, which have significantly improved coverage compared to whole-exome sequencing, and variants in non-coding regulatory regions. In addition to previously reported pathogenic regulatory variants, we have identified a previously unreported pathogenic intronic variant in CHM in two males with choroideremia. We have also identified 19 genes not previously known to be associated with inherited retinal disease, which harbor biallelic predicted protein-truncating variants in unsolved cases. Whole-genome sequencing is an increasingly important comprehensive method with which to investigate the genetic causes of inherited retinal disease.

Unique noncoding variants upstream of PRDM13 are associated with a spectrum of developmental retinal dystrophies including progressive bifocal chorioretinal atrophy.

The autosomal dominant progressive bifocal chorioretinal atrophy (PBCRA) disease locus has been mapped to chromosome 6q14-16.2 that overlaps the North Carolina macular dystrophy (NCMD) locus MCDR1. NCMD is a nonprogressive developmental macular dystrophy, in which variants upstream of PRDM13 have been implicated. Whole genome sequencing was performed to interrogate structural variants (SVs) and single nucleotide variants (SNVs) in eight individuals, six affected individuals from two families with PBCRA, and two individuals from an additional family with a related developmental macular dystrophy. A SNV (chr6:100,046,804T>C), located 7.8 kb upstream of the PRDM13 gene, was shared by all PBCRA-affected individuals in the disease locus. Haplotype analysis suggested that the variant arose independently in the two families. The two affected individuals from Family 3 were screened for rare variants in the PBCRA and NCMD loci. This revealed a de novo variant in the proband, 21 bp from the first SNV (chr6:100,046,783A>C). This study expands the noncoding variant spectrum upstream of PRDM13 and suggests altered spatio-temporal expression of PRDM13 as a candidate disease mechanism in the phenotypically distinct but related conditions, NCMD and PBCRA.

Mutations in REEP6 Cause Autosomal-Recessive Retinitis Pigmentosa.

Retinitis pigmentosa (RP) is the most frequent form of inherited retinal dystrophy. RP is genetically heterogeneous and the genes identified to date encode proteins involved in a wide range of functional pathways, including photoreceptor development, phototransduction, the retinoid cycle, cilia, and outer segment development. Here we report the identification of biallelic mutations in Receptor Expression Enhancer Protein 6 (REEP6) in seven individuals with autosomal-recessive RP from five unrelated families. REEP6 is a member of the REEP/Yop1 family of proteins that influence the structure of the endoplasmic reticulum but is relatively unstudied. The six variants identified include three frameshift variants, two missense variants, and a genomic rearrangement that disrupts exon 1. Human 3D organoid optic cups were used to investigate REEP6 expression and confirmed the expression of a retina-specific isoform REEP6.1, which is specifically affected by one of the frameshift mutations. Expression of the two missense variants (c.383C>T [p.Pro128Leu] and c.404T>C [p.Leu135Pro]) and the REEP6.1 frameshift mutant in cultured cells suggest that these changes destabilize the protein. Furthermore, CRISPR-Cas9-mediated gene editing was used to produce Reep6 knock-in mice with the p.Leu135Pro RP-associated variant identified in one RP-affected individual. The homozygous knock-in mice mimic the clinical phenotypes of RP, including progressive photoreceptor degeneration and dysfunction of the rod photoreceptors. Therefore, our study implicates REEP6 in retinal homeostasis and highlights a pathway previously uncharacterized in retinal dystrophy.

ABCA4-associated disease as a model for missing heritability in autosomal recessive disorders: novel noncoding splice, cis-regulatory, structural, and recurrent hypomorphic variants.

PURPOSE: ABCA4-associated disease, a recessive retinal dystrophy, is hallmarked by a large proportion of patients with only one pathogenic ABCA4 variant, suggestive for missing heritability. METHODS: By locus-specific analysis of ABCA4, combined with extensive functional studies, we aimed to unravel the missing alleles in a cohort of 67 patients (p), with one (p = 64) or no (p = 3) identified coding pathogenic variants of ABCA4. RESULTS: We identified eight pathogenic (deep-)intronic ABCA4 splice variants, of which five are novel and six structural variants, four of which are novel, including two duplications. Together, these variants account for the missing alleles in 40.3% of patients. Furthermore, two novel variants with a putative cis-regulatory effect were identified. The common hypomorphic variant c.5603A>T p.(Asn1868Ile) was found as a candidate second allele in 43.3% of patients. Overall, we have elucidated the missing heritability in 83.6% of our cohort. In addition, we successfully rescued three deep-intronic variants using antisense oligonucleotide (AON)-mediated treatment in HEK 293-T cells and in patient-derived fibroblast cells. CONCLUSION: Noncoding pathogenic variants, novel structural variants, and a common hypomorphic allele of the ABCA4 gene explain the majority of unsolved cases with ABCA4-associated disease, rendering this retinopathy a model for missing heritability in autosomal recessive disorders.

Deep-intronic ABCA4 variants explain missing heritability in Stargardt disease and allow correction of splice defects by antisense oligonucleotides.

PURPOSE: Using exome sequencing, the underlying variants in many persons with autosomal recessive diseases remain undetected. We explored autosomal recessive Stargardt disease (STGD1) as a model to identify the missing heritability. METHODS: Sequencing of ABCA4 was performed in 8 STGD1 cases with one variant and p.Asn1868Ile in trans, 25 cases with one variant, and 3 cases with no ABCA4 variant. The effect of intronic variants was analyzed using in vitro splice assays in HEK293T cells and patient-derived fibroblasts. Antisense oligonucleotides were used to correct splice defects. RESULTS: In 24 of the probands (67%), one known and five novel deep-intronic variants were found. The five novel variants resulted in messenger RNA pseudoexon inclusions, due to strengthening of cryptic splice sites or by disrupting a splicing silencer motif. Variant c.769-784C>T showed partial insertion of a pseudoexon and was found in cis with c.5603A>T (p.Asn1868Ile), so its causal role could not be fully established. Variant c.4253+43G>A resulted in partial skipping of exon 28. Remarkably, antisense oligonucleotides targeting the aberrant splice processes resulted in (partial) correction of all splicing defects. CONCLUSION: Our data demonstrate the importance of assessing noncoding variants in genetic diseases, and show the great potential of splice modulation therapy for deep-intronic variants.

Assessment of the incorporation of CNV surveillance into gene panel next-generation sequencing testing for inherited retinal diseases.

BACKGROUND: Diagnostic use of gene panel next-generation sequencing (NGS) techniques is commonplace for individuals with inherited retinal dystrophies (IRDs), a highly genetically heterogeneous group of disorders. However, these techniques have often failed to capture the complete spectrum of genomic variation causing IRD, including CNVs. This study assessed the applicability of introducing CNV surveillance into first-tier diagnostic gene panel NGS services for IRD. METHODS: Three read-depth algorithms were applied to gene panel NGS data sets for 550 referred individuals, and informatics strategies used for quality assurance and CNV filtering. CNV events were confirmed and reported to referring clinicians through an accredited diagnostic laboratory. RESULTS: We confirmed the presence of 33 deletions and 11 duplications, determining these findings to contribute to the confirmed or provisional molecular diagnosis of IRD for 25 individuals. We show that at least 7% of individuals referred for diagnostic testing for IRD have a CNV within genes relevant to their clinical diagnosis, and determined a positive predictive value of 79% for the employed CNV filtering techniques. CONCLUSION: Incorporation of CNV analysis increases diagnostic yield of gene panel NGS diagnostic tests for IRD, increases clarity in diagnostic reporting and expands the spectrum of known disease-causing mutations.