SpliceAI-visual: a free online tool to improve SpliceAI splicing variant interpretation

Genomic instability in ovarian cancer: Through the lens of single nucleotide polymorphisms

Ovarian cancer (OC) is the deadliest gynecological malignancy among all female reproductive cancers. It is characterized by high mortality rate and poor prognosis. Genomic instability caused by mutations, single nucleotide polymorphisms (SNPs), copy number variations (CNVs), microsatellite instability (MSI), and chromosomal instability (CIN) are associated with OC predisposition. SNPs, which are highly prevalent in the general population, show a greater relative risk contribution, particularly in sporadic cancers. Understanding OC etiology in terms of genetic basis can increase the use of molecular diagnostics and provide promising approaches for designing novel treatment modalities. This will help deliver personalized medicine to OC patients, which may soon be within reach. Given the pivotal impact of SNPs in cancers, the primary emphasis of this review is to shed light on their prevalence in key caretaker genes that closely monitor genomic integrity, viz., DNA damage response, repair, cell cycle checkpoints, telomerase maintenance, and apoptosis and their clinical implications in OC. We highlight the current challenges faced in different SNP-based studies. Various computational methods and bioinformatic tools employed to predict the functional impact of SNPs have also been comprehensively reviewed concerning OC research. Overall, this review identifies that variants in the DDR and HRR pathways are the most studied, implying their critical role in the disease. Conversely, variants in other pathways, such as NHEJ, MMR, cell cycle, apoptosis, telomere maintenance, and PARP genes, have been explored the least.

Classification and prediction of variants associated with hearing loss using sequence information in the vicinity of mutation sites

Hearing impairment is a major global health problem, affecting more than 5 % of the world's population at various ages, from neonates to the elderly. Among the common genetic variations in humans, single nucleotide variations and small insertions or deletions predominate. The study of hearing loss resulting from these variations is proving invaluable in the analysis and diagnosis of hearing disorders. The identification of pathogenic mutations is frequently a lengthy and laborious process. Existing computational prediction tools have been developed primarily for common diseases and genome-wide analyses, with less focus on deafness. This study proposes a novel approach that focuses on the regions surrounding mutation sites. Mutation sites associated with deafness and their flanking regions of different lengths were extracted from relevant databases and combined into seven distinct segments of different lengths. The information-theoretic features of these segments were computed. Five machine learning algorithms were then used for training, resulting in the construction of a model capable of classifying and predicting deafness-related mutations. For fragments encompassing the 250 bp regions upstream and downstream of the mutations, the average AUC of the five classifiers on the independent test set is 0.89 and the average ACC is 0.85, indicating that the model has a high recognition rate of the pathogenic deafness mutation site. An ensemble approach was also applied to predict variants of uncertain significance (VUS) that may be associated with deafness. These variants were then scored and ranked to assess their likelihood of contributing to the condition.

Congenital Titinopathies Linked to Mutations in TTN Metatranscript-Only Exons

Congenital titinopathies reported to date show autosomal recessive inheritance and are caused by a variety of genomic variants, most of them located in metatranscript (MTT)-only exons. The aim of this study was to describe additional patients and establish robust genotype-phenotype associations in titinopathies. This study involved analyzing molecular, clinical, pathological, and muscle imaging features in 20 patients who had at least one pathogenic or likely pathogenic TTN variant in MTT-only exons, with onset occurring antenatally or in the early postnatal stages. The 20 patients with recessive inheritance exhibited a heterogeneous range of phenotypes. These included fetal lethality, progressive weakness, cardiac or respiratory complications, hyper-CKemia, or dystrophic muscle biopsies. MRI revealed variable abnormalities in different muscles. All patients presented severe congenital myopathy at birth, characterized by arthrogryposis (either multiplex or axial-distal) or neonatal hypotonia in most cases. This study provides detailed genotype-phenotype correlations in congenital titinopathies caused by mutations in MTT-only exons. The findings highlight the variability in clinical presentation and the severity of phenotypes associated with these specific genetic alterations. RNA-seq analyses provided valuable insights into the molecular consequences of TTN variants, particularly in relation to splicing defects and nonsense-mediated RNA decay. In conclusion, this study reinforces the genotype-phenotype correlations between congenital myopathies and variants in TTN MTT-only exons, improves their molecular diagnosis, and provides a better understanding of their pathophysiology.

Novel Biallelic Synonymous Exonic Variant in VPS13A Affecting mRNA Splicing: Case Report

Objectives

Chorea-acanthocytosis is an autosomal recessively inherited condition caused by loss-of-function pathogenic variants in VPS13A. We identified a novel synonymous exonic variant leading to abnormal mRNA splicing in a patient with chorea-acanthocytosis.

Methods

A patient with focal epilepsy developed generalized chorea with orolingual dystonia, cognitive decline, and peripheral neuropathy, consistent with chorea-acanthocytosis. Her parents were first cousins, but there was otherwise no family history. Targeted gene sequencing for variants in VPS13A, mRNA splicing analysis, and Western blot for chorein were performed.

Results

A homozygous synonymous variant in exon 41 of VPS13A (NM_033305.3): c.5157C>T; p.Gly1719 = was identified; this was previously classified as a variant of uncertain significance. SpliceAI predicted a splice donor gain with a score of 0.75 2 base pairs upstream of the reported variant. RNA splicing analysis revealed the creation of a type III splice variant, resulting in a frameshift and a premature termination codon. Western blot showed absent chorein/VPS13A protein.

Discussion

The variant is reclassified as likely pathogenic based on the American College of Medical Genetics criteria. This is the first reported case of ChAc caused by a synonymous variant in VPS13A proven to affect splicing. Our report further expands the spectrum of variants known to cause ChAc.

In silico splicing analysis of the PMS2 gene: exploring alternative molecular mechanisms in PMS2-associated Lynch syndrome

Lynch syndrome (LS) is one of the most common hereditary cancer syndrome in human populations, associated with germline variants in MLH1, MSH2/EPCAM, MSH6 and PMS2 genes. The advent of next generation sequencing has proven a significant impact in germline variant detection in the causative genes; however, a large proportion of patients with clinical criteria still receive uncertain or negative results. PMS2 is the least frequent reported gene, associated with up to 15% of LS cases with late-onset disease and low penetrance phenotype; however, the proportion of PMS2-LS cases is considered to be highly underestimated. In this context, our analysis aimed to improve the current diagnostic yield by focusing on missense and intronic PMS2 variants available in public clinical databases (ClinVar, LOVD). We performed an in silico assessment of the wild-type DNA sequence and the reported genetic variants, employing splicing bioinformatics tools known for their effectiveness in other genes. Splicing variants were predicted in silico and using GTEx short-read RNA expression data. Out of the 2384 missense variants discovered, 90% were classified with uncertain significance (VUS). 4.9% of missense variants were shown to have a potential splicing consequence (DS > 0.2) using SpliceAI. As described in the original publication, SpliceAI-visual was proven effective in annotation of short intronic variants (< 50 bp). Four short intronic variants were identified using SpliceAI-visual as potentially splicing disturbing, in spite of using a lower threshold (DS > 0.1). Exons 2, 3, 4, 5, 6, 7, 8, 11, 12 and 14 were consistently predicted in at least three out of eight software with weak canonical splice sites. Additionally, we noted that both Exonic Splicing Enhancers (ESEs) and Exonic Splicing Silencers (ESSs) contribute significantly to alternative splicing and exonic selection in PMS2 gene. Specifically, ESE motifs were consistently more abundant in highly expressed exons 5, 11 and 14, while ESS motifs played a fundamental role in exons 6, 7 and 10. Computational analysis performed in our study serves as a valuable filtering step for guiding further RNA experiments. Additional functional data is necessary to validate our findings.

GoFCards: an integrated database and analytic platform for gain of function variants in humans

Gain-of-function (GOF) variants, which introduce new or amplify protein functions, are essential for understanding disease mechanisms. Despite advances in genomics and functional research, identifying and analyzing pathogenic GOF variants remains challenging owing to fragmented data and database limitations, underscoring the difficulty in accessing critical genetic information. To address this challenge, we manually reviewed the literature, pinpointing 3089 single-nucleotide variants and 72 insertions and deletions in 579 genes associated with 1299 diseases from 2069 studies, and integrated these with the 3.5 million predicted GOF variants. Our approach is complemented by a proprietary scoring system that prioritizes GOF variants on the basis of the evidence supporting their GOF effects and provides predictive scores for variants that lack existing documentation. We then developed a database named GoFCards for general geneticists and clinicians to easily obtain GOF variants in humans (http://www.genemed.tech/gofcards). This database also contains data from >150 sources and offers comprehensive variant-level and gene-level annotations, with the aim of providing users with convenient access to detailed and relevant genetic information. Furthermore, GoFCards empowers users with limited bioinformatic skills to analyze and annotate genetic data, and prioritize GOF variants. GoFCards offers an efficient platform for interpreting GOF variants and thereby advancing genetic research.

GPATCH11 variants cause mis-splicing and early-onset retinal dystrophy with neurological impairment

Here we conduct a study involving 12 individuals with retinal dystrophy, neurological impairment, and skeletal abnormalities, with special focus on GPATCH11, a lesser-known G-patch domain-containing protein, regulator of RNA metabolism. To elucidate its role, we study fibroblasts from unaffected individuals and patients carrying the recurring c.328+1 G > T mutation, which specifically removes the main part of the G-patch domain while preserving the other domains. Additionally, we generate a mouse model replicating the patients' phenotypic defects, including retinal dystrophy and behavioral abnormalities. Our results reveal a subcellular localization of GPATCH11 characterized by a diffuse presence in the nucleoplasm, as well as centrosomal localization, suggesting potential functions in RNA and cilia metabolism. Transcriptomic analysis performed on mouse retina detect dysregulation in both gene expression and splicing activity, impacting key processes such as photoreceptor light responses, RNA regulation, and primary cilia-associated metabolism. Proteomic analysis of mouse retina confirms the roles GPATCH11 plays in RNA processing, splicing, and transcription regulation, while also suggesting additional functions in synaptic plasticity and nuclear stress response. Our research provides insights into the diverse roles of GPATCH11 and identifies that the mutations affecting this protein are responsible for a recently characterized described syndrome.

Exome Sequencing Has a High Diagnostic Rate in Sporadic Congenital Hypopituitarism and Reveals Novel Candidate Genes

CONTEXT

The pituitary gland is key for childhood growth, puberty, and metabolism. Pituitary dysfunction is associated with a spectrum of phenotypes, from mild to severe. Congenital hypopituitarism (CH) is the most commonly reported pediatric endocrine dysfunction, with an incidence of 1:4000, yet low rates of genetic diagnosis have been reported.

OBJECTIVE

We aimed to unveil the genetic etiology of CH in a large cohort of patients from Argentina.

METHODS

We performed whole exome sequencing of 137 unrelated cases of CH, the largest cohort examined with this method to date.

RESULTS

Of the 137 cases, 19.1% and 16% carried pathogenic or likely pathogenic variants in known and new genes, respectively, while 28.2% carried variants of uncertain significance. This high yield was achieved through the integration of broad gene panels (genes described in animal models and/or other disorders), an unbiased candidate gene screen with a new bioinformatics pipeline (including genes with high loss-of-function intolerance), and analysis of copy number variants. Three novel findings emerged. First, the most prevalent affected gene encodes the cell adhesion factor ROBO1. Affected children had a spectrum of phenotypes, consistent with a role beyond pituitary stalk interruption syndrome. Second, we found that CHD7 mutations also produce a phenotypic spectrum, not always associated with full CHARGE syndrome. Third, we add new evidence of pathogenicity in the genes PIBF1 and TBC1D32, and report 13 novel candidate genes associated with CH (eg, PTPN6, ARID5B).

CONCLUSION

Overall, these results provide an unprecedented insight into the diverse genetic etiology of hypopituitarism.

Reclassification of Two MLH1 Variants of Uncertain Significance Utilizing Clinical and Functional Data

BACKGROUND

Pathogenic variants in the mismatch repair genes are associated with an elevated lifetime risk of colorectal cancer (CRC). We previously identified two variants of uncertain significance (VUS) in the MLH1 gene, c.696_698del, p.(Cys233del) and c.1919C > G, p.(Pro640Arg), in Danish families with numerous occurrences of CRC.

METHODS

To reclassify the variants we collected clinical data, initiated tumor and co-segregation analysis, and performed RNA splicing analysis, subcellular localization, and protein stability studies.

RESULTS

The functional analysis revealed that the c.696_698del, p.(Cys233del) variant had an effect at the RNA level, on subcellular localization, and on protein stability, while the c.1919C > G, p.(Pro640Arg) variant showed decreased expression in localization studies and decreased protein stability. These results suggest both variants disrupt DNA mismatch repair.

CONCLUSION

By applying all collected data and functional results we propose to reclassify the c.696_698del, p.(Cys233del) and the c.1919C > G, p.(Pro640Arg) variants as likely pathogenic (class 4) using MMR gene-specific ACMG/AMP guidelines. Consequently, the two MLH1 variants can now be used for risk assessment of variant carriers, while family members without the variants can be excluded from intensified cancer surveillance and follow population recommendations.

Biallelic variants in ERLIN1: a series of 13 individuals with spastic paraparesis

Biallelic variants in the ERLIN1 gene were recently reported as the cause of two motor neuron degeneration diseases, SPG62 and a recessive form of amyotrophic lateral sclerosis. However, only 12 individuals from five pedigrees have been identified so far. Thus, the description of the disease remains limited. Following the discovery of a homozygous pathogenic variant in a girl with SPG62, presenting with intellectual disability, and epilepsy, we gathered the largest series of SPG62 cases reported so far (13 individuals) to better understand the phenotype associated with ERLIN1. We collected molecular and clinical data for 13 individuals from six families with ERLIN1 biallelic variants. We performed RNA-seq analyses to characterize intronic variants and used Alphafold and a transcripts database to characterize the molecular consequences of the variants. We identified three new variants suspected to alter the bell-shaped ring formed by the ERLIN1/ERLIN2 complex. Affected individuals had childhood-onset paraparesis with slow progression. Six individuals presented with gait ataxia and three had superficial sensory loss. Aside from our proband, none had intellectual disability or epilepsy. Biallelic pathogenic ERLIN1 variants induce a rare, predominantly pure, spastic paraparesis, with possible cerebellar and peripheral nerve involvement.

GWAS and 3D chromatin mapping identifies multicancer risk genes associated with hormone-dependent cancers

Hormone-dependent cancers (HDCs) share several risk factors, suggesting a common aetiology. Using data from genome-wide association studies, we showed spatial clustering of risk variants across four HDCs (breast, endometrial, ovarian and prostate cancers), contrasting with genetically uncorrelated traits. We identified 44 multi-HDC risk regions across the genome, defined as overlapping risk regions for at least two HDCs: two regions contained risk variants for all four HDCs, 13 for three HDCs and 28 for two HDCs. Integrating GWAS data, epigenomic profiling and promoter capture HiC maps from diverse cell line models, we annotated 53 candidate risk genes at 22 multi-HDC risk regions. These targets were enriched for established genes from the COSMIC Cancer Gene Census, but many had no previously reported pleiotropic roles. Additionally, we pinpointed lncRNAs as potential HDC targets and identified risk alleles in several regions that altered transcription factors motifs, suggesting regulatory mechanisms. Known drug targets were over-represented among the candidate multi-HDC risk genes, implying that some may serve as targets for therapeutic development or facilitate the repurposing of existing treatments for HDC. Our approach provides a framework for identifying common target genes driving complex traits and enhances understanding of HDC susceptibility.

Genotypic spectrum of albinism in Mali

Albinism is a phenotypically and genetically heterogeneous condition characterized by a variable degree of hypopigmentation and by ocular features leading to reduced visual acuity. Whereas numerous genotypic studies have been conducted throughout the world, very little is known about the genotypic spectrum of albinism in Africa and especially in sub-Saharan Western Africa. Here we report the analysis of all known albinism genes in a series a 23 patients originating from Mali. Four were diagnosed with OCA 1 (oculocutaneous albinism type 1), 17 with OCA 2, and two with OCA 4. OCA2 variant NM_000275.3:c.819_822delinsGGTC was most frequently encountered. Four novel variants were identified (two in TYR, two in OCA2). A deep intronic variant was found to alter splicing of the OCA2 RNA by inclusion of a pseudo exon. Of note, the OCA2 exon 7 deletion commonly found in eastern, central, and southern Africa was absent from this series. African patients with OCA 1 and OCA 4 had only been reported twice and once, respectively, in previous publications. This study constitutes the first report of the genotypic spectrum of albinism in a western sub-Saharan country.

Genetic Variants at PRKCG Splice and UTR Sites Promote Cancer Susceptibility by Disrupting Epigenetic and miRNA Regulatory Network

The changes in the protein kinase C gamma gene (PRKCG) expression are associated with both coding and non-coding variants. No studies have specifically established the association between PRKCG 3'UTR, 5'UTR, donor and acceptor splice variants with post-transcriptional changes through utilizing in-silico tools. The current study intends to uncover this linkage. In total, 419 3' and 5'UTR variants were retrieved. 325 of these variant IDs were annotated as functionally significant. 18 variants impacted the transcription factors binding and therefore influenced the post-transcriptional regulatory activity while 7 variants affected regulatory mechanisms through histone modifications. 2 rsIDs (rs373228, rs446795) potentially impacted the interactions with RNA binding proteins. In addition to that, PRKCG showed high expression in brain cells and had variable expression in TCGA tumors, respectively. Furthermore, 5 3' UTR variants were identified to be targeted by miRNAs. In total, 5 of these miRNAs (hsa-miR-663a, hsa-miR-324-5p, hsa-miR-646, hsa-miR-1205 and hsa-miR-4270) that targeted 3'UTRs (rs57483118, rs181418157 and rs60891969) showed differential expressions in distinct cancer types. The presence of 3'UTR variants likely altered the secondary structure of mRNA. The 7 rsIDs at 3' UTR site caused the loss of function of authentic splice site at 10 positions was noted; at 1 position, gain of function was observed while at 2 positions no effect was identified. Moreover, the loss of donor and acceptor splice site was evident. Our results highlight the importance of non-coding regions that might boost our research capacity to predict and construct targeted therapeutic approaches.

CHRNE-related congenital myasthenic syndrome in Iran: Clinical and molecular insights

Variants in the CHRNE gene can lead to a condition called congenital myasthenic syndrome (CMS), which affects the neuromuscular junction (NMJ). CHRNE mutations are the most common cause of CMS. Seventy-seven patients with a possible diagnosis of CMS were referred to the neuromuscular clinic of Shariati Hospital affiliated with the Tehran University of Medical Sciences. We performed whole-exome sequencing (WES) to determine the underlying defect in a group of individuals with a possible diagnosis of CMS. Clinical features and morphological and molecular data on 33 patients with mutations in CHRNE were described. Age of onset, age at diagnosis, consanguinity, family history, motor milestone delay, ophthalmoparesis, generalized fatigue, dysphagia, neurophysiologic findings, and response to treatment of the patients were assessed. Nineteen CHRNE variants including 10 novel ones were identified. The most common mutations were c.1327del; (p.Glu443LysfsTer64) in four different families and c.1252-1267dup; (p.Cys423SerfsTer38) in three families. Clinical onset was mostly at birth or under one year with bilateral fatigable ptosis, ophthalmoplegia, bulbar weakness, and proximal muscle weakness. All patients were treated with pyridostigmine ± salbutamol, which resulted in improvement of motor function, dysphagia, and breathing.

Investigation of cryptic JAG1 splice variants as a cause of Alagille syndrome and performance evaluation of splice predictor tools

Haploinsufficiency of JAG1 is the primary cause of Alagille syndrome (ALGS), a rare, multisystem disorder. The identification of JAG1 intronic variants outside of the canonical splice region as well as missense variants, both of which lead to uncertain associations with disease, confuses diagnostics. Strategies to determine whether these variants affect splicing include the study of patient RNA or minigene constructs, which are not always available or can be laborious to design, as well as the utilization of computational splice prediction tools. These tools, including SpliceAI and Pangolin, use algorithms to calculate the probability that a variant results in a splice alteration, expressed as a Δ score, with higher Δ scores (>0.2 on a 0-1 scale) positively correlated with aberrant splicing. We studied the consequence of 10 putative splice variants in ALGS patient samples through RNA analysis and compared this to SpliceAI and Pangolin predictions. We identified eight variants with aberrant splicing, seven of which had not been previously validated. Combining these data with non-canonical and missense splice variants reported in the literature, we identified a predictive threshold for SpliceAI and Pangolin with high sensitivity (Δ score >0.6). Moreover, we showed reduced specificity for variants with low Δ scores (<0.2), highlighting a limitation of these tools that results in the misidentification of true splice variants. These results improve genomic diagnostics for ALGS by confirming splice effects for seven variants and suggest that the integration of splice prediction tools with RNA analysis is important to ensure accurate clinical variant classifications.

Further description of the phenotypic spectrum of neuronal ceroid lipofuscinosis type 11

PURPOSE

Ceroid lipofuscinosis type 11 (CLN11) is a very rare disease, being reported in only 13 unrelated families so far. Further reports are necessary to comprehend the clinical phenotype of this condition. This article aims to report 9 additional cases of CLN11 from 9 unrelated Latin American families presenting with relatively slow disease progression.

METHODS

This was a retrospective observational study including patients with CLN11. Patients were identified through an active search for granulin precursor gene (GRN) pathogenic variants across the entire database of next-generation sequencing of a commercial laboratory and by contacting attending physicians to check for clinical and radiologic findings compatible with a neuronal ceroid lipofuscinosis phenotype.

RESULTS

Nine CLN11 patients from unrelated families were evaluated. Age of onset varied between 3 to 17 years. The most common findings were visual impairment, cerebellar ataxia, seizures, myoclonus, and cognitive decline. One patient had a previously unreported finding of cervical, perioral, and tongue myoclonus. Most of the patients were able to walk unassisted after an average of 14.2 years (SD 4.76 y) from disease onset.

CONCLUSION

We describe 9 new cases of a very rare type of neuronal ceroid lipofuscinosis (CLN11) from Latin America with a recurrent p.(Gln257ProfsTer27) and a novel p.(Cys83Ter) nonsense variant. Our findings suggest that a slowly progressive neuronal ceroid lipofuscinosis might be a clue for the diagnosis of CLN11.

Dissecting the Reduced Penetrance of Putative Loss-of-Function Variants in Population-Scale Biobanks

Loss-of-function variants (LoFs) disrupt the activity of their impacted gene. They are often associated with clinical phenotypes, including autosomal dominant diseases driven by haploinsufficiency. Recent analyses using biobanks have suggested that LoF penetrance for some haploinsufficient disorders may be low, an observation that has important implications for population genomic screening. However, biobanks are also rife with missing data, and the reliability of these findings remains uncertain. Here, we examine the penetrance of putative LoFs (pLoFs) using a cohort of ≈24,000 carriers derived from two population-scale biobanks: the UK Biobank and the All of Us Research Program. We investigate several possible etiologies for reduced pLoF penetrance, including biobank recruitment biases, annotation artifacts, missed diagnoses, and incomplete clinical records. Systematically accounting for these factors increased penetrance, but widespread reduced penetrance remained. Therefore, we hypothesized that other factors must be driving this phenomenon. To test this, we trained machine learning models to identify pLoFs with high penetrance using the genomic features specific to each variant. These models were predictive of penetrance across a range of diseases and ploF types, including those with prior evidence for pathogenicity. This suggests that reduced ploF penetrance is in fact common, and care should be taken when counseling asymptomatic carriers.

Genetic variants associated with cardiac hypertrophy-related sudden cardiac death and cardiovascular outcomes in a Finnish population

BACKGROUND

Hypertrophic cardiomyopathy is a common cause of non-ischaemic sudden cardiac death (SCD). Left ventricular hypertrophy (LVH) without cardiomyopathy-related myocardial disarray is a common autopsy finding and is often associated with prior hypertension in SCD subjects. Our aim was to investigate novel rare gene variants among SCD subjects with presumably hypertension-related LVH and myocardial fibrosis at autopsy.

METHODS

Whole exome sequencing was used to study rare variants (minor allele frequency<0.005) estimated to be deleterious in 96 non-ischaemic SCD subjects with presumably hypertension-related LVH and myocardial fibrosis. Associations of the identified variants with cardiac disease endpoints were replicated in the Finnish national genetic study (FinnGen) dataset.

RESULTS

18 variants were estimated likely to affect protein function and 14 of these were associated with cardiomyopathies, heart failure, conduction abnormalities, hypertension and/or cardiac arrest in Finnish population (FinnGen). Three of the variants were classified as pathogenic or likely pathogenic. These include the splice site variant NM_000449.3:c.234-1G>A in regulatory factor X5 and frameshift variants NM_000449.3:c.234-1G>A in dehydrogenase/reductase 7C and NM_015873.3:c.1164del in villin like.

CONCLUSIONS

We identified rare deleterious variants associated with LVH in SCD subjects. Several of the identified rare variants associated with cardiovascular endpoints including heart failure, cardiomyopathies, cardiac arrest and hypertension in general population.

ParSE-seq: a calibrated multiplexed assay to facilitate the clinical classification of putative splice-altering variants

Interpreting the clinical significance of putative splice-altering variants outside canonical splice sites remains difficult without time-intensive experimental studies. To address this, we introduce Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed assay to quantify variant effects on RNA splicing. We first apply this technique to study hundreds of variants in the arrhythmia-associated gene SCN5A. Variants are studied in 'minigene' plasmids with molecular barcodes to allow pooled variant effect quantification. We perform experiments in two cell types, including disease-relevant induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The assay strongly separates known control variants from ClinVar, enabling quantitative calibration of the ParSE-seq assay. Using these evidence strengths and experimental data, we reclassify 29 of 34 variants with conflicting interpretations and 11 of 42 variants of uncertain significance. In addition to intronic variants, we show that many synonymous and missense variants disrupted RNA splicing. Two splice-altering variants in the assay also disrupt splicing and sodium current when introduced into iPSC-CMs by CRISPR-Cas9 editing. ParSE-seq provides high-throughput experimental data for RNA-splicing to support precision medicine efforts and can be readily adopted to study other loss-of-function genotype-phenotype relationships.

Case report: Two siblings with very late onset of holocarboxylase synthase deficiency and a mini-review

Holocarboxylase synthase (HCS) deficiency is an extremely rare metabolic disorder typically presenting as severe neonatal metabolic acidosis, lethargy, hypotonia, vomiting, and seizures. This report describes two siblings in a family with late-onset forms of HCS deficiency. The younger sister presented at the age of 11 years and manifested as acute metabolic acidosis, which promptly resolved following rehydration and biotin administration. The results of the organic urine profile confirmed multiple carboxylase deficiency, and genetic testing revealed a novel pathogenic variant in the HLCS gene (NM_000411.8) in the homozygous state: c.995A>G; p. (Gln332Arg). No further decompensation was observed for her during the 3-year follow-up period. His older brother was diagnosed at the age of 23 years-old through biochemical tests, without any history of acidotic decompensation. A mini-review of HCS deficiency with late onset (>1 year) or early onset (<1 month) revealed that splice variants are associated with late onset, while both variants p. (Leu216Arg) and p. (Leu237Pro) are associated with early onset. However, the majority of genotypes do not show a clear correlation with the timing of HCS deficiency onset. The most significant point here is the description of extremely late-onset cases of HCS deficiency. This can prompt metabolic investigations and raise suspicion of this rare disease in cases of unexplained metabolic acidosis, even beyond early childhood.

A Novel MAG Variant Causes Hereditary Spastic Paraplegia in a Consanguineous Pakistani Family

Background and objectives: Hereditary spastic paraplegia (HSP) is characterized by unsteady gait, motor incoordination, speech impairment, abnormal eye movement, progressive spasticity and lower limb weakness. Spastic paraplegia 75 (SPG75) results from a mutation in the gene that encodes myelin associated glycoprotein (MAG). Only a limited number of MAG variants associated with SPG75 in families of European, Middle Eastern, North African, Turkish and Palestinian ancestry have been documented so far. This study aims to provide further insight into the clinical and molecular manifestations of HSP. Methods: Using whole-exome sequencing, we investigated a consanguineous Pakistani family where three individuals presented with clinical signs of HSP. Sanger sequencing was used to carry out segregation analysis on available family members, and a minigene splicing assay was utilized to evaluate the effect of the splicing variant. Results: We identified a novel homozygous pathogenic splice donor variant in MAG (c.46 + 1G > T) associated with SPG75. RNA analysis revealed exon skipping that resulted in the loss of a start codon for ENST00000361922.8 isoform. Affected individuals exhibited variable combinations of nystagmus, developmental delay, cognitive impairments, spasticity, dysarthria, delayed gait and ataxia. The proband displayed a quadrupedal stride, and his siblings experienced frequent falls and ataxic gait as one of the prominent features that have not been previously reported in SPG75. Conclusions: Thus, the present study presents an uncommon manifestation of SPG75, the first from the Pakistani population, and broadens the spectrum of MAG variants.

Resolving unsolved whole-genome sequencing data in paediatric neurological disorders: a cohort study

OBJECTIVE

To resolve unsolved whole-genome sequencing (WGS) data in individuals with paediatric neurological disorders.

DESIGN

A cohort study method using updated bioinformatic tools, new analysis targets, clinical information and literature databases was employed to reanalyse existing unsolved genome data.

PARTICIPANTS

From January 2016 to September 2023, a total of 615 individuals who aged under 18 years old, exhibited neurological disorders and received singleton WGS were recruited. 364 cases were unsolved during initial WGS analysis, in which 102 consented to reanalyse existing singleton WGS data.

RESULTS

Median duration for reanalysis after initial negative WGS results was 2 years and 4 months. The diagnostic yield was 29 of 102 individuals (28.4%) through reanalysis. New disease gene discovery and new target acquisitions contributed to 13 of 29 solved cases (44.8%). The reasons of non-detected causative variants during initial WGS analysis were variant reclassification in 9 individuals (31%), analytical issue in 9 (31%), new emerging disease-gene association in 8 (27.6%) and clinical update in 3 (10.3%). The 29 new diagnoses increased the cumulative diagnostic yield of clinical WGS in the entire study cohort to 45.5% after reanalysis.

CONCLUSIONS

Unsolved paediatric WGS individuals with neurological disorders could obtain molecular diagnoses through reanalysis within a timeframe of 2-2.5 years. New disease gene, structural variations and deep intronic splice variants make a significant contribution to diagnostic yield. This approach can provide precise genetic counselling to positive reanalysis results and end a diagnostic odyssey.

Single-Base Substitution Causing Dual-Exon Skipping Event in PKD2 Gene: Unusual Molecular Finding from Exome Sequencing in a Patient with Autosomal Dominant Polycystic Kidney Disease

Background: Pathogenic variants in the Polycystic Kidney Disease 2 (PKD2) gene are associated with Autosomal Dominant Polycystic Kidney Disease (ADPKD) in approximately 30% of cases. In recent years, the high-throughput sequencing techniques have significantly increased the number of variants identified in affected patients. Here, we described the peculiar effect of a PKD2 splicing variant, the c.1717-2A>G, identified in an Italian male patient with ADPKD. This variant led to the unusual and rare skipping of two consecutive exons, causing a large in-frame deletion. Methods: The genetic evaluation of the patient was performed using the Next-Generation Sequencing (NGS) assay Clinical Exome Solution® (SOPHiA Genetics). Bioinformatics analysis was performed using the SOPHiA DDM platform (SOPHiA Genetics). Prediction of pathogenicity was carried out by integrating several in silico tools. RNA evaluation was performed to test the effect of the variant on the PKD2 splicing using a Reverse-Transcription PCR coupled with cDNA sequencing. Results: NGS revealed the presence of the PKD2 c.1717-2A>G variant that lies in the canonical splice site of intron 7. This rare variant was predicted to have a significant impact on the splicing, proved by the RNA-based analysis. We identified the presence of a transcript characterised by the simultaneous skipping of exons 8 and 9, with a retained reading frame and the merging of exons 7-10. Conclusions: We described for the first time a dual-exon skip event related to the presence of a single-base substitution in the PKD2 gene in an ADPKD-affected patient. We assumed that the molecular basis of such a rare mechanism lies in the specific order of intron removal. The finding represents novel evidence of an alternative and unusual splicing mechanism in the PKD2 gene, adding insights to the pathogenesis of the ADPKD.

Functional Characterization of Splice Variants in the Diagnosis of Albinism

Albinism is a genetically heterogeneous disease in which 21 genes are known so far. Its inheritance mode is autosomal recessive except for one X-linked form. The molecular analysis of exonic sequences of these genes allows for about a 70% diagnostic rate. About half (15%) of the unsolved cases are heterozygous for one pathogenic or probably pathogenic variant. Assuming that the missing variant may be located in non-coding regions, we performed sequencing for 122 such heterozygous patients of either the whole genome (27 patients) or our NGS panel (95 patients) that includes, in addition to all exons of the 21 genes, the introns and flanking sequences of five genes, TYR, OCA2, SLC45A2, GPR143 and HPS1. Rare variants (MAF < 0.01) in trans to the first variant were tested by RT-PCR and/or minigene assay. Of the 14 variants tested, nine caused either exon skipping or the inclusion of a pseudoexon, allowing for the diagnosis of 11 patients. This represents 9.8% (12/122) supplementary diagnosis for formerly unsolved patients and 75% (12/16) of those in whom the candidate variant was in trans to the first variant. Of note, one missense variant was demonstrated to cause skipping of the exon in which it is located, thus shedding new light on its pathogenic mechanism. Searching for non-coding variants and testing them for an effect on RNA splicing is warranted in order to increase the diagnostic rate.

Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

Whole genome sequencing identifies elusive variants in genetically unsolved Italian inherited retinal disease patients

Inherited retinal diseases (IRDs) are a group of rare monogenic diseases with high genetic heterogeneity (pathogenic variants identified in over 280 causative genes). The genetic diagnostic rate for IRDs is around 60%, mainly thanks to the routine application of next-generation sequencing (NGS) approaches such as extensive gene panels or whole exome analyses. Whole-genome sequencing (WGS) has been reported to improve this diagnostic rate by revealing elusive variants, such as structural variants (SVs) and deep intronic variants (DIVs). We performed WGS on 33 unsolved cases with suspected autosomal recessive IRD, aiming to identify causative genetic variants in non-coding regions or to detect SVs that were unexplored in the initial screening. Most of the selected cases (30 of 33, 90.9%) carried monoallelic pathogenic variants in genes associated with their clinical presentation, hence we first analyzed the non-coding regions of these candidate genes. Whenever additional pathogenic variants were not identified with this approach, we extended the search for SVs and DIVs to all IRD-associated genes. Overall, we identified the missing causative variants in 11 patients (11 of 33, 33.3%). These included three DIVs in ABCA4, CEP290 and RPGRIP1; one non-canonical splice site (NCSS) variant in PROM1 and three SVs (large deletions) in EYS, PCDH15 and USH2A. For the previously unreported DIV in CEP290 and for the NCCS variant in PROM1, we confirmed the effect on splicing by reverse transcription (RT)-PCR on patient-derived RNA. This study demonstrates the power and clinical utility of WGS as an all-in-one test to identify disease-causing variants missed by standard NGS diagnostic methodologies.

A systematic assessment of the impact of rare canonical splice site variants on splicing using functional and in silico methods

Canonical splice site variants (CSSVs) are often presumed to cause loss-of-function (LoF) and are assigned very strong evidence of pathogenicity (according to American College of Medical Genetics/Association for Molecular Pathology criterion PVS1). The exact nature and predictability of splicing effects of unselected rare CSSVs in blood-expressed genes are poorly understood. We identified 168 rare CSSVs in blood-expressed genes in 112 individuals using genome sequencing, and studied their impact on splicing using RNA sequencing (RNA-seq). There was no evidence of a frameshift, nor of reduced expression consistent with nonsense-mediated decay, for 25.6% of CSSVs: 17.9% had wildtype splicing only and normal junction depths, 3.6% resulted in cryptic splice site usage and in-frame insertions or deletions, 3.6% resulted in full exon skipping (in frame), and 0.6% resulted in full intron inclusion (in frame). Blind to these RNA-seq data, we attempted to predict the precise impact of CSSVs by applying in silico tools and the ClinGen Sequence Variant Interpretation Working Group 2018 guidelines for applying PVS1 criterion. The predicted impact on splicing using (1) SpliceAI, (2) MaxEntScan, and (3) AutoPVS1, an automatic classification tool for PVS1 interpretation of null variants that utilizes Ensembl Variant Effect Predictor and MaxEntScan, was concordant with RNA-seq analyses for 65%, 63%, and 61% of CSSVs, respectively. In summary, approximately one in four rare CSSVs did not show evidence for LoF based on analysis of RNA-seq data. Predictions from in silico methods were often discordant with findings from RNA-seq. More caution may be warranted in applying PVS1-level evidence to CSSVs in the absence of functional data.

Novel PRKAR1A mutation in Carney complex: a case report and literature review

Objective

Carney complex is a rare autosomal dominant syndrome that has been shown to be associated with inactivation due to PRKAR1A mutations. We revealed a novel PRKAR1A gene mutation in Chinese patient with Carney complex and review the literature to enhance understanding of Carney complex.

Case presentation

A 23-year-old Chinese male patient with a family history cardiac myxoma was admitted to our Department of Endocrinology because of central obesity and hyperpigmentation. Physical examination revealed a maximum blood pressure of 150/93mmHg, a waist circumference of 102cm, a weight of 70kg, a height of 170cm, and a BMI of 24.22kg/m2. Additionally, there was spotty skin pigmentation on the lip mucosa, purple striae on the abdomen, thin skin on both legs, and visible veins. Blood examination revealed hypercortisolemia, decreased adrenocorticotropic hormone (ACTH) levels and failure to suppress cortisol with low and high-dose dexamethasone suppression tests. Magnetic resonance imaging (MRI) scan revealed multiple small adrenal nodules and Retroperitoneal neurogenic tumor. Genetic testing showed a novel heterozygous mutation in exon 5 of PRKAR1A (c.500_502 + 8delAAGGTAAGGGC). The patient underwent resection of the right adrenal gland and retroperitoneal neoplasms in 2020. Postoperative pathology following the right adrenal gland resection showed nodular hyperplasia of the adrenal cortex. The pathology from the retroperitoneal tumor resection revealed spindle cell tumors rich in pigment and cells. The patient was diagnosed as Carney complex according to Stratakis CA in 2001 guidelines. After long-term follow-up, the patient's condition was stable, with weight loss, waist circumference reduction, significantly lower cortisol levels, and normal blood lipids.

Conclusion

This case reported a Carney complex in a Chinese patient, characterized clinically by non-ACTH-dependent Cushing's syndrome, familial recurrent cardiac myxomas, psammomatous melanotic schwannoma (PMS) and skin and mucosal pigmentation. A novel subtype of PRKAR1A mutation was discovered, which may affect the characteristics of the PRKAR1A protein and contribute to the development of Carney complex.

SpliceAPP: an interactive web server to predict splicing errors arising from human mutations

BACKGROUND

Splicing variants are a major class of pathogenic mutations, with their severity equivalent to nonsense mutations. However, redundant and degenerate splicing signals hinder functional assessments of sequence variations within introns, particularly at branch sites. We have established a massively parallel splicing assay to assess the impact on splicing of 11,191 disease-relevant variants. Based on the experimental results, we then applied regression-based methods to identify factors determining splicing decisions and their respective weights.

RESULTS

Our statistical modeling is highly sensitive, accurately annotating the splicing defects of near-exon intronic variants, outperforming state-of-the-art predictive tools. We have incorporated the algorithm and branchpoint information into a web-based tool, SpliceAPP, to provide an interactive application. This user-friendly website allows users to upload any genetic variants with genome coordinates (e.g., chr15 74,687,208 A G), and the tool will output predictions for splicing error scores and evaluate the impact on nearby splice sites. Additionally, users can query branch site information within the region of interest.

CONCLUSIONS

In summary, SpliceAPP represents a pioneering approach to screening pathogenic intronic variants, contributing to the development of precision medicine. It also facilitates the annotation of splicing motifs. SpliceAPP is freely accessible using the link https://bc.imb.sinica.edu.tw/SpliceAPP . Source code can be downloaded at https://github.com/hsinnan75/SpliceAPP .

Identification of a novel CFAP61 homozygous splicing variant associated with multiple morphological abnormalities of the flagella

In this study, we investigated the role of a newly identified homozygous variant (c.1245 + 6T > C) in the CFAP61 gene in the development of multiple morphologically abnormal flagella (MMAF) in an infertile patient. Using exome sequencing, we identified this variant, which led to exon 12 skipping and the production of a truncated CFAP61 protein. Transmission electron microscopy analysis of the patient's spermatozoa revealed various flagellar abnormalities, including defective nuclear chromatin condensation, axoneme disorganization, and mitochondria embedded in residual cytoplasmic droplets. Despite a fertilization rate of 83.3% through ICSI, there was no successful pregnancy due to poor embryo quality.Our findings suggest a link between the identified CFAP61 variant and MMAF, indicating potential disruption in radial spokes' assembly or function crucial for normal ciliary motility. Furthermore, nearly half of the observed sperm heads displayed chromatin condensation defects, possibly contributing to the low blastulation rate. This case underscores the significance of genetic counseling and testing, particularly for couples dealing with infertility and MMAF. Early identification of such genetic variants can guide appropriate interventions and improve reproductive outcomes.

Increasing the diagnostic yield of childhood glaucoma cases recruited into the 100,000 Genomes Project

Childhood glaucoma (CG) encompasses a heterogeneous group of genetic eye disorders that is responsible for approximately 5% of childhood blindness worldwide. Understanding the molecular aetiology is key to improving diagnosis, prognosis and unlocking the potential for optimising clinical management. In this study, we investigated 86 CG cases from 78 unrelated families of diverse ethnic backgrounds, recruited into the Genomics England 100,000 Genomes Project (GE100KGP) rare disease cohort, to improve the genetic diagnostic yield. Using the Genomics England/Genomic Medicine Centres (GE/GMC) diagnostic pipeline, 13 unrelated families were solved (13/78, 17%). Further interrogation using an expanded gene panel yielded a molecular diagnosis in 7 more unrelated families (7/78, 9%). This analysis effectively raises the total number of solved CG families in the GE100KGP to 26% (20/78 families). Twenty-five percent (5/20) of the solved families had primary congenital glaucoma (PCG), while 75% (15/20) had secondary CG; 53% of this group had non-acquired ocular anomalies (including iris hypoplasia, megalocornea, ectopia pupillae, retinal dystrophy, and refractive errors) and 47% had non-acquired systemic diseases such as cardiac abnormalities, hearing impairment, and developmental delay. CYP1B1 was the most frequently implicated gene, accounting for 55% (11/20) of the solved families. We identified two novel likely pathogenic variants in the TEK gene, in addition to one novel pathogenic copy number variant (CNV) in FOXC1. Variants that passed undetected in the GE100KGP diagnostic pipeline were likely due to limitations of the tiering process, the use of smaller gene panels during analysis, and the prioritisation of coding SNVs and indels over larger structural variants, CNVs, and non-coding variants.

Noncoding de novo mutations in SCN2A are associated with autism spectrum disorders

Coding de novo mutations (DNMs) contribute to the risk for autism spectrum disorders (ASD), but the contribution of noncoding DNMs remains relatively unexplored. Here we use whole genome sequencing (WGS) data of 12,411 individuals (including 3,508 probands and 2,218 unaffected siblings) from 3,357 families collected in Simons Foundation Powering Autism Research for Knowledge (SPARK) to detect DNMs associated with ASD, while examining Simons Simplex Collection (SSC) with 6383 individuals from 2274 families to replicate the results. For coding DNMs, SCN2A reached exome-wide significance (p=2.06×10-11) in SPARK. The 618 known dominant ASD genes as a group are strongly enriched for coding DNMs in cases than sibling controls (fold change=1.51, p =1.13×10-5 for SPARK; fold change=1.86, p =2.06×10-9 for SSC). For noncoding DNMs, we used two methods to assess statistical significance: a point-based test that analyzes sites with a Combined Annotation Dependent Depletion (CADD) score ≥15, and a segment-based test that analyzes 1kb genomic segments with segment-specific background mutation rates (inferred from expected rare mutations in Gnocchi genome constraint scores). The point-based test identified SCN2A as marginally significant (p=6.12×10-4) in SPARK, yet segment-based test identified CSMD1, RBFOX1 and CHD13 as exome-wide significant. We did not identify significant enrichment of noncoding DNMs (in all 1kb segments or those with Gnocchi>4) in the 618 known ASD genes as a group in cases than sibling controls. When combining evidence from both coding and noncoding DNMs, we found that SCN2A with 11 coding and 5 noncoding DNMs exhibited the strongest significance (p=4.15×10-13). In summary, we identified both coding and noncoding DNMs in SCN2A associated with ASD, while nominating additional candidates for further examination in future studies.

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

Artificial Intelligence and Computational Biology in Gene Therapy: A Review

One of the trending fields in almost all areas of science and technology is artificial intelligence. Computational biology and artificial intelligence can help gene therapy in many steps including: gene identification, gene editing, vector design, development of new macromolecules and modeling of gene delivery. There are various tools used by computational biology and artificial intelligence in this field, such as genomics, transcriptomic and proteomics data analysis, machine learning algorithms and molecular interaction studies. These tools can introduce new gene targets, novel vectors, optimized experiment conditions, predict the outcomes and suggest the best solutions to avoid undesired immune responses following gene therapy treatment.

Novel Synonymous Variant in IL7R Causes Preferential Expression of the Soluble Isoform

PURPOSE

The interleukin-7 receptor (IL-7R) is primarily expressed on lymphoid cells and plays a crucial role in the development, proliferation, and survival of T cells. Autosomal recessive mutations that disrupt IL-7Rα chain expression give rise to a severe combined immunodeficiency (SCID), which is characterized by lymphopenia and a T-B+NK+ phenotype. The objective here was to diagnose two siblings displaying the T-B+NK+ SCID phenotype as initial clinical genetic testing did not detect any variants in known SCID genes.

METHODS

Whole genome sequencing (WGS) was utilized to identify potential variants causing the SCID phenotype. Splicing prediction tools were employed to assess the deleterious impact of the mutation. Polymerase Chain Reaction (PCR), Sanger sequencing, flow cytometry, and ELISA were then used to validate the pathogenicity of the detected mutation.

RESULTS

We discovered a novel homozygous synonymous mutation in the IL7R gene. Our functional studies indicate that this variant is pathogenic, causing exon 6, which encodes the transmembrane domain, to be preferentially spliced out.

CONCLUSION

In this study, we identified a novel rare synonymous mutation causing a loss of IL-7Rα expression at the cellular membrane. This case demonstrates the value of reanalyzing genetic data based on the clinical phenotype and highlights the significance of functional studies in determining the pathogenicity of genetic variants.

CFAP47 is a novel causative gene implicated in X-linked polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is a well-described condition in which ~80% of cases have a genetic explanation, while the genetic basis of sporadic cystic kidney disease in adults remains unclear in ~30% of cases. This study aimed to identify novel genes associated with polycystic kidney disease (PKD) in patients with sporadic cystic kidney disease in which a clear genetic change was not identified in established genes. A next-generation sequencing panel analyzed known genes related to renal cysts in 118 sporadic cases, followed by whole-genome sequencing on 47 unrelated individuals without identified candidate variants. Three male patients were found to have rare missense variants in the X-linked gene Cilia And Flagella Associated Protein 47 (CFAP47). CFAP47 was expressed in primary cilia of human renal tubules, and knockout mice exhibited vacuolation of tubular cells and tubular dilation, providing evidence that CFAP47 is a causative gene involved in cyst formation. This discovery of CFAP47 as a newly identified gene associated with PKD, displaying X-linked inheritance, emphasizes the need for further cases to understand the role of CFAP47 in PKD.

Mutation analysis and clinical profile of South African patients with Neurofibromatosis type 1 (NF1) phenotype

Neurofibromatosis type 1 (NF1) is an autosomal dominant genetic condition with complete age-dependent penetrance, variable expressivity and a global prevalence of ∼1/3,000. It is characteriszed by numerous café-au-lait macules, skin freckling in the inguinal or axillary regions, Lisch nodules of the iris, optic gliomas, neurofibromas, and tumour predisposition. The diagnostic testing strategy for NF1 includes testing for DNA single nucleotide variants (SNVs), copy number variants (CNVs) as well as RNA analysis for deep intronic and splice variants, which can cumulatively identify the causative variant in 95% of patients. In the present study, NF1 patients were screened using a next-generation sequencing (NGS) assay targeting NF1 exons and intron/exon boundaries for SNV and NF1 multiple ligation-dependent probe amplification (MLPA) analysis for CNV detection. Twenty-six unrelated Southern African patients clinically suspected of having NF1, based on the clinical diagnostic criteria developed by the National Institute of Health (NIH), were included in the current study. A detection rate of 58% (15/26) was obtained, with SNVs identified in 80% (12/15) using a targeted gene panel and NF1 gene deletion in 20% (3/15) identified using MLPA. Ten patients (38%) had no variants identified, although they met NF1 diagnostic criteria. One VUS was identified in this study in a patient that met NF1 diagnostic criteria, however there was no sufficient information to classify variant as pathogenic. The clinical features of Southern African patients with NF1 are similar to that of the known NF1 phenotype, with the exception of a lower frequency of plexiform neurofibromas and a higher frequency of developmental/intellectual disability compared to other cohorts. This is the first clinical and molecular characterisation of a Southern African ancestry NF1 cohort using both next-generation sequencing and MLPA analysis. A significant number of patients remained without a diagnosis following DNA-level testing. The current study offers a potential molecular testing strategy for our low resource environment that could benefit a significant proportion of patients who previously only received a clinical diagnosis without molecular confirmation.

Preclinical evaluation of CRISPR-based therapies for Noonan syndrome caused by deep-intronic LZTR1 variants

Gene variants in LZTR1 are implicated to cause Noonan syndrome associated with a severe and early-onset hypertrophic cardiomyopathy. Mechanistically, LZTR1 deficiency results in accumulation of RAS GTPases and, as a consequence, in RAS-MAPK signaling hyperactivity, thereby causing the Noonan syndrome-associated phenotype. Despite its epidemiological relevance, pharmacological as well as invasive therapies remain limited. Here, personalized CRISPR-Cas9 gene therapies might offer a novel alternative for a curative treatment in this patient cohort. In this study, by utilizing a patient-specific screening platform based on iPSC-derived cardiomyocytes from two Noonan syndrome patients, we evaluated different clinically translatable therapeutic approaches using small Cas9 orthologs targeting a deep-intronic LZTR1 variant to cure the disease-associated molecular pathology. Despite high editing efficiencies in cardiomyocyte cultures transduced with lentivirus or all-in-one adeno-associated viruses, we observed crucial differences in editing outcomes in proliferative iPSCs vs. non-proliferative cardiomyocytes. While editing in iPSCs rescued the phenotype, the same editing approaches did not robustly restore LZTR1 function in cardiomyocytes, indicating critical differences in the activity of DNA double-strand break repair mechanisms between proliferative and non-proliferative cell types and highlighting the importance of cell type-specific screens for testing CRISPR-Cas9 gene therapies.

Emerging applications of machine learning in genomic medicine and healthcare

The integration of artificial intelligence technologies has propelled the progress of clinical and genomic medicine in recent years. The significant increase in computing power has facilitated the ability of artificial intelligence models to analyze and extract features from extensive medical data and images, thereby contributing to the advancement of intelligent diagnostic tools. Artificial intelligence (AI) models have been utilized in the field of personalized medicine to integrate clinical data and genomic information of patients. This integration allows for the identification of customized treatment recommendations, ultimately leading to enhanced patient outcomes. Notwithstanding the notable advancements, the application of artificial intelligence (AI) in the field of medicine is impeded by various obstacles such as the limited availability of clinical and genomic data, the diversity of datasets, ethical implications, and the inconclusive interpretation of AI models' results. In this review, a comprehensive evaluation of multiple machine learning algorithms utilized in the fields of clinical and genomic medicine is conducted. Furthermore, we present an overview of the implementation of artificial intelligence (AI) in the fields of clinical medicine, drug discovery, and genomic medicine. Finally, a number of constraints pertaining to the implementation of artificial intelligence within the healthcare industry are examined.

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants

BACKGROUND

Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis.

RESULTS

Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity.

CONCLUSIONS

By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

Bi-allelic truncating variants in CASP2 underlie a neurodevelopmental disorder with lissencephaly

Lissencephaly (LIS) is a malformation of cortical development due to deficient neuronal migration and abnormal formation of cerebral convolutions or gyri. Thirty-one LIS-associated genes have been previously described. Recently, biallelic pathogenic variants in CRADD and PIDD1, have associated with LIS impacting the previously established role of the PIDDosome in activating caspase-2. In this report, we describe biallelic truncating variants in CASP2, another subunit of PIDDosome complex. Seven patients from five independent families presenting with a neurodevelopmental phenotype were identified through GeneMatcher-facilitated international collaborations. Exome sequencing analysis was carried out and revealed two distinct novel homozygous (NM_032982.4:c.1156delT (p.Tyr386ThrfsTer25), and c.1174 C > T (p.Gln392Ter)) and compound heterozygous variants (c.[130 C > T];[876 + 1 G > T] p.[Arg44Ter];[?]) in CASP2 segregating within the families in a manner compatible with an autosomal recessive pattern. RNA studies of the c.876 + 1 G > T variant indicated usage of two cryptic splice donor sites, each introducing a premature stop codon. All patients from whom brain MRIs were available had a typical fronto-temporal LIS and pachygyria, remarkably resembling the CRADD and PIDD1-related neuroimaging findings. Other findings included developmental delay, attention deficit hyperactivity disorder, hypotonia, seizure, poor social skills, and autistic traits. In summary, we present patients with CASP2-related ID, anterior-predominant LIS, and pachygyria similar to previously reported patients with CRADD and PIDD1-related disorders, expanding the genetic spectrum of LIS and lending support that each component of the PIDDosome complex is critical for normal development of the human cerebral cortex and brain function.

Gene Distribution in Pediatric-Onset Inherited Peripheral Neuropathy: A Single Tertiary Center in Thailand

BACKGROUND

Inherited peripheral neuropathy presents a diagnostic and therapeutic challenge due to its association with mutations in over 100 genes. This condition leads to long-term disability and poses a substantial healthcare burden on society.

OBJECTIVE

This study aimed to investigate the distribution of genes and establish the genotype-phenotype correlations, focusing on pediatric-onset cases.

METHODS

Exome sequencing and other analytical techniques were employed to identify pathogenic variants, including duplication analysis of the PMP22 gene. Each patient underwent physical examination and electrophysiological studies. Genotypes were correlated with phenotypic features, such as age at disease onset and ulnar motor nerve conduction velocity.

RESULTS

We identified 35 patients with pediatric-onset inherited peripheral neuropathy. Pathogenic or likely pathogenic variants were confirmed in 24 out of 35 (68.6%) patients, with 4 of these variants being novel. A confirmed molecular diagnosis was achieved in 90.9% (10/11) of patients with demyelinating Charcot-Marie-Tooth disease (CMT) and 56.3% (9/16) of patients with axonal CMT. Among patients with infantile-onset CMT (≤2 years), the most common causative genes were MFN2 and NEFL, while GDAP1 and MFN2 were frequent causes among patients with childhood- or adolescent-onset CMT (3-9 years).

CONCLUSIONS

The MFN2 gene was the most commonly implicated gene, and the axonal type was predominant in this cohort of Thai patients with pediatric-onset inherited peripheral neuropathy.

[Genetic and phenotypic analysis of MYO15A rare variants associated with autosomal recessive hearing loss]

Objective:To analyze the phenotype and genotype characteristics of autosomal recessive hearing loss caused by MYO15A gene variants, and to provide genetic diagnosis and genetic counseling for patients and their families. Methods:Identification of MYO15A gene variants by next generation sequencing in two sporadic cases of hearing loss at Shanghai Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine. The sequence variants were verified by Sanger sequencing.The pathogenicity of these variants was determined according to the American College of Medical Genetics and Genomics（ACMG） variant classification guidelines, in conjuction with clinical data. Results:The probands of the two families have bilateral,severe or complete hearing loss.Four variants of MYO15A were identified, including one pathogenic variant that has been reported, two likely pathogenic variants,and one splicing variant of uncertain significance. Patient I carries c. 3524dupA（p. Ser1176Valfs*14）, a reported pathogenic variant, and a splicing variant c. 10082+3G>A of uncertain significance according to the ACMG guidelines. Patient I was treated with bilateral hearing aids with satisfactory effect, demonstrated average hearing thresholds of 37.5 dB in the right ear and 33.75 dB in the left ear. Patient Ⅱ carries c. 7441_7442del（p. Leu2481Glufs*86） and c. 10250_10252del（p. Ser3417del）,a pair of as likely pathogenic variants according to the ACMG guidelines. Patient Ⅱ, who underwent right cochlear implantation eight years ago, achieved scores of 9 on the Categorical Auditory Performance-Ⅱ（CAP-Ⅱ） and 5 on the Speech Intelligibility Rating（SIR）. Conclusion:This study's discovery of the rare c. 7441_7442del variant and the splicing variant c. 10082+3G>A in the MYO15A gene is closely associated with autosomal recessive hearing loss, expanding the MYO15A variant spectrum. Additionally, the pathogenicity assessment of the splicing variant facilitates classification of splicing variations.

Clinical and genetic analyses of a Swedish patient series diagnosed with ataxia

Hereditary ataxia is a heterogeneous group of complex neurological disorders. Next-generation sequencing methods have become a great help in clinical diagnostics, but it may remain challenging to determine if a genetic variant is the cause of the patient's disease. We compiled a consecutive single-center series of 87 patients from 76 families with progressive ataxia of known or unknown etiology. We investigated them clinically and genetically using whole exome or whole genome sequencing. Test methods were selected depending on family history, clinical phenotype, and availability. Genetic results were interpreted based on the American College of Medical Genetics criteria. For high-suspicion variants of uncertain significance, renewed bioinformatical and clinical evaluation was performed to assess the level of pathogenicity. Thirty (39.5%) of the 76 families had received a genetic diagnosis at the end of our study. We present the predominant etiologies of hereditary ataxia in a Swedish patient series. In two families, we established a clinical diagnosis, although the genetic variant was classified as "of uncertain significance" only, and in an additional three families, results are pending. We found a pathogenic variant in one family, but we suspect that it does not explain the complete clinical picture. We conclude that correctly interpreting genetic variants in complex neurogenetic diseases requires genetics and clinical expertise. The neurologist's careful phenotyping remains essential to confirm or reject a diagnosis, also by reassessing clinical findings after a candidate genetic variant is suggested. Collaboration between neurology and clinical genetics and combining clinical and research approaches optimizes diagnostic yield.

Waardenburg Syndrome: The Contribution of Next-Generation Sequencing to the Identification of Novel Causative Variants

Waardenburg syndrome (WS) is characterized by hearing loss and pigmentary abnormalities of the eyes, hair, and skin. The condition is genetically heterogeneous, and is classified into four clinical types differentiated by the presence of dystopia canthorum in type 1 and its absence in type 2. Additionally, limb musculoskeletal abnormalities and Hirschsprung disease differentiate types 3 and 4, respectively. Genes PAX3, MITF, SOX10, KITLG, EDNRB, and EDN3 are already known to be associated with WS. In WS, a certain degree of molecularly undetected patients remains, especially in type 2. This study aims to pinpoint causative variants using different NGS approaches in a cohort of 26 Brazilian probands with possible/probable diagnosis of WS1 (8) or WS2 (18). DNA from the patients was first analyzed by exome sequencing. Seven of these families were submitted to trio analysis. For inconclusive cases, we applied a targeted NGS panel targeting WS/neurocristopathies genes. Causative variants were detected in 20 of the 26 probands analyzed, these being five in PAX3, eight in MITF, two in SOX10, four in EDNRB, and one in ACTG1 (type 2 Baraitser-Winter syndrome, BWS2). In conclusion, in our cohort of patients, the detection rate of the causative variant was 77%, confirming the superior detection power of NGS in genetically heterogeneous diseases.

ODAD1 variants resulting from splice-site mutations retain partial function and cause primary ciliary dyskinesia with outer dynein arm defects

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous disorder caused by defects in motile ciliary function and/or structure. Outer dynein arm docking complex subunit 1 (ODAD1) is an important component of the outer dynein arm docking complex (ODA-DC). To date, 13 likely pathogenic mutations of ODAD1 have been reported. However, the pathogenesis of ODAD1 mutations remains elusive. To investigate the pathogenesis of splice-site mutations in ODAD1 discovered in this study and those reported previously, molecular and functional analyses were performed. Whole-exome sequencing revealed a compound mutation in ODAD1 (c.71-2A>C; c.598-2A>C) in a patient with PCD, with c.598-2A>C being a novel mutation that resulted in two mutant transcripts. The compound mutation in ODAD1 (c.71-2A>C; c.598-2A>C) led to aberrant splicing that resulted in the absence of the wild-type ODAD1 and defects of the outer dynein arm in ciliary axonemes, causing a decrease in ciliary beat frequency. Furthermore, we demonstrated that the truncated proteins resulting from splice-site mutations in ODAD1 could retain partial function and inhibit the interaction between wild-type ODAD1 and ODAD3. The results of this study expand the mutational and clinical spectrum of PCD, provide more evidence for genetic counseling, and offer new insights into gene-based therapeutic strategies for PCD.

O’Brien T, Perucca P, Ge Z, Anderson A. Applications for Deep Learning in Epilepsy Genetic Research

Epilepsy is a group of brain disorders characterised by an enduring predisposition to generate unprovoked seizures. Fuelled by advances in sequencing technologies and computational approaches, more than 900 genes have now been implicated in epilepsy. The development and optimisation of tools and methods for analysing the vast quantity of genomic data is a rapidly evolving area of research. Deep learning (DL) is a subset of machine learning (ML) that brings opportunity for novel investigative strategies that can be harnessed to gain new insights into the genomic risk of people with epilepsy. DL is being harnessed to address limitations in accuracy of long-read sequencing technologies, which improve on short-read methods. Tools that predict the functional consequence of genetic variation can represent breaking ground in addressing critical knowledge gaps, while methods that integrate independent but complimentary data enhance the predictive power of genetic data. We provide an overview of these DL tools and discuss how they may be applied to the analysis of genetic data for epilepsy research.

Novel Genetic and Phenotypic Expansion in GOSR2-Related Progressive Myoclonus Epilepsy

Biallelic variants in the Golgi SNAP receptor complex member 2 gene (GOSR2) have been reported in progressive myoclonus epilepsy with neurodegeneration. Typical clinical features include ataxia and areflexia during early childhood, followed by seizures, scoliosis, dysarthria, and myoclonus. Here, we report two novel patients from unrelated families with a GOSR2-related disorder and novel genetic and clinical findings. The first patient, a male compound heterozygous for the GOSR2 splice site variant c.336+1G>A and the novel c.364G>A,p.Glu122Lys missense variant showed global developmental delay and seizures at the age of 2 years, followed by myoclonus at the age of 8 years with partial response to clonazepam. The second patient, a female homozygous for the GOSR2 founder variant p.Gly144Trp, showed only mild fine motor developmental delay and generalized tonic-clonic seizures triggered by infections during adolescence, with seizure remission on levetiracetam. The associated movement disorder progressed atypically slowly during adolescence compared to its usual speed, from initial intention tremor and myoclonus to ataxia, hyporeflexia, dysmetria, and dystonia. These findings expand the genotype-phenotype spectrum of GOSR2-related disorders and suggest that GOSR2 should be included in the consideration of monogenetic causes of dystonia, global developmental delay, and seizures.