Minigene-Based Splice Assays Reveal the Effect of Non-Canonical Splice Site Variants in USH2A

Non-canonical splice site variants are increasingly recognized as a relevant cause of the USH2A-associated diseases, non-syndromic autosomal recessive retinitis pigmentosa and Usher syndrome type 2. Many non-canonical splice site variants have been reported in public databases, but an effect on pre-mRNA splicing has only been functionally verified for a subset of these variants. In this study, we aimed to extend the knowledge regarding splicing events by assessing a selected set of USH2A non-canonical splice site variants and to study their potential pathogenicity. Eleven non-canonical splice site variants were selected based on four splice prediction tools. Ten different USH2A constructs were generated and minigene splice assays were performed in HEK293T cells. An effect on pre-mRNA splicing was observed for all 11 variants. Various events, such as exon skipping, dual exon skipping and partial exon skipping were observed and eight of the tested variants had a full effect on splicing as no conventionally spliced mRNA was detected. We demonstrated that non-canonical splice site variants in USH2A are an important contributor to the genetic etiology of the associated disorders. This type of variant generally should not be neglected in genetic screening, both in USH2A-associated disease as well as other hereditary disorders. In addition, cases with these specific variants may now receive a conclusive genetic diagnosis.     

The Study of a 231 French Patient Cohort Significantly Extends the Mutational Spectrum of the Two Major Usher Genes MYO7A and USH2A

Usher syndrome is an autosomal recessive disorder characterized by congenital hearing loss combined with retinitis pigmentosa, and in some cases, vestibular areflexia. Three clinical subtypes are distinguished, and MYO7A and USH2A represent the two major causal genes involved in Usher type I, the most severe form, and type II, the most frequent form, respectively. Massively parallel sequencing was performed on a cohort of patients in the context of a molecular diagnosis to confirm clinical suspicion of Usher syndrome. We report here 231 pathogenic MYO7A and USH2A genotypes identified in 73 Usher type I and 158 Usher type II patients. Furthermore, we present the ACMG classification of the variants, which comprise all types. Among them, 68 have not been previously reported in the literature, including 12 missense and 16 splice variants. We also report a new deep intronic variant in USH2A. Despite the important number of molecular studies published on these two genes, we show that during the course of routine genetic diagnosis, undescribed variants continue to be identified at a high rate. This is particularly pertinent in the current era, where therapeutic strategies based on DNA or RNA technologies are being developed.     

Investigating Biomarkers for USH2A Retinopathy Using Multimodal Retinal Imaging

Pathogenic mutations in USH2A are a leading cause of visual loss secondary to non-syndromic or Usher syndrome-associated retinitis pigmentosa (RP). With an increasing number of RP-targeted clinical trials in progress, we sought to evaluate the photoreceptor topography underlying patterns of loss observed on clinical retinal imaging to guide surrogate endpoint selection in USH2A retinopathy. In this prospective cross-sectional study, twenty-five patients with molecularly confirmed USH2A-RP underwent fundus autofluorescence (FAF), spectral-domain optical coherence tomography (SD-OCT) and adaptive optics scanning laser ophthalmoscopy (AOSLO) retinal imaging. Analysis comprised measurement of FAF horizontal inner (IR) and outer (OR) hyperautofluorescent ring diameter; SD-OCT ellipsoid zone (EZ) and external limiting membrane (ELM) width, normalised EZ reflectance; AOSLO foveal cone density and intact macular photoreceptor mosaic (IMPM) diameter. Thirty-two eyes from 16 patients (mean age ± SD, 36.0 ± 14.2 years) with USH2A-associated Usher syndrome type 2 (n = 14) or non-syndromic RP (n = 2) met the inclusion criteria. Spatial alignment was observed between IR-EZ and OR-ELM diameters/widths (p < 0.001). The IMPM border occurred just lateral to EZ loss (p < 0.001), although sparser intact photoreceptor inner segments were detected until ELM disruption. EZ width and IR diameter displayed a biphasic relationship with cone density whereby slow cone loss occurred until retinal degeneration reached ~1350 μm from the fovea, beyond which greater reduction in cone density followed. Normalised EZ reflectance and cone density were significantly associated (p < 0.001). As the strongest correlate of cone density (p < 0.001) and best-corrected visual acuity (p < 0.001), EZ width is the most sensitive biomarker of structural and functional decline in USH2A retinopathy, rendering it a promising trial endpoint.     

Characteristics of Retinitis Pigmentosa Associated with ADGRV1 and Comparison with USH2A in Patients from a Multicentric Usher Syndrome Study Treatrush

In contrast to USH2A, variants in ADGRV1 are a minor cause of Usher syndrome type 2, and the associated phenotype is less known. The purpose of the study was to characterize the retinal phenotype of 18 ADGRV1 patients (9 male, 9 female; median age 52 years) and compare it with that of 204 USH2A patients (111 male, 93 female; median age 43 years) in terms of nyctalopia onset, best corrected visual acuity (BCVA), fundus autofluorescence (FAF), and optical coherence tomography (OCT) features. There was no statistical difference in the median age at onset (30 and 18 years; Mann–Whitney U test, p = 0.13); the mean age when 50% of the patients reached legal blindness (≥1.0 log MAR) based on visual acuity (64 years for both groups; log-rank, p = 0.3); the risk of developing advanced retinal degeneration (patch or atrophy) with age (multiple logistic regression, p = 0.8); or the frequency of cystoid macular edema (31% vs. 26%, Fisher’s exact test, p = 0.4). ADGRV1 and USH2A retinopathy were indistinguishable in all major functional and structural characteristics, suggesting that the loss of function of the corresponding proteins produces similar effects in the retina. The results are important for counseling ADGRV1 patients, who represent the minor patient subgroup.     

Single-Nucleotide Variations,Insertions/Deletions and Copy Number Variations in Myelodysplastic Syndrome during Disease Progression Revealed by a Single-Cell DNA Sequencing Platform

Simple SummaryMyelodysplastic syndrome (MDS) is a myeloid neoplasm associated with complex clonal architecture. The application of single-cell sequencing is capable of revealing the clonal dynamics of MDS during disease progression and treatment resistance. This has advantages over bulk-tumor sequencing which is limited by its resolution. In this study, we evaluated two patients with MDS for the clonal dynamics of pathogenic mutations at the single-cell level of disease progression and resistance to hypomethylating agents (HMAs). There were two key observations. First, changes in the clonal heterogeneity of the pathogenic FLT3-ITD, IDH2, EZH2, or GATA2 mutations was associated with disease progression and resistance to HMA. Secondly, disease progression and resistance to HMA was accompanied by the acquisition of copy number variations of DNMT3A, TET2, and GATA2.AbstractMyelodysplastic syndrome (MDS) is a clonal myeloid neoplasm characterized by ineffective hematopoiesis, cytopenia, dysplasia, and clonal instability, leading to leukemic transformation. Hypomethylating agents are the mainstay of treatment in higher-risk MDS. However, treatment resistance and disease transformation into acute myeloid leukemia (AML) is observed in the majority of patients and is indicative of a dismal outcome. The residual cell clones resistant to therapy or cell clones acquiring new genetic aberrations are two of the key events responsible for drug resistance. Bulk tumor sequencing often fails to detect these rare subclones that confer resistance to therapy. In this study, we employed a single-cell DNA (sc-DNA) sequencing approach to study the clonal heterogeneity and clonal evolution in two MDS patients refractory to HMA. In both patients, different single nucleotide variations (SNVs) or insertions and deletions (INDELs) were detected with bulk tumor sequencing. Rare cell clones with mutations that are undetectable by bulk tumor sequencing were detected by sc-DNA sequencing. In addition to SNVs and short INDELs, this study also revealed the presence of a clonal copy number loss of DNMT3A, TET2, and GATA2 as standalone events or in association with the small SNVs or INDELs detected during HMA resistance and disease progression.     

Clinical Utility of a Unique Genome-Wide DNA Methylation Signature for KMT2A-Related Syndrome

Wiedemann–Steiner syndrome (WDSTS) is a Mendelian syndromic intellectual disability (ID) condition associated with hypertrichosis cubiti, short stature, and characteristic facies caused by pathogenic variants in the KMT2A gene. Clinical features can be inconclusive in mild and unusual WDSTS presentations with variable ID (mild to severe), facies (typical or not) and other associated malformations (bone, cerebral, renal, cardiac and ophthalmological anomalies). Interpretation and classification of rare KMT2A variants can be challenging. A genome-wide DNA methylation episignature for KMT2A-related syndrome could allow functional classification of variants and provide insights into the pathophysiology of WDSTS. Therefore, we assessed genome-wide DNA methylation profiles in a cohort of 60 patients with clinical diagnosis for WDSTS or Kabuki and identified a unique highly sensitive and specific DNA methylation episignature as a molecular biomarker of WDSTS. WDSTS episignature enabled classification of variants of uncertain significance in the KMT2A gene as well as confirmation of diagnosis in patients with clinical presentation of WDSTS without known genetic variants. The changes in the methylation profile resulting from KMT2A mutations involve global reduction in methylation in various genes, including homeobox gene promoters. These findings provide novel insights into the molecular etiology of WDSTS and explain the broad phenotypic spectrum of the disease.     

Rare Titin (TTN) Variants in Diseases Associated with Sudden Cardiac Death

A leading cause of death in western countries is sudden cardiac death, and can be associated with genetic disease. Next-generation sequencing has allowed thorough analysis of genes associated with this entity, including, most recently, titin. We aimed to identify potentially pathogenic genetic variants in titin. A total of 1126 samples were analyzed using a custom sequencing panel including major genes related to sudden cardiac death. Our cohort was divided into three groups: 432 cases from patients with cardiomyopathies, 130 cases from patients with channelopathies, and 564 post-mortem samples from individuals showing anatomical healthy hearts and non-conclusive causes of death after comprehensive autopsy. None of the patients included had definite pathogenic variants in the genes analyzed by our custom cardio-panel. Retrospective analysis comparing the in-house database and available public databases also was performed. We identified 554 rare variants in titin, 282 of which were novel. Seven were previously reported as pathogenic. Of these 554 variants, 493 were missense variants, 233 of which were novel. Of all variants identified, 399 were unique and 155 were identified at least twice. No definite pathogenic variants were identified in any of genes analyzed. We identified rare, mostly novel, titin variants that seem to play a potentially pathogenic role in sudden cardiac death. Additional studies should be performed to clarify the role of these variants in sudden cardiac death.     

Novel TTLL5 Variants Associated with Cone-Rod Dystrophy and Early-Onset Severe Retinal Dystrophy

Variants of the TTLL5 gene, which encodes tubulin tyrosine ligase-like family member five, are a rare cause of cone dystrophy (COD) or cone-rod dystrophy (CORD). To date, only a few TTLL5 patients have been clinically and genetically described. In this study, we report five patients harbouring biallelic variants of TTLL5. Four adult patients presented either COD or CORD with onset in the late teenage years. The youngest patient had a phenotype of early onset severe retinal dystrophy (EOSRD). Genetic analysis was performed by targeted next generation sequencing of gene panels and assessment of copy number variants (CNV). We identified eight variants, of which six were novel, including two large multiexon deletions in patients with COD or CORD, while the EOSRD patient harboured the novel homozygous p.(Trp640*) variant and three distinct USH2A variants, which might explain the observed rod involvement. Our study highlights the role of TTLL5 in COD/CORD and the importance of large deletions. These findings suggest that COD or CORD patients lacking variants in known genes may harbour CNVs to be discovered in TTLL5, previously undetected by classical sequencing methods. In addition, variable phenotypes in TTLL5-associated patients might be due to the presence of additional gene defects.     

Genetic Variants of ABC and SLC Transporter Genes and Chronic Myeloid Leukaemia: Impact on Susceptibility and Prognosis

Solute carrier (SLC) and ATP-binding cassette (ABC) transporters comprise a variety of proteins expressed on cell membranes responsible for intrusion or extrusion of substrates, respectively, including nutrients, xenobiotics, and chemotherapeutic agents. These transporters mediate the cellular disposition of tyrosine kinase inhibitors (TKIs), and their genetic variants could affect its function, potentially predisposing patients to chronic myeloid leukaemia (CML) and modulating treatment response. We explored the impact of genetic variability (single nucleotide variants—SNVs) of drug transporter genes (ABCB1, ABCG2, SLC22A1, and SLC22A5) on CML susceptibility, drug response, and BCR-ABL1 mutation status. We genotyped 10 SNVs by tetra-primers-AMRS-PCR in 198 CML patients and 404 controls, and assessed their role in CML susceptibility and prognosis. We identified five SNVs associated with CML predisposition, with some variants increasing disease risk, including TT genotype ABCB1 (rs1045642), and others showing a protective effect (GG genotype SLC22A5 rs274558). We also observed different haplotypes and genotypic profiles associated with CML predisposition. Relating to drug response impact, we found that CML patients with the CC genotype (rs2231142 ABCG2) had an increased risk of TKI resistance (six-fold). Additionally, CML patients carrying the CG genotype (rs683369 SLC22A1) presented a 4.54-fold higher risk of BCR-ABL1 mutations. Our results suggest that drug transporters’ SNVs might be involved in CML susceptibility and TKI response, and predict the risk of BCR-ABL1 mutations, highlighting the impact that SNVs could have in therapeutic selection.     

Identification and Characterization of an Exonic Duplication in PALB2 in a Man with Synchronous Breast and Prostate Cancer

PALB2 (partner and localizer of BRCA2), as indicated by its name, is a BRCA2-interacting protein that plays an important role in homologous recombination (HR) and DNA double-strand break (DSB) repair. While pathogenic variants of PALB2 have been well proven to confer an increased risk of breast cancer, data on its involvement in prostate cancer (PrC) have not been clearly demonstrated. We investigated, using targeted next generation sequencing (NGS), a 59-year-old Caucasian man who developed synchronous breast and prostate cancers. This genetic investigation allowed to identify an intragenic germline heterozygous duplication in PALB2, implicating intronic repetitive sequences spanning exon 11. This variant was confirmed by multiplex ligation probe amplification (MLPA), and genomic breakpoints have been identified and characterized at the nucleotide level (c.3114-811_3202-1756dup) using an approach based on walking PCR, long range PCR, and Sanger sequencing. RT-PCR using mRNA extracted from lymphocytes and followed by Sanger sequencing revealed a tandem duplication r.3114_3201dup; p.(Gly1068Glufs * 14). This duplication results in the synthesis of a truncated, and most-likely, non-functional protein. These findings expand the phenotypic spectrum of PALB2 variants and may improve the yield of genetic diagnoses in this field.     

Thrombocytopenia-Absent Radius Syndrome: Descriptions of Three New Cases and a Novel Splicing Variant in RBM8A That Expands the Spectrum of Null Alleles

Thrombocytopenia-absent radius (TAR) syndrome is a rare congenital disorder characterized by the bilateral absence of the radius and thrombocytopenia, and sometimes by other skeletal, gastrointestinal, cardiac, and renal abnormalities. The underlying genetic defect is usually the compound inheritance of a microdeletion in 1q21.1 (null allele) and a low-frequency, non-coding single nucleotide variant (SNV) in the RBM8A gene (hypomorphic allele). We report three new cases from two unrelated families. The two siblings presented the common genotype, namely the compound heterozygosity for a 1q21.1 microdeletion and the hypomorphic SNV c.-21G>A in RBM8A, whereas the third, unrelated patient presented a rare genotype comprised by two RBM8A variants: c.-21G>A (hypomorphic allele) and a novel pathogenic variant, c.343-2A>G (null allele). Of the eight documented RBM8A variants identified in TAR syndrome patients, four have hypomorphic expression and four behave as null alleles. The present report expands the RBM8A null allele spectrum and corroborates the particularities of RBM8A involvement in TAR syndrome pathogenesis.     

New Candidates for Autism/Intellectual Disability Identified by Whole-Exome Sequencing

Intellectual disability (ID) is characterized by impairments in the cognitive processes and in the tasks of daily life. It encompasses a clinically and genetically heterogeneous group of neurodevelopmental disorders often associated with autism spectrum disorder (ASD). Social and communication abilities are strongly compromised in ASD. The prevalence of ID/ASD is 1–3%, and approximately 30% of the patients remain without a molecular diagnosis. Considering the extreme genetic locus heterogeneity, next-generation sequencing approaches have provided powerful tools for candidate gene identification. Molecular diagnosis is crucial to improve outcome, prevent complications, and hopefully start a therapeutic approach. Here, we performed parent–offspring trio whole-exome sequencing (WES) in a cohort of 60 mostly syndromic ID/ASD patients and we detected 8 pathogenic variants in genes already known to be associated with ID/ASD (SYNGAP1, SMAD6, PACS1, SHANK3, KMT2A, KCNQ2, ACTB, and POGZ). We found four de novo disruptive variants of four novel candidate ASD/ID genes: MBP, PCDHA1, PCDH15, PDPR. We additionally selected via bioinformatic tools many variants in unknown genes that alone or in combination can contribute to the phenotype. In conclusion, our data confirm the efficacy of WES in detecting pathogenic variants of known and novel ID/ASD genes.     

Peripheral Ion Channel Genes Screening in Painful Small Fiber Neuropathy

Neuropathic pain is a characteristic feature of small fiber neuropathy (SFN), which in 18% of the cases is caused by genetic variants in voltage-gated sodium ion channels. In this study, we assessed the role of fifteen other ion channels in neuropathic pain. Patients with SFN (n = 414) were analyzed for ANO1, ANO3, HCN1, KCNA2, KCNA4, KCNK18, KCNN1, KCNQ3, KCNQ5, KCNS1, TRPA1, TRPM8, TRPV1, TRPV3 and TRPV4 variants by single-molecule molecular inversion probes–next-generation sequencing. These patients did not have genetic variants in SCN3A, SCN7A-SCN11A and SCN1B-SCN4B. In twenty patients (20/414, 4.8%), a potentially pathogenic heterozygous variant was identified in an ion-channel gene (ICG). Variants were present in seven genes, for two patients (0.5%) in ANO3, one (0.2%) in KCNK18, two (0.5%) in KCNQ3, seven (1.7%) in TRPA1, three (0.7%) in TRPM8, three (0.7%) in TRPV1 and two (0.5%) in TRPV3. Variants in the TRP genes were the most frequent (n = 15, 3.6%), partly in patients with high mean maximal pain scores VAS = 9.65 ± 0.7 (n = 4). Patients with ICG variants reported more severe pain compared to patients without such variants (VAS = 9.36 ± 0.72 vs. VAS = 7.47 ± 2.37). This cohort study identified ICG variants in neuropathic pain in SFN, complementing previous findings of ICG variants in diabetic neuropathy. These data show that ICG variants are central in neuropathic pain of different etiologies and provides promising gene candidates for future research.     

Fanconi Anemia Patients from an Indigenous Community in Mexico Carry a New Founder Pathogenic Variant in FANCG

Fanconi anemia (FA) is a rare genetic disorder caused by pathogenic variants (PV) in at least 22 genes, which cooperate in the Fanconi anemia/Breast Cancer (FA/BRCA) pathway to maintain genome stability. PV in FANCA, FANCC, and FANCG account for most cases (~90%). This study evaluated the chromosomal, molecular, and physical phenotypic findings of a novel founder FANCG PV, identified in three patients with FA from the Mixe community of Oaxaca, Mexico. All patients presented chromosomal instability and a homozygous PV, FANCG: c.511-3_511-2delCA, identified by next-generation sequencing analysis. Bioinformatic predictions suggest that this deletion disrupts a splice acceptor site promoting the exon 5 skipping. Analysis of Cytoscan 750 K arrays for haplotyping and global ancestry supported the Mexican origin and founder effect of the variant, reaffirming the high frequency of founder PV in FANCG. The degree of bone marrow failure and physical findings (described through the acronyms VACTERL-H and PHENOS) were used to depict the phenotype of the patients. Despite having a similar frequency of chromosomal aberrations and genetic constitution, the phenotype showed a wide spectrum of severity. The identification of a founder PV could help for a systematic and accurate genetic screening of patients with FA suspicion in this population.     

SCN5A Variants as Genetic Arrhythmias Triggers for Familial Bileaflet Mitral Valve Prolapse

Mitral valve prolapse (MVP) is a common valvular heart defect with variable outcomes. Several studies reported MVP as an underestimated cause of life-threatening arrhythmias and sudden cardiac death (SCD), mostly in young adult women. Herein, we report a clinical and genetic investigation of a family with bileaflet MVP and a history of syncopes and resuscitated sudden cardiac death. Using family based whole exome sequencing, we identified two missense variants in the SCN5A gene. A rare variant SCN5A:p.Ala572Asp and the well-known functional SCN5A:p.His558Arg polymorphism. Both variants are shared between the mother and her daughter with a history of resuscitated SCD and syncopes, respectively. The second daughter with prodromal MVP as well as her healthy father and sister carried only the SCN5A:p.His558Arg polymorphism. Our study is highly suggestive of the contribution of SCN5A mutations as the potential genetic cause of the electric instability leading to ventricular arrhythmias in familial MVP cases with syncope and/or SCD history.     

Classification of MSH6 Variants of Uncertain Significance Using Functional Assays

Lynch syndrome (LS) is one of the most common hereditary cancer predisposition syndromes worldwide. Individuals with LS have a high risk of developing colorectal or endometrial cancer, as well as several other cancers. LS is caused by autosomal dominant pathogenic variants in one of the DNA mismatch repair (MMR) genes MLH1, MSH2, PMS2 or MSH6, and typically include truncating variants, such as frameshift, nonsense or splicing variants. However, a significant number of missense, intronic, or silent variants, or small in-frame insertions/deletions, are detected during genetic screening of the MMR genes. The clinical effects of these variants are often more difficult to predict, and a large fraction of these variants are classified as variants of uncertain significance (VUS). It is pivotal for the clinical management of LS patients to have a clear genetic diagnosis, since patients benefit widely from screening, preventive and personal therapeutic measures. Moreover, in families where a pathogenic variant is identified, testing can be offered to family members, where non-carriers can be spared frequent surveillance, while carriers can be included in cancer surveillance programs. It is therefore important to reclassify VUSs, and, in this regard, functional assays can provide insight into the effect of a variant on the protein or mRNA level. Here, we briefly describe the disorders that are related to MMR deficiency, as well as the structure and function of MSH6. Moreover, we review the functional assays that are used to examine VUS identified in MSH6 and discuss the results obtained in relation to the ACMG/AMP PS3/BS3 criterion. We also provide a compiled list of the MSH6 variants examined by these assays. Finally, we provide a future perspective on high-throughput functional analyses with specific emphasis on the MMR genes.     

Molecular Genetics of GLUT1DS Italian Pediatric Cohort: 10 Novel Disease-Related Variants and Structural Analysis

GLUT1 deficiency syndrome (GLUT1DS1; OMIM #606777) is a rare genetic metabolic disease, characterized by infantile-onset epileptic encephalopathy, global developmental delay, progressive microcephaly, and movement disorders (e.g., spasticity and dystonia). It is caused by heterozygous mutations in the SLC2A1 gene, which encodes the GLUT1 protein, a glucose transporter across the blood-brain barrier (BBB). Most commonly, these variants arise de novo resulting in sporadic cases, although several familial cases with AD inheritance pattern have been described. Twenty-seven Italian pediatric patients, clinically suspect of GLUT1DS from both sporadic and familial cases, have been enrolled. We detected by trios sequencing analysis 25 different variants causing GLUT1DS. Of these, 40% of the identified variants (10 out of 25) had never been reported before, including missense, frameshift, and splice site variants. Their structural mapping on the X-ray structure of GLUT1 strongly suggested the potential pathogenic effects of these novel disease-related mutations, broadening the genotypic spectrum heterogeneity found in the SLC2A1 gene. Moreover, 24% is located in a vulnerable region of the GLUT1 protein that involves transmembrane 4 and 5 helices encoded by exon 4, confirming a mutational hotspot in the SLC2A1 gene. Lastly, we investigated possible correlations between mutation type and clinical and biochemical data observed in our GLUT1DS cohort, revealing that splice site and frameshift variants are related to a more severe phenotype and low CSF parameters.