A deep intronic TCTN2 variant activating a cryptic exon predicted by SpliceRover in a patient with Joubert syndrome

Expanding the Phenotypic Spectrum of Pathogenic KIAA0586 Variants: From Joubert Syndrome to Hydrolethalus Syndrome

KIAA0586 variants have been associated with a wide range of ciliopathies, mainly Joubert syndrome (JS, OMIM #616490) and short-rib thoracic dysplasia syndrome (SRTD, OMIM #616546). However, the hypothesis that this gene is involved with hydrolethalus syndrome (HSL, OMIM #614120) and orofaciodigital syndrome IV (OMIM #258860) has already been raised. Ciliopathies' clinical features are often overlapped despite differing in phenotype severity. Besides KIAA0586, HYLS1 and KIF7 are also known for being causative of ciliopathies, indicating that all three genes may have similar or converging genomic pathways. Overall, the genotypic and phenotypic spectrum of ciliopathies becomes wider and conflicting while more and more new variants are added to this group of disorders' molecular pot. In this case report we discuss the first Brazilian individual clinically diagnosed with hydrolethalus syndrome and molecular findings that demonstrate the role of KIAA0586 as a causative gene of a group of genetic disorders. Also, recent reports on individuals with intronic and exonic variants combined leading to ciliopathies support our patient's molecular diagnosis. At the same time, we discuss variable expressivity and overlapping features in ciliopathies.

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.

RNA sequencing and target long-read sequencing reveal an intronic transposon insertion causing aberrant splicing

More than half of cases with suspected genetic disorders remain unsolved by genetic analysis using short-read sequencing such as exome sequencing (ES) and genome sequencing (GS). RNA sequencing (RNA-seq) and long-read sequencing (LRS) are useful for interpretation of candidate variants and detection of structural variants containing repeat sequences, respectively. Recently, adaptive sampling on nanopore sequencers enables target LRS more easily. Here, we present a Japanese girl with premature chromatid separation (PCS)/mosaic variegated aneuploidy (MVA) syndrome. ES detected a known pathogenic maternal heterozygous variant (c.1402-5A>G) in intron 10 of BUB1B (NM_001211.6), a known responsive gene for PCS/MVA syndrome with autosomal recessive inheritance. Minigene splicing assay revealed that almost all transcripts from the c.1402-5G allele have mis-splicing with 4-bp insertion. GS could not detect another pathogenic variant, while RNA-seq revealed abnormal reads in intron 2. To extensively explore variants in intron 2, we performed adaptive sampling and identified a paternal 3.0 kb insertion. Consensus sequence of 16 reads spanning the insertion showed that the insertion consists of Alu and SVA elements. Realignment of RNA-seq reads to the new reference sequence containing the insertion revealed that 16 reads have 5' splice site within the insertion and 3' splice site at exon 3, demonstrating causal relationship between the insertion and aberrant splicing. In addition, immunoblotting showed severely diminished BUB1B protein level in patient derived cells. These data suggest that detection of transcriptomic abnormalities by RNA-seq can be a clue for identifying pathogenic variants, and determination of insert sequences is one of merits of LRS.