Germline DNA testing using the next-gene-ration sequencing (NGS) technology has become the analytical standard for the diagnostics of hereditary diseases, including cancer. Its increasing use places high demands on correct sample identification, independent confirmation of prioritized variants, and their functional and clinical interpretation. To streamline these processes, we introduced parallel DNA and RNA capture-based NGS using identical capture panel CZECANCA, which is routinely used for DNA analysis of hereditary cancer predisposition. Here, we present the analytical workflow for RNA sample processing and its analytical and diagnostic performance. Parallel DNA/RNA analysis allowed credible sample identification by calculating the kinship coefficient. The RNA capture-based approach enriched transcriptional targets for the majority of clinically relevant cancer predisposition genes to a degree that allowed analysis of the effect of identified DNA variants on mRNA processing. By comparing the panel and whole-exome RNA enrichment, we demonstrated that the tissue-specific gene expression pattern is independent of the capture panel. Moreover, technical replicates confirmed high reproducibility of the tested RNA analysis. We concluded that parallel DNA/RNA NGS using the identical gene panel is a robust and cost-effective diagnostic strategy. In our setting, it allows routine analysis of 48 DNA/RNA pairs using NextSeq 500/550 Mid Output Kit v2.5 (150 cycles) in a single run with sufficient coverage to analyse 226 cancer predisposition and candidate ge-nes. This approach can replace laborious Sanger confirmatory sequencing, increase testing turnaround, reduce analysis costs, and improve interpretation of the impact of variants by analysing their effect on mRNA processing.

Most pathogenic DMD variants are detectable and interpretable by standard genetic testing for dystrophinopthies. However, approximately 1∼3% of dystrophinopthies patients still do not have a detectable DMD variant after standard genetic testing, most likely due to structural chromosome rearrangements and/or deep intronic pseudoexon-activating variants. Here, we report on a boy with a suspected diagnosis of Becker muscular dystrophy (BMD) who remained without a detectable DMD variant after exonic DNA-based standard genetic testing. Dystrophin mRNA studies and genomic Sanger sequencing were performed in the boy, followed by in silico splicing analyses. We successfully detected a novel deep intronic disease-causing variant in the DMD gene (c.2380 + 3317A > T), which consequently resulting in a new dystrophin pseudoexon activation through the enhancement of a cryptic donor splice site. The patient was therefore genetically diagnosed with BMD. Our case report further emphasizes the significant role of disease-causing splicing variants within deep intronic regions in genetically undiagnosed dystrophinopathies.

Genetic variants in ATP7A are associated with a spectrum of X-linked disorders. In descending order of severity, these are Menkes disease, occipital horn syndrome, and X-linked distal spinal muscular atrophy. After 30 years of diagnostic investigation, we identified a deep intronic ATP7A variant in four males from a family affected to variable degrees by a predominantly skeletal phenotype, featuring bowing of long bones, elbow joints with restricted mobility which dislocate frequently, coarse curly hair, chronic diarrhoea, and motor coordination difficulties. Analysis of whole genome sequencing data from the Genomics England 100,000 Genomes Project following clinical re-evaluation identified a deep intronic ATP7A variant, which was predicted by SpliceAI to have a modest splicing effect. Using a mini-gene splicing assay, we determined that the intronic variant results in aberrant splicing. Sanger sequencing of patient cDNA revealed ATP7A transcripts with exon 5 skipping, or inclusion of a novel intron 4 pseudoexon. In both instances, frameshift leading to premature termination are predicted. Quantification of ATP7A mRNA transcripts using a qPCR assay indicated that the majority of transcripts (86.1 %) have non-canonical splicing, with 68.0 % featuring exon 5 skipping, and 18.1 % featuring the novel pseudoexon. We suggest that the variability of the phenotypes within the affected males results from the stochastic effects of splicing. This deep intronic variant, resulting in aberrant ATP7A splicing, expands the understanding of intronic variation on the ATP7A-related disease spectrum.

Propionic acidemia (PA) is an autosomal recessive metabolic disorder caused by variants in PCCA or PCCB, both sub-units of the propionyl-CoA carboxylase (PCC) enzyme. PCC is required for the catabolism of certain amino acids and odd-chain fatty acids. In its absence, the accumulated toxic metabolites cause metabolic acidosis, neurologic symptoms, multi-organ dysfunction and possible death. The clinical presentation of PA is highly variable, with typical onset in the neonatal or early infantile period. We encountered two families, whose children were diagnosed with PA. Exome sequencing (ES) failed to identify a pathogenic variant, and we proceeded with genome sequencing (GS), demonstrating homozygosity to a deep intronic PCCB variant. RNA analysis established that this variant creates a pseudoexon with a premature stop codon. The parents are variant carriers, though three of them display pseudo-homozygosity due to a common large benign intronic deletion on the second allele. The parental presumed homozygosity merits special attention, as it masked the causative variant at first, which was resolved only by RNA studies. Arriving at a rapid diagnosis, whether biochemical or genetic, can be crucial in directing lifesaving care, concluding the diagnostic odyssey, and allowing the family prenatal testing in subsequent pregnancies. This study demonstrates the power of integrative genetic studies in reaching a diagnosis, utilizing GS and RNA analysis to overcome ES limitations and define pathogenicity. Importantly, it highlights that intronic deletions should be taken into consideration when analyzing genomic data, so that pseudo-homozygosity would not be misinterpreted as true homozygosity, and pathogenic variants will not be mislabeled as benign.

Anophthalmia and microphthalmia (A/M) are among the most severe congenital developmental eye disorders. Despite the advancements in genome screening technologies, more than half of A/M patients do not receive a molecular diagnosis. We included seven consanguineous families affected with A/M from Pakistani cohort and an unknown molecular basis. Single gene testing of FOXE3 was performed, followed by genome sequencing for unsolved probands in order to establish a genetic diagnosis for these families. All seven families were provided with a genetic diagnosis. The identified variants were all homozygous, classified as (likely) pathogenic and present in an A/M-associated gene. Targeted FOXE3 sequencing revealed two previously reported pathogenic FOXE3 variants in four families. In the remaining families, genome sequencing revealed a known pathogenic PXDN variant, a novel 13bp deletion in VSX2, and one novel deep intronic splice variant in PXDN. An in vitro splice assay was performed for the PXDN splice variant which revealed a severe splicing defect. Our study confirmed the utility of genome sequencing as a diagnostic tool for A/M-affected individuals. Furthermore, the identification of a novel deep intronic pathogenic variant in PXDN highlights the role of non-coding variants in A/M-disorders and the value of genome sequencing for the identification of this type of variants.

The rare autosomal dominant brain disorder DLG4-related synaptopathy is caused by de novo variants in DLG4 (encoding PSD-95), the majority of which are predicted to be protein-truncating. In addition to splice site variants, a number of synonymous and missense DLG4 variants are predicted to exert their effect through altered RNA splicing, although the pathogenicity of these variants is uncertain without functional RNA studies. Here, we describe a young boy with a deep intronic DLG4 variant (c.2105+235C>T) identified using whole genome sequencing. By using reverse-transcription PCR on RNA derived from peripheral blood, we demonstrate that DLG4 mRNA expression is detectable in blood and the deep intronic variant gives rise to two alternative DLG4 transcripts, one of which includes a pseudoexon. Both alternative transcripts are out-of-frame and predicted to result in protein-truncation, thereby establishing the genetic diagnosis for the proband. This adds to the evidence concerning the pathogenic potential of deep intronic variants and underlines the importance of functional studies, even in cases where reported tissue-specific gene expression might suggest otherwise.

BACKGROUND: Fabry disease is an X-linked lysosomal storage disorder caused by insufficient α-galactosidase A (GLA) activity resulting from variants in the GLA gene, which leads to glycosphingolipid accumulation and life-threatening, multi-organ complications. Approximately 50 variants have been reported that cause splicing abnormalities in GLA. Most were found within canonical splice sites, which are highly conserved GT and AG splice acceptor and donor dinucleotides, whereas one-third were located outside canonical splice sites, making it difficult to interpret their pathogenicity. In this study, we aimed to investigate the genetic pathogenicity of variants located in non-canonical splice sites within the GLA gene.
METHODS: 13 variants, including four deep intronic variants, were selected from the Human Gene Variant Database Professional. We performed an in vitro splicing assay to identify splicing abnormalities in the variants.
RESULTS: All candidate non-canonical splice site variants in GLA caused aberrant splicing. Additionally, all but one variant was protein-truncating. The four deep intronic variants generated abnormal transcripts, including a cryptic exon, as well as normal transcripts, with the proportion of each differing in a cell-specific manner.
CONCLUSIONS: Validation of splicing effects using an in vitro splicing assay is useful for confirming pathogenicity and determining associations with clinical phenotypes.

Phenylketonuria (PKU) is an autosomal recessive congenital metabolic disorder caused by PAH variants. Previously, approximately 5% of PKU patients remained undiagnosed after Sanger sequencing and multiplex ligation-dependent probe amplification. To date, increasing numbers of pathogenic deep intronic variants have been reported in more than 100 disease-associated genes.
In this study, we performed full-length sequencing of PAH to investigate the deep intronic variants in PAH of PKU patients without definite genetic diagnosis.
We identified five deep intronic variants (c.1199+502A>T, c.1065+241C>A, c.706+368T>C, c.706+531>C, and c.706+608A>C). Of these, the c.1199+502A>T variant was found at high frequency and may be a hotspot PAH variant in Chinese PKU. c.706+531T>C and c.706+608A>C are two novel variants that extend the deep intronic variant spectrum of PAH.
Deep intronic variant pathogenicity analysis can further improve the genetic diagnosis of PKU patients. In silico prediction and minigene analysis are powerful approaches for studying the functions and effects of deep intronic variants. Targeted sequencing after full-length gene amplification is an economical and effective tool for the detection of deep intron variation in genes with small fragments.

EMC1 encodes subunit 1 of the endoplasmic reticulum (ER) membrane protein complex (EMC), a transmembrane domain insertase involved in membrane protein biosynthesis. Variants in EMC1 are described as a cause of global developmental delay, hypotonia, cortical visual impairment, and commonly, cerebral atrophy on MRI scan. We report an individual with severe global developmental delay and progressive cerebellar atrophy in whom exome sequencing identified a heterozygous essential splice-site variant in intron-3 of EMC1 (NM_015047.3:c.287-1G>A). Whole genome sequencing (WGS) identified a deep intronic variant in intron-20 of EMC1 (NM_015047.3:c.2588-771C>G) that was poorly predicted by in silico programs to disrupt pre-mRNA splicing. Reverse Transcription-PCR (RT-PCR) revealed stochastic activation of a pseudo-exon associated with the c.2588-771C>G variant and mis-splicing arising from the c.287-1G>A variant. This case highlights the utility of WGS and RNA studies to identify and assess likely pathogenicity of deep intronic variants and expands the genotypic and phenotypic spectrum of EMC1-related disorders.

The purpose of this study was to detect the missing heritability of patients with KIF11-related retinopathy and to describe their clinical and genetic characteristics. We enrolled 10 individuals from 7 unrelated families harboring a pathogenic monoallelic variant in KIF11. All subjects underwent ophthalmic assessment and extraocular phenotype evaluations, as well as comprehensive molecular genetic analyses using next-generation sequencing. Minigene assays were performed to observe the effects of one novel deep intron variant (DIV) and one novel synonymous variant on pre-mRNA splicing. We detected 6 novel different disease-causing variants of KIF11 in the seven pedigrees. Co-segregation analysis and ultra-deep sequencing results indicated that 5 variants arose de novo in 5 families (71%). Functional validation revealed that the synonymous variant leads to an exon skip, while the DIV causes a pseudoexon (PE) inclusion. The patients presented with high variations in their phenotype, and two families exhibited incomplete penetrance. Ocular manifestations and characteristic facial features were observed in all patients, as well as microcephaly in seven patients, intellectual disability in five patients, and lymphedema in one patient. The key retinal features for KIF11-related retinopathy were retinal folds, tractional retinal detachment, and chorioretinal dysplasia. All seven probands had more severe visual detects than other affected family members. Our findings widen the genetic spectrum of KIF11 variants. DIV explained rare unresolved cases with KIF11-related retinopathy. The patients displayed a variable phenotype expressivity and incomplete penetrance, indicating the importance of genetic analysis for patients with KIF11-related retinopathy.