Examination of Subbasal Nerve Plexus and Central Corneal Stromal Microstructure in Subjects With Congenital Aniridia, Using in Vivo Confocal Laser Scanning Microscopy

PURPOSE

During life up to 70% of aniridia subjects develop aniridia-associated keratopathy (AAK). AAK is characterized by limbal stem cell insufficiency, impaired corneal epithelial cell differentiation and abnormal cell adhesion, which leads to centripetal spreading vascularization, conjunctivalization, and thickening of the cornea. Our aim was to examine the subbasal nerve plexus and central corneal stromal microstructure in subjects with congenital aniridia, using in vivo confocal laser scanning microscopy CLSM.

METHODS

31 eyes of 18 patients (55.6% males, mean age: 25.22 ± 16.35 years) with congenital aniridia and 46 eyes of 29 healthy subjects (41.4% males, mean age 30 ± 14.82 years) were examined using the Rostock Cornea Module of Heidelberg Retina Tomograph-III. At the subbasal nerve plexus, corneal nerve fiber density (CNFD), corneal nerve fiber length (CNFL), corneal total branch density (CTBD), and corneal nerve fiber width (CNFW) were analyzed using ACCMetrics software. Keratocyte density in the anterior, middle and posterior stroma was assessed manually.

RESULTS

The CNFD (2.02 ± 4.08 vs 13.99 ± 6.34/mm2), CNFL (5.78 ± 2.68 vs 10.56 ± 2.82 mm/mm2) and CTBD (15.08 ± 15.62 vs 27.44 ± 15.05/mm2) were significantly lower in congenital aniridia subjects than in controls (p < 0.001 for all). CNFW was significantly higher in aniridia subjects than in controls (0.03 ± 0.004 vs 0.02 ± 0.003 mm/mm2) (p = 0.003). Keratocyte density was significantly lower in all stromal layers of aniridia subjects than in controls (p < 0.001 for all). Stromal alterations included confluent keratocytes, keratocytes with long extensions and hyperreflective dots between keratocytes in aniridia.

CONCLUSIONS

Decrease in CNFD, CNFL, and CTBD, as well as increase in CNFW well refer to the congenital aniridia-associated neuropathy. The decreased keratocyte density and the stromal alterations may be related to an increased cell death in congenital aniridia, nevertheless, stromal changes in different stages of AAK have to be further analyzed in detail.

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.

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.

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.

[Congenital aniridia patients' experience on their visual impairment in Hungary.]

INTRODUCTION

Aniridia is a rare congenital panocular disease associated with varying degrees of visual acuity impairment.

OBJECTIVE

To assess the experiences of congenital aniridia patients in Hungary, with visual impairment using a questionnaire developed by the ANIRIDIA-NET.

PATIENTS AND METHOD

Patients completed the Hungarian version of the 20-item ANIRIDIA-NET questionnaire with our assistance. The questionnaire covered demographic data, the most common complaints caused by the disease, the difficulties caused by low vision in different life situations and the frequency of low vision aids used in daily life.

RESULTS

33 subjects (17 female [51.51%] and 16 male [48.48%]), 16 (48.5%) children and 17 (51.5%) adults completed the questionnaire, with an age of 25.69 ± 17.49 years (5-59 years). Daily photosensitivity was reported by 27 (81.8%), dry eyes by 5 (15.2%), tearing by 4 (12.1%), fluctuating vision by 3 (9.1%), and eye pain by 2 (6.1%) subjects. The majority of respondents said that personal communication with schoolmates (16 [48.5%]) or colleagues at work (11 [33.3%]) never caused difficulties because of their visual impairment. 29 people (87.9%) never needed help with daily routines at home, 24 (72.7%) with getting to school/work and 17 (51.5%) with various activities. 29 people (87.8%) never used low vision aids for communication, 23 (69.7%) for travelling, 20 (60.6%) for participating in social activities, 18 (54.5%) for studying/work.

CONCLUSION

Although aniridia is associated with reduced visual acuity, the majority of people with congenital aniridia, especially in childhood, manage to cope with personal communication and various life situations without difficulty, despite their eye complaints. Low vision aids can be an important aid for them as they grow into adulthood and as they age. Orv Hetil. 2023; 164(34): 1342-1349.

[Staging of aniridia-associated keratopathy]

INTRODUCTION

Congenital aniridia is a rare panocular disease that affects almost all eye structures leading in most patients to reduced visual acuity. Ophthalmological signs include aniridia-associated keratopathy, secondary glaucoma, cataract, macular and optic nerve head hypoplasia, nystagmus. Although the term aniridia-associated keratopathy has long been used in the literature, various staging proposals have been described.

OBJECTIVE

To analyze aniridia-associated keratopathy stages, using available literature classifications, in patients with aniridia in Hungary.

PATIENTS AND METHODS

We examined 65 eyes of 33 patients with congenital aniridia (age: 25.69 ± 17.49 [5-59] years, 17 females [51.51%]). We recorded the corneal status by slit-lamp examination and classified the corneal abnormalities according to the Mackman, Mayer, López-García and Lagali staging.

RESULTS

According to Mackman's classification, 8 eyes (12.3%) were in stage 0, 0 eye in stage 1A, 38 eyes (58.46%) in stage 1B and 19 eyes (29.23%) in stage 2. According to Mayer, stage I included 8 eyes (12.3%), stage II 38 eyes (58.46%), stage III 5 eyes (7.7%), stage IV 7 eyes (10.77%) and stage V 7 eyes (10.77%). In López-García's classification, 8 eyes (12.3%) could not be grouped, 20 eyes (30.77%) were in stage 1, 18 eyes (27.7%) in stage 2 and 19 eyes (29.3%) in stage 3. Lagali's classification included 8 eyes (12.3%) in stage 0, 20 eyes (30.77%) in stage 1, 18 eyes (27.7%) in stage 2, 5 eyes (7.7%) in stage 3 and 14 eyes (21.54%) in stage 4.

CONCLUSION

We recommend using Lagali's staging scheme for aniridia-associated keratoptahy due to its ease of use, detailed progression assessment, and treatment planning. In stage 1 according to Lagali, blood vessels cross the limbus by up to 1 mm, in stage 2 the central 2-3 mm of the corneal area is spared of blood vessels. When the blood vessels reach the center of the cornea, it is stage 3, followed by opaque, uneven corneal pannus in stage 4. Orv Hetil. 2023; 164(27): 1063-1069.

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.

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.

New Cerebellar Ataxia, Neuropathy, Vestibular Areflexia Syndrome cases are caused by the presence of a nonsense variant in compound heterozygosity with the pathogenic repeat expansion in the RFC1 gene

The biallelic pathogenic repeat (AAGGG)_400-2000 intronic expansion in the RFC1 gene has been recently described as the cause of cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) and as a major cause of late-onset ataxia. Since then, many heterozygous carriers have been identified, with an estimated allele frequency of 0.7% to 4% in the healthy population. Here, we describe in two affected CANVAS sisters the presence of the nonsense c.724C > T p.(Arg242*) variant in compound heterozygosity with the pathogenic repeat expansion in the RFC1 gene. Further RNA analysis demonstrated a reduced expression of the p.Arg242* allele in patients confirming an efficient nonsense-mediated mRNA decay. We also highlight the importance of considering the sequencing of the RFC1 gene for the diagnosis, especially in patients with CANVAS diagnosis carriers of the AAGGG repeat expansion.

[Congenital aniridia - Hungarian data of a spectrum disease]

INTRODUCTION

Congenital aniridia is a rare disease, characterised by the complete or partial absence of the iris, but lesions may be present in all structures of the eye.

OBJECTIVE

To determine the prevalence of ocular diseases in congenital aniridia by analyzing patients from a Hungarian centre.

PATIENTS AND METHODS

Patients at the Department of Ophthalmology of Semmelweis University, examined between October 2005 and May 2022, have been included. After taking the patients' medical history, a detailed ophthalmological examination has been performed.

RESULTS

Of the 82 patients in the database, 33 (age 25.69 ± 17.49 [5-59] years, 17 females [51.51%]) presented for examination and 65 eyes were examined. Nystagmus was found in 45 eyes of 23 patients (69.23%), and the patients' uncorrected distance visual acuity was 0.14 ± 0.128 (0.9 logMAR; 0.63-0.005). The aniridia-associated keratopathy was Grade 0 in 8 eyes (12.3%), Grade 1 in 10 eyes (15.38%), Grade 2 in 16 eyes (24.62%), Grade 3 in 4 eyes (6.15%) and Grade 4 in 25 eyes (38.46%). 30 eyes (46.15%) of 15 patients had secondary glaucoma, 6 eyes (9.2%) of 3 patients were glaucoma suspect. 8 eyes (12.3%) had a clear lens, 44 eyes (67.69%) had cataract, of which 22 (33.84%) were anterior cortical polar cataracts. 13 eyes (20%) were pseudophakic (PCL) and 7 eyes (10.77%) had lens dislocation or zonular insufficiency. Macular hypoplasia was found in 6 eyes of 3 patients (4.6%) and optic nerve head malformation in 2 eyes of 1 patient (3.03%).

CONCLUSION

The ocular signs of congenital aniridia are aniridia-associated keratopathy, secondary glaucoma, cataract, macular and optic nerve head hypoplasia. Systematic collaboration of different ophthalmological specialties is required for the management and care of all these ocular abnormalities. Orv Hetil. 2023; 164(4): 148-155.

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.

Individuals with heterozygous variants in the Wnt-signalling pathway gene FZD5 delineate a phenotype characterized by isolated coloboma and variable expressivity

BACKGROUND

Anophthalmia, microphthalmia and coloboma are a genetically heterogenous spectrum of developmental eye disorders. Recently, variants in the Wnt-pathway gene Frizzled Class Receptor 5 (FZD5) have been identified in individuals with coloboma and rarely microphthalmia, sometimes with additional phenotypes and variable penetrance.

MATERIALS AND METHODS

We identified variants in FZD5 in individuals with developmental eye disorders from the UK (including the DDD Study [www.ddduk.org/access.html]), France and Spain using whole genome/exome sequencing or customized NGS panels of ocular development genes.

RESULTS

We report eight new families with FZD5 variants and ocular coloboma. Three individuals presented with additional syndromic features, two explicable by additional variants in other genes (SLC12A2 and DDX3X). In two families initially showing incomplete penetrance, re-examination of apparently unaffected carrier individuals revealed subtle ocular colobomatous phenotypes. Finally, we report two families with microphthalmia in addition to coloboma, representing the second and third reported cases of this phenotype in conjunction with FZD5 variants.

CONCLUSIONS

Our findings indicate FZD5 variants are typically associated with isolated ocular coloboma, occasionally microphthalmia, and that extraocular phenotypes are likely to be explained by other gene alterations.

Allelic overload and its clinical modifier effect in Bardet-Biedl syndrome

Bardet–Biedl syndrome (BBS) is an autosomal recessive ciliopathy characterized by extensive inter- and intra-familial variability, in which oligogenic interactions have been also reported. Our main goal is to elucidate the role of mutational load in the clinical variability of BBS. A cohort of 99 patients from 77 different families with biallelic pathogenic variants in a BBS-associated gene was retrospectively recruited. Human Phenotype Ontology terms were used in the annotation of clinical symptoms. The mutational load in 39 BBS-related genes was studied in index cases using different molecular and next-generation sequencing (NGS) approaches. Candidate allele combinations were analysed using the in silico tools ORVAL and DiGePred. After clinical annotation, 76 out of the 99 cases a priori fulfilled established criteria for diagnosis of BBS or BBS-like. BBS1 alleles, found in 42% of families, were the most represented in our cohort. An increased mutational load was excluded in 41% of the index cases (22/54). Oligogenic inheritance was suspected in 52% of the screened families (23/45), being 40 tested by means of NGS data and 5 only by traditional methods. Together, ORVAL and DiGePred platforms predicted an oligogenic effect in 44% of the triallelic families (10/23). Intrafamilial variable severity could be clinically confirmed in six of the families. Our findings show that the presence of more than two alleles in BBS-associated genes correlated in six families with a more severe phenotype and associated with specific findings, highlighting the role of the mutational load in the management of BBS cases.

Presence of rare potential pathogenic variants in subjects under 65 years old with very severe or fatal COVID-19

Rare variants affecting host defense against pathogens could be involved in COVID-19 severity and may help explain fatal outcomes in young and middle-aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID-19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization-based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID-19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID-19 severity.

Novel genes and sex differences in COVID-19 severity

Here we describe the results of a genome-wide study conducted in 11 939 COVID-19 positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (p < 5x10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (p = 1.3x10-22 and p = 8.1x10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (p = 4.4x10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (p = 2.7x10-8) and ARHGAP33 (p = 1.3x10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, p = 4.1x10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥ 60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

Novel homozygous mutations in the transcription factor NRL cause non-syndromic retinitis pigmentosa

Purpose

To describe the clinical phenotype and genetic basis of non-syndromic retinitis pigmentosa (RP) in one family and two sporadic cases with biallelic mutations in the transcription factor neural retina leucine zipper (NRL).

Methods

Exome sequencing was performed in one affected family member. Microsatellite genotyping was used for haplotype analysis. PCR and Sanger sequencing were used to confirm mutations in and screen other family members where they were available. The SMART tool for domain prediction helped us build the protein schematic diagram.

Results

For family MM1 of Pakistani origin, whole-exome sequencing and microsatellite genotyping revealed homozygosity on chromosome 14 and identified a homozygous stop-loss mutation in NRL, NM_006177.5: c.713G>T, p.*238Lext57, which is predicted to add an extra 57 amino acids to the normal protein chain. The variant segregated with disease symptoms in the family. For case RP-3051 of Spanish ancestry, clinical exome sequencing focusing on the morbid genome highlighted a homozygous nonsense mutation in NRL, c.238C>T, p.Gln80*, as the most likely disease candidate. For case RP-1553 of Romanian ethnicity, targeted-exome sequencing of 73 RP/LCA genes identified a homozygous nonsense mutation in NRL, c.544C>T, p.Gln182*. The variants were either rare or absent in the gnomAD database.

Conclusions

NRL mutations predominantly cause dominant retinal disease, but there have been five published reports of mutations causing recessive disease. Here, we present three further examples of recessive RP due to NRL mutations. The phenotypes observed are consistent with those in the previous reports, and the observed mutation types and distribution further confirm distinct patterns for variants in NRL causing recessive and dominant diseases.

Impact of Next Generation Sequencing in Unraveling the Genetics of 1036 Spanish Families With Inherited Macular Dystrophies

Purpose

To assess the potential of next-generation sequencing (NGS) technologies to characterize cases diagnosed with autosomal recessive (ar) or sporadic (s) macular dystrophies (ar/sMD) and describe their mutational spectrum.

Methods

A cohort of 1036 families was classified according to their suspected clinical diagnosis—Stargardt disease (STGD), cone and cone-rod dystrophy (CCRD) or other maculopathies (otherMD). Molecular studies included genotyping microarrays, Sanger sequencing, NGS, and sequencing of intronic regions of the ABCA4 gene. Clinical reclassification was done after the genetic study.

Results

At the end of the study, 677 patients (65%) had a confirmed genetic diagnosis, representing 78%, 63%, and 38% of STGD, CCRD, and otherMD groups of patients, respectively. ABCA4 is the most mutated gene in all groups, and a second pathogenic variant was found in 76% of STGD patients with one previously identified mutated ABCA4 allele. Autosomal dominant or X-linked mutations were found in 5% of cases together with not-MD genes (CHM, EYS, RHO, RPGR, RLBP1, OPA1, and USH2A among others) leading to their reclassification. Novel variants in the very rare genes PLA2G5 and TTLL5 revealed additional phenotypic associations.

Conclusions

This study provides for the first time a genetic landscape of 1036 ar/sMD families according to their suspected diagnosis. The analysis of >200 genes associated with retinal dystrophies and the entire locus of ABCA4 increase the rate of characterization, even regardless of available clinical and familiar data. The use of the suspected a priori diagnosis referred by the clinicians, especially in the past, could lead to clinical reclassifications to other inherited retinal dystrophies.

CIBERER: Spanish National Network for Research on Rare Diseases: a highly productive collaborative initiative

CIBER (Center for Biomedical Network Research; Centro de Investigación Biomédica En Red) is a public national consortium created in 2006 under the umbrella of the Spanish National Institute of Health Carlos III (ISCIII). This innovative research structure comprises 11 different specific areas dedicated to the main public health priorities in the National Health System. CIBERER, the thematic area of CIBER focused on rare diseases (RDs) currently consists of 75 research groups belonging to universities, research centers, and hospitals of the entire country. CIBERER's mission is to be a center prioritizing and favoring collaboration and cooperation between biomedical and clinical research groups, with special emphasis on the aspects of genetic, molecular, biochemical, and cellular research of RDs. This research is the basis for providing new tools for the diagnosis and therapy of low‐prevalence diseases, in line with the International Rare Diseases Research Consortium (IRDiRC) objectives, thus favoring translational research between the scientific environment of the laboratory and the clinical setting of health centers. In this article, we intend to review CIBERER's 15‐year journey and summarize the main results obtained in terms of internationalization, scientific production, contributions toward the discovery of new therapies and novel genes associated to diseases, cooperation with patients' associations and many other topics related to RD research.

RPE65-related retinal dystrophy: Mutational and phenotypic spectrum in 45 affected patients

INTRODUCTION

Biallelic pathogenic RPE65 variants are related to a spectrum of clinically overlapping inherited retinal dystrophies (IRD). Most affected individuals progress to severe disease, with 50% of patients becoming legally blind by 20 years of age. Deeper knowledge of the mutational spectrum and the phenotype-genotype correlation in RPE65-related IRD is needed.

PATIENTS AND METHODS

Forty-five affected subjects from 27 unrelated families with a clinical diagnosis of RPE65-related IRD were included. Clinical evaluation consisted of self-reported ophthalmological history and objective ophthalmological examination. Patients' genotype was classified according to variant class (truncating or missense) or to variant location at different protein domains. The main phenotypic outcome measure was age at onset (AAO) of symptomatic disease and a Kaplan-Meier analysis of disease symptom event-free survival was performed.

RESULTS

Twenty-nine different RPE65 variants were identified in our cohort, 7 of them novel. Patients carrying two missense alleles showed a later disease onset than those with 1 or 2 truncating variants (log-rank test p <0.05). While 60% of patients carrying a missense/missense genotype presented symptoms before or during the first year of life, almost all patients with at least 1 truncating allele (91%) had an AAO ≤1 year (p <0.05).

CONCLUSION

Our findings suggest an association between the type of RPE65 variant carried and AAO. These findings provide useful data on RPE65-associated IRD phenotypes and may help improve clinical and therapeutic management of these patients.

NGS and phenotypic ontology-based approaches increase the diagnostic yield in syndromic retinal diseases

Syndromic retinal diseases (SRDs) are a group of complex inherited systemic disorders, with challenging molecular underpinnings and clinical management. Our main goal is to improve clinical and molecular SRDs diagnosis, by applying a structured phenotypic ontology and next-generation sequencing (NGS)-based pipelines. A prospective and retrospective cohort study was performed on 100 probands with an a priori diagnosis of non-Usher SRDs, using available clinical data, including Human Phenotype Ontology annotation, and further classification into seven clinical categories (ciliopathies, specific syndromes and five others). Retrospective molecular diagnosis was assessed using different molecular and bioinformatic methods depending on availability. Subsequently, uncharacterized probands were prospectively screened using other NGS approaches to extend the number of analyzed genes. After phenotypic classification, ciliopathies were the most common SRD (35%). A global characterization rate of 52% was obtained, with six cases incompletely characterized for a gene that partially explained the phenotype. An improved characterization rate was achieved addressing prospective cases (83%) and well-recognizable syndrome (62%) subgroups. The 27% of the fully characterized cases were reclassified into a different clinical category after identification of the disease-causing gene. Clinical-exome sequencing is the most appropriate first-tier approach for prospective cases, whereas whole-exome sequencing and bioinformatic reanalysis increases the diagnosis of uncharacterized retrospective cases to 45%, mostly those with unspecific symptoms. Our study describes a comprehensive approach to SRDs in daily clinical practice and the importance of thorough clinical assessment and selection of the most appropriate molecular test to be used to solve these complex cases and elucidate novel associations.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

Pathogenic variants in IMPG1 cause autosomal dominant and autosomal recessive retinitis pigmentosa

BACKGROUND

Inherited retinal disorders are a clinically and genetically heterogeneous group of conditions and a major cause of visual impairment. Common disease subtypes include vitelliform macular dystrophy (VMD) and retinitis pigmentosa (RP). Despite the identification of over 90 genes associated with RP, conventional genetic testing fails to detect a molecular diagnosis in about one third of patients with RP.

METHODS

Exome sequencing was carried out for identifying the disease-causing gene in a family with autosomal dominant RP. Gene panel testing and exome sequencing were performed in 596 RP and VMD families to identified additional IMPG1 variants. In vivo analysis in the medaka fish system by knockdown assays was performed to screen IMPG1 possible pathogenic role.

RESULTS

Exome sequencing of a family with RP revealed a splice variant in IMPG1. Subsequently, the same variant was identified in individuals from two families with either RP or VMD. A retrospective study of patients with RP or VMD revealed eight additional families with different missense or nonsense variants in IMPG1. In addition, the clinical diagnosis of the IMPG1 retinopathy-associated variant, originally described as benign concentric annular macular dystrophy, was also revised to RP with early macular involvement. Using morpholino-mediated ablation of Impg1 and its paralog Impg2 in medaka fish, we confirmed a phenotype consistent with that observed in the families, including a decreased length of rod and cone photoreceptor outer segments.

CONCLUSION

This study discusses a previously unreported association between monoallelic or biallelic IMPG1 variants and RP. Notably, similar observations have been reported for IMPG2.

Author Correction: Genetic landscape of 6089 inherited retinal dystrophies affected cases in Spain and their therapeutic and extended epidemiological implications

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Activation of cryptic donor splice sites by non-coding and coding PAX6 variants contributes to congenital aniridia

BACKGROUND

The paired-domain transcription factor paired box gene 6 (PAX6) causes a wide spectrum of ocular developmental anomalies, including congenital aniridia, Peters anomaly and microphthalmia. Here, we aimed to functionally assess the involvement of seven potentially non-canonical splicing variants on missplicing of exon 6, which represents the main hotspot region for loss-of-function PAX6 variants.

METHODS

By locus-specific analysis of PAX6 using Sanger and/or targeted next-generation sequencing, we screened a Spanish cohort of 106 patients with PAX6-related diseases. Functional splicing assays were performed by in vitro minigene approaches or directly in RNA from patient-derived lymphocytes cell line, when available.

RESULTS

Five out seven variants, including three synonymous changes, one small exonic deletion and one non-canonical splice variant, showed anomalous splicing patterns yielding partial exon skipping and/or elongation.

CONCLUSION

We describe new spliceogenic mechanisms for PAX6 variants mediated by creating or strengthening five different cryptic donor sites at exon 6. Our work revealed that the activation of cryptic PAX6 splicing sites seems to be a recurrent and underestimated cause of aniridia. Our findings pointed out the importance of functional assessment of apparently silent PAX6 variants to uncover hidden genetic alterations and to improve variant interpretation for genetic counselling in aniridia.

Prioritizing variants of uncertain significance for reclassification using a rule-based algorithm in inherited retinal dystrophies

Inherited retinal dystrophies (IRD) are a highly heterogeneous group of rare diseases with a molecular diagnostic rate of >50%. Reclassification of variants of uncertain significance (VUS) poses a challenge for IRD diagnosis. We collected 668 IRD cases analyzed by our geneticists using two different clinical exome-sequencing tests. We identified 114 unsolved cases pending reclassification of 125 VUS and studied their genomic, functional, and laboratory-specific features, comparing them to pathogenic and likely pathogenic variants from the same cohort (N = 390). While the clinical exome used did not show differences in diagnostic rate, the more IRD-experienced geneticist reported more VUS (p = 4.07e-04). Significantly fewer VUS were reported in recessive cases (p = 2.14e-04) compared to other inheritance patterns, and of all the genes analyzed, ABCA4 and IMPG2 had the lowest and highest VUS frequencies, respectively (p = 3.89e-04, p = 6.93e-03). Moreover, few frameshift and stop-gain variants were found to be informed VUS (p = 6.73e-08 and p = 2.93e-06). Last, we applied five pathogenicity predictors and found there is a significant proof of deleteriousness when all score for pathogenicity in missense variants. Altogether, these results provided input for a set of rules that correctly reclassified ~70% of VUS as pathogenic in validation datasets. Disease- and setting-specific features influence VUS reporting. Comparison with pathogenic and likely pathogenic variants can prioritize VUS more likely to be reclassified as causal.

Genetic landscape of 6089 inherited retinal dystrophies affected cases in Spain and their therapeutic and extended epidemiological implications

Inherited retinal diseases (IRDs), defined by dysfunction or progressive loss of photoreceptors, are disorders characterized by elevated heterogeneity, both at the clinical and genetic levels. Our main goal was to address the genetic landscape of IRD in the largest cohort of Spanish patients reported to date. A retrospective hospital-based cross-sectional study was carried out on 6089 IRD affected individuals (from 4403 unrelated families), referred for genetic testing from all the Spanish autonomous communities. Clinical, demographic and familiar data were collected from each patient, including family pedigree, age of appearance of visual symptoms, presence of any systemic findings and geographical origin. Genetic studies were performed to the 3951 families with available DNA using different molecular techniques. Overall, 53.2% (2100/3951) of the studied families were genetically characterized, and 1549 different likely causative variants in 142 genes were identified. The most common phenotype encountered is retinitis pigmentosa (RP) (55.6% of families, 2447/4403). The most recurrently mutated genes were PRPH2, ABCA4 and RS1 in autosomal dominant (AD), autosomal recessive (AR) and X-linked (XL) NON-RP cases, respectively; RHO, USH2A and RPGR in AD, AR and XL for non-syndromic RP; and USH2A and MYO7A in syndromic IRD. Pathogenic variants c.3386G > T (p.Arg1129Leu) in ABCA4 and c.2276G > T (p.Cys759Phe) in USH2A were the most frequent variants identified. Our study provides the general landscape for IRD in Spain, reporting the largest cohort ever presented. Our results have important implications for genetic diagnosis, counselling and new therapeutic strategies to both the Spanish population and other related populations.

CSVS, a crowdsourcing database of the Spanish population genetic variability

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

Genotype-Phenotype Correlations in a Spanish Cohort of 506 Families With Biallelic ABCA4 Pathogenic Variants

PURPOSE

To define genotype-phenotype correlations in the largest cohort study worldwide of patients with biallelic ABCA4 variants, including 434 patients with Stargardt disease (STGD1) and 72 with cone-rod dystrophy (CRD).

DESIGN

Cohort study.

METHODS

We characterized 506 patients with ABCA4 variants using conventional genetic tools and next-generation sequencing technologies. Medical history and ophthalmologic data were obtained from 372 patients. Genotype-phenotype correlation studies were carried out for the following variables: variant type, age at symptom onset (AO), and clinical phenotype.

RESULTS

A total of 228 different pathogenic variants were identified in 506 ABCA4 patients, 50 of which were novel. Genotype-phenotype correlations showed that most of the patients with biallelic truncating variants presented with CRD and that these cases had a significantly earlier AO than patients with STGD1. Three missense variants are associated with CRD for the first time (c.1804C>T; p.[Arg602Trp], c.3056C>T; p.[Thr1019Met], and c.6320G>C; p.[Arg2107Pro]). Analysis of the most prevalent ABCA4 variant in Spain, c.3386G>T; p.(Arg1129Leu), revealed that is correlated to STGD1, later AO, and foveal sparing.

CONCLUSIONS

Our study, conducted in the largest ABCA4-associated disease cohort reported to date, updates the genotype-phenotype model established for ABCA4 variants and broadens the mutational spectrum of the gene. According to our observations, patients with ABCA4 presenting with 2 truncating variants may first present features of STGD1 but eventually develop rod dysfunction, and specific missense variants may be associated with a different phenotype, underscoring the importance of an accurate genetic diagnosis. Also, it is a prerequisite for enrollment in clinical trials, and to date, no other treatment has been approved for STGD1.

Expanding the phenotype of CRYAA nucleotide variants to a complex presentation of anterior segment dysgenesis

Background

Mutations in CRYAA, which encodes the α-crystallin protein, are associated with a spectrum of congenital cataract–microcornea syndromes.

Results

In this study, we performed clinical examination and subsequent genetic analysis in two unrelated sporadic cases of different geographical origins presenting with a complex phenotype of ocular malformation. Both cases manifested bilateral microphthalmia and severe anterior segment dysgenesis, primarily characterized by congenital aphakia, microcornea, and iris hypoplasia/aniridia. NGS-based analysis revealed two novel single nucleotide variants occurring de novo and affecting the translation termination codon of the CRYAA gene, c.520T > C and c.521A > C. Both variants are predicted to elongate the C-terminal protein domain by one-third of the original length.

Conclusions

Our report not only expands the mutational spectrum of CRYAA but also identifies the genetic cause of the unusual ocular phenotype described in this report.

Exome sequencing identifies PEX6 mutations in three cases diagnosed with Retinitis Pigmentosa and hearing impairment

Purpose

The aim of the present work is the molecular diagnosis of three patients with deafness and retinal degeneration.

Methods

Three patients from two unrelated families were initially analyzed with custom gene panels for Usher genes, non-syndromic hearing loss, or inherited syndromic retinopathies and further investigated by means of clinical or whole exome sequencing.

Results

The study allowed us to detect likely pathogenic variants in PEX6, a gene typically involved in peroxisomal biogenesis disorders (PBDs). Beside deaf-blindness, both families showed additional features: Siblings from Family 1 showed enamel alteration and abnormal peroxisome. In addition, the brother had mild neurodevelopmental delay and nephrolithiasis. The case II:1 from Family 2 showed intellectual disability, enamel alteration, and dysmorphism.

Conclusions

We have reported three new cases with pathogenic variants in PEX6 presenting with milder forms of the Zellweger spectrum disorders (ZSD). The three cases showed distinct clinical features. Thus, expanding the phenotypic spectrum of PBDs and ascertaining exome sequencing is an effective strategy for an accurate diagnosis of clinically overlapping and genetically heterogeneous disorders such as deafness-blindness association.

The Spectrum of PAX6 Mutations and Genotype-Phenotype Correlations in the Eye

The transcription factor PAX6 is essential in ocular development in vertebrates, being considered the master regulator of the eye. During eye development, it is essential for the correct patterning and formation of the multi-layered optic cup and it is involved in the developing lens and corneal epithelium. In adulthood, it is mostly expressed in cornea, iris, and lens. PAX6 is a dosage-sensitive gene and it is highly regulated by several elements located upstream, downstream, and within the gene. There are more than 500 different mutations described to affect PAX6 and its regulatory regions, the majority of which lead to PAX6 haploinsufficiency, causing several ocular and systemic abnormalities. Aniridia is an autosomal dominant disorder that is marked by the complete or partial absence of the iris, foveal hypoplasia, and nystagmus, and is caused by heterozygous PAX6 mutations. Other ocular abnormalities have also been associated with PAX6 changes, and genotype-phenotype correlations are emerging. This review will cover recent advancements in PAX6 regulation, particularly the role of several enhancers that are known to regulate PAX6 during eye development and disease. We will also present an updated overview of the mutation spectrum, where an increasing number of mutations in the non-coding regions have been reported. Novel genotype-phenotype correlations will also be discussed.

Expanded Phenotypic Spectrum of Retinopathies Associated with Autosomal Recessive and Dominant Mutations in PROM1

PURPOSE

To describe the genetic and phenotypic characteristics of a cohort of patients with PROM1 variants.

DESIGN

Case-case study.

METHODS

We screened a cohort of 2216 families with inherited retinal dystrophies using classical molecular techniques and next-generation sequencing approaches. The clinical histories of 25 patients were reviewed to determine age of onset of symptoms and the results of ophthalmoscopy, best-corrected visual acuity, full-field electroretinography, and visual field studies. Fundus autofluorescence and spectral-domain optical coherence tomography were further assessed in 7 patients.

RESULTS

PROM1 variants were identified in 32 families. Disease-causing variants were found in 18 autosomal recessive and 4 autosomal dominant families. Monoallelic pathogenic variants or variants of unknown significance were identified in the remaining 10 families. Comprehensive phenotyping of 25 patients from 22 families carrying likely disease-causing variants revealed clinical heterogeneity associated with the PROM1 gene. Most of these patients presented cone-rod dystrophy and some exhibited macular dystrophy or retinitis pigmentosa, while all presented with macular damage. Phenotypic association of a dominant splicing variant with late-onset mild maculopathy was established. This variant is one of the 3 likely founder variants identified in our Spanish cohort.

CONCLUSIONS

We report the largest cohort of patients with PROM1 variants, describing in detail the phenotype in 25 of them. Interestingly, within the variability of phenotypes related to this gene, macular involvement is a common feature in all patients.

Diabetes-mediated promotion of colon mucosa carcinogenesis is associated with mitochondrial dysfunction

Type 2 diabetes mellitus (T2DM) has been associated with an increased risk of cancer, including colon cancer (CC). However, we recently reported no influence of T2DM on CC prognosis, suggesting that any effect might be at the early stages of tumor development. We hypothesized that T2DM may create an environment in the healthy tissue, which acts as a carcinogenesis driver in agreement with the field of cancerization concept. Here, we focused on early carcinogenesis by analyzing paired tumor and normal colonic mucosa samples from the same patients. The proteome of CC and paired mucosa was quantitatively analyzed in 28 individuals (12 diabetics and 16 nondiabetics) by mass spectrometry with isobaric labeling. Out of 3076 identified proteins, 425 were differentially expressed at the tumor in diabetics compared with nondiabetics. In the adjacent mucosa, 143 proteins were differentially expressed in diabetics and nondiabetics. An enrichment analysis of this signature pointed to mitochondria, ribosome, and translation. Only six proteins were upregulated by diabetes both in tumor and mucosa, of which five were mitochondrial proteins. Differential expression in diabetic versus nondiabetic mucosa was confirmed for MRPL53, MRPL18, and TIMM8B. Higher levels of MRPL18, TIMM8B, and EIF1A were also found in normal colon epithelial cells exposed to high‐glucose conditions. We conclude that T2DM is associated with specific molecular changes in the normal mucosa of CC patients, consistent with field of cancerization in a diabetic environment. The mitochondrial protein signature identifies a potential therapeutic target that could underlie the higher risk of CC in diabetics.

Genomic Landscape of Sporadic Retinitis Pigmentosa: Findings from 877 Spanish Cases

PURPOSE

We aimed to unravel the molecular basis of sporadic retinitis pigmentosa (sRP) in the largest cohort reported to date.

DESIGN

Case series.

PARTICIPANTS

A cohort of 877 unrelated Spanish sporadic cases with a clinical diagnosis of retinitis pigmentosa (RP) and negative family history.

METHODS

The cohort was studied by classic genotyping or targeted next-generation sequencing (NGS). Multiplex ligation-dependent probe amplification (MLPA) and array-based comparative genomic hybridization were performed to confirm copy number variations detected by NGS. Quantitative fluorescent polymerase chain reaction was assessed in sRP cases carrying de novo variants to confirm paternity.

MAIN OUTCOME MEASURES

The study of the sRP cohort showed a high proportion of causal autosomal dominant (AD) and X-linked (XL) variants, most of them being de novo.

RESULTS

Causative variants were identified in 38% of the patients studied, segregating recessively in 84.5% of the solved cases. Biallelic variants detected in only 6 different autosomal recessive genes explained 50% of the cases characterized. Causal AD and XL variants were found in 7.6% and 7.9% of cases, respectively. Remarkably, 20 de novo variants were confirmed after trio analysis, explaining 6% of the cases. In addition, 17% of the solved sRP cases were reclassified to a different retinopathy phenotype.

CONCLUSIONS

This study highlights the clinical utility of NGS testing for sRP cases, expands the mutational spectrum, and provides accurate prevalence of mutated genes. Our findings evidence the underestimated role of de novo variants in the etiology of RP, emphasizing the importance of segregation analysis as well as comprehensive screening of genes carrying XL and AD variants in sporadic cases. Such in-depth study is essential for accurate family counseling and future enrollment in gene therapy-based treatments.

Molecular evidence of field cancerization initiated by diabetes in colon cancer patients

The potential involvement of type 2 diabetes mellitus (T2DM) as a risk factor for colon cancer (CC) has been previously reported. While several clinical studies show a higher incidence of CC and a lower survival rate in diabetics, others report no association. Our own experience indicates that diabetes does not seem to worsen the prognosis once the tumor is present. Despite this controversy, there are no wide‐spectrum molecular studies that delve into the impact of T2DM‐related mechanisms in colon carcinogenesis. Here, we present a transcriptomic and proteomic profiling of paired tumor and normal colon mucosa samples in a cohort of 42 CC patients, 23 of which have T2DM. We used gene set enrichment and network approaches to extract relevant pathways in diabetics, referenced them to current knowledge, and tested them using in vitro techniques. Through our transcriptomics approach, we identified an unexpected overlap of pathways overrepresented in diabetics compared to nondiabetics, in both tumor and normal mucosa, including diabetes‐related metabolic and signaling processes. Proteomic approaches highlighted several cancer‐related signaling routes in diabetics found only in normal mucosa, not in tumors. An integration of the transcriptome and proteome analyses suggested the deregulation of key pathways related to colon carcinogenesis which converged on tumor initiation axis TEAD/YAP‐TAZ as a potential initiator of the process. In vitro studies confirmed upregulation of this pathway in nontumor colon cells under high‐glucose conditions. In conclusion, T2DM associates with deregulation of cancer‐related processes in normal colon mucosa adjacent to tissue which has undergone a malignant transformation. These data support that in diabetic patients, the local microenvironment in normal colon mucosa may be a factor driving field cancerization promoting carcinogenesis. Our results set a new framework to study links between diabetes and colon cancer, including a new role of the TEAD/YAP‐TAZ complex as a potential driver.

Toward the Mutational Landscape of Autosomal Dominant Retinitis Pigmentosa: A Comprehensive Analysis of 258 Spanish Families

Purpose

To provide a comprehensive overview of the molecular basis of autosomal dominant retinitis pigmentosa (adRP) in Spanish families. Thus, we established the molecular characterization rate, gene prevalence, and mutational spectrum in the largest European cohort reported to date.

Methods

A total of 258 unrelated Spanish families with a clinical diagnosis of RP and suspected autosomal dominant inheritance were included. Clinical diagnosis was based on complete ophthalmologic examination and family history. Retrospective and prospective analysis of Spanish adRP families was carried out using a combined strategy consisting of classic genetic techniques and next-generation sequencing (NGS) for single-nucleotide variants and copy number variation (CNV) screening.

Results

Overall, 60% of our families were genetically solved. Interestingly, 3.1% of the cohort carried pathogenic CNVs. Disease-causing variants were found in an autosomal dominant gene in 55% of the families; however, X-linked and autosomal recessive forms were also identified in 3% and 2%, respectively. Four genes (RHO, PRPF31, RP1, and PRPH2) explained up to 62% of the solved families. Missense changes were most frequently found in adRP-associated genes; however, CNVs represented a relevant disease cause in PRPF31- and CRX-associated forms.

Conclusions

Implementation of NGS technologies in the adRP study clearly increased the diagnostic yield compared with classic approaches. Our study outcome expands the spectrum of disease-causing variants, provides accurate data on mutation gene prevalence, and highlights the implication of CNVs as important contributors to adRP etiology.

Parental Mosaicism in PAX6 Causes Intra-Familial Variability: Implications for Genetic Counseling of Congenital Aniridia and Microphthalmia

Mutations in PAX6 are involved in several developmental eye disorders. These disorders have considerable phenotypic variability, ranging from panocular forms of congenital aniridia and microphthalmia to isolated anomalies of the anterior or posterior segment. Here, we describe 3 families with variable inter-generational ocular expression of aniridia, iris coloboma, or microphthalmia, and an unusual transmission of PAX6 mutations from an unaffected or mildly affected parent; all of which raised suspicion of gonosomal mosaicism. We first identified two previously known nonsense mutations and one novel likely pathogenic missense variant in PAX6 in probands by means of targeted NGS. The subsequent segregation analysis by Sanger sequencing evidenced the presence of highly probable mosaic events in paternal blood samples. Mosaicism was further confirmed by droplet digital PCR analysis in several somatic tissues of mosaic fathers. Quantification of the mutant allele fraction in parental samples showed a marked deviation from 50%, with a range between 12 and 29% depending on cell type. Gonosomal mosaicsm was definitively confirmed in one of the families thanks to the availability of a sperm sample from the mosaic father. Thus, the recurrence risk in this family was estimated to be about one-third. This is the first report confirming parental PAX6 mosaicism as a cause of disease recurrence in aniridia and other related phenotypes. In addition, we demonstrated that post-zygotic mosaicism is a frequent and underestimated pathogenic mechanism in aniridia, explaining intra-familial phenotypic variability in many cases. Our findings may have substantial implications for genetic counseling in congenital aniridia. Thus, we also highlight the importance of comprehensive genetic screening of parents for new sporadic cases with aniridia or related developmental eye disease to more accurately assess recurrence risk. In conclusion, somatic and/or gonosomal mosaicism should be taken into consideration as a genetic factor to explain not only families with unaffected parents despite multiple affected children but also variable expressivity, apparent de novo cases, and even uncharacterized cases of aniridia and related developmental eye disorders, apparently lacking PAX6 mutations.

Implication of non-coding PAX6 mutations in aniridia

There is an increasing implication of non-coding regions in pathological processes of genetic origin. This is partly due to the emergence of sophisticated techniques that have transformed research into gene expression by allowing a more global understanding of the genome, both at the genomic, epigenomic and chromatin levels. Here, we implemented the analysis of PAX6, whose coding loss-of-function variants are mainly implied in aniridia, by studying its non-coding regions (untranslated regions, introns and cis-regulatory sequences). In particular, we have taken advantage of the development of high-throughput approaches to screen the upstream and downstream regulatory regions of PAX6 in 47 aniridia patients without identified mutation in the coding sequence. This was made possible through the use of custom targeted resequencing and/or CGH array to analyze the entire PAX6 locus on 11p13. We found candidate variants in 30 of the 47 patients. 9/30 correspond to the well-known described 3' deletions encompassing SIMO and other enhancer elements. In addition, we identified numerous different variants in various non-coding regions, in particular untranslated regions. Among these latter, most of them demonstrated an in vitro functional effect using a minigene strategy, and 12/21 are thus considered as causative mutations or very likely to explain the phenotypes. This new analysis strategy brings molecular diagnosis to more than 90% of our aniridia patients. This study revealed an outstanding mutation pattern in non-coding PAX6 regions confirming that PAX6 remains the major gene for aniridia.

Expanding the phenotype of the X-linked BCOR microphthalmia syndromes

Two distinct syndromes arise from pathogenic variants in the X-linked gene BCOR (BCL-6 corepressor): oculofaciocardiodental (OFCD) syndrome, which affects females, and a severe microphthalmia ('Lenz'-type) syndrome affecting males. OFCD is an X-linked dominant syndrome caused by a variety of BCOR null mutations. As it manifests only in females, it is presumed to be lethal in males. The severe male X-linked recessive microphthalmia syndrome ('Lenz') usually includes developmental delay in addition to the eye findings and is caused by hypomorphic BCOR variants, mainly by a specific missense variant c.254C > T, p.(Pro85Leu). Here, we detail 16 new cases (11 females with 4 additional, genetically confirmed, affected female relatives; 5 male cases each with unaffected carrier mothers). We describe new variants and broaden the phenotypic description for OFCD to include neuropathy, muscle hypotonia, pituitary underdevelopment, brain atrophy, lipoma and the first description of childhood lymphoma in an OFCD case. Our male X-linked recessive cases show significant new phenotypes: developmental delay (without eye anomalies) in two affected half-brothers with a novel BCOR variant, and one male with high myopia, megalophthalmos, posterior embryotoxon, developmental delay, and heart and bony anomalies with a previously undescribed BCOR splice site variant. Our female OFCD cases and their affected female relatives showed variable features, but consistently had early onset cataracts. We show that a mosaic carrier mother manifested early cataract and dental anomalies. All female carriers of the male X-linked recessive cases for whom genetic confirmation was available showed skewed X-inactivation and were unaffected. In view of the extended phenotype, we suggest a new term of X-linked BCOR-related syndrome.

Unravelling the pathogenic role and genotype-phenotype correlation of the USH2A p.(Cys759Phe) variant among Spanish families

INTRODUCTION

Mutations in USH2A cause both isolated Retinitis Pigmentosa (RP) and Usher syndrome (that implies RP and hearing impairment). One of the most frequent variants identified in this gene and among these patients is the p.(Cys759Phe) change. However, the pathogenic role of this allele has been questioned since it was found in homozygosity in two healthy siblings of a Spanish family. To assess the causative role of USH2A p.(Cys759Phe) in autosomal recessive RP (ARRP) and Usher syndrome type II (USH2) and to establish possible genotype-phenotype correlations associated with p.(Cys759Phe), we performed a comprehensive genetic and clinical study in patients suffering from any of the two above-mentioned diseases and carrying at least one p.(Cys759Phe) allele.

MATERIALS AND METHODS

Diagnosis was set according to previously reported protocols. Genetic analyses were performed by using classical molecular and Next-Generation Sequencing approaches. Probands of 57 unrelated families were molecularly studied and 63 patients belonging to these families were phenotypically evaluated.

RESULTS

Molecular analysis characterized 100% of the cases, identifying: 11 homozygous patients for USH2A p.(Cys759Phe), 42 compound heterozygous patients (12 of them with another missense USH2A pathogenic variant and 30 with a truncating USH2A variant), and 4 patients carrying the p.(Cys759Phe) allele and a pathogenic variant in another RP gene (PROM1, CNGB1 or RP1). No additional causative variants were identified in symptomatic homozygous patients. Statistical analysis of clinical differences between zygosity states yielded differences (p≤0.05) in age at diagnosis of RP and hypoacusis, and progression of visual field loss. Homozygosity of p.(Cys759Phe) and compound heterozygosity with another USH2A missense variant is associated with ARRP or ARRP plus late onset hypoacusis (OR = 20.62, CI = 95%, p = 0.041).

CONCLUSIONS

The present study supports the role of USH2A p.(Cys759Phe) in ARRP and USH2 pathogenesis, and demonstrates the clinical differences between different zygosity states. Phenotype-genotype correlations may guide the genetic characterization based upon specific clinical signs and may advise on the clinical management and prognosis based upon a specific genotype.

Combining targeted panel-based resequencing and copy-number variation analysis for the diagnosis of inherited syndromic retinopathies and associated ciliopathies

Inherited syndromic retinopathies are a highly heterogeneous group of diseases that involve retinal anomalies and systemic manifestations. They include retinal ciliopathies, other well-defined clinical syndromes presenting with retinal alterations and cases of non-specific multisystemic diseases. The heterogeneity of these conditions makes molecular and clinical characterization of patients challenging in daily clinical practice. We explored the capacity of targeted resequencing and copy-number variation analysis to improve diagnosis of a heterogeneous cohort of 47 patients mainly comprising atypical cases that did not clearly fit a specific clinical diagnosis. Thirty-three likely pathogenic variants were identified in 18 genes (ABCC6, ALMS1, BBS1, BBS2, BBS12, CEP41, CEP290, IFT172, IFT27, MKKS, MYO7A, OTX2, PDZD7, PEX1, RPGRIP1, USH2A, VPS13B, and WDPCP). Molecular findings and additional clinical reassessments made it possible to accurately characterize 14 probands (30% of the total). Notably, clinical refinement of complex phenotypes was achieved in 4 cases, including 2 de novo OTX2-related syndromes, a novel phenotypic association for the ciliary CEP41 gene, and the co-existence of biallelic USH2A variants and a Koolen-de-Vries syndrome-related 17q21.31 microdeletion. We demonstrate that combining next-generation sequencing and CNV analysis is a comprehensive and useful approach to unravel the extensive phenotypic and genotypic complexity of inherited syndromic retinopathies.

New GJA8 variants and phenotypes highlight its critical role in a broad spectrum of eye anomalies

GJA8 encodes connexin 50 (Cx50), a transmembrane protein involved in the formation of lens gap junctions. GJA8 mutations have been linked to early onset cataracts in humans and animal models. In mice, missense mutations and homozygous Gja8 deletions lead to smaller lenses and microphthalmia in addition to cataract, suggesting that Gja8 may play a role in both lens development and ocular growth. Following screening of GJA8 in a cohort of 426 individuals with severe congenital eye anomalies, primarily anophthalmia, microphthalmia and coloboma, we identified four known [p.(Thr39Arg), p.(Trp45Leu), p.(Asp51Asn), and p.(Gly94Arg)] and two novel [p.(Phe70Leu) and p.(Val97Gly)] likely pathogenic variants in seven families. Five of these co-segregated with cataracts and microphthalmia, whereas the variant p.(Gly94Arg) was identified in an individual with congenital aphakia, sclerocornea, microphthalmia and coloboma. Four missense variants of unknown or unlikely clinical significance were also identified. Furthermore, the screening of GJA8 structural variants in a subgroup of 188 individuals identified heterozygous 1q21 microdeletions in five families with coloboma and other ocular and/or extraocular findings. However, the exact genotype-phenotype correlation of these structural variants remains to be established. Our data expand the spectrum of GJA8 variants and associated phenotypes, confirming the importance of this gene in early eye development.

Identification of PITX3 mutations in individuals with various ocular developmental defects

BACKGROUND

Congenital cataract displays large phenotypic (syndromic and isolated cataracts) and genetic heterogeneity. Mutations in several transcription factors involved in eye development, like PITX3, have been associated with congenital cataracts and anterior segment mesenchymal disorders.

MATERIALS AND METHODS

Targeted sequencing of 187 genes involved in ocular development was performed in 96 patients with mainly anophthalmia and microphthalmia. Additionally, Sanger sequencing analysis of PITX3 was performed on a second cohort of 32 index cases with congenital cataract and Peters anomaly and/or sclereocornea.

RESULTS

We described five families with four different PITX3 mutations, two of which were novel. In Family 1, the heterozygous recurrent c.640_656dup (p.Gly220Profs*95) mutation cosegregated with eye anomalies ranging from congenital cataract to Peters anomaly. In Family 2, the novel c.669del [p.(Leu225Trpfs*84)] mutation cosegregated with dominantly inherited eye anomalies ranging from posterior embryotoxon to congenital cataract in heterozygous carriers and congenital sclereocornea and cataract in a patient homozygous for this mutation. In Family 3, we identified the recurrent heterozygous c.640_656dup (p.Gly220Profs*95) mutation segregating with congenital cataract. In Family 4, the de novo c.582del [p.(Ile194Metfs*115)] mutation was identified in a patient with congenital cataract, microphthalmia, developmental delay and autism. In Family 5, the c.38G>A (p.Ser13Asn) mutation segregated dominantly in a family with Peters anomaly, which is a novel phenotype associated with the c.38G>A variant compared with the previously reported isolated congenital cataract.

CONCLUSIONS

Our study unveils different phenotypes associated with known and novel mutations in PITX3, which will improve the genetic counselling of patients and their families.

Establishment of a human DOA 'plus' iPSC line, IISHDOi003-A, with the mutation in the OPA1 gene: c.1635C>A; p.Ser545Arg

We have generated a human iPSC line IISHDOi003-A from fibroblasts of a patient with a dominant optic atrophy 'plus' phenotype, harbouring a heterozygous mutation, c.1635C>A; p.Ser545Arg, in the OPA1 gene. Reprogramming factors Oct3/4, Sox2, Klf4, and c-Myc were delivered using Sendai virus.

Analysis of the PRPF31 Gene in Spanish Autosomal Dominant Retinitis Pigmentosa Patients: A Novel Genomic Rearrangement

Purpose

The aim was to determine the prevalence of PRPF31 mutations in a cohort of Spanish autosomal dominant retinitis pigmentosa (adRP) families to deepen knowledge of the pathogenic mechanisms underlying the disease and to assess genotype-phenotype correlations.

Methods

A cohort of 211 adRP patients was screened for variants in PRPF31 by using a combined strategy comprising next-generation sequencing approaches and copy-number variation (CNV) analysis. Quantitative RT-PCR and CNV analysis of the regulatory MSR1 element were also performed to assess PRPF31 gene expression. Phenotype was assessed by using ophthalmologic examination protocols.

Results

Fifteen different causative mutations and genomic rearrangements were identified, revealing five novel mutations. Prevalence of PRPF31 mutations, genomic rearrangements, and lack of penetrance were 7.6%, 1.9%, and 66.7%, respectively. Interestingly, we identified a tandem duplication and a partial PRPF31 deletion in different affected individuals from the same family. PRPF31 gene expression was significantly decreased in symptomatic cases carrying either PRPF31 duplication or deletion as compared to controls. The 4 MSR1 allele in cis with the PRPF31 wild-type allele was apparently a protective factor. The mutated phenotype varied from no symptoms to typical retinitis pigmentosa with variable onset and course depending on the kind of mutation, with the duplication case the most severe.

Conclusions

In view of the high genetic heterogeneity of PRPF31 mutations, the screening must include the entire gene, as well as CNV assays, to detect large rearrangements. This is the first report of a variable phenotype correlation as well as a gross duplication and deletion within the same family.

2017 update on the relationship between diabetes and colorectal cancer: epidemiology, potential molecular mechanisms and therapeutic implications

Worldwide deaths from diabetes mellitus (DM) and colorectal cancer increased by 90% and 57%, respectively, over the past 20 years. The risk of colorectal cancer was estimated to be 27% higher in patients with type 2 DM than in non-diabetic controls. However, there are potential confounders, information from lower income countries is scarce, across the globe there is no correlation between DM prevalence and colorectal cancer incidence and the association has evolved over time, suggesting the impact of additional environmental factors. The clinical relevance of these associations depends on understanding the mechanism involved. Although evidence is limited, insulin use has been associated with increased and metformin with decreased incidence of colorectal cancer. In addition, colorectal cancer shares some cellular and molecular pathways with diabetes target organ damage, exemplified by diabetic kidney disease. These include epithelial cell injury, activation of inflammation and Wnt/β-catenin pathways and iron homeostasis defects, among others. Indeed, some drugs have undergone clinical trials for both cancer and diabetic kidney disease. Genome-wide association studies have identified diabetes-associated genes (e.g. TCF7L2) that may also contribute to colorectal cancer. We review the epidemiological evidence, potential pathophysiological mechanisms and therapeutic implications of the association between DM and colorectal cancer. Further studies should clarify the worldwide association between DM and colorectal cancer, strengthen the biological plausibility of a cause-and-effect relationship through characterization of the molecular pathways involved, search for specific molecular signatures of colorectal cancer under diabetic conditions, and eventually explore DM-specific strategies to prevent or treat colorectal cancer.