Deep intronic GPR143 mutation in a Japanese family with ocular albinism

A progeroid syndrome caused by a deep intronic variant in TAPT1 is revealed by RNA/SI-NET sequencing

Exome sequencing has introduced a paradigm shift for the identification of germline variations responsible for Mendelian diseases. However, non-coding regions, which make up 98% of the genome, cannot be captured. The lack of functional annotation for intronic and intergenic variants makes RNA-seq a powerful companion diagnostic. Here, we illustrate this point by identifying six patients with a recessive Osteogenesis Imperfecta (OI) and neonatal progeria syndrome. By integrating homozygosity mapping and RNA-seq, we delineated a deep intronic TAPT1 mutation (c.1237-52 G>A) that segregated with the disease. Using SI-NET-seq, we document that TAPT1's nascent transcription was not affected in patients' fibroblasts, indicating instead that this variant leads to an alteration of pre-mRNA processing. Predicted to serve as an alternative splicing branchpoint, this mutation enhances TAPT1 exon 12 skipping, creating a protein-null allele. Additionally, our study reveals dysregulation of pathways involved in collagen and extracellular matrix biology in disease-relevant cells. Overall, our work highlights the power of transcriptomic approaches in deciphering the repercussions of non-coding variants, as well as in illuminating the molecular mechanisms of human diseases.

Identifying gene variants underlying the pathogenesis of diabetic retinopathy based on integrated genomic and transcriptomic analysis of clinical extreme phenotypes

Diabetic retinopathy (DR) is a common complication and the leading cause of blindness in patients with type 2 diabetes. DR has been shown to be closely correlated with blood glucose levels and the duration of diabetes. However, the onset and progression of DR also display clinical heterogeneity. We applied whole-exome sequencing and RNA-seq approaches to study the gene mutation and transcription profiles in three groups of diabetic patients with extreme clinical phenotypes in DR onset, timing, and disease progression, aiming to identify genetic variants that may play roles in the pathogenesis of DR. We identified 23 putatively pathogenic genes, and ingenuity pathway analysis of these mutated genes reveals their functional association with glucose metabolism, diabetic complications, neural system activity, and dysregulated immune responses. In addition, ten potentially protective genes were also proposed. These findings shed light on the mechanisms underlying the pathogenesis of DR and may provide potential targets for developing new strategies to combat DR.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Siblings with MAN1B1-CDG Showing Novel Biochemical Profiles

Congenital disorders of glycosylation (CDG), inherited metabolic diseases caused by defects in glycosylation, are characterized by a high frequency of intellectual disability (ID) and various clinical manifestations. Two siblings with ID, dysmorphic features, and epilepsy were examined using mass spectrometry of serum transferrin, which revealed a CDG type 2 pattern. Whole-exome sequencing showed that both patients were homozygous for a novel pathogenic variant of MAN1B1 (NM_016219.4:c.1837del) inherited from their healthy parents. We conducted a HPLC analysis of sialylated N-linked glycans released from total plasma proteins and characterized the α1,2-mannosidase I activity of the lymphocyte microsome fraction. The accumulation of monosialoglycans was observed in MAN1B1-deficient patients, indicating N-glycan-processing defects. The enzymatic activity of MAN1B1 was compromised in patient-derived lymphocytes. The present patients exhibited unique manifestations including early-onset epileptic encephalopathy and cerebral infarction. They also showed coagulation abnormalities and hypertransaminasemia. Neither sibling had truncal obesity, which is one of the characteristic features of MAN1B1-CDG.

In silico identification of pseudo-exon activation events in personal genome and transcriptome data

Causative mutations for human genetic disorders have mainly been identified in exonic regions that code for amino acid sequences. Recently, however, it has been reported that mutations in deep intronic regions can also cause certain human genetic disorders by creating novel splice sites, leading to pseudo-exon activation. To investigate how frequently pseudo-exon activation events occur in normal individuals, we conducted in silico identification of such events using personal genome data and corresponding high-quality transcriptome data. With rather stringent conditions, on average, 2.6 pseudo-exon activation events per individual were identified. More pseudo-exon activation events were found in 5' donor splice sites than in 3' acceptor splice sites. Although pseudo-exon activation events have sporadically been reported as causative mutations in genetic disorders, it is revealed in this study that such events can be observed in normal individuals at a certain frequency. We estimate that human genomes typically contain on average at least 10 pseudo-exon activation events. The actual number should be higher than this, because we used stringent criteria to identify pseudo-exon activation events. This suggests that it is worth considering the possibility of pseudo-exon activation when searching for causative mutations of genetic disorders if candidate mutations are not identified in coding regions or RNA splice sites.

Transposon clusters as substrates for aberrant splice-site activation

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5ʹ splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5ʹ splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

Autosomal recessive complete STAT1 deficiency caused by compound heterozygous intronic mutations

Autosomal recessive (AR) complete signal transducer and activator of transcription 1 (STAT1) deficiency is an extremely rare primary immunodeficiency that causes life-threatening mycobacterial and viral infections. Only seven patients from five unrelated families with this disorder have been so far reported. All causal STAT1 mutations reported are exonic and homozygous. We studied a patient with susceptibility to mycobacteria and virus infections, resulting in identification of AR complete STAT1 deficiency due to compound heterozygous mutations, both located in introns: c.128+2 T>G and c.542-8 A>G. Both mutations were the first intronic STAT1 mutations to cause AR complete STAT1 deficiency. Targeted RNA-seq documented the impairment of STAT1 mRNA expression and contributed to the identification of the intronic mutations. The patient's cells showed a lack of STAT1 expression and phosphorylation, and severe impairment of the cellular response to IFN-γ and IFN-α. The case reflects the importance of accurate clinical diagnosis and precise evaluation, to include intronic mutations, in the comprehensive genomic study when the patient lacks molecular pathogenesis. In conclusion, AR complete STAT1 deficiency can be caused by compound heterozygous and intronic mutations. Targeted RNA-seq-based systemic gene expression assay may help to increase diagnostic yield in inconclusive cases after comprehensive genomic study.

The identification of a RNA splice variant in TULP1 in two siblings with early-onset photoreceptor dystrophy

Background

Early‐onset photoreceptor dystrophies are a major cause of irreversible visual impairment in children and young adults. This clinically heterogeneous group of disorders can be caused by mutations in many genes. Nevertheless, to date, 30%–40% of cases remain genetically unexplained. In view of expanding therapeutic options, it is essential to obtain a molecular diagnosis in these patients as well. In this study, we aimed to identify the genetic cause in two siblings with genetically unexplained retinal disease.

Methods

Whole exome sequencing was performed to identify the causative variants in two siblings in whom a single pathogenic variant in TULP1 was found previously. Patients were clinically evaluated, including assessment of the medical history, slit‐lamp biomicroscopy, and ophthalmoscopy. In addition, a functional analysis of the putative splice variant in TULP1 was performed using a midigene assay.

Results

Clinical assessment showed a typical early‐onset photoreceptor dystrophy in both the patients. Whole exome sequencing identified two pathogenic variants in TULP1, a c.1445G>A (p.(Arg482Gln)) missense mutation and an intronic c.718+23G>A variant. Segregation analysis confirmed that both siblings were compound heterozygous for the TULP1 c.718+23G>A and c.1445G>A variants, while the unaffected parents were heterozygous. The midigene assay for the c.718+23G>A variant confirmed an elongation of exon 7 leading to a frameshift.

Conclusion

Here, we report the first near‐exon RNA splice variant that is not present in a consensus splice site sequence in TULP1, which was found in a compound heterozygous manner with a previously described pathogenic TULP1 variant in two patients with an early‐onset photoreceptor dystrophy. We provide proof of pathogenicity for this splice variant by performing an in vitro midigene splice assay, and highlight the importance of analysis of noncoding regions beyond the noncanonical splice sites in patients with inherited retinal diseases.

Understanding human DNA variants affecting pre-mRNA splicing in the NGS era

Pre-mRNA splicing, an essential step in eukaryotic gene expression, relies on recognition of short sequences on the primary transcript intron ends and takes place along transcription by RNA polymerase II. Exonic and intronic auxiliary elements may modify the strength of exon definition and intron recognition. Splicing DNA variants (SV) have been associated with human genetic diseases at canonical intron sites, as well as exonic substitutions putatively classified as nonsense, missense or synonymous variants. Their effects on mRNA may be modulated by cryptic splice sites associated to the SV allele, comprehending exon skipping or shortening, and partial or complete intron retention. As splicing mRNA outputs result from combinatorial effects of both intrinsic and extrinsic factors, in vitro functional assays supported by computational analyses are recommended to assist SV pathogenicity assessment for human Mendelian inheritance diseases. The increasing use of next-generating sequencing (NGS) targeting full genomic gene sequence has raised awareness of the relevance of deep intronic SV in genetic diseases and inclusion of pseudo-exons into mRNA. Finally, we take advantage of recent advances in sequencing and computational technologies to analyze alternative splicing in cancer. We explore the Catalog of Somatic Mutations in Cancer (COSMIC) to describe the proportion of splice-site mutations in cis and trans regulatory elements. Genomic data from large cohorts of different cancer types are increasingly available, in addition to repositories of normal and somatic genetic variations. These are likely to bring new insights to understanding the genetic control of alternative splicing by mapping splicing quantitative trait loci in tumors.

A rare deep intronic mutation of PKHD1 gene, c.8798-459 C > A, causes autosomal recessive polycystic kidney disease by pseudoexon activation

Autosomal recessive polycystic kidney disease (ARPKD), is a rare hepatorenal fibrocystic disorder primarily associated with progressive growth of multiple cysts in the kidneys causing progressive loss of renal function. The disease is linked to mutations in the PKHD1 gene. In this study, we describe the gene diagnosis and prenatal diagnosis for a consanguineous family with two fetuses diagnosed with polycystic kidney disease by fetal sonography during the pregnancy. Sequence analysis of cDNA synthesized from the PKHD1 mRNA of the second induced fetus identified a 111-nucleotide insert at the junction of exon 56 and 57 that originated from intervening sequence (IVS) 56. Further genomic sequencing of IVS 56 of the PKHD1 gene identified a rare homozygous deep intronic mutation (c.8798-459 C > A), which was inherited from the parents and not detectable in 100 unrelated control subjects. Moreover, we explored the pathogenicity of this deep intronic mutation by conducting a minigene splicing assay experiment, which demonstrated that the mutation causes a pseudoexon insertion, which results in a frameshift followed by a premature termination codon in exon 57. Eventually, the parents had a healthy baby by undergoing prenatal genetic diagnosis based on the targeted detection of the intron mutation. The newly identified deep intronic mutation is associated with a rare mechanism of abnormal splicing that expands the spectrum of known PKHD1 gene mutations. It can be used in evidence-based genetic and reproductive counseling for families with ARPKD.

RNA sequencing solved the most common but unrecognized NEB pathogenic variant in Japanese nemaline myopathy

PURPOSE

The diagnostic rate for Mendelian diseases by exome sequencing (ES) is typically 20-40%. The low rate is partly because ES misses deep-intronic or synonymous variants leading to aberrant splicing. In this study, we aimed to apply RNA sequencing (RNA-seq) to efficiently detect the aberrant splicings and their related variants.

METHODS

Aberrant splicing in biopsied muscles from six nemaline myopathy (NM) cases unresolved by ES were analyzed with RNA-seq. Variants related to detected aberrant splicing events were analyzed with Sanger sequencing. Detected variants were screened in NM patients unresolved by ES.

RESULTS

We identified a novel deep-intronic NEB pathogenic variant, c.1569+339A>G in one case, and another novel synonymous NEB pathogenic variant, c.24684G>C (p.Ser8228Ser) in three cases. The c.24684G>C variant was observed to be the most frequent among all NEB pathogenic variants in normal Japanese populations with a frequency of 1 in 178 (20 alleles in 3552 individuals), but was previously unrecognized. Expanded screening of the variant identified it in a further four previously unsolved nemaline myopathy cases.

CONCLUSION

These results indicated that RNA-seq may be able to solve a large proportion of previously undiagnosed muscle diseases.

Identification and Rescue of Splice Defects Caused by Two Neighboring Deep-Intronic ABCA4 Mutations Underlying Stargardt Disease

Sequence analysis of the coding regions and splice site sequences in inherited retinal diseases is not able to uncover ∼40% of the causal variants. Whole-genome sequencing can identify most of the non-coding variants, but their interpretation is still very challenging, in particular when the relevant gene is expressed in a tissue-specific manner. Deep-intronic variants in ABCA4 have been associated with autosomal-recessive Stargardt disease (STGD1), but the exact pathogenic mechanism is unknown. By generating photoreceptor precursor cells (PPCs) from fibroblasts obtained from individuals with STGD1, we demonstrated that two neighboring deep-intronic ABCA4 variants (c.4539+2001G>A and c.4539+2028C>T) result in a retina-specific 345-nt pseudoexon insertion (predicted protein change: p.Arg1514Leufs^∗36), likely due to the creation of exonic enhancers. Administration of antisense oligonucleotides (AONs) targeting the 345-nt pseudoexon can significantly rescue the splicing defect observed in PPCs of two individuals with these mutations. Intriguingly, an AON that is complementary to c.4539+2001G>A rescued the splicing defect only in PPCs derived from an individual with STGD1 with this but not the other mutation, demonstrating the high specificity of AONs. In addition, a single AON molecule rescued splicing defects associated with different neighboring mutations, thereby providing new strategies for the treatment of persons with STGD1. As many genes associated with human genetic conditions are expressed in specific tissues and pre-mRNA splicing may also rely on organ-specific factors, our approach to investigate and treat splicing variants using differentiated cells derived from individuals with STGD1 can be applied to any tissue of interest.

Accuracy of Next-Generation Sequencing for Molecular Diagnosis in Patients With Infantile Nystagmus Syndrome

Importance

Infantile nystagmus syndrome (INS) is a group of disorders presenting with genetic and clinical heterogeneities that have challenged the genetic and clinical diagnoses of INS. Precise molecular diagnosis in early infancy may result in more accurate genetic counseling and improved patient management.

Objective

To assess the accuracy of genomic data from next-generation sequencing (NGS) and phenotypic data to enhance the definitive diagnosis of INS.

Design, Setting, and Participants

A single-center retrospective case series was conducted in 48 unrelated, consecutive patients with INS, with or without associated ocular or systemic conditions, who underwent genetic testing between June 1, 2015, and January 31, 2017. Next-generation sequencing analysis was performed using a target panel that included 113 genes associated with INS (n = 47) or a TruSight One sequencing panel that included 4813 genes associated with known human phenotypes (n = 1). Variants were filtered and prioritized by in-depth clinical review, and finally classified according to the American College of Medical Genetics and Genomics guidelines. Patients underwent a detailed ophthalmic examination, including electroretinography and optical coherence tomography, if feasible.

Main Outcomes and Measures

Diagnostic yield of targeted NGS testing.

Results

Among the 48 patients (21 female and 27 male; mean [SD] age at genetic testing, 9.2 [10.3] years), 8 had a family history of nystagmus and 40 were simplex. All patients were of a single ethnicity (Korean). Genetic variants that were highly likely to be causative were identified in 28 of the 48 patients, corresponding to a molecular diagnostic yield of 58.3% (95% CI, 44.4%-72.2%). FRMD7, GPR143, and PAX6 mutations appeared to be the major genetic causes of familial INS. A total of 10 patients (21%) were reclassified to a different diagnosis based on results of NGS testing, enabling accurate clinical management.

Conclusions and Relevance

These findings suggest that NGS is an accurate diagnostic tool to differentiate causes of INS because diagnostic tests, such as electroretinography and optical coherence tomography, are not easily applicable in young infants. Accurate application of NGS using a standardized, stepwise, team-based approach in early childhood not only facilitated early molecular diagnosis but also led to improved personalized management in patients with INS.

Novel human mutation and CRISPR/Cas genome-edited mice reveal the importance of C-terminal domain of MSX1 in tooth and palate development

Several mutations, located mainly in the MSX1 homeodomain, have been identified in non-syndromic tooth agenesis predominantly affecting premolars and third molars. We identified a novel frameshift mutation of the highly conserved C-terminal domain of MSX1, known as Msx homology domain 6 (MH6), in a Japanese family with non-syndromic tooth agenesis. To investigate the importance of MH6 in tooth development, Msx1 was targeted in mice with CRISPR/Cas system. Although heterozygous MH6 disruption did not alter craniofacial development, homozygous mice exhibited agenesis of lower incisors with or without cleft palate at E16.5. In addition, agenesis of the upper third molars and the lower second and third molars were observed in 4-week-old mutant mice. Although the upper second molars were present, they were abnormally small. These results suggest that the C-terminal domain of MSX1 is important for tooth and palate development, and demonstrate that that CRISPR/Cas system can be used as a tool to assess causality of human disorders in vivo and to study the importance of conserved domains in genes.

A novel missense mutation of COL5A2 in a patient with Ehlers-Danlos syndrome

Ehlers–Danlos syndrome (EDS) is a group of inherited connective tissue disorders characterized by hyperextensible skin, joint hypermobility and soft tissue fragility. For molecular diagnosis, targeted exome sequencing was performed on a 9-year-old male patient who was clinically suspected to have EDS. The patient presented with progressive kyphoscoliosis, joint hypermobility and hyperextensible skin without scars. Ultimately, classical EDS was diagnosed by identifying a novel, mono-allelic mutation in COL5A2 [NM_000393.3(COL5A2_v001):c.682G>A, p.Gly228Arg].

A novel COL11A1 missense mutation in siblings with non-ocular Stickler syndrome

Stickler syndrome (STL) is an autosomal, dominantly inherited, clinically variable and genetically heterogeneous connective tissue disorder characterized by ocular, auditory, orofacial and skeletal abnormalities. We conducted targeted resequencing using a next-generation sequencer for molecular diagnosis of a 2-year-old girl who was clinically suspected of having STL with Pierre Robin sequence. We detected a novel heterozygous missense mutation, NM_001854.3:n.4838G>A [NM_001854.3 (COL11A1_v001):c.4520G>A], in COL11A1, resulting in a Gly to Asp substitution at position 1507 [NM_001854.3(COL11A1_i001)] within one of the collagen-like domains of the triple helical region. The same mutation was detected in her 4-year-old brother with cleft palate and high-frequency sensorineural hearing loss.

A novel frameshift mutation of CHD7 in a Japanese patient with CHARGE syndrome

CHARGE syndrome is a rare autosomal dominant developmental disorder involving multiple organs. CHD7 is a major causative gene of CHARGE syndrome. We performed targeted-exome sequencing using a next-generation sequencer for molecular diagnosis of a 4-month-old male patient who was clinically suspected to have CHARGE syndrome, and report a novel monoallelic mutation in CHD7, NM_017780.3(CHD7_v001):c.2966del causing a reading frameshift [p.(Cys989Serfs*3)].

Mutation Detection in Patients with Retinal Dystrophies Using Targeted Next Generation Sequencing

Retinal dystrophies (RD) constitute a group of blinding diseases that are characterized by clinical variability and pronounced genetic heterogeneity. The different nonsyndromic and syndromic forms of RD can be attributed to mutations in more than 200 genes. Consequently, next generation sequencing (NGS) technologies are among the most promising approaches to identify mutations in RD. We screened a large cohort of patients comprising 89 independent cases and families with various subforms of RD applying different NGS platforms. While mutation screening in 50 cases was performed using a RD gene capture panel, 47 cases were analyzed using whole exome sequencing. One family was analyzed using whole genome sequencing. A detection rate of 61% was achieved including mutations in 34 known and two novel RD genes. A total of 69 distinct mutations were identified, including 39 novel mutations. Notably, genetic findings in several families were not consistent with the initial clinical diagnosis. Clinical reassessment resulted in refinement of the clinical diagnosis in some of these families and confirmed the broad clinical spectrum associated with mutations in RD genes.

A novel COL11A1 mutation affecting splicing in a patient with Stickler syndrome

Stickler syndrome is a clinically and genetically heterogeneous collagenopathy characterized by ocular, auditory, skeletal and orofacial abnormalities, commonly occurring as an autosomal dominant trait. We conducted target resequencing to analyze candidate genes associated with known clinical phenotypes from a 4-year-old girl with Stickler syndrome. We detected a novel heterozygous intronic mutation (NM_001854.3:c.3168+5G>A) in COL11A1 that may impair splicing, which was suggested by in silico prediction and a minigene assay.