Splicing in action: assessing disease causing sequence changes

Functional analysis of the GPAT4 gene mutation predicted to affect splicing

The GPAT4 gene is considered as a potential functional candidate for single nucleotide polymorphism (SNP) studies in dairy cattle breeding due to its association with dairy performance in cattle by encoding an enzyme responsible for the presence of diacylglycerols and triacylglycerols in milk. Using the example of the GPAT4 gene, we applied the minigene splicing assay to analyze the functional consequences of its variant that was predicted to affect normal splicing. The results of functional analysis revealed the sequence variations (rs442541537), transfection experiments in a wild type and mutant cell line model system demonstrated that the investigated mutation in the second intron of the GPAT4 gene was responsible for the presence of a second exon in mature messenger RNA (mRNA). The cases of its absence in the spliced mature mRNA transcript resulted in a truncated dysfunctional protein due to the appearance of a stop codon. Thus, the discovered SNP led to alternative splicing in pre-mRNA by the 'cassette exon' ('exon skipping') mechanism. The studied mutation can potentially be a molecular genetic marker for alternative splicing for the GPAT4 gene and, therefore contributes to economic benefits in cattle breeding programs.

Antisense oligonucleotides enhance SLC20A2 expression and suppress brain calcification in a humanized mouse model

Primary familial brain calcification (PFBC) is a genetic neurological disease, yet no effective treatment is currently available. Here, we identified five novel intronic variants in SLC20A2 gene from six PFBC families. Three of these variants increased aberrant SLC20A2 pre-mRNA splicing by altering the binding affinity of splicing machineries to newly characterized cryptic exons, ultimately causing premature termination of SLC20A2 translation. Inhibiting the cryptic-exon incorporation with splice-switching ASOs increased the expression levels of functional SLC20A2 in cells carrying SLC20A2 mutations. Moreover, by knocking in a humanized SLC20A2 intron 2 sequence carrying a PFBC-associated intronic variant, the SLC20A2-KI mice exhibited increased inorganic phosphate (Pi) levels in cerebrospinal fluid (CSF) and progressive brain calcification. Intracerebroventricular administration of ASOs to these SLC20A2-KI mice reduced CSF Pi levels and suppressed brain calcification. Together, our findings expand the genetic etiology of PFBC and demonstrate ASO-mediated splice modulation as a potential therapy for PFBC patients with SLC20A2 haploinsufficiency.

GOAT rs10096097 and CREB1 rs6740584 single nucleotide polymorphisms are associated with type 2 diabetes mellitus in Egyptians

Diabetes mellitus (DM) is a chronic disorder that affects nearly half a billion people around the world and causes millions of deaths annually. Treatment of diabetes or related complications represents an economic burden not only for developing countries but also for the developed ones. Hence, new efficient therapeutic and preventive strategies and screening tools are necessary. The current work aimed to assess the potential association of single nucleotide polymorphisms (SNPs) in ghrelin O-acyltransferase (GOAT) rs10096097, cyclic adenosine monophosphate (cAMP) response element binding protein (CREB) rs6740584, and v-maf musculoaponeurotic fibrosarcoma oncogene homolog A (MafA) rs62521874 genes with type 2 DM susceptibility in Egyptians. A total of 96 patients with type 2 DM along with 72 healthy individuals participated in this study. Genotyping was executed via real-time polymerase chain reaction (PCR), and the serum protein levels of GOAT, CREB, and MafA were measured by enzyme-linked immunosorbent assay (ELISA). Genotyping revealed a significant association of GOAT rs10096097 and CREB1 rs6740584 SNPs with type 2 diabetes risk, with significantly higher GOAT rs10096097 G allele and CREB1 rs6740584 T allele frequencies in diabetic patients than in controls. However, insignificant association was identified between the MafA rs62521874 SNP and diabetes in the examined sample of the Egyptian residents. Serum GOAT, CREB1, and MafA protein levels did not vary significantly between diabetic and control individuals. Yet, significant variation in serum GOAT and CREB1 levels was detected between CREB1 rs6740584 genotypes within the diabetic group, with CT and TT genotype carriers showing higher levels than AA genotype patients. GOAT rs10096097 and CREB1 rs6740584, but not MafA rs62521874, SNPs are associated with type 2 diabetes risk in the studied Egyptians.

RBM7 deficiency promotes breast cancer metastasis by coordinating MFGE8 splicing switch and NF-kB pathway

Aberrant alternative splicing is well-known to be closely associated with tumorigenesis of various cancers. However, the intricate mechanisms underlying breast cancer metastasis driven by deregulated splicing events remain largely unexplored. Here, we unveiled that RBM7 is decreased in lymph node and distant organ metastases of breast cancer as compared to primary lesions and low expression of RBM7 is correlated with the reduced disease-free survival of breast cancer patients. Breast cancer cells with RBM7 depletion exhibited an increased potential for lung metastasis compared to scramble control cells. The absence of RBM7 stimulated breast cancer cell migration, invasion, and angiogenesis. Mechanistically, RBM7 controlled the splicing switch of MFGE8, favoring the production of the predominant isoform of MFGE8, MFGE8-L. This resulted in the attenuation of STAT1 phosphorylation and alterations in cell adhesion molecules. MFGE8-L exerted an inhibitory effect on the migratory and invasive capability of breast cancer cells, while the truncated isoform MFGE8-S, which lack the second F5/8 type C domain had the opposite effect. In addition, RBM7 negatively regulates the NF-κB cascade and an NF-κB inhibitor could obstruct the increase in HUVEC tube formation caused by RBM7 silencing. Clinically, we noticed a positive correlation between RBM7 expression and MFGE8 exon7 inclusion in breast cancer tissues, providing new mechanistic insights for molecular-targeted therapy in combating breast cancer.

Activation of Cryptic Donor Splice Sites within the UDP-Glucuronosyltransferase (UGT)1A First-Exon Region Generates Variant Transcripts That Encode UGT1A Proteins with Truncated Aglycone-Binding Domains

The human UDP-glucuronosyltransferases (UGTs) have crucial roles in metabolizing and clearing numerous small lipophilic compounds. The UGT1A locus generates nine UGT1A mRNAs, 65 spliced transcripts, and 34 circular RNAs. In this study, our analysis of published UGT-RNA capture sequencing (CaptureSeq) datasets identified novel splice junctions that predict 24 variant UGT1A transcripts derived from ligation of exon 2 to unique sequences within the UGT1A first-exon region using cryptic donor splice sites. Of these variants, seven (1A1_n1, 1A3_n3, 1A4_n4, 1A5_n1, 1A8_n2, 1A9_n2, 1A10_n7) are predicted to encode UGT1A proteins with truncated aglycone-binding domains. We assessed their expression profiles and deregulation in cancer using four RNA sequencing (RNA-Seq) datasets of paired normal and cancerous drug-metabolizing tissues from large patient cohorts. Variants were generally coexpressed with their canonical counterparts with a higher relative abundance in tumor than in normal tissues. Variants showed tissue-specific expression with high interindividual variability but overall low abundance. However, 1A8_n2 showed high abundance in normal and cancerous colorectal tissues, with levels that approached or surpassed canonical 1A8 mRNA levels in many samples. We cloned 1A8_n2 and showed expression of the predicted protein (1A8_i3) in human embryonic kidney (HEK)293T cells. Glucuronidation assays with 4-methylumbelliferone (4MU) showed that 1A8_i3 had no activity and was unable to inhibit the activity of 1A8_i1 protein. In summary, the activation of cryptic donor splice sites within the UGT1A first-exon region expands the UGT1A transcriptome and proteome. The 1A8_n2 cryptic donor splice site is highly active in colorectal tissues, representing an important cis-regulatory element that negatively regulates the function of the UGT1A8 gene through pre-mRNA splicing. SIGNIFICANT STATEMENT: The UGT1A locus generates nine canonical mRNAs, 65 alternately spliced transcripts, and 34 different circular RNAs. The present study reports a series of novel UDP-glucuronosyltransferase (UGT)1A variants resulting from use of cryptic donor splice sites in both normal and cancerous tissues, several of which are predicted to encode variant UGT1A proteins with truncated aglycone-binding domains. Of these, 1A8_n2 shows exceptionally high abundance in colorectal tissues, highlighting its potential role in the first-pass metabolism in gut through the glucuronidation pathway.

Tumor Antigens beyond the Human Exome

With the advent of immunotherapeutics, a new era in the combat against cancer has begun. Particularly promising are neo-epitope-targeted therapies as the expression of neo-antigens is tumor-specific. In turn, this allows the selective targeting and killing of cancer cells whilst healthy cells remain largely unaffected. So far, many advances have been made in the development of treatment options which are tailored to the individual neo-epitope repertoire. The next big step is the achievement of efficacious "off-the-shelf" immunotherapies. For this, shared neo-epitopes propose an optimal target. Given the tremendous potential, a thorough understanding of the underlying mechanisms which lead to the formation of neo-antigens is of fundamental importance. Here, we review the various processes which result in the formation of neo-epitopes. Broadly, the origin of neo-epitopes can be categorized into three groups: canonical, noncanonical, and viral neo-epitopes. For the canonical neo-antigens that arise in direct consequence of somatic mutations, we summarize past and recent findings. Beyond that, our main focus is put on the discussion of noncanonical and viral neo-epitopes as we believe that targeting those provides an encouraging perspective to shape the future of cancer immunotherapeutics.

Molecular and phenotypic characteristics of Bardet-Biedl syndrome in Chinese patients

BACKGROUND

Bardet-Biedl syndrome (BBS) is a type of non-motile ciliopathy. To date, 26 genes have been reported to be associated with BBS. However, BBS is genetically heterogeneous, with significant clinical overlap with other ciliopathies, which complicates diagnosis. Disability and mortality rates are high in BBS patients; therefore, it is urgent to improve our understanding of BBS. Thus, our study aimed to describe the genotypic and phenotypic spectra of BBS in China and to elucidate genotype-phenotype correlations.

METHODS

Twenty Chinese patients diagnosed with BBS were enrolled in this study. We compared the phenotypes of Chinese BBS patients in this study with those from other countries to analyze the phenotypic differences across patients worldwide. In addition, genotype-phenotype correlations were described for our cohort. We also summarized all previously reported cases of BBS in Chinese patients (71 patients) and identified common and specific genetic variants in the Chinese population.

RESULTS

Twenty-eight variants, of which 10 are novel, in 5 different BBS-associated genes were identified in 20 Chinese BBS patients. By comparing the phenotypes of BBSome-coding genes (BBS2,7,9) with those of chaperonin-coding genes (BBS10,12), we found that patients with mutations in BBS10 and 12 had an earlier age of onset (1.10 Vs. 2.20, p < 0.01) and diagnosis (4.64 Vs. 13.17, p < 0.01), whereas patients with mutations in BBS2, 7, and 9 had a higher body mass index (28.35 Vs. 24.21, p < 0.05) and more vision problems (p < 0.05). Furthermore, in 91 Chinese BBS patients, mutations were predominant in BBS2 (28.89%) and BBS7 (15.56%), and the most frequent variants were in BBS2: c.534 + 1G > T (10/182 alleles) and BBS7: c.1002delT (7/182 alleles), marking a difference from the genotypic spectra of BBS reported abroad.

CONCLUSIONS

We recruited 20 Chinese patients with BBS for genetic and phenotypic analyses, and identified common clinical manifestations, pathogenic genes, and variants. We also described the phenotypic differences across patients worldwide and among different BBS-associated genes. This study involved the largest cohort of Chinese patients with BBS, and provides new insights into the distinctive clinical features of specific pathogenic variants.

The first pineoblastoma case report of a patient with Sotos syndrome harboring NSD1 germline mutation

Germline mutations of NSD1 are associated with Sotos syndrome, characterized by distinctive facial features, overgrowth, and developmental delay. Approximately 3% of individuals with Sotos syndrome develop tumors. In this study, we describe an infant in pineoblastoma with facial anomalies, learning disability and mild autism at 1 years diagnosed as Sotos syndrome owing to carrying a novel mutation de novo germline NSD1 likely pathogenic variant. This patient expands both the mutation and phenotype spectrum of the Sotos Syndrome and provides new clinical insights into the potential mechanism of underlying pinealoblastoma pathology.

Digenic inheritance accounts for phenotypic variability in amelogenesis imperfecta

Amelogenesis imperfecta (AI) represents a group of clinically and genetically heterogeneous disorders that affect enamel formation and mineralization. Although AI is commonly considered a monogenic disorder, digenic inheritance is rarely reported. In this study, we recruited two nonconsanguineous Chinese families exhibiting diverse phenotypes of enamel defects among affected family members. Digenic variants were discovered in both probands. In family 1, the proband inherited a paternal frameshift variant in LAMA3 (NM_198129.4:c.3712dup) and a maternal deletion encompassing the entire AMELX gene. This resulted in a combined hypoplastic and hypomineralized AI phenotype, which was distinct from the parents' manifestations. In family 2, whole-exome sequencing analysis revealed the proband carried a maternal heterozygous splicing variant in COL17A1 (NC_000010.11 (NM_000494.3): c.4156 + 2dup) and compound heterozygous variants in RELT (paternal: NM_032871.4:c.260A > T; maternal: NM_032871.4:c.521 T > G). These genetic changes caused the abundant irregular enamel defects observed in the proband, whereas other affected family members carrying heterozygous variants in both COL17A1 and RELT displayed only horizontal grooves as their phenotype. The pathogenicity of the novel COL17A1 splice site variant was confirmed through RT-PCR and minigene assay. This study enhances our understanding by highlighting the potential association between the co-occurrence of variants in two genes and variable phenotypes observed in AI patients.

Nonsense-mediated mRNA decay affects hyperactive root formation in oculo-facio-cardio-dental syndrome via up-frameshift protein 1

OBJECTIVES

Oculo-facio-cardio-dental (OFCD) syndrome is a rare X-linked genetic disorder caused by mutations in the BCL6 co-repressor (BCOR) and is mainly characterized by radiculomegaly (elongated dental roots). All BCOR mutations reported to date have been associated with premature termination codons, indicating that nonsense-mediated mRNA decay (NMD) might play a vital role in the pathogenesis of OFCD syndrome. However, the molecular mechanisms underlying NMD remain unclear. In this study, we investigated the involvement of up-frameshift protein 1 (UPF1), which plays a central role in NMD, in the hyperactive root formation caused by BCOR mutations.

METHODS

Periodontal ligament cells, isolated from a Japanese woman with a c.3668delC frameshift mutation in BCOR, and primary human periodontal ligament fibroblasts (HPdLFs) were used for an RNA immunoprecipitation assay to confirm the binding of UPF1 to mutated BCOR. Additionally, the effects of UPF1 on the BCOR transcription levels and corresponding gene expression were determined by performing relative quantitative real-time polymerase chain reactions.

RESULTS

RNA immunoprecipitation revealed that UPF1 binds to exon 9 of mutated BCOR. Additionally, UPF1 knockdown via siRNA upregulated the transcription of BCOR, whereas overexpression of wild-type and mutated BCOR with the same frameshift mutation in HPdLFs altered bone morphogenetic protein 2 (BMP2) expression.

CONCLUSIONS

Our findings indicate that BCOR mutations regulate the transcription of BCOR via UPF1, which may in turn regulate the expression of BMP2. NMD, caused by a c.3668delC mutation, potentially leads to an OFCD syndrome phenotype, including elongated dental roots.

Spectrum of variants associated with inherited retinal dystrophies in Northeast Mexico

BACKGROUND

Inherited retinal dystrophies are hereditary diseases which have in common the progressive degeneration of photoreceptors. They are a group of diseases with clinical, genetic, and allelic heterogeneity. There is limited information regarding the genetic landscape of inherited retinal diseases in Mexico, therefore, the present study was conducted in the northeast region of the country.

METHODS

Patients with inherited retinal dystrophies were included. A complete history, full ophthalmological and medical genetics evaluations, and genetic analysis through a targeted NGS panel for inherited retinal dystrophies comprising at least 293 genes were undertaken.

RESULTS

A total of 126 patients were included. Cases were solved in 74.6% of the study's population. Retinitis pigmentosa accounted for the most found inherited retinal disease. Ninety-nine causal variants were found, being USH2A and ABCA4 the most affected genes (26 and 15 cases, respectively).

CONCLUSIONS

The present study documents the most prevalent causative genes in IRDs, as USH2A, in northeastern Mexico. This contrasts with previous reports of IRDs in other zones of the country. Further studies, targeting previously unstudied populations in Mexico are important to document the genetic background of inherited retinal dystrophies in the country.

Relevance of Extending FGFR3 Gene Analysis in Osteochondrodysplasia to Non-Coding Sequences: A Case Report

Skeletal dysplasia, also called osteochondrodysplasia, is a category of disorders affecting bone development and children's growth. Up to 552 genes, including fibroblast growth factor receptor 3 (FGFR3), have been implicated by pathogenic variations in its genesis. Frequently identified causal mutations in osteochondrodysplasia arise in the coding sequences of the FGFR3 gene: c.1138G>A and c.1138G>C in achondroplasia and c.1620C>A and c.1620C>G in hypochondroplasia. However, in some cases, the diagnostic investigations undertaken thus far have failed to identify the causal anomaly, which strengthens the relevance of the diagnostic strategies being further refined. We observed a Caucasian adult with clinical and radiographic features of achondroplasia, with no common pathogenic variant. Exome sequencing detected an FGFR3(NM_000142.4):c.1075+95C>G heterozygous intronic variation. In vitro studies showed that this variant results in the aberrant exonization of a 90-nucleotide 5' segment of intron 8, resulting in the substitution of the alanine (Ala359) for a glycine (Gly) and the in-frame insertion of 30 amino acids. This change may alter FGFR3's function. Our report provides the first clinical description of an adult carrying this variant, which completes the phenotype description previously provided in children and confirms the recurrence, the autosomal-dominant pathogenicity, and the diagnostic relevance of this FGFR3 intronic variant. We support its inclusion in routinely used diagnostic tests for osteochondrodysplasia. This may increase the detection rate of causal variants and therefore could have a positive impact on patient management. Finally, FGFR3 alteration via non-coding sequence exonization should be considered a recurrent disease mechanism to be taken into account for new drug design and clinical trial strategies.

A novel homozygous splice-site mutation of JK gene leads to Jk(a-b-) phenotype

OBJECTIVES

This study aimed to investigate the molecular mechanism of the Jk(a-b-) phenotype in a Chinese transfusion patient.

BACKGROUND

Many different mutation types relating to Jk(a-b-) phenotype have been reported. However, the splice-site mutation is relatively rare and the related functional verification is lacking.

MATERIALS AND METHODS

In this study, the blood sample was collected from a transfusion patient with the Jk(a-b-) phenotype. Serotyping was performed using routine serological methods. The exons sequences and coding regions of the JK gene were amplified using polymerase chain reaction and directly sequenced. To perform a minigene splicing assay, the intronic mutation sequences were cloned into a pSPL3 splice reporting vector. The splicing reporter minigene assay was performed in HEK 293T cells.

RESULTS

The Jk(a-b-) phenotype of the blood sample was identified through serological testing. Sequencing results revealed that the sample had a novel homozygous splice-site mutation JK*02N (NM_015865.7: c.663+3A>C). Further analysis, including cDNA sequencing and minigene splicing assay, confirmed that the novel splice-site mutation resulted in exon skipping. Interestingly, different numbers of exons being skipped were obtained by the two methods.

CONCLUSION

This study revealed a novel homozygous splicing-site mutation associated with the Jk(a-b-) phenotype in Chinese population. Our results emphasise the importance of the in vitro functional method minigene splicing assay, while also acknowledging its potential limitations when compared to cDNA sequencing.

Germline mutational variants of Turkish ovarian cancer patients suspected of Hereditary Breast and Ovarian Cancer (HBOC) by next-generation sequencing

Hereditary Breast and Ovarian Cancer (HBOC) syndrome is characterized by an increased risk of developing breast cancer (BC) and ovarian cancer (OC) due to inherited genetic mutations. Understanding the genetic variants associated with HBOC is crucial for identifying individuals at high risk and implementing appropriate preventive measures. The study included 630 Turkish OC patients with confirmed diagnostic criteria of The National Comprehensive Cancer Network (NCCN) concerning HBOC. Genomic DNA was extracted from peripheral blood samples, and targeted Next-generation sequencing (NGS) was performed. Bioinformatics analysis and variant interpretation were conducted to identify pathogenic variants (PVs). Our analysis revealed a spectrum of germline pathogenic variants associated with HBOC in Turkish OC patients. Notably, several pathogenic variants in BRCA1, BRCA2, and other DNA repair genes were identified. Specifically, we observed germline PVs in 130 individuals, accounting for 20.63% of the total cohort. 76 distinct PVs in genes, BRCA1 (40 PVs), BRCA2 (29 PVs), ATM (1 PV), CHEK2 (2 PVs), ERCC2 (1 PV), MUTYH (1 PV), RAD51C (1 PV), and TP53 (1PV) and also, two different PVs (i.e., c.135-2 A>G p.? in BRCA1 and c.6466_6469delTCTC in BRCA2) were detected in a 34-year-old OC patient. In conclusion, our study contributes to a better understanding of the genetic variants underlying HBOC in Turkish OC patients. These findings provide valuable insights into the genetic architecture of HBOC in the Turkish population and shed light on the potential contribution of specific germline PVs to the increased risk of OC.

Clinical-genetic analysis of selected genes involved in the development of the human skeleton in 128 Czech patients with suspected congenital skeletal abnormalities

BACKGROUND

Congenital skeletal abnormalities are a heterogeneous group of diseases most commonly associated with small or disproportionate growth, cranial and facial dysmorphisms, delayed bone maturation, etc. Nonetheless, no detailed genotype-phenotype correlation in patients with specific genetic variants is readily available. Ergo, this study focuses on the analysis of patient phenotypes with candidate variants in genes involved in bone growth as detected by molecular genetic analysis.

METHODS

In this study we used molecular genetic methods to analyse the ACAN, COL2A1, FGFR3, IGFALS, IGF1, IGF1R, GHR, NPR2, STAT5B and SHOX genes in 128 Czech children with suspected congenital skeletal abnormalities. Pathogenic variants and variants of unclear clinical significance were identified and we compared their frequency in this study cohort to the European non-Finnish population. Furthermore, a prediction tool was utilised to determine their possible impact on the final protein. All clinical patient data was obtained during pre-test genetic counselling.

RESULTS

Pathogenic variants were identified in the FGFR3, GHR, COL2A1 and SHOX genes in a total of six patients. Furthermore, we identified 23 variants with unclear clinical significance and high allelic frequency in this cohort of patients with skeletal abnormalities. Five of them have not yet been reported in the scientific literature.

CONCLUSION

Congenital skeletal abnormalities may lead to a number of musculoskeletal, neurological, cardiovascular problems. Knowledge of specific pathogenic variants may help us in therapeutic procedures.

Opportunistic genetic screening increases the diagnostic yield and is medically valuable for care of patients and their relatives with hereditary cancer

BACKGROUND

Multigene panel testing by next-generation sequencing (MGP-NGS) enables the detection of germline pathogenic or likely pathogenic variants (PVs/LPVs) in genes beyond those associated with a certain cancer phenotype. Opportunistic genetic screening based on MGP-NGS in patients with suspicion of hereditary cancer reveals these incidental findings (IFs).

METHODS

MGP-NGS was performed in patients who fulfilled the clinical criteria to undergo genetic testing according to the Catalan Health Service guidelines. Variants were classified following the American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines and the Cancer Variant Interpretation Group UK guidelines.

RESULTS

IFs were identified in 10 (1.22%) of the 817 patients who underwent MGP-NGS. The mean age at cancer diagnosis was 49.4±9.5 years. Three IFs (30.0%) were detected in PMS2, two (20.0%) in ATM and TP53 and one (10.0%) in MSH6, NTHL1 and VHL. Seven (70.0%) IFs were single-nucleotide substitutions, two (20.0%) were deletions and one (10.0%) was a duplication. Three (30.0) IFs were located in intronic regions, three (30.3%) were nonsense, two (20.0%) were frameshift and two (20.0%) were missense variations. Six (60.0%) IFs were classified as PVs and four (40.0%) as LPVs.

CONCLUSIONS

Opportunistic genetic screening increased the diagnostic yield by 1.22% in our cohort. Most of the identified IFs were present in clinically actionable genes (n=7; 70.0%), providing these families with an opportunity to join cancer early detection programmes, as well as secondary cancer prevention. IFs might facilitate the diagnosis of asymptomatic individuals and the early management of cancer once it develops.

Genome editing of patient-derived iPSCs identifies a deep intronic variant causing aberrant splicing in hemophilia A

The importance of genetic diagnosis for patients with hemophilia has been recently demonstrated. However, the pathological variant cannot be identified in some patients. Here, we aimed to identify the pathogenic intronic variant causing hemophilia A using induced pluripotent stem cells (iPSCs) from patients and genome editing. We analyzed siblings with moderate hemophilia A and without abnormalities in the F8 exon. Next-generation sequencing of the entire F8 revealed 23 common intron variants. Variant effect predictor software indicated that the deep intronic variant at c.5220-8563A>G (intron 14) might act as a splicing acceptor. We developed iPSCs from patients and used genome editing to insert the elongation factor 1α promoter to express F8 messenger RNA (mRNA). Then, we confirmed the existence of abnormal F8 mRNA derived from aberrant splicing, resulting in a premature terminal codon as well as a significant reduction in F8 mRNA in iPSCs due to nonsense-mediated RNA decay. Gene repair by genome editing recovered whole F8 mRNA expression. Introduction of the intron variant into human B-domain-deleted F8 complementary DNA suppressed factor VIII (FVIII) activity and produced abnormal FVIII lacking the light chain in HEK293 cells. Furthermore, genome editing of the intron variant restored FVIII production. In summary, we have directly proven that the deep intronic variant in F8 results in aberrant splicing, leading to abnormal mRNA and nonsense-mediated RNA decay. Additionally, genome editing targeting the variant restored F8 mRNA and FVIII production. Our approach could be useful not only for identifying causal variants but also for verifying the therapeutic effect of personalized genome editing.

Long-term follow-up of children with carbamoyl phosphate synthase 1 deficiency detected in newborn screening

OBJECTIVES

To investigate genotype-phenotype characteristics and long-term prognosis of neonatal carbamoyl phosphate synthetase 1 (CPS1) deficiency among children through newborn screening in Zhejiang province.

METHODS

The clinical and follow-up data of children with CPS1 deficiency detected through neonatal screening and confirmed by tandem mass spectrometry and genetic testing in Zhejiang Province Newborn Disease Screening Center from September 2013 to August 2023 were retrospectively analyzed.

RESULTS

A total of 4 056 755 newborns were screened and 6 cases of CPS1 deficiency were diagnosed through phenotypic and genetic testing. Ten different variations of CPS1 genewere identified in genetic testing, including 2 known pathogenic variations (c.2359C>T and c.1549+1G>T) and 8 unreported variations (c.3405-1G>T, c.2372C>T, c.1436C>T, c.2228T>C, c.2441G>A, c.3031G>A, c.3075T>C and c.390-403del). All patients had decreased citrulline levels (2.72-6.21 μmol/L), and varying degrees of elevated blood ammonia. The patients received restricted natural protein intake (special formula), arginine and supportive therapy after diagnosis, and were followed-up for a period ranging from 9 months to 10 years. Three patients experienced hyperammonemia, and one patient each had attention deficit hyperactivity disorder, transient facial twitching and increased muscle tone. One patient died, while the other five surviving patients had normal scores of the Ages & Stages Questionnaires (ASQ) and Griffiths Development Scales up to the present time; 4 cases had combined height or weight lag and one case was normal in height and weight.

CONCLUSIONS

Low citrulline levels and hyperammonemia are common in CPS1 deficiency patients in Zhejiang. Most gene variants identified were specific to individual families, and no hotspot mutations were found. Early diagnosis through newborn screening and following standardized treatment can significantly improve the prognosis of the patients.

Dual Role of Dysfunctional Asc-1 Transporter in Distinct Human Pathologies, Human Startle Disease, and Developmental Delay

Human startle disease is associated with mutations in distinct genes encoding glycine receptors, transporters or interacting proteins at glycinergic synapses in spinal cord and brainstem. However, a significant number of diagnosed patients does not carry a mutation in the common genes GLRA1, GLRB, and SLC6A5 Recently, studies on solute carrier 7 subfamily 10 (SLC7A10; Asc-1, alanine-serine-cysteine transporter) knock-out (KO) mice displaying a startle disease-like phenotype hypothesized that this transporter might represent a novel candidate for human startle disease. Here, we screened 51 patients from our patient cohort negative for the common genes and found three exonic (one missense, two synonymous), seven intronic, and single nucleotide changes in the 5' and 3' untranslated regions (UTRs) in Asc-1. The identified missense mutation Asc-1G307R from a patient with startle disease and developmental delay was investigated in functional studies. At the molecular level, the mutation Asc-1G307R did not interfere with cell-surface expression, but disrupted glycine uptake. Substitution of glycine at position 307 to other amino acids, e.g., to alanine or tryptophan did not affect trafficking or glycine transport. By contrast, G307K disrupted glycine transport similar to the G307R mutation found in the patient. Structurally, the disrupted function in variants carrying positively charged residues can be explained by local structural rearrangements because of the large positively charged side chain. Thus, our data suggest that SLC7A10 may represent a rare but novel gene associated with human startle disease and developmental delay.

A novel splicing variation in L1CAM is responsible for recurrent fetal hydrocephalus

BACKGROUND

The L1 cell adhesion molecule (L1CAM, OMIM 308840) gene is primarily expressed in the nervous system and encodes the L1 adhesion molecule protein. Variations in L1CAM cause a wide spectrum of X-linked neurological disorders summarized as the L1 syndrome.

METHODS

We report a 29-year-old pregnant woman who experienced multiple adverse pregnancy outcomes due to recurrent fetal hydrocephalus with an X-linked recessive inheritance. Genomic DNA was extracted from the third aborted male fetus and analyzed via trio whole-exome sequencing (WES). Total RNA was isolated from the pregnant woman to assess splicing variation at the mRNA level, and amniotic fluid was extracted from the woman for prenatal diagnosis on her fourth fetus.

RESULTS

All four male fetuses were affected by severe hydrocephalus. We identified a maternally derived hemizygous splicing variation NM_000425.5:[c.3046 + 1G > A]; NP_000416.1 p.(Gly1016AspfsTer6) (chrX:153130275) in Intron 22 of the L1CAM. This variation disrupts the donor splice site and causes the retention of Intron 22, which results in frame shift and a premature termination codon at position 1021 with the ability to produce a truncated protein without the fifth fibronectin-repeat III, transmembrane, and cytoplasmic domains or to induce the degradation of mRNAs by nonsense-mediated mRNA decay. The same hemizygous variant was also detected in the amniocytes.

CONCLUSION

This report enhances our knowledge of genetic and phenotypic characteristics of X-linked fetal hydrocephalus, providing a new genetic basis for prenatal diagnosis and pre-implantation prenatal diagnosis.

Investigation of a Fused in Sarcoma Splicing Mutation in a Chinese Amyotrophic Lateral Sclerosis Patient

OBJECTIVE

Genetic mutations of fused in sarcoma (FUS) causing amyotrophic lateral sclerosis (ALS) may disrupt mRNA splicing events. For example, the FUS c.1394-2delA variant was reported in two western ALS patients, but its molecular mechanism is unclear. In this study, we aim to investigate FUS splice site mutations in Chinese ALS patients.

METHODS

Sanger sequencing was used to identify FUS splicing mutations in Chinese ALS patients. We combined a deep learning tool (SpliceAI), RNA sequencing, and RT-PCR/RT-qPCR to analyze the effect of FUS c.1394-2delA mutation on RNA splicing and expression. AlphaFold was used to predict the protein structure of mutant FUS. In transfected cell lines, we used immunofluorescence to assess cytoplasmic mislocalization of mutant FUS protein.

RESULTS

We identified a de novo FUS splice acceptor site mutation (c.1394-2delA, p. Gly466Valfs*14) in one Chinese sporadic ALS patient, which is linked to exon 14 skipping, and upregulated total FUS mRNA expression. The FUS splice site mutation was predicted to be translated into a truncated protein product at C-terminal. In vitro studies revealed that the FUS mutation increased cytoplasmic mislocalization in both HEK293T and SH-SY5Y cells.

CONCLUSIONS

We identified a de novo FUS splicing mutation (c.1394-2delA, p. Gly466Valfs*14) in 1 out of 233 Chinese ALS patients. It caused abnormal RNA splicing, upregulated gene expression, truncated FUS translation, and cytosolic mislocalization. Our findings suggested that FUS splice site mutation is rare in Chinese ALS patients and extended our knowledge of molecular mechanisms of the FUS c.1394-2delA mutation.

Case report: Analysis of a gene variant and prenatal diagnosis in a family with megalencephalic leukoencephalopathy with subcortical cysts

Megalencephalic leukoencephalopathy with subcortical cysts (MLC) is a rare inherited cerebral white matter disorder in children. Pathogenic variations in the causative gene MLC1 are found in approximately 76% of patients and are inherited in an autosomal recessive manner. In this study, we identified an IVS2 + 1delG variant in MLC1 in the firstborn girl of a pregnant woman who has the clinical features of MLC, including macrocephaly, motor development delay, progressive functional deterioration, and myelinopathy, whereas no obvious subcortical cysts were observed by magnetic resonance imaging of the brain. The proband is homozygous for the IVS2 + 1delG mutation, which was inherited from the parents. This variant disrupts the donor splice site, causing an abnormal transcript that results in a premature termination codon and produces a truncated protein, which was confirmed to affect splicing by MLC1 cDNA analysis. This variant was also detected in family members, and a prenatal diagnosis for the fetus was undertaken. Eventually, the couple gave birth to an unaffected baby. Furthermore, we conducted a long-term follow-up of the proband's clinical course. This report improves our understanding of the genetic and phenotypic characteristics of MLC and provides a new genetic basis for prenatal diagnosis and genetic counseling.

The spectrum of hereditary neuromuscular disorders in the Pakistani population

Hereditary neuromuscular disorders (NMDs) are a broad group of clinically heterogeneous disorders with varying inheritance patterns, that are associated with over 500 implicated genes. In the context of a highly consanguineous Pakistani population, we expect that autosomal recessive NMDs may have a higher prevalence compared with patients of European descent. This is the first study to offer a detailed description of the spectrum of genes causing hereditary NMDs in the Pakistani population using NGS testing. To study the clinical and genetic profiles of patients presenting for evaluation of a hereditary neuromuscular disorder. This is a retrospective chart review of patients seen in the Neuromuscular Disorders Clinic and referred to the Genetics Clinic with a suspected hereditary neuromuscular disorder, between 2016 and 2020 at the Aga Khan University Hospital, Karachi and Mukhtiar A. Sheikh Hospital, Multan, Pakistan. The genetic testing for these patients included NGS-based single gene sequencing, NGS-based multi-gene panel and whole exome sequencing. In a total of 112 patients studied, 35 (31.3%) were female. The mean age of onset in all patients was 14.6 years (SD ±12.1 years), with the average age at presentation to the clinic of 22.4 years (SD ±14.10 years). Forty-seven (41.9%) patients had a positive genetic test result, 53 (47.3%) had one or more variants of uncertain significance (VUS), and 12 (10.7%) had a negative result. Upon further genotype-phenotype correlation and family segregation analysis, the diagnostic yield improved, with 59 (52.7%) patients reaching a diagnosis of a hereditary NMD. We also report probable founder variants in COL6A2, FKTN, GNE, and SGCB, previously reported in populations that have possible shared ancestry with the Pakistani population. Our findings reemphasizes that the rate of VUSs can be reduced by clinical correlation and family segregation studies.

Genomic variants exclusively identified in children with birth defects and concurrent malignant tumors predispose to cancer development

Children with birth defects (BD) express distinct clinical features that often have various medical consequences, one of which is predisposition to the development of cancers. Identification of the underlying genetic mechanisms related to the development of cancer in BD patients would allow for preventive measures. We performed a whole genome sequencing (WGS) study on blood-derived DNA samples from 1566 individuals without chromosomal anomalies, including 454 BD probands with at least one type of malignant tumors, 767 cancer-free BD probands, and 345 healthy individuals. Exclusive recurrent variants were identified in BD-cancer and BD-only patients and mapped to their corresponding genomic regions. We observed statistically significant overlaps for protein-coding/ncRNA with exclusive variants in exons, introns, ncRNAs, and 3'UTR regions. Exclusive exonic variants, especially synonymous variants, tend to occur in prior exons locus in BD-cancer children. Intronic variants close to splicing site (< 500 bp from exon) have little overlaps in BD-cancer and BD-only patients. Exonic variants in non-coding RNA (ncRNA) tend to occur in different ncRNAs exons regardless of the overlaps. Notably, genes with 5' UTR variants are almost mutually exclusive between the two phenotypes. In conclusion, we conducted the first genomic study to explore the impact of recurrent variants exclusive to the two distinguished clinical phenotypes under study, BD with or without cancer, demonstrating enrichment of selective protein-coding/ncRNAs differentially expressed between these two phenotypes, suggesting that selective genetic factors may underlie the molecular processes of pediatric cancer development in BD children.

ICOS gene polymorphisms in systemic lupus erythematosus: A case-control study

The inducible T-cell costimulator (ICOS) may play an important role in adaptive immunity by regulating the interaction between T cells and antigen-presenting cells. Disruption of this molecule can lead to autoimmune diseases, in particular systemic lupus erythematosus (SLE). In this study, we aimed to explore the possible association between ICOS gene polymorphisms and SLE as well as their influence on disease susceptibility and clinical outcomes. A further objective was to assess the potential impact of these polymorphisms on RNA expression. A case-control study, including 151 patients with SLE, and 291 unrelated healthy controls (HC) matched in gender, and geographical origin, was performed to genotype two polymorphisms located in the ICOS gene: rs11889031 (-693 G/A) and rs10932029 (IVS1 + 173 T/C); using the polymerase chain reaction (PCR)-restriction fragment length polymorphism method. The different genotypes were validated by direct sequencing. The expression level of ICOS mRNA was assessed by quantitative PCR in peripheral blood mononuclear cells of SLE patients and HC. The results were analysed using Shesis and spss.20. Our results revealed a significant association between ICOS gene rs11889031 > CC genotype and SLE disease (codominant genetic model 1, (C/C vs. C/T), p = .001, odds ratio [OR] = 2.18 IC [1.36-3.49]); codominant genetic model 2, (C/C vs. T/T) p = .007, OR = 15.29 IC [1.97-118.5]); dominant genetic model, (C/C vs. C/T + T/T) p = .0001, OR = 2.44 IC [1.53-3.9]). Besides, there was a marginal association between rs11889031 > TT genotype and T allele with a protective role from SLE (recessive genetic model, p = .016, OR = 0.08 IC [0.01-0.63] and p = 7.6904E - 05, OR = 0.43 IC = [0.28-0.66], respectively). Moreover, statistical analysis indicated that the rs11889031 > CC genotype was linked with clinical and serological manifestations of SLE, including blood pressure, and anti-SSA antibodies production in SLE patients. However, the ICOS gene rs10932029 polymorphism was not associated with susceptibility to SLE. On the other side, we did not note any effect of the two selected polymorphisms on the level of ICOS mRNA gene expression. The study showed a significant predisposing association of the ICOS rs11889031 > CC genotype with SLE, in contrast to a protective effect of rs11889031 > TT genotype in Tunisian patients. Our results suggest that ICOS rs11889031 may act as a risk factor for SLE and could be used as a genetic susceptibility biomarker.

Cardiac splicing as a diagnostic and therapeutic target

Despite advances in therapeutics for heart failure and arrhythmias, a substantial proportion of patients with cardiomyopathy do not respond to interventions, indicating a need to identify novel modifiable myocardial pathobiology. Human genetic variation associated with severe forms of cardiomyopathy and arrhythmias has highlighted the crucial role of alternative splicing in myocardial health and disease, given that it determines which mature RNA transcripts drive the mechanical, structural, signalling and metabolic properties of the heart. In this Review, we discuss how the analysis of cardiac isoform expression has been facilitated by technical advances in multiomics and long-read and single-cell sequencing technologies. The resulting insights into the regulation of alternative splicing - including the identification of cardiac splice regulators as therapeutic targets and the development of a translational pipeline to evaluate splice modulators in human engineered heart tissue, animal models and clinical trials - provide a basis for improved diagnosis and therapy. Finally, we consider how the medical and scientific communities can benefit from facilitated acquisition and interpretation of splicing data towards improved clinical decision-making and patient care.

Identification of a Novel de novo Splicing Mutation in Duchenne Muscular Dystrophy Gene in an Iranian Family

Introduction

Duchenne muscular dystrophy (DMD) (NM_004006.3) is one of the most notable neuromuscular disorders of early years. The majority of DMD cases are caused by deletions or duplications in dystrophin, while point mutations are less prevalent in dystrophin abnormalities. It is a common knowledge that the severity of the disease depends on the effect of the mutation on the translational reading frame of the dystrophin mRNA.

Case Report

We studied an 8-year-old boy with relevant clinical presentations for DMD. Deletion/duplication screening was performed by using multiplex ligation-dependent probe amplification, and whole-exome sequencing was conducted in order to identify potential variants. A novel de novo splice site variant was identified in the DMD gene (DMD: c.8548-2A>G). To explore the effect of a novel variant in DMD, various in silico analyses were carried out to investigate the pathogenicity of the causative variant. To study the structure of a DMD protein and information on how the genetic variant impacts splicing site in models of wild-type and mutated DMD, we carried out different computational studies. Sanger sequencing was performed for the purpose of variant confirmation and familial segregation analysis.

Discussion

This novel de novo variant was predicted to have an effect on splicing, which leads to DMD due to its significant impacts on dystrophin functionality. The novel mutation would be expected to disrupt the protein structure.

Using a Dual CRISPR/Cas9 Approach to Gain Insight into the Role of LRP1B in Glioblastoma

LRP1B remains one of the most altered genes in cancer, although its relevance in cancer biology is still unclear. Recent advances in gene editing techniques, particularly CRISPR/Cas9 systems, offer new opportunities to evaluate the function of large genes, such as LRP1B. Using a dual sgRNA CRISPR/Cas9 gene editing approach, this study aimed to assess the impact of disrupting LRP1B in glioblastoma cell biology. Four sgRNAs were designed for the dual targeting of two LRP1B exons (1 and 85). The U87 glioblastoma (GB) cell line was transfected with CRISPR/Cas9 PX459 vectors. To assess LRP1B-gene-induced alterations and expression, PCR, Sanger DNA sequencing, and qRT-PCR were carried out. Three clones (clones B9, E6, and H7) were further evaluated. All clones presented altered cellular morphology, increased cellular and nuclear size, and changes in ploidy. Two clones (E6 and H7) showed a significant decrease in cell growth, both in vitro and in the in vivo CAM assay. Proteomic analysis of the clones' secretome identified differentially expressed proteins that had not been previously associated with LRP1B alterations. This study demonstrates that the dual sgRNA CRISPR/Cas9 strategy can effectively edit LRP1B in GB cells, providing new insights into the impact of LRP1B deletions in GBM biology.

Homozygous splice variant (c.1741-6G>A) of the COL6A1 gene in three patients with Ullrich congenital muscular dystrophy

The three major collagen VI genes: COL6A1, COL6A2, and COL6A3 encode microfibrillar components of extracellular matrices in multiple tissues including muscles and tendons. Pathogenic variants in the collagen VI genes cause collagen VI-related dystrophies representing a continuum of conditions from Bethlem myopathy at the milder end to Ullrich congenital muscular dystrophy at the more severe end. Here we describe a pathogenic variant in the COL6A1 gene (NM_001848.3; c.1741-6G>A) found in homozygosity in three patients with Ullrich congenital muscular dystrophy. The patients suffered from severe muscle impairment characterised by proximal weakness, distal hyperlaxity, joint contractures, wheelchair-dependency, and use of nocturnal non-invasive ventilation. The pathogenicity was verified by RNA analyses showing that the variant induced aberrant splicing leading to a frameshift and loss of function. The analyses were in line with immunocytochemistry studies of patient-derived skin fibroblasts and muscle tissue demonstrating impaired secretion of collagen VI into the extracellular matrix. Thereby, we add the variant c.1741-6G>A to the list of pathogenic, recessive, splice variants in COL6A1 causing Ullrich congenital muscular dystrophy. The variant is listed in ClinVar as of "uncertain significance" and "likely benign" and may presumably have been overlooked in other patients.

Intronic position +9 and -9 are potentially splicing sites boundary from intronic variants analysis of whole exome sequencing data

Whole exome sequencing (WES) can also detect some intronic variants, which may affect splicing and gene expression, but how to use these intronic variants, and the characteristics about them has not been reported. This study aims to reveal the characteristics of intronic variant in WES data, to further improve the clinical diagnostic value of WES. A total of 269 WES data was analyzed, 688,778 raw variants were called, among these 367,469 intronic variants were in intronic regions flanking exons which was upstream/downstream region of the exon (default is 200 bps). Contrary to expectation, the number of intronic variants with quality control (QC) passed was the lowest at the +2 and -2 positions but not at the +1 and -1 positions. The plausible explanation was that the former had the worst effect on trans-splicing, whereas the latter did not completely abolish splicing. And surprisingly, the number of intronic variants that passed QC was the highest at the +9 and -9 positions, indicating a potential splicing site boundary. The proportion of variants which could not pass QC filtering (false variants) in the intronic regions flanking exons generally accord with "S"-shaped curve. At +5 and -5 positions, the number of variants predicted damaging by software was most. This was also the position at which many pathogenic variants had been reported in recent years. Our study revealed the characteristics of intronic variant in WES data for the first time, we found the +9 and -9 positions might be a potentially splicing sites boundary and +5 and -5 positions were potentially important sites affecting splicing or gene expression, the +2 and -2 positions seem more important splicing site than +1 and -1 positions, and we found variants in intronic regions flanking exons over ± 50 bps may be unreliable. This result can help researchers find more useful variants and demonstrate that WES data is valuable for intronic variants analysis.

Role of non-coding variants in cardiovascular disease

Cardiovascular diseases (CVDs) constitute one of the significant causes of death worldwide. Different pathological states are linked to CVDs, which despite interventions and treatments, still have poor prognoses. The genetic component, as a beneficial tool in the risk stratification of CVD development, plays a role in the pathogenesis of this group of diseases. The emergence of genome-wide association studies (GWAS) have led to the identification of non-coding parts associated with cardiovascular traits and disorders. Variants located in functional non-coding regions, including promoters/enhancers, introns, miRNAs and 5'/3' UTRs, account for 90% of all identified single-nucleotide polymorphisms associated with CVDs. Here, for the first time, we conducted a comprehensive review on the reported non-coding variants for different CVDs, including hypercholesterolemia, cardiomyopathies, congenital heart diseases, thoracic aortic aneurysms/dissections and coronary artery diseases. Additionally, we present the most commonly reported genes involved in each CVD. In total, 1469 non-coding variants constitute most reports on familial hypercholesterolemia, hypertrophic cardiomyopathy and dilated cardiomyopathy. The application and identification of non-coding variants are beneficial for the genetic diagnosis and better therapeutic management of CVDs.

The remarkable journey of one female individual with ornithine transcarbamylase deficiency diagnosed post-mortem

Urea cycle disorders (UCDs) comprise a group of inborn errors of metabolism with impaired ammonia clearance and an incidence of ~1:35 000 individuals. First described in the 1970s, the diagnosis and management of these disorders has evolved dramatically. We report on a 59-year-old woman with a UCD who contributed to advances in the understanding and treatment of this group of disorders. This individual was diagnosed with carbamoyl phosphate synthetase 1 deficiency based on a biochemical assay under a research context predating genetic sequencing, treated longitudinally as having this metabolic disorder, and was among the first participants to trial UCD pharmaceutical therapies. She ultimately succumbed to a SARS-CoV-2 infection while maintaining unexpectedly normal ammonium levels. Postmortem genetic testing revealed ornithine transcarbamylase deficiency. This individual's contributions to the field of UCDs is discussed herein.

Analysis of BRCA1 and BRCA2 alternative splicing in predisposition to ovarian cancer

BACKGROUND

The mRNA splicing is regulated on multiple levels, resulting in the proper distribution of genes' transcripts in each cell and maintaining cell homeostasis. At the same time, the expression of alternative transcripts can change in response to underlying genetic variants, often missed during routine diagnostics.

AIM

The main aim of this study was to define the frequency of aberrant splicing in BRCA1 and BRCA2 genes in blood RNA extracted from ovarian cancer patients who were previously found negative for the presence of pathogenic alterations in the 25 most commonly analysed ovarian cancer genes, including BRCA1 and BRCA2.

MATERIAL AND METHODS

Frequency and spectrum of splicing alterations in BRCA1 and BRCA2 genes were analysed in blood RNA from 101 ovarian cancer patients and healthy controls (80 healthy women) using PCR followed by gel electrophoresis and Sanger sequencing. The expression of splicing events was examined using RT-qPCR.

RESULTS

We did not identify any novel, potentially pathogenic splicing alterations. Nevertheless, we detected six naturally occurring transcripts, named BRCA1ΔE9-10, BRCA1ΔE11, BRCA1ΔE11q, and BRCA2ΔE3, BRCA2ΔE12 and BRCA2ΔE17-18 of which three (BRCA1ΔE11q, BRCA1ΔE11 and BRCA2ΔE3) were significantly higher expressed in the ovarian cancer cohort than in healthy controls (p ≤ 0.0001).

CONCLUSIONS

This observation indicates that the upregulation of selected naturally occurring transcripts can be stimulated by non-genetic mechanisms and be a potential systemic response to disease progression and/or treatment. However, this hypothesis requires further examination.

Molecular Basis of Unequal Alternative Splicing of Human SCD5 and Its Alteration by Natural Genetic Variations

Alternative splicing (AS) is a major means of post-transcriptional control of gene expression, and provides a dynamic versatility of protein isoforms. Cancer-related AS disorders have diagnostic, prognostic and therapeutic values. Changes in the expression and AS of human stearoyl-CoA desaturase-5 (SCD5) are promising specific tumor markers, although the transcript variants (TVs) of the gene have not yet been confirmed. Our in silico, in vitro and in vivo study focuses on the distribution of SCD5 TVs (A and B) in human tissues, the functionality of the relevant splice sites, and their modulation by certain single-nucleotide variations (SNVs). An order of magnitude higher SCD5A expression was found compared with SCD5B. This unequal splicing is attributed to a weaker recognition of the SCD5B-specific splicing acceptor site, based on predictions confirmed by an optimized minigene assay. The pronounced dominance of SCD5A was largely modified (rs1430176385_A, rs1011850309_A) or even inverted (rs1011850309_C) by natural SNVs at the TV-specific splice sites. Our results provide long missing data on the proportion of SCD5 TVs in human tissues and reveal mutation-driven changes in SCD5 AS, potentially affecting tumor-associated reprogramming of lipid metabolism, thus having prognostic significance, which may be utilized for novel and personalized therapeutic approaches.

Improving Hereditary Hemorrhagic Telangiectasia Molecular Diagnosis: A Referral Center Experience

BACKGROUND

Hereditary hemorrhagic telangiectasia (HHT) is a rare vascular disease inherited in an autosomal dominant manner. Disease-causing variants in endoglin (ENG) and activin A receptor type II-like 1 (ACVRL1) genes are detected in more than 90% of the patients undergoing molecular testing. The identification of variants of unknown significance is often seen as a challenge in clinical practice that makes family screening and genetic counseling difficult. Here, we show that the implementation of cDNA analysis to assess the effect of splice site variants on mRNA splicing is a powerful tool.

METHODS

Gene panel sequencing of genes associated with HHT and other arteriovenous malformation-related syndromes was performed. To evaluate the effect of the splice site variants, cDNA analysis of ENG and ACVRL1 genes was carried out.

RESULTS

three novel splice site variants were identified in ENG (c.68-2A > T and c.1311+4_1311+8del) and ACVLR1 (c.526-6C > G) genes correspondingly in three individuals with HHT that met ≥ 3 Curaçao criteria. All three variants led to an aberrant splicing inducing exon skipping (ENG:c.68-2A > T and ACVRL1:c.526-6C > G) or intron retention (ENG:c.1311+4_1311+8del) allowing the confirmation of the predicted effect on splicing and the reclassification from unknown significance to pathogenic/likely pathogenic of two of them.

CONCLUSIONS

RNA analysis should be performed to assess and/or confirm the impact of variants on splicing. The molecular diagnosis of HHT patients is crucial to allow family screening and accurate genetic counseling. A multidisciplinary approach including clinicians and geneticists is crucial when dealing with patients with rare diseases.

Dysregulated RNA processing and metabolism: a new hallmark of ageing and provocation for cellular senescence

The human genome is capable of producing hundreds of thousands of different proteins and non-coding RNAs from <20 000 genes, in a co-ordinated and regulated fashion. This is achieved by a collection of phenomena known as mRNA processing and metabolism, and encompasses events in the life cycle of an RNA from synthesis to degradation. These factors are critical determinants of cellular adaptability and plasticity, which allows the cell to adjust its transcriptomic output in response to its internal and external environment. Evidence is building that dysfunctional RNA processing and metabolism may be a key contributor to the development of cellular senescence. Senescent cells by definition have exited cell cycle, but have gained functional features such as the secretion of the senescence-associated secretory phenotype (SASP), a known driver of chronic disease and perhaps even ageing itself. In this review, I will outline the impact of dysregulated mRNA processing and metabolism on senescence and ageing at the level of genes, cells and systems, and describe the mechanisms by which progressive deterioration in these processes may impact senescence and organismal ageing. Finally, I will present the evidence implicating this important process as a new hallmark of ageing, which could be harnessed in the future to develop new senotherapeutic interventions for chronic disease.

Genetic variations in OLR1 gene associated with PCOS and atherosclerotic risk factors

Polycystic ovary syndrome (PCOS) is the most common endocrinopathy in women of reproductive age. The aim of this study was to investigate the association of oxidized low-density lipoprotein receptor 1 (OLR1) gene variations with the susceptibility of PCOS and to examine the relationship between the frequencies of OLR1 gene variations and atherosclerotic risk factors. Genomic DNA was extracted from blood samples collected from 49 patients with PCOS and 43 healthy controls. The variants in the OLR1 gene were identified using next-generation sequencing (NGS). Heterozygous rs11053646 (K167N), rs11611438, rs11611453, and rs35688880 genotype frequencies were significantly higher in the PCOS group than that of control group. Single nucleotide polymorphism (SNP) rs34163097 minor A allele increased the PCOS risk by ∼10-fold (p = 0.03). SNPs rs11053646, rs11611438, rs11611453, rs34163097, and rs35688880 were positively correlated with body mass index (BMI). The logistic regression model (area under the curve: 0.770, p = 0.000) further revealed a combination of 2-h plasma glucose (PG-2 h), dehydroepiandrosterone sulfate (DHEAS), and rs11053646 as predictors of PCOS phenotype. This is the first study reporting the NGS data of OLR1 gene variants which might be associated with the pathogenesis of PCOS and several atherosclerotic risk factors, particularly higher BMI and DHEAS. To fully understand the genetic basis of PCOS and the contribution of OLR1 gene variants to PCOS pathogenesis, additional large-scale studies are warranted.

Minigene Assay as an Effective Molecular Diagnostic Strategy in Determining the Pathogenicity of Noncanonical Splice-Site Variants in FLCN

Primary spontaneous pneumothorax (PSP) or pulmonary cyst is one of the manifestations of Birt-Hogg-Dubé syndrome, which is caused by pathogenic variants in FLCN gene. Genetic testing in patients with PSP identifies a certain number of missense or intronic variants. These variants are usually considered as variants of uncertain significance, whose functional interpretations pose a challenge in clinical genetics. To improve recognition of pathogenic splice-altering variants in FLCN gene, computational tools are used to prioritize potential splice-altering variants and then a hybrid minigene assay is performed to verify the RNA splicing pattern. Herein, variants in FLCN exon 11 and its flanking sequence are focused. Eight variants detected in 11 patients with PSP are evaluated, and six variants are prioritized by in silico tools as potential splice-altering variants of uncertain significance. Four variants (c.1177-5_1177-3delCTC, c.1292_1300+4del, c.1300+4C>T, and c.1300+5G>A) are demonstrated by minigene assay to alter RNA splicing of FLCN, and the last three of them are novel. RT-PCR of patient-derived RNA gives consistent results. Genotype-phenotype correlation analysis in patients with PSP with these variants demonstrates good concordance. Our results underline the importance of RNA analysis, which could provide molecular evidence for pathogenicity of a variant, and provide essential information for the clinical interpretation of variants. Combining the clinical information, a definitive diagnosis could be made.

Novel homozygous leptin receptor mutation in an infant with monogenic obesity

Monogenic obesity can be caused by a mutation in one of the single genes involved in hunger and satiety. The most common mutations affect melanocortin 4 (MC4) followed by the leptin gene and its receptor. Leptin receptor (LEPR) gene mutation is an extremely rare endocrine disease characterized by early-onset obesity, hyperphagia in addition to pituitary hormone deficiency, and metabolic abnormalities. We report the case of a 12-month-old male infant born of a non-consanguineous marriage. He presented to us with rapid weight gain from 2 months of age along with hyperphagia. Biochemistry revealed a deranged lipid profile, elevated transaminases, and markedly raised serum leptin levels. On genetic analysis, a novel mutation was detected, which was a homozygous variation In exon 12 of the LEPR gene (chr1:g.65608901G>A) that resulted in the synonymous amino acid change of lysine at codon 584 proximal to donor splice site (p.Lys584). The in silico prediction of the variant was 'damaging' by MutationTaster2. The mutation was classified as a 'variant of uncertain significance' due to a lack of published literature and had to be correlated carefully with the clinical symptoms. It was recommended to do Sanger sequencing of the parents and other family members. However, due to financial constraints, the family could not afford the same. At the time of writing, funds were being arranged for procuring setmelanotide, which is a novel and effective therapy for monogenic obesity due to LepR mutation.

CD5 B-Cell Predominant Primary Immunodeficiency: Part of the Spectrum of MAGT1 Deficiency

Background

Selective anti-polysaccharide antibody deficiency (SPAD) with CD5 B-cell predominance and autoimmune phenomena was identified in a male cohort first reported by Antall et al in 1999. The phenotypically likewise and genotypically identical X-linked immunodeficiency with magnesium defect, Epstein-Barr Virus infection, and neoplasia (XMEN) disease was defined as a novel primary immunodeficiency (PID) in 2011. Recent studies of the magnesium transporter 1 (MAGT1) gene mutation reveal glycosylation defects contributing to more phenotypic variance than the "XMEN" title pathologies. The updated title, "X-linked MAGT1 deficiency with increased susceptibility to EBV-infection and N-linked glycosylation defect," was proposed in 2020.

Objectives

To reflect the patient population more accurately, a prospective classification update may consider MAGT1 glycobiological errors contributing to phenotypic variance but also pre-genetic testing era reports with CD5 B-cell predominance.

Methods

Patient 1 from Antall et al presented at 28 years of age for further immunological evaluation of his CD5/CD19 B-cell predominance diagnosed at 5 years old.

Design

Immune re-evaluation done through flow cytometry and next-generation sequencing.

Results

Flow cytometry B-cell phenotyping revealed persistent CD5+CD19+ (93%). Flow cytometric histogram quantified reduced activator CD16+CD56+ natural killer and CD8+ T-cell receptor, Group 2, Member D (NKG2D) glycoprotein expression. A c.923-1_934 deletion loss of function mutation was identified in the MAGT1 gene.

Conclusion

We suggest the novel PID XMEN, based on its CD5 B-cell predominance, had been discovered and reported over a decade earlier as CD5+ PID based on the MAGT1 mutation found in the same. We encourage consideration of combining these labels and recent findings to offer the most accurate classification of this disease.

Clinical and genetic characteristics of Tunisian children with infantile nephropathic cystinosis

BACKGROUND

Nephropathic cystinosis is an autosomal recessive disease caused by a mutation in the CTNS gene which encodes cystinosin, a lysosomal cystine transporter. The spectrum of mutations in the CTNS gene is not well defined in the North African population. Here, we investigated twelve patients with nephropathic cystinosis belonging to eight Tunisian families in order to analyze the clinical and genetic characteristics of Tunisian children with infantile nephropathic cystinosis.

METHODS

Clinical data were collected retrospectively. Molecular analysis of the CTNS gene was performed by Sanger sequencing.

RESULTS

We describe a new splicing mutation c.971-1G > C in the homozygous state in 6/12 patients which seems to be a founder mutation. The reported deletion of 23nt c.771_793 Del (p.Gly258Serfs*30) was detected in a homozygous state in one patient and in a heterozygous compound state with the c.971-1G > C mutation in 3/12 patients. Two of 12 patients have a deletion of exons 4 and 5 of the CTNS gene. None of our patients had the most common 57-kb deletion.

CONCLUSIONS

The mutational spectrum in the Tunisian population is different from previously described populations. Thus, a molecular diagnostic strategy must be implemented in Tunisia, by targeting as a priority the common mutations described in this country. Such a strategy will allow a cost-effective diagnosis confirmation as well as early administration of treatment with oral cysteamine. A higher resolution version of the Graphical abstract is available as Supplementary information.

The regulatory role of alternative splicing in inflammatory bowel disease

Inflammatory bowel disease (IBD) mainly includes Crohn's disease and ulcerative colitis. These diseases have a progressive course of chronic relapse and remission and affect a large number of children and adults worldwide. The burden of IBD is rising worldwide, with levels and trends varying greatly in countries and regions. Like most chronic diseases, the costs associated with IBD are high, including hospitalizations, outpatient and emergency visits, surgeries, and pharmacotherapies. However, there is no radical cure for it yet, and its therapeutic targets still need further study. Currently, the pathogenesis of IBD remains unclear. It is generally assumed that the occurrence and development of IBD are related to the environmental factors, gut microbiota, immune imbalance, and genetic susceptibility. Alternative splicing contributes to a various diseases, such as spinal muscular atrophy, liver diseases, and cancers. In the past, it has been reported that alternative splicing events, splicing factors, and splicing mutations were associated with IBD, but there were no reports on the practical application for clinical diagnosis and treatment of IBD using splicing-related methods. Therefore, this article reviews research progress on alternative splicing events, splicing factors, and splicing mutations associated with IBD.

The multi-faceted nature of 15 CFTR exonic variations: Impact on their functional classification and perspectives for therapy

BACKGROUND

The majority of variants of unknown clinical significance (VUCS) in the CFTR gene are missense variants. While change on the CFTR protein structure or function is often suspected, impact on splicing may be neglected. Such undetected splicing default of variants may complicate the interpretation of genetic analyses and the use of an appropriate pharmacotherapy.

METHODS

We selected 15 variants suspected to impact CFTR splicing after in silico predictions on 319 missense variants (214 VUCS), reported in the CFTR-France database. Six specialized laboratories assessed the impact of nucleotide substitutions on splicing (minigenes), mRNA expression levels (quantitative PCR), synthesis and maturation (western blot), cellular localization (immunofluorescence) and channel function (patch clamp) of the CFTR protein. We also studied maturation and function of the truncated protein, consecutive to in-frame aberrant splicing, on additional plasmid constructs.

RESULTS

Six of the 15 variants had a major impact on CFTR splicing by in-frame (n = 3) or out-of-frame (n = 3) exon skipping. We reclassified variants into: splicing variants; variants causing a splicing defect and the impairment of CFTR folding and/or function related to the amino acid substitution; deleterious missense variants that impair CFTR folding and/or function; and variants with no consequence on the different processes tested.

CONCLUSION

The 15 variants have been reclassified by our comprehensive approach of in vitro experiments that should be used to properly interpret very rare exonic variants of the CFTR gene. Targeted therapies may thus be adapted to the molecular defects regarding the results of laboratory experiments.

Genetic association between the PTPN22, IRF5 and TYK2 gene variants and susceptibility to juvenile idiopathic arthritis

Juvenile idiopathic arthritis (JIA) refers to a group of chronic childhood arthropathies of unknown etiology. In the present study, the genetic association between the variants in PTPN22, IRF5 and TYK2 genes and susceptibility to JIA was investigated. The distributions of 16 variants in PTPN22, IRF5 and TYK2 genes were analyzed by direct sequencing in 378 patients with JIA and 378 healthy controls. Odds ratios and 95% confidence intervals were used to evaluate the association between the gene variants and JIA. The gene-gene interactions were investigated using multifactor dimensionality reduction. All allelic and dominant models of PTPN22 rs1214414, rs1214418, rs1746853, rs3765598 and rs3811021 were significantly associated with JIA risk (P<0.05). IRF5 rs10954213 in both allelic and dominant models, as well as the allelic model of rs2004640, was significantly related to JIA risk (P<0.05). In addition, the allelic, recessive and dominant models of TYK2 rs280500, rs280519, rs2304256 and rs12720270 were significantly related to JIA risk (P<0.05). In addition, three haplotypes (H_{C A G T C C}, H_{C A G T T C} and H_{C G T T C T} ) in PTPN22 gene, three haplotypes (H_{D T A A}, H_{I T A C} and H_{D T G C}) in IRF5 gene and two haplotypes (H_{A G G A T} and H_{G A G G T}) in TYK2 gene were associated with the risk of JIA (P<0.05). Furthermore, a three-way interaction between IRF5 rs10954213, rs2004640 and PTPN22 rs1214414 was shown to be associated with JIA risk. In conclusion, PTPN22 rs1214418, rs1746853, rs3765598, IRF5 rs2004640, TYK2 rs280500, rs2304256 and a three-way interaction between IRF5 rs10954213, rs2004640 and PTPN22 rs1214414 may be risk factors for JIA.

Adenosine Deaminase 2 Deficiency Caused by Biallele Variants Including Splicing Variant: The First Case in Korea

Deficiency of adenosine deaminase 2 (DADA2) is an autoinflammatory disease caused by pathogenic variants of the ADA2 gene and has similar clinical features to polyarteritis nodosa (PAN). We, herein, report a case of DADA2 in Korea that was diagnosed in a patient with childhood-onset PAN. The patient had a truncal ataxia and facial palsy caused by thalamic infarction at 34 months of age. Livedo reticularis with Raynaud phenomenon and abdominal pain with fever were followed. Radiologic examination showed multiple infarctions in brain and kidney. She was diagnosed with PAN using skin biopsy and angiography. She had severe hemorrhagic strokes despite medical treatments. Her disease activity was controlled after adding a tumor necrosis factor-α inhibitor. Molecular analysis revealed compound heterozygous pathogenic variants of ADA2 gene. This is the first case of DADA2 in Korea. Genetic analysis for ADA2 gene should be considered in patients with childhood-onset PAN.

Integrative Functional Genomic Analysis in Multiplex Autism Families from Kazakhstan

The study of extended pedigrees containing autism spectrum disorder- (ASD-) related broader autism phenotypes (BAP) offers a promising approach to the search for ASD candidate variants. Here, a total of 650,000 genetic markers were tested in four Kazakhstani multiplex families with ASD and BAP to obtain data on de novo mutations (DNMs), common, and rare inherited variants that may contribute to the genetic risk for developing autistic traits. The variants were analyzed in the context of gene networks and pathways. Several previously well-described enriched pathways were identified, including ion channel activity, regulation of synaptic function, and membrane depolarization. Perhaps these pathways are crucial not only for the development of ASD but also for ВАР. The results also point to several additional biological pathways (circadian entrainment, NCAM and BTN family interactions, and interaction between L1 and Ankyrins) and hub genes (CFTR, NOD2, PPP2R2B, and TTR). The obtained results suggest that further exploration of PPI networks combining ASD and BAP risk genes can be used to identify novel or overlooked ASD molecular mechanisms.

Frequency of actionable Exomic secondary findings in 160 Colombian patients: Impact in the healthcare system

INTRODUCTION

By 2021, the American College of Medical Genetics and Genomics (ACMG) published the last version of their secondary findings (SF) reporting recommendations for cases in which a person receives a genetic test.

OBJECTIVE

To determine in a sample of the Colombian population the prevalence of SF for the 59 genes on the ACMG SF v2.0 list associated with 27 genetic diseases.

MATERIALS AND METHODS

An analytical cross-sectional study was developed by examining the sequences of 160 exomes. Based on the ACMG guidelines, a variant classification algorithm was designed to filter and select reportable SF.

RESULTS

Eleven pathogenic variants were identified in 13/160 (8.13%) patients in genes APOB, BRCA2, CACNA1S, COL3A1, LDLR, MYBPC3, PCSK9, PKP2, PMS2 and RYR2. No association was found between the sociodemographic variables and the SF to report (P > 0,05).

CONCLUSION

We reported the first approach of actionable pathogenic variants spectrum in the Colombian population. Given the frequency found in this study and the clinical impact of genomic variants on health, it is essential to actively search for SF having the opportunity to receive genetic counselling, prevention and clinical management.

Associations between liver X receptor polymorphisms and blood lipids: A systematic review and meta-analysis

Genetic susceptibility to dyslipidaemia remains incompletely understood. The liver X receptors (LXRs), members of the nuclear receptor superfamily of ligand dependent transcription factors, are homeostatic regulators of lipid metabolism. Multiple single nucleotide polymorphisms (SNPs)have been identified previously in the coding and regulatory regions of the LXRs. The aim of this systematic review and meta-analysis was to summarise associations between SNPs of LXRs (α and β isoforms) with blood lipid and lipoprotein traits. Five databases (PubMed, Ovid Embase, Scopus, Web of Science, and the Cochrane Library) were systematically searched for population-based studies that assessed associations between one or more blood lipid/lipoprotein traits and LXR SNPs. Of seventeen articles included in the qualitative synthesis, ten were eligible for meta-analysis. Nine LXRα SNPs and five LXRβ SNPs were identified, and the three most studied LXRα SNPs were quantitatively summarised. Carriers of the minor allele A of LXRα rs12221497 (-115G>A) had higher triglyceride levels than GG homozygotes (0.13 mmol/L; 95%CI: [0.03, 0.23], P = 0.01). Heterozygote carriers of LXRα rs2279238 (297C/T) had higher total cholesterol levels (0.12 mmol/L; (95%CI: [0.01, 0.23], P = 0.04) than either CC or TT homozygotes. For LXRα rs11039155 (-6G>A), no significant differences in blood levels of either triglyceride (P = 0.39) or HDL-C (P = 0.98) were detected between genotypes in meta-analyses. In addition, there were no strong associations for other SNPs of LXRα and LXRβ. This study provides the evidence of an association between LXRα, but not LXRβ, SNPs and blood-lipid traits. Systematic review registration: PROSPERO No. CRD42021246158.