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.

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.

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.

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.

CLAPO syndrome: identification of somatic activating PIK3CA mutations and delineation of the natural history and phenotype.

PURPOSE: CLAPO syndrome is a rare vascular disorder characterized by capillary malformation of the lower lip, lymphatic malformation predominant on the face and neck, asymmetry, and partial/generalized overgrowth. Here we tested the hypothesis that, although the genetic cause is not known, the tissue distribution of the clinical manifestations in CLAPO seems to follow a pattern of somatic mosaicism. METHODS: We clinically evaluated a cohort of 13 patients with CLAPO and screened 20 DNA blood/tissue samples from 9 patients using high-throughput, deep sequencing. RESULTS: We identified five activating mutations in the PIK3CA gene in affected tissues from 6 of the 9 patients studied; one of the variants (NM_006218.2:c.248T>C; p.Phe83Ser) has not been previously described in developmental disorders. CONCLUSION: We describe for the first time the presence of somatic activating PIK3CA mutations in patients with CLAPO. We also report an update of the phenotype and natural history of the syndrome.

Age-dependent association of clonal hematopoiesis with COVID-19 mortality in patients over 60 years.

Clonal hematopoiesis, especially that of indeterminate potential (CHIP), has been associated with age-related diseases, such as those contributing to a more severe COVID-19. Four studies have attempted to associate CHIP with COVID-19 severity without conclusive findings. In the present work, we explore the association between CHIP and COVID-19 mortality. Genomic DNA extracted from peripheral blood of COVID-19 patients (n = 241 deceased, n = 239 survivors) was sequenced with the Myeloid Solutions™ panel of SOPHiA Genetics. The association between clonality and age and clonality and mortality was studied using logistic regression models adjusted for sex, ethnicity, and comorbidities. The association with mortality was performed with patients stratified into four groups of age according to the quartiles of the distribution: 60-74 years, 75-84 years, 85-91 years, and 92-101 years. Clonality was found in 38% of the cohort. The presence of CHIP variants, but not the number, significantly increased with age in the entire cohort of COVID-19 patients, as well as in the group of survivors (p < 0.001). When patients were stratified by age and the analysis adjusted, CHIP classified as pathogenic/likely pathogenic was significantly more represented in deceased patients compared with survivors in the group of 75-84 years (34.6% vs 13.7%, p = 0.020). We confirmed the well-established linear relationship between age and clonality in the cohort of COVID-19 patients and found a significant association between pathogenic/likely pathogenic CHIP and mortality in patients from 75 to 84 years that needs to be further validated.

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.

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.

Derivation of a human DOA iPSC line, IISHDOi006-A, with a mutation in the ACO2 gene: c.1999G>A; p.Glu667Lys.

Human iPSC line, IISHDOi006-A, was obtained from fibroblasts of a patient with Dominant Optic Atrophy (DOA) carrying a heterozygous mutation in the gene ACO2: c.1999G>A; p.Glu667Lys. Reprogramming factors Oct3/4, Sox2, Klf4, and c-Myc were delivered using a non-integrative methodology that involves the use of Sendai virus.