Novel genes and sex differences in COVID-19 severity.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (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 < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-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.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) 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.1 × 10-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.

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.

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.

VIsoQLR: an interactive tool for the detection, quantification and fine-tuning of isoforms in selected genes using long-read sequencing.

DNA variants altering the pre-mRNA splicing process represent an underestimated cause of human genetic diseases. Their association with disease traits should be confirmed using functional assays from patient cell lines or alternative models to detect aberrant mRNAs. Long-read sequencing is a suitable technique to identify and quantify mRNA isoforms. Available isoform detection and/or quantification tools are generally designed for the whole transcriptome analysis. However experiments focusing on genes of interest need more precise data fine-tuning and visualization tools.Here we describe VIsoQLR, an interactive analyzer, viewer and editor for the semi-automated identification and quantification of known and novel isoforms using long-read sequencing data. VIsoQLR is tailored to thoroughly analyze mRNA expression in splicing assays of selected genes. Our tool takes sequences aligned to a reference, and for each gene, it defines consensus splice sites and quantifies isoforms. VIsoQLR introduces features to edit the splice sites through dynamic and interactive graphics and tables, allowing accurate manual curation. Known isoforms detected by other methods can also be imported as references for comparison. A benchmark against two other popular transcriptome-based tools shows VIsoQLR accurate performance on both detection and quantification of isoforms. Here, we present VIsoQLR principles and features and its applicability in a case study example using nanopore-based long-read sequencing. VIsoQLR is available at https://github.com/TBLabFJD/VIsoQLR .

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.

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.

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/.

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

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

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

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

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.

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.

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.

IL-6-based mortality prediction model for COVID-19: Validation and update in multicenter and second wave cohorts.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is a highly variable condition. Validated tools to assist in the early detection of patients at high risk of mortality can help guide medical decisions. OBJECTIVE: We sought to validate externally, as well as in patients from the second pandemic wave in Europe, our previously developed mortality prediction model for hospitalized COVID-19 patients. METHODS: Three validation cohorts were generated: 2 external with 185 and 730 patients from the first wave and 1 internal with 119 patients from the second wave. The probability of death was calculated for all subjects using our prediction model, which includes peripheral blood oxygen saturation/fraction of inspired oxygen ratio, neutrophil-to-lymphocyte ratio, lactate dehydrogenase, IL-6, and age. Discrimination and calibration were evaluated in the validation cohorts. The prediction model was updated by reestimating individual risk factor effects in the overall cohort (N = 1477). RESULTS: The mortality prediction model showed good performance in the external validation cohorts 1 and 2, and in the second wave validation cohort 3 (area under the receiver-operating characteristic curve, 0.94, 0.86, and 0.86, respectively), with excellent calibration (calibration slope, 0.86, 0.94, and 0.79; intercept, 0.05, 0.03, and 0.10, respectively). The updated model accurately predicted mortality in the overall cohort (area under the receiver-operating characteristic curve, 0.91), which included patients from both the first and second COVID-19 waves. The updated model was also useful to predict fatal outcome in patients without respiratory distress at the time of evaluation. CONCLUSIONS: This is the first COVID-19 mortality prediction model validated in patients from the first and second pandemic waves. The COR+12 online calculator is freely available to facilitate its implementation (https://utrero-rico.shinyapps.io/COR12_Score/).

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies.

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort's IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

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.

Gene Correction Recovers Phagocytosis in Retinal Pigment Epithelium Derived from Retinitis Pigmentosa-Human-Induced Pluripotent Stem Cells.

Hereditary retinal dystrophies (HRD) represent a significant cause of blindness, affecting mostly retinal pigment epithelium (RPE) and photoreceptors (PRs), and currently suffer from a lack of effective treatments. Highly specialized RPE and PR cells interact mutually in the functional retina, therefore primary HRD affecting one cell type leading to a secondary HRD in the other cells. Phagocytosis is one of the primary functions of the RPE and studies have discovered that mutations in the phagocytosis-associated gene Mer tyrosine kinase receptor (MERTK) lead to primary RPE dystrophy. Treatment strategies for this rare disease include the replacement of diseased RPE with healthy autologous RPE to prevent PR degeneration. The generation and directed differentiation of patient-derived human-induced pluripotent stem cells (hiPSCs) may provide a means to generate autologous therapeutically-relevant adult cells, including RPE and PR. However, the continued presence of the MERTK gene mutation in patient-derived hiPSCs represents a significant drawback. Recently, we reported the generation of a hiPSC model of MERTK-associated Retinitis Pigmentosa (RP) that recapitulates disease phenotype and the subsequent creation of gene-corrected RP-hiPSCs using Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9. In this study, we differentiated gene-corrected RP-hiPSCs into RPE and found that these cells had recovered both wild-type MERTK protein expression and the lost phagocytosis of fluorescently-labeled photoreceptor outer segments observed in uncorrected RP-hiPSC-RPE. These findings provide proof-of-principle for the utility of gene-corrected hiPSCs as an unlimited cell source for personalized cell therapy of rare vision disorders.

Fine Breakpoint Mapping by Genome Sequencing Reveals the First Large X Inversion Disrupting the NHS Gene in a Patient with Syndromic Cataracts.

Inversions are structural variants that are generally balanced. However, they could lead to gene disruptions or have positional effects leading to diseases. Mutations in the NHS gene cause Nance-Horan syndrome, an X-linked disorder characterised by congenital cataracts and dental anomalies. Here, we aimed to characterise a balanced pericentric inversion X(p22q27), maternally inherited, in a child with syndromic bilateral cataracts by breakpoint mapping using whole-genome sequencing (WGS). 30× Illumina paired-end WGS was performed in the proband, and breakpoints were confirmed by Sanger sequencing. EdU assays and FISH analysis were used to assess skewed X-inactivation patterns. RNA expression of involved genes in the breakpoint boundaries was evaluated by droplet-digital PCR. We defined the breakpoint position of the inversion at Xp22.13, with a 15 bp deletion, disrupting the unusually large intron 1 of the canonical NHS isoform, and also perturbing topologically-associated domains (TADs). Moreover, a microhomology region of 5 bp was found on both sides. RNA analysis confirmed null and reduced NHS expression in the proband and his unaffected mother, respectively. In conclusion, we report the first chromosomal inversion disrupting NHS, fine-mapped by WGS. Our data expand the clinical spectrum and the pathogenic mechanisms underlying the NHS defects.

CPAMD8 loss-of-function underlies non-dominant congenital glaucoma with variable anterior segment dysgenesis and abnormal extracellular matrix.

Abnormal development of the ocular anterior segment may lead to a spectrum of clinical phenotypes ranging from primary congenital glaucoma (PCG) to variable anterior segment dysgenesis (ASD). The main objective of this study was to identify the genetic alterations underlying recessive congenital glaucoma with ASD (CG-ASD). Next-generation DNA sequencing identified rare biallelic CPAMD8 variants in four patients with CG-ASD and in one case with PCG. CPAMD8 is a gene of unknown function and recently associated with ASD. Bioinformatic and in vitro functional evaluation of the variants using quantitative reverse transcription PCR and minigene analysis supported a loss-of-function pathogenic mechanism. Optical and electron microscopy of the trabeculectomy specimen from one of the CG-ASD cases revealed an abnormal anterior chamber angle, with altered extracellular matrix, and apoptotic trabecular meshwork cells. The CPAMD8 protein was immunodetected in adult human ocular fluids and anterior segment tissues involved in glaucoma and ASD (i.e., aqueous humor, non-pigmented ciliary epithelium, and iris muscles), as well as in periocular mesenchyme-like cells of zebrafish embryos. CRISPR/Cas9 disruption of this gene in F0 zebrafish embryos (96 hpf) resulted in varying degrees of gross developmental abnormalities, including microphthalmia, pharyngeal maldevelopment, and pericardial and periocular edemas. Optical and electron microscopy examination of these embryos showed iridocorneal angle hypoplasia (characterized by altered iris stroma cells, reduced anterior chamber, and collagen disorganized corneal stroma extracellular matrix), recapitulating some patients' features. Our data support the notion that CPAMD8 loss-of-function underlies a spectrum of recessive CG-ASD phenotypes associated with extracellular matrix disorganization and provide new insights into the normal and disease roles of this gene.

Novel PXDN biallelic variants in patients with microphthalmia and anterior segment dysgenesis.

Microphthalmia, anophthalmia, and anterior segment dysgenesis are severe ocular developmental defects. There is a wide genetic heterogeneity leading to these ocular malformations. By using whole genome, exome and targeted sequencing in patients with ocular developmental anomalies, six biallelic pathogenic variants (including five novel variants) were identified in the PXDN gene in four families with microphthalmia and anterior segment dysgenesis. Only 11 different mutations (11 families) have been described in this gene to date. The phenotype of these patients is variable in severity, ranging from cataract and developmental glaucoma to complex microphthalmia. Interestingly, two unrelated patients of our series presented with an ocular phenotype including aniridia and microspherophakia. However, despite various phenotypic presentations and types of mutations, no genotype-phenotype correlation could be made. Thus, this work improves our knowledge of the recessive phenotype associated with biallelic variants in this gene and highlights the importance of screening PXDN in patients with anterior segment dysgenesis with or without microphthalmia.