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.

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.

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.

De novo small deletion affecting transcription start site of short isoform of AUTS2 gene in a patient with syndromic neurodevelopmental defects.

Disruption of the autism susceptibility candidate 2 (AUTS2) gene through genomic rearrangements, copy number variations (CNVs), and intragenic deletions and mutations, has been recurrently involved in syndromic forms of developmental delay and intellectual disability, known as AUTS2 syndrome. The AUTS2 gene plays an important role in regulation of neuronal migration, and when altered, associates with a variable phenotype from severely to mildly affected patients. The more severe phenotypes significantly correlate with the presence of defects affecting the C-terminus part of the gene. This article reports a new patient with a syndromic neurodevelopmental disorder, who presents a deletion of 30 nucleotides in the exon 9 of the AUTS2 gene. Importantly, this deletion includes the transcription start site for the AUTS2 short transcript isoform, which has an important role in brain development. Gene expression analysis of AUTS2 full-length and short isoforms revealed that the deletion found in this patient causes a remarkable reduction in the expression level, not only of the short isoform, but also of the full AUTS2 transcripts. This report adds more evidence for the role of mutated AUTS2 short transcripts in the development of a severe phenotype in the AUTS2 syndrome.

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.

Mitotic Catastrophe Induced in HeLa Tumor Cells by Photodynamic Therapy with Methyl-aminolevulinate.

Photodynamic therapy (PDT) constitutes a cancer treatment modality based on the administration of a photosensitizer, which accumulates in tumor cells. The subsequent irradiation of the tumoral area triggers the formation of reactive oxygen species responsible for cancer cell death. One of the compounds approved in clinical practice is methyl-aminolevulinate (MAL), a protoporphyrin IX (PpIX) precursor intermediate of heme synthesis. We have identified the mitotic catastrophe (MC) process after MAL-PDT in HeLa human carcinoma cells. The fluorescence microscopy revealed that PpIX was located mainly at plasma membrane and lysosomes of HeLa cells, although some fluorescence was also detected at endoplasmic reticulum and Golgi apparatus. Cell blockage at metaphase-anaphase transition was observed 24 h after PDT by phase contrast microscopy and flow cytometry. Mitotic apparatus components evaluation by immunofluorescence and Western blot indicated: multipolar spindles and disorganized chromosomes in the equatorial plate accompanied with dispersion of centromeres and alterations in aurora kinase proteins. The mitotic blockage induced by MAL-PDT resembled that induced by two compounds used in chemotherapy, taxol and nocodazole, both targeting microtubules. The alterations in tumoral cells provided evidence of MC induced by MAL-PDT, resolving mainly by apoptosis, directly or through the formation of multinucleate cells.

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

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

An evaluation of pipelines for DNA variant detection can guide a reanalysis protocol to increase the diagnostic ratio of genetic diseases.

Clinical exome (CE) sequencing has become a first-tier diagnostic test for hereditary diseases; however, its diagnostic rate is around 30-50%. In this study, we aimed to increase the diagnostic yield of CE using a custom reanalysis algorithm. Sequencing data were available for three cohorts using two commercial protocols applied as part of the diagnostic process. Using these cohorts, we compared the performance of general and clinically relevant variant calling and the efficacy of an in-house bioinformatic protocol (FJD-pipeline) in detecting causal variants as compared to commercial protocols. On the whole, the FJD-pipeline detected 99.74% of the causal variants identified by the commercial protocol in previously solved cases. In the unsolved cases, FJD-pipeline detects more INDELs and non-exonic variants, and is able to increase the diagnostic yield in 2.5% and 3.2% in the re-analysis of 78 cancer and 62 cardiovascular cases. These results were considered to design a reanalysis, filtering and prioritization algorithm that was tested by reassessing 68 inconclusive cases of monoallelic autosomal recessive retinal dystrophies increasing the diagnosis by 4.4%. In conclusion, a guided NGS reanalysis of unsolved cases increases the diagnostic yield in genetic disorders, making it a useful diagnostic tool in medical genetics.

Cellular intrinsic factors involved in the resistance of squamous cell carcinoma to photodynamic therapy.

Photodynamic therapy (PDT) is widely used to treat non-melanoma skin cancer. However, some patients affected with squamous cell carcinoma (SCC) do not respond adequately to PDT with methyl-δ-aminolevulinic acid (MAL-PDT) and the tumors acquire an infiltrative phenotype and became histologically more aggressive, less differentiated, and more fibroblastic. To search for potential factors implicated in SCC resistance to PDT, we have used the SCC-13 cell line (parental) and resistant SCC-13 cells obtained by repeated MAL-PDT treatments (5th and 10th PDT-resistant generations). Xenografts assays in immunodeficient mice showed that the tumors generated by resistant cells were bigger than those induced by parental cells. Comparative genomic hybridization array (aCGH) showed that the three cell types presented amplicons in 3p12.1 CADM2, 7p11.2 EFGR, and 11q13.3 CCND1 genes. The 5th and 10th PDT-resistant cells showed an amplicon in 5q11.2 MAP3K1, which was not present in parental cells. The changes detected by aCGH on CCND1, EFGR, and MAP3K1 were confirmed in extracts of SCC-13 cells by reverse-transcriptase PCR and by western blot, and by immunohistochemistry in human biopsies from persistent tumors after MAL-PDT. Our data suggest that genomic imbalances related to CCND1, EFGR, and particularly MAP3K1 seem to be involved in the development of the resistance of SCC to PDT.