Deciphering novel TCF4-driven mechanisms underlying a common triplet repeat expansion-mediated disease.

Fuchs endothelial corneal dystrophy (FECD) is an age-related cause of vision loss, and the most common repeat expansion-mediated disease in humans characterised to date. Up to 80% of European FECD cases have been attributed to expansion of a non-coding CTG repeat element (termed CTG18.1) located within the ubiquitously expressed transcription factor encoding gene, TCF4. The non-coding nature of the repeat and the transcriptomic complexity of TCF4 have made it extremely challenging to experimentally decipher the molecular mechanisms underlying this disease. Here we comprehensively describe CTG18.1 expansion-driven molecular components of disease within primary patient-derived corneal endothelial cells (CECs), generated from a large cohort of individuals with CTG18.1-expanded (Exp+) and CTG 18.1-independent (Exp-) FECD. We employ long-read, short-read, and spatial transcriptomic techniques to interrogate expansion-specific transcriptomic biomarkers. Interrogation of long-read sequencing and alternative splicing analysis of short-read transcriptomic data together reveals the global extent of altered splicing occurring within Exp+ FECD, and unique transcripts associated with CTG18.1-expansions. Similarly, differential gene expression analysis highlights the total transcriptomic consequences of Exp+ FECD within CECs. Furthermore, differential exon usage, pathway enrichment and spatial transcriptomics reveal TCF4 isoform ratio skewing solely in Exp+ FECD with potential downstream functional consequences. Lastly, exome data from 134 Exp- FECD cases identified rare (minor allele frequency <0.005) and potentially deleterious (CADD>15) TCF4 variants in 7/134 FECD Exp- cases, suggesting that TCF4 variants independent of CTG18.1 may increase FECD risk. In summary, our study supports the hypothesis that at least two distinct pathogenic mechanisms, RNA toxicity and TCF4 isoform-specific dysregulation, both underpin the pathophysiology of FECD. We anticipate these data will inform and guide the development of translational interventions for this common triplet-repeat mediated disease.

Disease-Causing TIMP3 Variants and Deep Phenotyping of Two Czech Families with Sorsby Fundus Dystrophy Associated with Novel p.(Tyr152Cys) Mutation.

We aim to report the ocular phenotype and molecular genetic findings in two Czech families with Sorsby fundus dystrophy and to review all the reported TIMP3 pathogenic variants. Two probands with Sorsby fundus dystrophy and three first-degree relatives underwent ocular examination and retinal imaging, including optical coherence tomography angiography. The DNA of the first proband was screened using a targeted ocular gene panel, while, in the second proband, direct sequencing of the TIMP3 coding region was performed. Sanger sequencing was also used for segregation analysis within the families. All the previously reported TIMP3 variants were reviewed using the American College of Medical Genetics and the Association for Molecular Pathology interpretation framework. A novel heterozygous variant, c.455A>G p.(Tyr152Cys), in TIMP3 was identified in both families and potentially de novo in one. Optical coherence tomography angiography documented in one patient the development of a choroidal neovascular membrane at 54 years. Including this study, 23 heterozygous variants in TIMP3 have been reported as disease-causing. Application of gene-specific criteria denoted eleven variants as pathogenic, eleven as likely pathogenic, and one as a variant of unknown significance. Our study expands the spectrum of TIMP3 pathogenic variants and highlights the importance of optical coherence tomography angiography for early detection of choroidal neovascular membranes.

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.

Current methods for monitoring Pseudomonas syringae biofilm development.

This work reviews biofilm investigation techniques and highlights the benefits and drawbacks of each approach focusing especially on Pseudomonas syringae and may serve as a comprehensive guide for any early-career researchers starting with the topic of biofilm. Each approach with applications of individual microscopy and spectroscopy techniques is summarized together with characterization of Pseudomonas syringae and its role in pathogenesis.

LEGO-lipophosphonoxins: length of hydrophobic module affects permeabilizing activity in target membranes of different phospholipid composition.

In the past few decades, society has faced rapid development and spreading of antimicrobial resistance due to antibiotic misuse and overuse and the immense adaptability of bacteria. Difficulties in obtaining effective antimicrobial molecules from natural sources challenged scientists to develop synthetic molecules with antimicrobial effect. We developed modular molecules named LEGO-Lipophosphonoxins (LEGO-LPPO) capable of inducing cytoplasmic membrane perforation. In this structure-activity relationship study we focused on the role of the LEGO-LPPO hydrophobic module directing the molecule insertion into the cytoplasmic membrane. We selected three LEGO-LPPO molecules named C9, C8 and C7 differing in the length of their hydrophobic chain and consisting of an alkenyl group containing one double bond. The molecule with the long hydrophobic chain (C9) was shown to be the most effective with the lowest MIC and highest perforation rate both in vivo and in vitro. We observed high antimicrobial activity against both G+ and G- bacteria with significant differences in LEGO-LPPOs mechanism of action on these two cell types. We observed a highly cooperative mechanism of LEGO-LPPO action on G- bacteria as well as on liposomes resembling G- bacteria. LEGO-LPPO action on G- bacteria was significantly slower compared to G+ bacteria suggesting the role of the outer membrane in affecting the LEGO-LPPOs perforation rate. This notion was supported by the higher sensitivity of the E. coli strain with a compromised outer membrane. Finally, we noted that the composition of the cytoplasmic membrane affects the activity of LEGO-LPPOs since the presence of phosphatidylethanolamine increases their membrane disrupting activity.

Lisch Epithelial Corneal Dystrophy Is Caused by Heterozygous Loss-of-Function Variants in MCOLN1.

PURPOSE: To report the genetic etiology of Lisch epithelial corneal dystrophy (LECD). DESIGN: Multicenter cohort study. METHODS: A discovery cohort of 27 individuals with LECD from 17 families, including 7 affected members from the original LECD family, 6 patients from 2 new families and 14 simplex cases, was recruited. A cohort of 6 individuals carrying a pathogenic MCOLN1 (mucolipin 1) variant was reviewed for signs of LECD. Next-generation sequencing or targeted Sanger sequencing were used in all patients to identify pathogenic or likely pathogenic variants and penetrance of variants. RESULTS: Nine rare heterozygous MCOLN1 variants were identified in 23 of 27 affected individuals from 13 families. The truncating nature of 7 variants and functional testing of 1 missense variant indicated that they result in MCOLN1 haploinsufficiency. Importantly, in the homozygous and compound-heterozygous state, 4 of 9 LECD-associated variants cause the rare lysosomal storage disorder mucolipidosis IV (MLIV). Autosomal recessive MLIV is a systemic disease and comprises neurodegeneration as well as corneal opacity of infantile-onset with epithelial autofluorescent lysosomal inclusions. However, the 6 parents of 3 patients with MLIV confirmed to carry pathogenic MCOLN1 variants did not have the LECD phenotype, suggesting MCOLN1 haploinsufficiency may be associated with reduced penetrance and variable expressivity. CONCLUSIONS: MCOLN1 haploinsufficiency is the major cause of LECD. Based on the overlapping clinical features of corneal epithelial cells with autofluorescent inclusions reported in both LECD and MLIV, it is concluded that some carriers of MCOLN1 haploinsufficiency-causing variants present with LECD.

Ovol2 promoter mutations in mice and human illuminate species-specific phenotypic divergence.

Pathogenic variants in the highly conserved OVOL2 promoter region cause posterior polymorphous corneal dystrophy (PPCD) 1 by inducing an ectopic expression of the endothelial OVOL2 mRNA. Here we produced an allelic series of Ovol2 promoter mutations in the mouse model including the heterozygous c.-307T>C variant (RefSeq NM_021220.4) causing PPCD1 in humans. Despite the high evolutionary conservation of the Ovol2 promoter, only some alterations of its sequence had phenotypic consequences in mice. Four independent sequence variants in the distal part of the Ovol2 promoter had no significant effect on endothelial Ovol2 mRNA level or caused any ocular phenotype. In contrast, the mutation c.-307T>C resulted in increased Ovol2 expression in the corneal endothelium. However, only a small fraction of adult mice c.-307T>C heterozygotes developed ocular phenotypes such as irido-corneal adhesions, and corneal opacity. Interestingly, phenotypic penetrance was increased at embryonic stages. Notably, c.-307T>C mutation is located next to the Ovol1/Ovol2 transcription factor binding site. Mice carrying an allele with a deletion encompassing the Ovol2 binding site c.-307_-320del showed significant Ovol2 gene upregulation in the cornea endothelium and exhibited phenotypes similar to the c.-307T>C mutation. In conclusion, although the mutations c.-307T>C and -307_-320del lead to a comparably strong increase in endothelial Ovol2 expression as seen in PPCD1 patients, endothelial dystrophy was not observed in the mouse model, implicating species-specific differences in endothelial cell biology. Nonetheless, the emergence of dominant ocular phenotypes associated with Ovol2 promoter variants in mice implies a potential role of this gene in eye development and disease.

Dysfunction of peripheral somatic and autonomic nervous system in patients with severe forms of Crohn's disease on biological therapy with TNFα inhibitors-A single center study.

OBJECTIVE: Crohn's disease (CD) can be associated with a wide range of extraintestinal manifestations (EIMs), including neurological ones. Published studies differ in their conclusions about the epidemiology and etiopathogenesis of neurological EIMs. The aims of this study were to demonstrate the presence and find risk factors of peripheral (somatic and autonomic) neuropathy patients with severe CD on anti-TNFα biological therapy. MATERIAL AND METHODS: A clinical examination focusing on detection of peripheral sensor-motor nervous dysfunction (including Sudoscan) and examination of autonomic nervous system dysfunction (using Ewing´s battery tests and spectral analysis) together with laboratory tests and collection of demographic data followed by administration of questionnaires were performed on a total of 30 neurologically asymptomatic outpatients with severe CD on anti-TNFα biological therapy. RESULTS: Peripheral sensor-motor nervous function via clinical neurological examination was pathological in 36.7% and Sudoscan in 33.3% of cases. Statistically significant associations between vibration perception test and age, CD and biological therapy duration, body mass index and Crohn's Disease Activity Index were proved while statistically significant associations between temperature perception test and age and BMI were proved as well. Additionally, a decrease of total protein in a patient´s serum below the physiological cut-off in the 6 months prior to measurement was associated with a pathological result of a Sudoscan. Cardiovascular autonomic neuropathy based on Ewing´s battery tests was present in 56.7% of patients, no statistically significant risk factors were found. Our peripheral neuropathy questionnaire correlated with the results of the Sudoscan test and some tests of the clinical examination of peripheral sensor-motor nervous function (discriminatory contact perception test, temperature perception test). CONCLUSIONS: This study demonstrated a relatively high prevalence of peripheral (especially autonomic) neuropathy and verified some risk factors for the development of peripheral somatic neuropathy in asymptomatic patients with severe form of CD on anti-TNFα biological therapy.

Targeted sequencing and in vitro splice assays shed light on ABCA4-associated retinopathies missing heritability.

The ABCA4 gene is the most frequently mutated Mendelian retinopathy-associated gene. Biallelic variants lead to a variety of phenotypes, however, for thousands of cases the underlying variants remain unknown. Here, we aim to shed further light on the missing heritability of ABCA4-associated retinopathy by analyzing a large cohort of macular dystrophy probands. A total of 858 probands were collected from 26 centers, of whom 722 carried no or one pathogenic ABCA4 variant, while 136 cases carried two ABCA4 alleles, one of which was a frequent mild variant, suggesting that deep-intronic variants (DIVs) or other cis-modifiers might have been missed. After single molecule molecular inversion probes (smMIPs)-based sequencing of the complete 128-kb ABCA4 locus, the effect of putative splice variants was assessed in vitro by midigene splice assays in HEK293T cells. The breakpoints of copy number variants (CNVs) were determined by junction PCR and Sanger sequencing. ABCA4 sequence analysis solved 207 of 520 (39.8%) naive or unsolved cases and 70 of 202 (34.7%) monoallelic cases, while additional causal variants were identified in 54 of 136 (39.7%) probands carrying two variants. Seven novel DIVs and six novel non-canonical splice site variants were detected in a total of 35 alleles and characterized, including the c.6283-321C>G variant leading to a complex splicing defect. Additionally, four novel CNVs were identified and characterized in five alleles. These results confirm that smMIPs-based sequencing of the complete ABCA4 gene provides a cost-effective method to genetically solve retinopathy cases and that several rare structural and splice altering defects remain undiscovered in Stargardt disease cases.

Light-responsive microRNA molecules in human retinal organoids are differentially regulated by distinct wavelengths of light.

Cells in the human retina must rapidly adapt to constantly changing visual stimuli. This fast adaptation to varying levels and wavelengths of light helps to regulate circadian rhythms and allows for adaptation to high levels of illumination, thereby enabling the rest of the visual system to remain responsive. It has been shown that retinal microRNA (miRNA) molecules play a key role in regulating these processes. However, despite extensive research using various model organisms, light-regulated miRNAs in human retinal cells remain unknown. Here, we aim to characterize these miRNAs. We generated light-responsive human retinal organoids that express miRNA families and clusters typically found in the retina. Using an in-house developed photostimulation device, we identified a subset of light-regulated miRNAs. Importantly, we found that these miRNAs are differentially regulated by distinct wavelengths of light and have a rapid turnover, highlighting the dynamic and adaptive nature of the human retina.

Exome sequencing reveals IFT172 variants in patients with non-syndromic cholestatic liver disease.

BACKGROUND AND AIM: Gene defects contribute to the aetiology of intrahepatic cholestasis. We aimed to explore the outcome of whole-exome sequencing (WES) in a cohort of 51 patients with this diagnosis. PATIENTS AND METHODS: Both paediatric (n = 33) and adult (n = 18) patients with cholestatic liver disease of unknown aetiology were eligible. WES was used for reassessment of 34 patients (23 children) without diagnostic genotypes in ABCB11, ATP8B1, ABCB4 or JAG1 demonstrable by previous Sanger sequencing, and for primary assessment of additional 17 patients (10 children). Nasopharyngeal swab mRNA was analysed to address variant pathogenicity in two families. RESULTS: WES revealed biallelic variation in 3 ciliopathy genes (PKHD1, TMEM67 and IFT172) in 4 clinically unrelated index subjects (3 children and 1 adult), heterozygosity for a known variant in PPOX in one adult index subject, and homozygosity for an unreported splice-site variation in F11R in one child. Whereas phenotypes of the index patients with mutated PKHD1, TMEM67, and PPOX corresponded with those elsewhere reported, how F11R variation underlies liver disease remains unclear. Two unrelated patients harboured different novel biallelic variants in IFT172, a gene implicated in short-rib thoracic dysplasia 10 and Bardet-Biedl syndrome 20. One patient, a homozygote for IFT172 rs780205001 c.167A>C p.(Lys56Thr) born to first cousins, had liver disease, interpreted on biopsy aged 4y as glycogen storage disease, followed by adult-onset nephronophthisis at 25y. The other, a compound heterozygote for novel frameshift variant IFT172 NM_015662.3 c.2070del p.(Met690Ilefs*11) and 2 syntenic missense variants IFT172 rs776310391 c.157T>A p.(Phe53Ile) and rs746462745 c.164C>G p.(Thr55Ser), had a severe 8mo cholestatic episode in early infancy, with persisting hyperbilirubinemia and fibrosis on imaging studies at 17y. No patient had skeletal malformations. CONCLUSION: Our findings suggest association of IFT172 variants with non-syndromic cholestatic liver disease.

MIR204 n.37C>T variant as a cause of chorioretinal dystrophy variably associated with iris coloboma, early-onset cataracts and congenital glaucoma.

Four members of a three-generation Czech family with early-onset chorioretinal dystrophy were shown to be heterozygous carriers of the n.37C>T in MIR204. The identification of this previously reported pathogenic variant confirms the existence of a distinct clinical entity caused by a sequence change in MIR204. Chorioretinal dystrophy was variably associated with iris coloboma, congenital glaucoma, and premature cataracts extending the phenotypic range of the condition. In silico analysis of the n.37C>T variant revealed 713 novel targets. Additionally, four family members were shown to be affected by albinism resulting from biallelic pathogenic OCA2 variants. Haplotype analysis excluded relatedness with the original family reported to harbour the n.37C>T variant in MIR204. Identification of a second independent family confirms the existence of a distinct MIR204-associated clinical entity and suggests that the phenotype may also involve congenital glaucoma.

Case report: A rare variant m.4135T>C in the MT-ND1 gene leads to Leber hereditary optic neuropathy and altered respiratory chain supercomplexes.

Leber hereditary optic neuropathy is a primary mitochondrial disease characterized by acute visual loss due to the degeneration of retinal ganglion cells. In this study, we describe a patient carrying a rare missense heteroplasmic variant in MT-ND1, NC_012920.1:m.4135T>C (p.Tyr277His) manifesting with a typical bilateral painless decrease of the visual function, triggered by physical exercise or higher ambient temperature. Functional studies in muscle and fibroblasts show that amino acid substitution Tyr277 with His leads to only a negligibly decreased level of respiratory chain complex I (CI), but the formation of supercomplexes and the activity of the enzyme are disturbed noticeably. Our data indicate that although CI is successfully assembled in the patient's mitochondria, its function is hampered by the m.4135T>C variant, probably by stabilizing CI in its inactive form. We conclude that the m.4135T>C variant together with a combination of external factors is necessary to manifest the phenotype.

Development of the Pseudomonas syringae pv. morsprunorum Biofilm Monitored in Real Time Using Attenuated Total Reflection Fourier Transform Infrared Measurements in a Flow Cell Chamber.

Biofilms of sessile Pseudomonas syringae cells formed on top of plant host's leaves or fruits allow surviving harsh environmental conditions (desiccation) and improve their resistance to antibacterial treatments of crops. A better understanding of these biofilms can help minimize their effect on harvests. In the present study, infrared attenuated total reflection spectroscopy coupled with optical and confocal laser scanning microscopy has been applied for the first time to analyze Pseudomonas syringae pathovar morsprunorum biofilm development in real time. The biofilm development was observed within a spectral window 4000-800 cm-1 under constant flow conditions for 72 h. The kinetics of representative integrated band areas (nucleic acids with polysaccharides at 1141-1006 cm-1, amino acid side chains with free fatty acids at 1420-1380 cm-1, proteins at 1580-1490 cm-1, and lipids with proteins at 2935-2915 cm-1) were analyzed with regard to the observed biofilm structure and the following P. syringae biofilm developmental stages were attributed: The inoculation phase, washing of weakly attached bacteria closely followed by recolonization of the vacated surface, the restructuration phase, and finally the maturation phase.

Phenotype and genotype of concurrent keratoconus and Fuchs endothelial corneal dystrophy.

PURPOSE: To characterise the phenotype and genotype of concurrent keratoconus and Fuchs endothelial corneal dystrophy (KC + FECD). METHODS: We recruited 20 patients with concurrent KC + FECD for a retrospective observational case series from the United Kingdom and the Czech Republic. We compared eight parameters of corneal shape (Pentacam, Oculus) with two groups of age-matched controls who had either isolated keratoconus (KC) or isolated FECD. We genotyped probands for an intronic triplet TCF4 repeat expansion (CTG18.1) and the ZEB1 variant c.1920G >T p.(Gln640His). RESULTS: The median age at diagnosis of patients with KC + FECD was 54 (interquartile range 46 to 66) years, with no evidence of KC progression (median follow-up 84 months, range 12 to 120 months). The mean (standard deviation (SD)) of the minimum corneal thickness, 493 (62.7) µm, was greater than eyes with KC, 458 (51.1) µm, but less than eyes with FECD, 590 (55.6) µm. Seven other parameters of corneal shape were more like KC than FECD. Seven (35%) probands with KC + FECD had a TCF4 repeat expansion of ≥50 compared to five controls with isolated FECD. The average of the largest TCF4 expansion in cases with KC + FECD (46 repeats, SD 36 repeats) was similar to the age-matched controls with isolated FECD (36 repeats, SD 28 repeats; p = 0.299). No patient with KC + FECD harboured the ZEB1 variant. CONCLUSIONS: The KC + FECD phenotype is consistent with KC but with superimposed stromal swelling from endothelial disease. The proportion of cases with a TCF4 expansion is similar in concurrent KC + FECD and age-matched controls with isolated FECD.

Whole genome sequencing for USH2A-associated disease reveals several pathogenic deep-intronic variants that are amenable to splice correction.

A significant number of individuals with a rare disorder such as Usher syndrome (USH) and (non-)syndromic autosomal recessive retinitis pigmentosa (arRP) remain genetically unexplained. Therefore, we assessed subjects suspected of USH2A-associated disease and no or mono-allelic USH2A variants using whole genome sequencing (WGS) followed by an improved pipeline for variant interpretation to provide a conclusive diagnosis. One hundred subjects were screened using WGS to identify causative variants in USH2A or other USH/arRP-associated genes. In addition to the existing variant interpretation pipeline, a particular focus was put on assessing splice-affecting properties of variants, both in silico and in vitro. Also structural variants were extensively addressed. For variants resulting in pseudoexon inclusion, we designed and evaluated antisense oligonucleotides (AONs) using minigene splice assays and patient-derived photoreceptor precursor cells. Biallelic variants were identified in 49 of 100 subjects, including novel splice-affecting variants and structural variants, in USH2A or arRP/USH-associated genes. Thirteen variants were shown to affect USH2A pre-mRNA splicing, including four deep-intronic USH2A variants resulting in pseudoexon inclusion, which could be corrected upon AON treatment. We have shown that WGS, combined with a thorough variant interpretation pipeline focused on assessing pre-mRNA splicing defects and structural variants, is a powerful method to provide subjects with a rare genetic condition, a (likely) conclusive genetic diagnosis. This is essential for the development of future personalized treatments and for patients to be eligible for such treatments.

Axenfeld-Rieger syndrome: more than meets the eye.

BACKGROUND: Axenfeld-Rieger syndrome (ARS) is characterised by typical anterior segment anomalies, with or without systemic features. The discovery of causative genes identified ARS subtypes with distinct phenotypes, but our understanding is incomplete, complicated by the rarity of the condition. METHODS: Genetic and phenotypic characterisation of the largest reported ARS cohort through comprehensive genetic and clinical data analyses. RESULTS: 128 individuals with causative variants in PITX2 or FOXC1, including 81 new cases, were investigated. Ocular anomalies showed significant overlap but with broader variability and earlier onset of glaucoma for FOXC1-related ARS. Systemic anomalies were seen in all individuals with PITX2-related ARS and the majority of those with FOXC1-related ARS. PITX2-related ARS demonstrated typical umbilical anomalies and dental microdontia/hypodontia/oligodontia, along with a novel high rate of Meckel diverticulum. FOXC1-related ARS exhibited characteristic hearing loss and congenital heart defects as well as previously unrecognised phenotypes of dental enamel hypoplasia and/or crowding, a range of skeletal and joint anomalies, hypotonia/early delay and feeding disorders with structural oesophageal anomalies in some. Brain imaging revealed highly penetrant white matter hyperintensities, colpocephaly/ventriculomegaly and frequent arachnoid cysts. The expanded phenotype of FOXC1-related ARS identified here was found to fully overlap features of De Hauwere syndrome. The results were used to generate gene-specific management plans for the two types of ARS. CONCLUSION: Since clinical features of ARS vary significantly based on the affected gene, it is critical that families are provided with a gene-specific diagnosis, PITX2-related ARS or FOXC1-related ARS. De Hauwere syndrome is proposed to be a FOXC1opathy.

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

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

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy.

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Snail Track Lesion with Flat Keratometry in Anterior Segment Dysgenesis Caused by a Novel FOXC1 Variant.

We report the phenotype of a 15-year-old female patient with anterior segment dysgenesis (ASD) caused by a novel heterozygous loss-of-function FOXC1 variant. The proband underwent an ophthalmic examination as well as a molecular genetic investigation comprising exome sequencing, a single nucleotide polymorphism array to access copy number and Sanger sequencing to exclude non-coding causal variants. There was bilateral mild iris hypoplasia with pupil deformation and iridocorneal adhesions. In addition to these features of ASD, the corneas were flat, with mean keratometry readings of 38.8 diopters in the right eye and 39.5 diopters in the left eye. There was a snail track lesion of the left cornea at the level of the Descemet membrane. The central corneal endothelial cell density was reduced bilaterally at 1964 and 1373 cells/mm2 in the right and left eyes, respectively. Molecular genetic analysis revealed that the proband was a carrier of a novel heterozygous frameshifting variant in FOXC1, c.605del p.(Pro202Argfs*113). Neither parent had this change, suggesting a de novo origin which was supported by paternity testing. We found no possibly pathogenic variants in the other genes associated with posterior corneal dystrophies or ASD. Further studies are warranted to verify whether there is a true association between snail track lesions, corneal flattening, and pathogenic variants in FOXC1.