Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads.

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

Clinical and Genetic Analysis of Patients With TK2 Deficiency.

Objectives: Thymidine kinase 2 deficiency (TK2d) is a rare autosomal recessive disorder that stems from a perturbation of the mitochondrial DNA maintenance. Nucleoside treatment has recently shown promise as a disease-modifying therapy. TK2d was initially associated with rapidly progressive fatal myopathy in children featuring mitochondrial DNA depletion. Subsequently, less severe variants of the disease were described, with onset of symptoms during adolescence or adulthood and associated with the presence of multiple mtDNA deletions. These less severe phenotypes have been reported in only 15% of the approximately 120 patients described worldwide. However, some reports suggest that these juvenile and adult-onset presentations may be more common. The objective of this study was to describe the clinical phenotype in a sample of patients from Spain. Methods: This study includes 53 patients harboring biallelic TK2 pathogenic variants, compiling data retrospectively from 7 Spanish centers. We analyzed allele frequency, investigated the most recent common ancestor of core haplotypes, and used the Runs of Homozygosity approach to investigate variant coalescence. Results: Symptom onset distribution revealed that 32 patients (60%) experienced symptoms beyond 12 years of age. Approximately 30% of patients died of respiratory insufficiency, while 56% of surviving patients needed mechanical ventilation. Genetic analysis identified 16 distinct variants in TK2. Two variants, p.Lys202del and p.Thr108Met, exhibited significantly higher prevalence in the Spanish population than that reported in gnomAD database (86-fold and 13-fold, respectively). These variants are estimated to have originated approximately 16.8 generations ago for p.Thr108Met and 95.2 generations ago for p.Lys202del within the Spanish population, with the increase in frequency attributed to various forms of inbreeding. In late-onset cases, 46.9% carried the p.Lys202del variant. Discussion: The higher frequency of TK2d in Spain can be partially attributed to the increased prevalence of 2 variants and consanguinity. Notably, in 60% of the cohort, the disease was late-onset, emphasizing the potential underdiagnosis of this subgroup of patients in other regions. Raising awareness of this potentially treatable disorder is of utmost importance because early interventions can significantly affect the quality of life and survival of affected individuals.

PRPH2-Related Retinal Dystrophies: Mutational Spectrum in 103 Families from a Spanish Cohort.

PRPH2, one of the most frequently inherited retinal dystrophy (IRD)-causing genes, implies a high phenotypic variability. This study aims to analyze the PRPH2 mutational spectrum in one of the largest cohorts worldwide, and to describe novel pathogenic variants and genotype-phenotype correlations. A study of 220 patients from 103 families recruited from a database of 5000 families. A molecular diagnosis was performed using classical molecular approaches and next-generation sequencing. Common haplotypes were ascertained by analyzing single-nucleotide polymorphisms. We identified 56 variants, including 11 novel variants. Most of them were missense variants (64%) and were located in the D2-loop protein domain (77%). The most frequently occurring variants were p.Gly167Ser, p.Gly208Asp and p.Pro221_Cys222del. Haplotype analysis revealed a shared region in families carrying p.Leu41Pro or p.Pro221_Cys222del. Patients with retinitis pigmentosa presented an earlier disease onset. We describe the largest cohort of IRD families associated with PRPH2 from a single center. Most variants were located in the D2-loop domain, highlighting its importance in interacting with other proteins. Our work suggests a likely founder effect for the variants p.Leu41Pro and p.Pro221_Cys222del in our Spanish cohort. Phenotypes with a primary rod alteration presented more severe affectation. Finally, the high phenotypic variability in PRPH2 hinders the possibility of drawing genotype-phenotype correlations.

Towards Uncovering the Role of Incomplete Penetrance in Maculopathies through Sequencing of 105 Disease-Associated Genes.

Inherited macular dystrophies (iMDs) are a group of genetic disorders, which affect the central region of the retina. To investigate the genetic basis of iMDs, we used single-molecule Molecular Inversion Probes to sequence 105 maculopathy-associated genes in 1352 patients diagnosed with iMDs. Within this cohort, 39.8% of patients were considered genetically explained by 460 different variants in 49 distinct genes of which 73 were novel variants, with some affecting splicing. The top five most frequent causative genes were ABCA4 (37.2%), PRPH2 (6.7%), CDHR1 (6.1%), PROM1 (4.3%) and RP1L1 (3.1%). Interestingly, variants with incomplete penetrance were revealed in almost one-third of patients considered solved (28.1%), and therefore, a proportion of patients may not be explained solely by the variants reported. This includes eight previously reported variants with incomplete penetrance in addition to CDHR1:c.783G>A and CNGB3:c.1208G>A. Notably, segregation analysis was not routinely performed for variant phasing-a limitation, which may also impact the overall diagnostic yield. The relatively high proportion of probands without any putative causal variant (60.2%) highlights the need to explore variants with incomplete penetrance, the potential modifiers of disease and the genetic overlap between iMDs and age-related macular degeneration. Our results provide valuable insights into the genetic landscape of iMDs and warrant future exploration to determine the involvement of other maculopathy genes.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease.

BACKGROUND: 5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs). METHODS: We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches. RESULTS: Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels. CONCLUSIONS: This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

Apolipoprotein E-ε2 and Resistance to Atherosclerosis in Midlife: The PESA Observational Study.

BACKGROUND: APOE is a known genetic contributor to cardiovascular disease, but the differential role APOE alleles play in subclinical atherosclerosis remains unclear. METHODS: The PESA (Progression of Early Subclinical Atherosclerosis) is an observational cohort study that recruited 4184 middle-aged asymptomatic individuals to be screened for cardiovascular risk and multiterritorial subclinical atherosclerosis. Participants were APOE-genotyped, and omics data were additionally evaluated. RESULTS: In the PESA study, the frequencies for APOE -ε2, -ε3, and -ε4 alleles were 0.060, 0.844, and 0.096, respectively. This study included a subcohort of 3887 participants (45.8±4.3 years of age; 62% males). As expected, APOE-ε4 carriers were at the highest risk for cardiovascular disease and had significantly greater odds of having subclinical atherosclerosis compared with ε3/ε3 carriers, which was mainly explained by their higher levels of low-density lipoprotein (LDL)-cholesterol. In turn, APOE-ε2 carriers were at the lowest risk for cardiovascular disease and had significantly lower odds of having subclinical atherosclerosis in several vascular territories (carotids: 0.62 [95% CI, 0.47-0.81]; P=0.00043; femorals: 0.60 [0.47-0.78]; P=9.96×10-5; coronaries: 0.53 [0.39-0.74]; P=0.00013; and increased PESA score: 0.58 [0.48-0.71]; P=3.16×10-8). This APOE-ε2 atheroprotective effect was mostly independent of the associated lower LDL-cholesterol levels and other cardiovascular risk factors. The protection conferred by the ε2 allele was greater with age (50-54 years: 0.49 [95% CI, 0.32-0.73]; P=0.00045), and normal (<150 mg/dL) levels of triglycerides (0.54 [0.44-0.66]; P=4.70×10-9 versus 0.90 [0.57-1.43]; P=0.67 if ≥150 mg/dL). Omics analysis revealed an enrichment of several canonical pathways associated with anti-inflammatory mechanisms together with the modulation of erythrocyte homeostasis, coagulation, and complement activation in ε2 carriers that might play a relevant role in the ε2's atheroprotective effect. CONCLUSIONS: This work sheds light on the role of APOE in cardiovascular disease development with important therapeutic and prevention implications on cardiovascular health, especially in early midlife. REGISTRATION: URL: https://www.clinicaltrials.gov: NCT01410318.

The p.C759F Variant in USH2A Is a Pathogenic Mutation: Systematic Literature Review and Meta-Analysis of 667 Genotypes.

BACKGROUND: Although the p.C759F (c.2276G>T, p.Cys759Phe) variant in the USH2A gene has been identified in association with retinal degeneration by several authors, its pathogenicity has been questioned once by the publication of two unaffected homozygotes from a single family. OBJECTIVES: The objective of the study was to ascertain the role of p.C759F in hereditary retinal disease. METHODS: We examined 87 research articles reporting on patients carrying this variant and then used this information as primary data for a series of meta-analytical tests. RESULTS: Independent statistical analyses showed that p.C759F (i) is highly enriched in patients with respect to healthy individuals, (ii) represents a clear-cut recessive allele causing disease when it is in trans with other mutations, (iii) is pathogenic in homozygotes. CONCLUSIONS: Our results confirm that p.C759F is a bona fide mutation, leading to retinal blindness according to a recessive pattern of inheritance.

ABCA4 c.6480-35A>G, a novel branchpoint variant associated with Stargardt disease.

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

KCNV2-associated retinopathy: genotype-phenotype correlations - KCNV2 study group report 3.

BACKGROUND/AIMS: To investigate genotype-phenotype associations in patients with KCNV2 retinopathy. METHODS: Review of clinical notes, best-corrected visual acuity (BCVA), molecular variants, electroretinography (ERG) and retinal imaging. Subjects were grouped according to the combination of KCNV2 variants-two loss-of-function (TLOF), two missense (TM) or one of each (MLOF)-and parameters were compared. RESULTS: Ninety-two patients were included. The mean age of onset (mean±SD) in TLOF (n=55), TM (n=23) and MLOF (n=14) groups was 3.51±0.58, 4.07±2.76 and 5.54±3.38 years, respectively. The mean LogMAR BCVA (±SD) at baseline in TLOF, TM and MLOF groups was 0.89±0.25, 0.67±0.38 and 0.81±0.35 for right, and 0.88±0.26, 0.69±0.33 and 0.78±0.33 for left eyes, respectively. The difference in BCVA between groups at baseline was significant in right (p=0.03) and left eyes (p=0.035). Mean outer nuclear layer thickness (±SD) at baseline in TLOF, MLOF and TM groups was 37.07±15.20 µm, 40.67±12.53 and 40.38±18.67, respectively, which was not significantly different (p=0.85). The mean ellipsoid zone width (EZW) loss (±SD) was 2051 µm (±1318) for patients in the TLOF, and 1314 µm (±965) for MLOF. Only one patient in the TM group had EZW loss at presentation. There was considerable overlap in ERG findings, although the largest DA 10 ERG b-waves were associated with TLOF and the smallest with TM variants. CONCLUSIONS: Patients with missense alterations had better BCVA and greater structural integrity. This is important for patient prognostication and counselling, as well as stratification for future gene therapy trials.

Snijders Blok-Campeau Syndrome: Description of 20 Additional Individuals with Variants in CHD3 and Literature Review.

Snijders Blok-Campeau syndrome (SNIBCPS, OMIM# 618205) is an extremely infrequent disease with only approximately 60 cases reported so far. SNIBCPS belongs to the group of neurodevelopmental disorders (NDDs). Clinical features of patients with SNIBCPS include global developmental delay, intellectual disability, speech and language difficulties and behavioral disorders like autism spectrum disorder. In addition, patients with SNIBCPS exhibit typical dysmorphic features including macrocephaly, hypertelorism, sparse eyebrows, broad forehead, prominent nose and pointed chin. The severity of the neurological effects as well as the presence of other features is variable among subjects. SNIBCPS is caused likely by pathogenic and pathogenic variants in CHD3 (Chromodomain Helicase DNA Binding Protein 3), which seems to be involved in chromatin remodeling by deacetylating histones. Here, we report 20 additional patients with clinical features compatible with SNIBCPS from 17 unrelated families with confirmed likely pathogenic/pathogenic variants in CHD3. Patients were analyzed by whole exome sequencing and segregation studies were performed by Sanger sequencing. Patients in this study showed different pathogenic variants affecting several functional domains of the protein. Additionally, none of the variants described here were reported in control population databases, and most computational predictors suggest that they are deleterious. The most common clinical features of the whole cohort of patients are global developmental delay (98%) and speech disorder/delay (92%). Other frequent features (51-74%) include intellectual disability, hypotonia, hypertelorism, abnormality of vision, macrocephaly and prominent forehead, among others. This study expands the number of individuals with confirmed SNIBCPS due to pathogenic or likely pathogenic variants in CHD3. Furthermore, we add evidence of the importance of the application of massive parallel sequencing for NDD patients for whom the clinical diagnosis might be challenging and where deep phenotyping is extremely useful to accurately manage and follow up the patients.

TBC1D32 variants disrupt retinal ciliogenesis and cause retinitis pigmentosa.

Retinitis pigmentosa (RP) is the most common inherited retinal disease (IRD) and is characterized by photoreceptor degeneration and progressive vision loss. We report 4 patients presenting with RP from 3 unrelated families with variants in TBC1D32, which to date has never been associated with an IRD. To validate TBC1D32 as a putative RP causative gene, we combined Xenopus in vivo approaches and human induced pluripotent stem cell-derived (iPSC-derived) retinal models. Our data showed that TBC1D32 was expressed during retinal development and that it played an important role in retinal pigment epithelium (RPE) differentiation. Furthermore, we identified a role for TBC1D32 in ciliogenesis of the RPE. We demonstrated elongated ciliary defects that resulted in disrupted apical tight junctions, loss of functionality (delayed retinoid cycling and altered secretion balance), and the onset of an epithelial-mesenchymal transition-like phenotype. Last, our results suggested photoreceptor differentiation defects, including connecting cilium anomalies, that resulted in impaired trafficking to the outer segment in cones and rods in TBC1D32 iPSC-derived retinal organoids. Overall, our data highlight a critical role for TBC1D32 in the retina and demonstrate that TBC1D32 mutations lead to RP. We thus identify TBC1D32 as an IRD-causative gene.

Risk of Treatment Failure and Death after Ablation in Hepatocellular Carcinoma Patients-A Multiparametric Prediction.

BACKGROUND: Ablation is a first-line treatment for Barcelona Clinic Liver Cancer (BCLC)-0/A hepatocellular carcinoma (HCC). However, there are scarce data about patients' outcomes after recurrence. The present study evaluates the impact of patient and tumor characteristics at baseline and at recurrence on the Clinical Decision-Making process. METHODS: We evaluated BCLC-0/A patients treated with percutaneous ablation from January 2010 to November 2018. Clinical and radiological data such as age, tumor location at ablation, pattern of recurrence/progression, and comorbidities during follow-up were registered. Tumor location was divided into 'suboptimal' vs. 'optimal' locations for ablation. The Clinical Decision-Making was based on tumor burden, liver dysfunction, or comorbidities. The statistical analysis included the time-to-recurrence/progression, censoring at time of death, date of last follow-up or liver transplantation, and time-to-event was estimated by the Kaplan-Meier method and Cox regression models to evaluate the risk of an event of death and change of treatment strategy. RESULTS: A total of 225 patients [39.1% BCLC-0 and 60.9% BCLC-A] were included, 190 had unifocal HCC and 82.6% were ≤3 cm. The complete response rate and median overall survival were 96% and 60.7 months. The HCC nodules number (Hazard Ratio-HR 3.1), Child-Pugh (HR 2.4), and Albumin-Bilirubin score (HR 3.2) were associated with increased risk of death during follow-up. HCC in 'suboptimal location' presented a shorter time to recurrence. When comorbidities prevented further loco-regional or systemic treatment, the risk of death was significantly increased (HR 2.0, p = 0.0369) in comparison to those who received treatment. CONCLUSIONS: These results expose the impact of non-liver comorbidities when considering treatment for recurrence after ablation in the real-world setting and in research trials. Ultimately, we identified an orphan population for which effective interventions are needed.

Long-read genome sequencing identifies cryptic structural variants in congenital aniridia cases.

BACKGROUND: Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved "PAX6-negative" cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches. RESULTS: Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1. CONCLUSIONS: In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

Comprehensive Genotyping and Phenotyping Analysis of GUCY2D-Associated Rod- and Cone-Dominated Dystrophies.

PURPOSE: To describe the genetic and clinical spectrum of GUCY2D-associated retinopathies and to accurately establish their prevalence in a large cohort of patients. DESIGN: Retrospective case series. METHODS: Institutional study of 47 patients from 27 unrelated families with retinal dystrophies carrying disease-causing GUCY2D variants from the Fundación Jiménez Díaz hospital dataset of 8000 patients. Patients underwent ophthalmological examination and molecular testing by Sanger or exome sequencing approaches. Statistical and principal component analyses were performed to determine genotype-phenotype correlations. RESULTS: Four clinically different associated phenotypes were identified: 66.7% of families with cone/cone-rod dystrophy, 22.2% with Leber congenital amaurosis, 7.4% with early-onset retinitis pigmentosa, and 3.7% with congenital night blindness. Twenty-three disease-causing GUCY2D variants were identified, including 6 novel variants. Biallelic variants accounted for 28% of patients, whereas most carried dominant alleles associated with cone/cone-rod dystrophy. The disease onset had statistically significant differences according to the functional variant effect. Patients carrying GUCY2D variants were projected into 3 subgroups by allelic combination, disease onset, and presence of nystagmus or night blindness. In contrast to patients with the most severe phenotype of Leber congenital amaurosis, 7 patients with biallelic GUCY2D had a later and milder rod form with night blindness in infancy as the first symptom. CONCLUSIONS: This study represents the largest GUCY2D cohort in which 4 distinctly different phenotypes were identified, including rare intermediate presentations of rod-dominated retinopathies. We established that GUCY2D is linked to about 1% of approximately 3000 molecularly characterized families of our cohort. All of these findings are critical for defining cohorts for inclusion in future clinical trials.

Clinical and genetic analysis further delineates the phenotypic spectrum of ALDH1A3-related anophthalmia and microphthalmia.

Biallelic pathogenic variants in ALDH1A3 are responsible for approximately 11% of recessively inherited cases of severe developmental eye anomalies. Some individuals can display variable neurodevelopmental features, but the relationship to the ALDH1A3 variants remains unclear. Here, we describe seven unrelated families with biallelic pathogenic ALDH1A3 variants: four compound heterozygous and three homozygous. All affected individuals had bilateral anophthalmia/microphthalmia (A/M), three with additional intellectual or developmental delay, one with autism and seizures and three with facial dysmorphic features. This study confirms that individuals with biallelic pathogenic ALDH1A3 variants consistently manifest A/M, but additionally display neurodevelopmental features with significant intra- and interfamilial variability. Furthermore, we describe the first case with cataract and highlight the importance of screening ALDH1A3 variants in nonconsanguineous families with A/M.

Minigene Splicing Assays and Long-Read Sequencing to Unravel Pathogenic Deep-Intronic Variants in PAX6 in Congenital Aniridia.

PAX6 haploinsufficiency causes aniridia, a congenital eye disorder that involves the iris, and foveal hypoplasia. Comprehensive screening of the PAX6 locus, including the non-coding regions, by next-generation sequencing revealed four deep-intronic variants with potential effects on pre-RNA splicing. Nevertheless, without a functional analysis, their pathogenicity could not be established. We aimed to decipher their impact on the canonical PAX6 splicing using in vitro minigene splicing assays and nanopore-based long-read sequencing. Two multi-exonic PAX6 constructs were generated, and minigene assays were carried out. An aberrant splicing pattern was observed for two variants in intron 6, c.357+136G>A and c.357+334G>A. In both cases, several exonization events, such as pseudoexon inclusions and partial intronic retention, were observed due to the creation or activation of new/cryptic non-canonical splicing sites, including a shared intronic donor site. In contrast, two variants identified in intron 11, c.1032+170A>T and c.1033-275A>C, seemed not to affect splicing processes. We confirmed the high complexity of alternative splicing of PAX6 exon 6, which also involves unreported cryptic intronic sites. Our study highlights the importance of integrating functional studies into diagnostic algorithms to decipher the potential implication of non-coding variants, usually classified as variants of unknown significance, thus allowing variant reclassification to achieve a conclusive genetic diagnosis.

Prioritization of New Candidate Genes for Rare Genetic Diseases by a Disease-Aware Evaluation of Heterogeneous Molecular Networks.

Screening for pathogenic variants in the diagnosis of rare genetic diseases can now be performed on all genes thanks to the application of whole exome and genome sequencing (WES, WGS). Yet the repertoire of gene-disease associations is not complete. Several computer-based algorithms and databases integrate distinct gene-gene functional networks to accelerate the discovery of gene-disease associations. We hypothesize that the ability of every type of information to extract relevant insights is disease-dependent. We compiled 33 functional networks classified into 13 knowledge categories (KCs) and observed large variability in their ability to recover genes associated with 91 genetic diseases, as measured using efficiency and exclusivity. We developed GLOWgenes, a network-based algorithm that applies random walk with restart to evaluate KCs' ability to recover genes from a given list associated with a phenotype and modulates the prediction of new candidates accordingly. Comparison with other integration strategies and tools shows that our disease-aware approach can boost the discovery of new gene-disease associations, especially for the less obvious ones. KC contribution also varies if obtained using recently discovered genes. Applied to 15 unsolved WES, GLOWgenes proposed three new genes to be involved in the phenotypes of patients with syndromic inherited retinal dystrophies.

Biallelic intragenic tandem duplication of CPLANE1 in Joubert syndrome: A case report.

Joubert syndrome (JS) is a clinically and genetically heterogeneous genetic disorder. To date, 40 JS-causing genes have been reported and CPLANE1 is one of the most frequently mutated, with biallelic pathogenic missense and truncating variants explaining up to 14% of JS cases. We present a case of JS diagnosed after the identification of a novel biallelic intragenic duplication of exons 20-46 of CPLANE1. The quadruplication was identified by short-read sequencing and copy number variant analysis and confirmed in tandem by long PCR with the breakpoints defined by a nanopore-based long-read sequencing approach. Based on the genetic findings and the clinical presentation of the patient, a brain MRI was ordered, evidencing the molar tooth sign, which confirmed the diagnosis of JS in the patient. This is, to the best of our knowledge, the first report of an intragenic duplication in this gene as the potential molecular mechanism of JS.