A Deep Intronic PKHD1 Variant Identified by SpliceAI in a Deceased Neonate With Autosomal Recessive Polycystic Kidney Disease

The etiologies of newborn deaths in neonatal intensive care units usually remain unknown, even after genetic testing. Whole-genome sequencing, combined with artificial intelligence-based methods for predicting the effects of non-coding variants, provide an avenue for resolving these deaths. Using one such method, SpliceAI, we identified a maternally inherited deep intronic PKHD1 splice variant (chr6:52030169T>C), in trans with a pathogenic missense variant (p.Thr36Met), in a newborn who died of autosomal recessive polycystic kidney disease at age 2 days. We validated the deep intronic variant's impact in maternal urine-derived cells expressing PKHD1. Reverse transcription polymerase chain reaction followed by Sanger sequencing showed that the variant causes inclusion of 147bp of the canonical intron between exons 29 and 30 of PKHD1 into the mRNA, including a premature stop codon. Allele-specific expression analysis at a heterozygous site in the mother showed that the mutant allele completely suppresses canonical splicing. In an unrelated healthy control, there was no evidence of transcripts including the novel splice junction. We returned a diagnostic report to the parents, who underwent in vitro embryo selection.

Combined Bioinformatic and Splicing Analysis of Likely Benign Intronic and Synonymous Variants Reveals Evidence for Pathogenicity

Background

Current clinical variant analysis pipelines focus on coding variants and intronic variants within 10-20 bases of an exon-intron boundary that may affect splicing. The impact of newer splicing prediction algorithms combined with in vitro splicing assays on rare variants currently considered Benign/Likely Benign (B/LB) is unknown.

Methods

Exome sequencing data from 576 pediatric cancer patients enrolled in the Texas KidsCanSeq study were filtered for intronic or synonymous variants absent from population databases, predicted to alter splicing via SpliceAI (>0.20), and scored as potentially deleterious by CADD (>10.0). Total cellular RNA was extracted from monocytes and RT-PCR products analyzed. Subsequently, rare synonymous or intronic B/LB variants in a subset of genes submitted to ClinVar were similarly evaluated. Variants predicted to lead to a frameshifted splicing product were functionally assessed using an in vitro splicing reporter assay in HEK-293T cells.

Results

KidsCanSeq exome data analysis revealed a rare, heterozygous, intronic variant (NM_177438.3(DICER1):c.574-26A>G) predicted by SpliceAI to result in gain of a secondary splice acceptor site. The proband had a personal and family history of pleuropulmonary blastoma consistent with DICER1 syndrome but negative clinical sequencing reports. Proband RNA analysis revealed alternative DICER1 transcripts including the SpliceAI-predicted transcript.Similar bioinformatic analysis of synonymous or intronic B/LB variants (n=31,715) in ClinVar from 61 Mendelian disease genes yielded 18 variants, none of which could be scored by MaxEntScan. Eight of these variants were assessed (DICER1 n=4, CDH1 n=2, PALB2 n=2) using in vitro splice reporter assay and demonstrated abnormal splice products (mean 66%; range 6% to 100%). Available phenotypic information from submitting laboratories demonstrated DICER1 phenotypes in 2 families (1 variant) and breast cancer phenotypes for PALB2 in 3 families (2 variants).

Conclusions

Our results demonstrate the power of newer predictive splicing algorithms to highlight rare variants previously considered B/LB in patients with features of hereditary conditions. Incorporation of SpliceAI annotation of existing variant data combined with either direct RNA analysis or in vitro assays has the potential to identify disease-associated variants in patients without a molecular diagnosis.

Engineered U1 snRNAs to modulate alternatively spliced exons

Alternative splicing accounts for a considerable portion of transcriptomic diversity, as most protein-coding genes are spliced into multiple mRNA isoforms. However, errors in splicing patterns can give rise to mis-splicing with pathological consequences, such as the congenital diseases familial dysautonomia, Duchenne muscular dystrophy, and spinal muscular atrophy. Small nuclear RNA (snRNA) components of the U snRNP family have been proposed as a therapeutic modality for the treatment of mis-splicing. U1 snRNAs offer great promise, with prior studies demonstrating in vivo efficacy, suggesting additional preclinical development is merited. Improvements in enabling technologies, including screening methodologies, gene delivery vectors, and relevant considerations from gene editing approaches justify further advancement of U1 snRNA as a therapeutic and research tool. With the goal of providing a user-friendly protocol, we compile and demonstrate a methodological toolkit for sequence-specific targeted perturbation of alternatively spliced pre-mRNA with engineered U1 snRNAs. We observe robust modulation of endogenous pre-mRNA transcripts targeted in two contrasting splicing contexts, SMN2 exon 7 and FAS exon 6, exhibiting the utility and applicability of engineered U1 snRNA to both inclusion and exclusion of targeted exons. We anticipate that these demonstrations will contribute to the usability of U1 snRNA in investigating splicing modulation in eukaryotic cells, increasing accessibility to the broader research community.

Two intronic cis-acting variants in both alleles of the POLR3A gene cause progressive spastic ataxia with hypodontia

POLR3A encodes the largest subunit of the DNA-dependent RNA polymerase III. Pathogenic variants in this gene are associated with dysregulation of tRNA production and other non-coding RNAs. POLR3A-related disorders include variable phenotypes. The genotype-phenotype correlation is still unclear. Phenotypic analysis and exome sequencing were performed in four affected siblings diagnosed clinically with hereditary spastic ataxia, two healthy siblings and their unaffected mother. All four affected siblings (ages 46-55) had similar clinical features of early childhood-onset hypodontia and adolescent-onset progressive spastic ataxia. None had progeria, gonadal dysfunction or dysmorphism. All affected individuals had biallelic POLR3A pathogenic variants composed by two cis-acting intronic splicing-altering variants, c.1909 + 22G > A and c.3337-11 T > C. The two healthy siblings had wild-type alleles. The mother and another unaffected sibling were heterozygous for the allele containing both variants. This is the first report addressing the clinical consequence associated with homozygosity for a unique pathogenic intronic allele in the POLR3A gene. This allele was previously reported in compound heterozygous combinations in patients with Wiedemann-Rautenstrauch syndrome, a severe progeroid POLR3A-associated phenotype. We show that homozygosity for this allele is associated with spastic ataxia with hypodontia, and not with progeroid features. These findings contribute to the characterization of genotype-phenotype correlation in POLR3A-related disorders.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics

PURPOSE

Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies.

METHODS

Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations.

RESULTS

Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases.

CONCLUSION

NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Novel mutations causing biotinidase deficiency in individuals identified by newborn screening in Michigan including an unique intronic mutation that alters mRNA expression of the biotinidase gene

Biotinidase deficiency (BD) is an autosomal recessive disorder resulting in the inability to recycle the vitamin biotin. Individuals with biotinidase deficiency can develop neurological and cutaneous symptoms if they are not treated with biotin. To date, more than 165 mutations in the biotinidase gene (BTD) have been reported. Essentially all the mutations result in enzymatic activities with less than 10% of mean normal serum enzyme activity (profound biotinidase deficiency) with the exception of the c.1330G>C (p.D444H) mutation, which results in an enzyme having 50% of mean normal serum activity and causes partial biotinidase deficiency (10-30% of mean normal serum biotinidase activity) if there is a mutation for profound biotinidase deficiency on the second allele. We now reported eight novel mutations in ten children identified by newborn screening in Michigan from 1988 to the end of 2012. Interestingly, one intronic mutation, c.310-15delT, results in an approximately two-fold down-regulation of BTD mRNA expression by Quantitative real-time reverse-transcription PCR (qRT-PCR). This is the first report of an intronic mutation in the BTD gene with demonstration of its effect on enzymatic activity by altering mRNA expression. This study identified three other mutations likely to cause partial biotinidase deficiency. These results emphasize the importance of full gene sequencing of BTD on patients with biotinidase deficiency to better understand the genotype and phenotype correlation in the future.