ESEfinder: A web resource to identify exonic splicing enhancers

Whole-transcriptome sequencing-based concomitant detection of viral and human genetic determinants of cutaneous lesions

Severe viral infections of the skin can occur in patients with inborn errors of immunity (IEI). We report an all-in-one whole-transcriptome sequencing–based method by RNA-Seq on a single skin biopsy for concomitantly identifying the cutaneous virome and the underlying IEI. Skin biopsies were obtained from healthy and lesional skin from patients with cutaneous infections suspected to be of viral origin. RNA-Seq was utilized as the first-tier strategy for unbiased human genome-wide rare variant detection. Reads unaligned to the human genome were utilized for the exploration of 926 viruses in a viral genome catalog. In 9 families studied, the patients carried pathogenic variants in 6 human IEI genes, including IL2RG, WAS, CIB1, STK4, GATA2, and DOCK8. Gene expression profiling also confirmed pathogenicity of the human variants and permitted genome-wide homozygosity mapping, which assisted in identification of candidate genes in consanguineous families. This automated, online, all-in-one computational pipeline, called VirPy, enables simultaneous detection of the viral triggers and the human genetic variants underlying skin lesions in patients with suspected IEI and viral dermatosis.

Silent mutations reveal therapeutic vulnerability in RAS Q61 cancers

RAS family members are the most frequently mutated oncogenes in human cancers. Although KRAS(G12C)-specific inhibitors show clinical activity in patients with cancer^1-3, there are no direct inhibitors of NRAS, HRAS or non-G12C KRAS variants. Here we uncover the requirement of the silent KRAS^G60G mutation for cells to produce a functional KRAS(Q61K). In the absence of this G60G mutation in KRAS^Q61K, a cryptic splice donor site is formed, promoting alternative splicing and premature protein termination. A G60G silent mutation eliminates the splice donor site, yielding a functional KRAS(Q61K) variant. We detected a concordance of KRAS^Q61K and a G60G/A59A silent mutation in three independent pan-cancer cohorts. The region around RAS Q61 is enriched in exonic splicing enhancer (ESE) motifs and we designed mutant-specific oligonucleotides to interfere with ESE-mediated splicing, rendering the RAS(Q61) protein non-functional in a mutant-selective manner. The induction of aberrant splicing by antisense oligonucleotides demonstrated therapeutic effects in vitro and in vivo. By studying the splicing necessary for a functional KRAS(Q61K), we uncover a mutant-selective treatment strategy for RAS^Q61 cancer and expose a mutant-specific vulnerability, which could potentially be exploited for therapy in other genetic contexts.

Novel Intronic Mutations of TBK1 Promote Aberrant Splicing Modes in Amyotrophic Lateral Sclerosis

TANK-binding kinase 1 (TBK1) has been identified as a causative gene of amyotrophic lateral sclerosis (ALS) in the Caucasian population in 2015. Here, we sequenced for TBK1 variants in a cohort of 15 familial ALS (fALS) and 275 sporadic ALS (sALS) of Chinese origin by targeted next-generation sequencing. We identified one likely benign missense variant (p. Ser398Pro), two missense variants of uncertain significance (p. Ile37Leu and p. Tyr677Asn), and two novel heterozygous variants in introns of TBK1, c.1522-3T > G and c.2066 + 4A > G. We performed splicing assays through minigene plasmids and RNA pull-down assay to determine that the two substitutions of nucleotides disrupted the binding of the important splicing regulator hnRNPA1 and promoted aberrant pre-mRNA splicing modes. The c.1522-3T > G variant promoted nearly 50.0% of abnormal transcripts (3 different types of insertions and deletions (indels) in junction of intron 13-exon 14) and the c.2066 + 4A > G variant inhibited about 75.0% inclusion of exon 19, both causing premature stop codon and producing TBK1 protein without CCD2. Immunofluorescence analysis showed that the expression of TBK1 with intronic variants was lower since less TBK1 distribution was observed in HEK293T cells. Both patients carrying TBK1 c.1522-3T > G and c.2066 + 4A > G variants developed a rapidly progressive ALS, with a survival of 31 and 10 months, respectively. The frequency of loss of function (LoF) variants in TBK1 was 0.73% in sALS in our cohort. We emphasize that intronic sequencing and pre-mRNA splicing analysis cannot be ignored to demonstrate the complex mutational spectrum and pathogenesis of ALS.

Exon Skipping Through Chimeric Antisense U1 snRNAs to Correct Retinitis Pigmentosa GTPase-Regulator (RPGR) Splice Defect

Inherited retinal dystrophies are caused by mutations in more than 250 genes, each of them carrying several types of mutations that can lead to different clinical phenotypes. Mutations in Retinitis Pigmentosa GTPase-Regulator (RPGR) cause X-linked Retinitis pigmentosa (RP). A nucleotide substitution in intron 9 of RPGR causes the increase of an alternatively spliced isoform of the mature mRNA, bearing exon 9a (E9a). This introduces a stop codon, leading to truncation of the protein. Aiming at restoring impaired gene expression, we developed an antisense RNA-based therapeutic approach for the skipping of RPGR E9a. We designed a set of specific U1 antisense snRNAs (U1_asRNAs) and tested their efficacy in vitro, upon transient cotransfection with RPGR minigene reporter systems in HEK-293T, 661W, and PC-12 cell lines.

Cytoprotective Activity of Polyamines Is Associated with the Alternative Splicing of RAD51A Pre-mRNA in Normal Human CD4+ T Lymphocytes

Physiological polyamines are ubiquitous polycations with pleiotropic biochemical activities, including regulation of gene expression and cell proliferation as well as modulation of cell signaling. They can also decrease DNA damage and promote cell survival. In the present study, we demonstrated that polyamines have cytoprotective effects on normal human CD4+ T lymphocytes but not on cancer Jurkat or K562 cells. Pretreatment of lymphocytes with polyamines resulted in a significant reduction in cells with DNA damage induced by doxorubicin, cisplatin, or irinotecan, leading to an increase in cell survival and viability. The induction of RAD51A expression was in response to DNA damage in both cancer and normal cells. However, in normal cells, putrescin pretreatment resulted in alternative splicing of RAD51A and the switch of the predominant expression from the splice variant with the deletion of exon 4 to the full-length variant. Induction of RAD51A alternative splicing by splice-switching oligonucleotides resulted in a decrease in DNA damage and cell protection against cisplatin-induced apoptosis. The results of this study suggest that the cytoprotective activity of polyamines is associated with the alternative splicing of RAD51A pre-mRNA in normal human CD4+ T lymphocytes. The difference in the sensitivity of normal and cancer cells to polyamines may become the basis for the use of these compounds to protect normal lymphocytes during lymphoblastic chemotherapy.

PDE4B Proposed as a High Myopia Susceptibility Gene in Chinese Population

Myopia is the most common cause of refractive error worldwide. High myopia is a severe type of myopia, which usually accompanies pathological changes in the fundus. To identify high myopia susceptibility genes, DNA-pooling based genome-wide association analysis was used to search for a correlation between single nucleotide polymorphisms and high myopia in a Han Chinese cohort (cases vs. controls in discovery stage: 507 vs. 294; replication stage 1: 991 vs. 1,025; replication stage 2: 1,021 vs. 52,708). Three variants (rs10889602T/G, rs2193015T/C, rs9676191A/C) were identified as being significantly associated with high myopia in the discovery, and replication stage. rs10889602T/G is located at the third intron of phosphodiesterase 4B (PDE4B), whose functional assays were performed by comparing the effects of rs10889602T/T deletion of this risk allele on PDE4B and COL1A1 gene and protein expression levels in the rs10889602T/T^del/del, rs10889602T/T^del/wt, and normal control A549 cell lines. The declines in the PDE4B and COL1A1 gene expression levels were larger in the rs10889602T/T deleted A549 cells than in the normal control A549 cells (one-way ANOVA, p < 0.001). The knockdown of PDE4B by siRNA in human scleral fibroblasts led to downregulation of COL1A1. This correspondence between the declines in rs10889602 of the PDE4B gene, PDE4B knockdown, and COL1A1 protein expression levels suggest that PDE4B may be a novel high myopia susceptibility gene, which regulates myopia progression through controlling scleral collagen I expression levels. More studies are needed to determine if there is a correlation between PDE4B and high myopia in other larger sample sized cohorts.

Genetic Control of Splicing at SIRPG Modulates Risk of Type 1 Diabetes

Signal regulatory protein SIRPγ (CD172G) is expressed on the surface of lymphocytes, where it acts by engaging its ligand, CD47. SIRPG, which encodes SIRPγ, contains a nonsynonymous coding variant, rs6043409, which is significantly associated with risk for type 1 diabetes. SIRPG produces multiple transcript isoforms via alternative splicing, all encoding potentially functional proteins. We show that rs6043409 alters a predicted exonic splicing enhancer, resulting in significant shifts in the distribution of SIRPG transcript isoforms. All of these transcript isoforms produced protein upon transient expression in vitro. However, CRISPR/Cas9 targeting of one of the alternatively spliced exons in SIRPG eliminated all SIRPγ expression in Jurkat T cells. These targeted cells formed fewer cell-cell conjugates with each other than with wild-type Jurkat cells, expressed reduced levels of genes associated with CD47 signaling, and had significantly increased levels of cell-surface CD47. In primary CD4+ and CD8+ T cells, cell-surface SIRPγ levels in response to anti-CD3 stimulation varied quantitatively by rs6043409 genotype. Our results suggest that SIRPG is the most likely causative gene for type 1 diabetes risk in the 20p13 region and highlight the role of alternative splicing in lymphocytes in mediating the genetic risk for autoimmunity.

Oligonucleotide correction of an intronic TIMMDC1 variant in cells of patients with severe neurodegenerative disorder

TIMMDC1 encodes the Translocase of Inner Mitochondrial Membrane Domain-Containing protein 1 (TIMMDC1) subunit of complex I of the electron transport chain responsible for ATP production. We studied a consanguineous family with two affected children, now deceased, who presented with failure to thrive in the early postnatal period, poor feeding, hypotonia, peripheral neuropathy and drug-resistant epilepsy. Genome sequencing data revealed a known, deep intronic pathogenic variant TIMMDC1 c.597-1340A>G, also present in gnomAD (~1/5000 frequency), that enhances aberrant splicing. Using RNA and protein analysis we show almost complete loss of TIMMDC1 protein and compromised mitochondrial complex I function. We have designed and applied two different splice-switching antisense oligonucleotides (SSO) to restore normal TIMMDC1 mRNA processing and protein levels in patients' cells. Quantitative proteomics and real-time metabolic analysis of mitochondrial function on patient fibroblasts treated with SSOs showed restoration of complex I subunit abundance and function. SSO-mediated therapy of this inevitably fatal TIMMDC1 neurologic disorder is an attractive possibility.

Open reading frame correction using splice-switching antisense oligonucleotides for the treatment of cystic fibrosis

CFTR gene mutations that result in the introduction of premature termination codons (PTCs) are common in cystic fibrosis (CF). This mutation type causes a severe form of the disease, likely because of low CFTR messenger RNA (mRNA) expression as a result of nonsense-mediated mRNA decay, as well as the production of a nonfunctional, truncated CFTR protein. Current therapeutics for CF, which target residual protein function, are less effective in patients with these types of mutations due in part to low CFTR protein levels. Splice-switching antisense oligonucleotides (ASOs), designed to induce skipping of exons in order to restore the mRNA open reading frame, have shown therapeutic promise preclinically and clinically for a number of diseases. We hypothesized that ASO-mediated skipping of CFTR exon 23 would recover CFTR activity associated with terminating mutations in the exon, including CFTR p.W1282X, the fifth most common mutation in CF. Here, we show that CFTR lacking the amino acids encoding exon 23 is partially functional and responsive to corrector and modulator drugs currently in clinical use. ASO-induced exon 23 skipping rescued CFTR expression and chloride current in primary human bronchial epithelial cells isolated from a homozygote CFTR-W1282X patient. These results support the use of ASOs in treating CF patients with CFTR class I mutations in exon 23 that result in unstable CFTR mRNA and truncations of the CFTR protein.

Exon-skipping antisense oligonucleotides for cystic fibrosis therapy

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), and the CFTR-W1282X nonsense mutation causes a severe form of CF. Although Trikafta and other CFTR-modulation therapies benefit most CF patients, targeted therapy for patients with the W1282X mutation is lacking. The CFTR-W1282X protein has residual activity but is expressed at a very low level due to nonsense-mediated messenger RNA (mRNA) decay (NMD). NMD-suppression therapy and read-through therapy are actively being researched for CFTR nonsense mutants. NMD suppression could increase the mutant CFTR mRNA, and read-through therapies may increase the levels of full-length CFTR protein. However, these approaches have limitations and potential side effects: because the NMD machinery also regulates the expression of many normal mRNAs, broad inhibition of the pathway is not desirable, and read-through drugs are inefficient partly because the mutant mRNA template is subject to NMD. To bypass these issues, we pursued an exon-skipping antisense oligonucleotide (ASO) strategy to achieve gene-specific NMD evasion. A cocktail of two splice-site-targeting ASOs induced the expression of CFTR mRNA without the premature-termination-codon-containing exon 23 (CFTR-Δex23), which is an in-frame exon. Treatment of human bronchial epithelial cells with this cocktail of ASOs that target the splice sites flanking exon 23 results in efficient skipping of exon 23 and an increase in CFTR-Δex23 protein. The splice-switching ASO cocktail increases the CFTR-mediated chloride current in human bronchial epithelial cells. Our results set the stage for developing an allele-specific therapy for CF caused by the W1282X mutation.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

Reclassification of a TMC1 synonymous substitution as a variant disrupting splicing regulatory elements associated with recessive hearing loss

Alterations of the transmembrane channel-like 1 gene (TMC1) are involved in autosomal recessive and dominant nonsyndromic hearing loss (NSHL). To date, up to 117 causal variants including substitutions, insertions and splice variants have been reported in families from different populations. In a patient suffering from severe prelingual NSHL, we identified, in the homozygous state, the previously considered likely benign synonymous c.627C>T; p.(Leu209=) substitution. We used in silico tools predicting variant-induced alterations of splicing regulatory elements (SREs) and pinpointed this transition as a candidate splice-altering variation. Functional splicing analysis, using a minigene assay, confirmed that the variant altered a critical regulatory sequence which is essential for the exon 11 inclusion in the TMC1 transcripts. This result was reinforced by the analysis of orthologous TMC1 mammalian sequences for which the deleterious effect on the mRNA processing of a native thymidine was always counteracted by the presence of a stronger donor splice site or additional enhancer motifs. This study demonstrates, for the first time, the pathogenicity of the c.627C>T alteration leading to its reclassification as a causal variant impacting SREs and highlights the major importance of exhaustive studies to accurately evaluate the pathogenicity of a variant, regardless of the variation type.

OUP accepted manuscript

Alternative splicing is frequently involved in the diversification of protein function and can also be modulated for therapeutic purposes. Here we develop a predictive model, called Exon ByPASS (predicting Exon skipping Based on Protein amino acid SequenceS), to assess the criticality of exon inclusion based solely on information contained in the amino acid sequence upstream and downstream of the exon junctions. By focusing on protein sequence, Exon ByPASS predicts exon skipping independent of tissue and species in the absence of any intronic information. We validate model predictions using transcriptomic and proteomic data and show that the model can capture exon skipping in different tissues and species. Additionally, we reveal potential therapeutic opportunities by predicting synthetically skippable exons and neo-junctions arising in cancer cells.

Silent variant in F8:c.222G>T (p.Thr74Thr) causes a partial exon skipping in a patient with mild hemophilia A

One of the challenges of genetic testing in patients with hemophilia A is the interpretation of sequence variants. Here we report a silent variant found in exon 2 in the F8 gene in a 47‐year‐old patient with a previous von Willebrand disease (VWD) type 1 diagnosis. Clinically he had mild bleeding symptoms restricted to prolonged bleeding from minor wounds. Sanger sequencing of F8 gene using genomic DNA showed a hemizygous silent variant in exon 2: c.222G>T, p.Thr74Thr. When applying ACMG criteria, the variant was predicted to be “likely benign” in the analyzing software or VUS after curating. Sanger sequencing of the patient's cDNA after nested polymerase chain reaction showed that the patient had both a normal transcript containing exons 1–4 and a defect transcript lacking exon 2. These findings explain the patient's low FVIII:C level and led to the diagnosis of mild hemophilia A instead of VWD type 1. This case illustrates that mRNA work‐up may be needed to clarify a patient's phenotype–genotype.

Genetic association of ARID5B with the risk of colorectal cancer within Jammu and Kashmir, India

Colorectal cancer (CRC), which includes the development of cancer from the colon or rectum, is one of the highly prevalent cancers in the populations of Jammu and Kashmir (J&K) in India. However, case-control genetic association studies on CRC are lacking in this population. Various genome-wide association studies have previously shown that single-nucleotide polymorphisms (SNPs) of the AT-rich interaction domain 5B (ARID5B) gene located on chromosome 10q21.2 contribute substantially to the development of colorectal cancer. The association between ARID5B and CRC risk in north Indian population groups is still unknown. To understand the role of ARID5B SNPs in CRC in the population of J&K, we designed a case-control study to investigate the association of the cancer susceptibility variant rs10740055 of ARID5B with CRC in the population of J&K. The study included 180 cases and 390 healthy controls. Genotyping of the rs10740055 variant was performed by RT-PCR using the TaqMan assay technique. Hardy-Weinberg equilibrium of the variant was assessed using the chi-squared test. The allele- and genotype-specific risks were estimated by odds ratios (ORs) with 95% confidence intervals (CIs). The rs10740055 variant showed a higher risk for colorectal cancer with an OR of 3.35 (1.99-5.65 at 95% CI) and P = 0.000005 corrected for age, gender, ethnicity, BMI, alcohol intake and smoking. Our results indicate that the A allele of rs10740055 imparts risk to the population and also that a larger sample size is needed for further statistical validation. The association of other variants in other ARID family genes should also be tested as their role cannot be ruled out.

Associations of MCM8 rs3761873 and rs16991617 variants with abnormal uterine bleeding induced by copper intrauterine device

AIM

Intrauterine device (IUD) is a commonly used contraceptive method worldwide. Abnormal uterine bleeding (AUB) is one of the most common side effects of Cu-IUDs. Since AUB varies among Cu-IUD users, changes in the bleeding-related genetic factors may contribute to AUB. This study aimed to determine the genetic risk factors of AUB after Cu-IUD insertion.

METHODS

We conducted a case-control study on women who experienced AUB after Cu-IUD insertion (case:control = 62:59). Six candidate variants were genotyped using the Sequenom MassARRAY. Genotype and allele frequencies were analyzed using SHEsisPlus. We performed Pearson's Chi-squared test to analyze categorical data, and ESEfinder to predict the impact on splicing regulation.

RESULTS

MCM8 coding sequence variants: rs3761873-A>C was in Exon 7 and rs16991617 A>G was in Exon 12 of all 19 exons, both of which were significantly different between cases and controls (p_allele  = 0.039 and p_genotype  = 0.092). rs6022 and rs6029 in F5 gene and rs3761873 and rs16991617 in the MCM8 gene showed strong linkage disequilibrium (R²  > 0.8). ESEfinder indicated that the variants of MCM8 may affect the splicing regulation.

CONCLUSIONS

MCM8 rs376187 and rs16991617 were associated with AUB in Cu-IUDs users. MCM8 may play a role in AUB by regulating functions of reproductive organs and primary ovarian insufficiency. Our findings may improve the understanding of the genetic basis of AUB caused by Cu-IUDs.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Cas9-specific immune responses compromise local and systemic AAV CRISPR therapy in multiple dystrophic canine models

Adeno-associated virus (AAV)-mediated CRISPR-Cas9 editing holds promise to treat many diseases. The immune response to bacterial-derived Cas9 has been speculated as a hurdle for AAV-CRISPR therapy. However, immunological consequences of AAV-mediated Cas9 expression have thus far not been thoroughly investigated in large mammals. We evaluate Cas9-specific immune responses in canine models of Duchenne muscular dystrophy (DMD) following intramuscular and intravenous AAV-CRISPR therapy. Treatment results initially in robust dystrophin restoration in affected dogs but also induces muscle inflammation, and Cas9-specific humoral and cytotoxic T-lymphocyte (CTL) responses that are not prevented by the muscle-specific promoter and transient prednisolone immune suppression. In normal dogs, AAV-mediated Cas9 expression induces similar, though milder, immune responses. In contrast, other therapeutic (micro-dystrophin and SERCA2a) and reporter (alkaline phosphatase, AP) vectors result in persistent expression without inducing muscle inflammation. Our results suggest Cas9 immunity may represent a critical barrier for AAV-CRISPR therapy in large mammals.

Six splice site variations, three of them novel, in the ABO gene occurring in nine individuals with ABO subtypes

Background

Nucleotide mutations in the ABO gene may reduce the activity of glycosyltransferase, resulting in lower levels of A or B antigen expression in red blood cells. Six known splice sites have been identified according to the database of red cell immunogenetics and the blood group terminology of the International Society of Blood Transfusion. Here, we describe six distinct splice site variants in individuals with ABO subtypes.

Methods

The ABO phenotype was examined using a conventional serological method. A polymerase chain reaction sequence-based typing method was used to examine the whole coding sequence of the ABO gene. The ABO gene haplotypes were studied using allele-specific primer amplification or cloning technology. In silico analytic tools were used to assess the functional effect of splice site variations.

Results

Six distinct variants in the ABO gene splice sites were identified in nine individuals with ABO subtypes, including c.28 + 1_2delGT, c.28 + 5G > A, c.28 + 5G > C, c.155 + 5G > A, c.204-1G > A and c.374 + 5G > A. c.28 + 1_2delGT was detected in an A_w individual, while c.28 + 5G > A, c.28 + 5G > C, and c.204-1G > A were detected in B_el individuals. c.155 + 5G > A was detected in one B₃ and two AB₃ individuals, whereas c.374 + 5G > A was identified in two A_el individuals. Three novel splice site variants (c.28 + 1_2delGT, c.28 + 5G > A and c.28 + 5G > C) in the ABO gene were discovered, all of which resulted in low antigen expression. In silico analysis revealed that all variants had the potential to alter splice transcripts.

Conclusions

Three novel splice site variations in the ABO gene were identified in Chinese individuals, resulting in decreased A or B antigen expression and the formation of ABO subtypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03141-5.

Functional pre-therapeutic evaluation by genome editing of variants of uncertain significance of essential tumor suppressor genes

BACKGROUND

Targeted therapies in oncology are promising but variants of uncertain significance (VUS) limit their use for clinical management and necessitate functional testing in vitro. Using BRCA1 and BRCA2 variants, which have consequences on PARP inhibitor sensitivity, and POLE variants, potential biomarkers of immunotherapy response, we developed a rapid functional assay based on CRISPR-Cas9 genome editing to determine the functional consequences of these variants having potentially direct implications on patients' access to targeted therapies.

METHODS

We first evaluated the functional impact of 26 BRCA1 and 7 BRCA2 variants by editing and comparing NGS results between the variant of interest and a silent control variant. Ten of these variants had already been classified as benign or pathogenic and were used as controls. Finally, we extended this method to the characterization of POLE VUS.

RESULTS

For the 23 variants that were unclassified or for which conflicting interpretations had been reported, 15 were classified as functionally normal and 6 as functionally abnormal. Another two variants were found to have intermediate consequences, both with potential impacts on splicing. We then compared these scores to the patients' responses to PARP inhibitors when possible. Finally, to prove the application of our method to the classification of variants from other tumor suppressor genes, we exemplified with three POLE VUS. Among them, two were classified with an intermediate functional impact and one was functionally abnormal. Eventually, four POLE variants previously classified in databases were also evaluated. However, we found evidence of a discordance with the classification, variant p.Leu424Val being found here functionally normal.

CONCLUSIONS

Our new rapid functional assay can be used to characterize the functional implication of BRCA1 and BRCA2 variants, giving patients whose variants were evaluated as functionally abnormal access to PARP inhibitor treatment. Retrospective analysis of patients' responses to PARP inhibitors, where accessible, was consistent with our functional score evaluation and confirmed the accuracy of our protocol. This method could potentially be extended to the classification of VUS from all essential tumor suppressor genes and can be performed within a timeframe compatible with clinical applications, thereby having a direct theranostic impact.

Identification of novel deep intronic PAH gene variants in patients diagnosed with phenylketonuria

Phenylketonuria (PKU) is caused by phenylalanine hydroxylase (PAH) gene variants. Previously, 94.21% of variants were identified using Sanger sequencing and multiplex ligation-dependent probe amplification. To investigate the remaining variants, we performed whole-genome sequencing for four patients with PKU and unknown genotypes to identify deep intronic or structural variants. We identified three novel heterozygous variants (c.706+368T>C, c.1065+241C>A, and c.1199+502A>T) in a deep PAH gene intron. We detected a c.1199+502A>T variant in 60% (6/10) of PKU patients with genetically undetermined PKU. In silico predictions indicated that the three deep variants may impact splice site selection and result in the inclusion of a pseudo-exon. A c.1199+502A>T PAH minigene and reverse transcription PCR (RT-PCR) on blood RNA from a PKU patient with biallelic variants c.1199+502A>T and c.1199G>A confirmed that the c.1199+502A>T variant may strengthen the predicted branch point and leads to the inclusion of a 25-nt pseudo-exon in the PAH mRNA. Reverse transcription polymerase chain reaction (RT-PCR) on the minigene revealed that c.706+368T>C may create an SRSF2 (SC35) binding site via a 313-nt pseudo-exon, whereas c.1065+241C>A may produce an 81-nt pseudo-exon that strengthens the predicted SRSF1 (SF2/ASF) binding site. These results augment current knowledge of PAH genotypes and show that deep intronic analysis of PAH can genetically diagnose PKU.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.

SLC4A11 mutations causative of congenital hereditary endothelial dystrophy (CHED) progressing to Harboyan syndrome in consanguineous Pakistani families

BACKGROUND

Autosomal recessive corneal hereditary endothelial dystrophy (CHED) is a rare congenital disorder of cornea. Mutations in SLC4A11 gene are associated with CHED phenotype. CHED is also an early feature of Harboyan syndrome. The aim of the present study was to identify genetic mutations in the SLC4A11 gene in CHED cases belonging to inbred Pakistani families. Furthermore, all homozygous mutation carriers were investigated for hearing deficit.

METHODS AND RESULTS

This study included consanguineous CHED families presented at Al-Shifa Trust Eye Hospital, Rawalpindi, Pakistan from June 2018 to September 2018. DNA was extracted from blood samples. Direct sequencing of SLC4A11 gene was performed. All identified variants were evaluated by in silico programs i.e., SIFT, PolyPhen-2, and MutationTaster. Pathogenicity of the two identified splice site variants was analyzed by Human Splicing Finder and MaxEnt Scan. Screening of five CHED families revealed a total of three previously un reported (p.Arg128Gly, c.2241-2A > T and c.1898-2A > C in family CHED19, CHED22 and CHED26 respectively) and two already reported homozygous disease causing variants (p.Arg869Cys and p.Val824Met in family CHED24 and CHED25 respectively) as predicted by mutation taster. All of these variants segregated with disease phenotype and were not detected in controls.

CONCLUSION

Affected individuals of the five CHED families screened in this study had the disease due to SLC4A11 mutations and progressing to Harboyan syndrome. Identification of previously unreported mutations aid to heterogeneity of SLC4A11 and CHED pathogenesis as well as helped to provide genetic counseling to affected families.

A novel mouse model of PMS2 founder mutation that causes mismatch repair defect due to aberrant splicing

Hereditary non-polyposis colorectal cancer, now known as Lynch syndrome (LS) is one of the most common cancer predisposition syndromes and is caused by germline pathogenic variants (GPVs) in DNA mismatch repair (MMR) genes. A common founder GPV in PMS2 in the Canadian Inuit population, NM_000535.5: c.2002A>G, leads to a benign missense (p.I668V) but also acts as a de novo splice site that creates a 5 bp deletion resulting in a truncated protein (p.I668*). Individuals homozygous for this GPV are predisposed to atypical constitutional MMR deficiency with a delayed onset of first primary malignancy. We have generated mice with an equivalent germline mutation (Pms2c.1993A>G) and demonstrate that it results in a splicing defect similar to those observed in humans. Homozygous mutant mice are viable like the Pms2 null mice. However, unlike the Pms2 null mice, these mutant mice are fertile, like humans homozygous for this variant. Furthermore, these mice exhibit a significant increase in microsatellite instability and intestinal adenomas on an Apc mutant background. Rectification of the splicing defect in human and murine fibroblasts using antisense morpholinos suggests that this novel mouse model can be valuable in evaluating the efficacy aimed at targeting the splicing defect in PMS2 that is highly prevalent among the Canadian Inuits.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

A multistep approach to the genotype-phenotype analysis of Polish patients with tuberous sclerosis complex

The aim of this study was to evaluate a cost-effective diagnostic strategy for identification of casual variants for tuberous sclerosis complex (TSC) in the Polish population and to correlate the genetic results with selected phenotypic features. Fifty-five patients, aged 3-44 years, with a clinical diagnosis of TSC were enrolled into the study. All patients received a three-step analysis: next generation sequencing screening (NGS), multiplex ligation-dependent probe amplification (MLPA) and deep sequencing. This multistep approach obtained positive results in 51/55 (93%) patients: of the 51 positives TSC1 variants were observed in 16 (31%) and TSC2 variants in 35 (69%); these included 13 novel variants and two patients with mosaicism. Four patients (7%) had no mutation identified (NMI). Among the TSC1 gene variants, there were five nonsense, four frameshift, three large deletions, two missense and two splicing variants. For the TSC2 gene, 11 were missense, eight splicing, six frameshift, four large deletions, two in-frame deletions and four nonsense variants. The patients with TSC2 changes had their clinical diagnosis of TSC at a younger age than those with TSC1 changes (one year vs three years, p = 0.041). The TSC2 group demonstrated a higher number of major symptoms per patient (p = 0.04). Subependymal giant cell astrocytoma with concomitance of other brain lesions was more common in patients with missense mutations in either gene (23% vs 0%, p = 0.02). Such a multistep molecular diagnostic strategy could increase the possibility of detecting causal variants for TSC and may allow detection of mosaicism at low levels. Missense pathogenic variants in TSC1 or TSC2 gene might be associated with a higher risk of brain lesions.

Induction of cryptic pre-mRNA splice-switching by antisense oligonucleotides

Antisense oligomers (AOs) are increasingly being used to modulate RNA splicing in live cells, both for research and for the development of therapeutics. While the most common intended effect of these AOs is to induce skipping of whole exons, rare examples are emerging of AOs that induce skipping of only part of an exon, through activation of an internal cryptic splice site. In this report, we examined seven AO-induced cryptic splice sites in six genes. Five of these cryptic splice sites were discovered through our own experiments, and two originated from other published reports. We modelled the predicted effects of AO binding on the secondary structure of each of the RNA targets, and how these alterations would in turn affect the accessibility of the RNA to splice factors. We observed that a common predicted effect of AO binding was disruption of the exon definition signal within the exon's excluded segment.

Therapeutic manipulation of IKBKAP mis-splicing with a small molecule to cure familial dysautonomia

Approximately half of genetic disease-associated mutations cause aberrant splicing. However, a widely applicable therapeutic strategy to splicing diseases is yet to be developed. Here, we analyze the mechanism whereby IKBKAP-familial dysautonomia (FD) exon 20 inclusion is specifically promoted by a small molecule splice modulator, RECTAS, even though IKBKAP-FD exon 20 has a suboptimal 5' splice site due to the IVS20 + 6 T > C mutation. Knockdown experiments reveal that exon 20 inclusion is suppressed in the absence of serine/arginine-rich splicing factor 6 (SRSF6) binding to an intronic splicing enhancer in intron 20. We show that RECTAS directly interacts with CDC-like kinases (CLKs) and enhances SRSF6 phosphorylation. Consistently, exon 20 splicing is bidirectionally manipulated by targeting cellular CLK activity with RECTAS versus CLK inhibitors. The therapeutic potential of RECTAS is validated in multiple FD disease models. Our study indicates that small synthetic molecules affecting phosphorylation state of SRSFs is available as a new therapeutic modality for mechanism-oriented precision medicine of splicing diseases.

IntSplice2: Prediction of the Splicing Effects of Intronic Single-Nucleotide Variants Using LightGBM Modeling

Prediction of the effect of a single-nucleotide variant (SNV) in an intronic region on aberrant pre-mRNA splicing is challenging except for an SNV affecting the canonical GU/AG splice sites (ss). To predict pathogenicity of SNVs at intronic positions −50 (Int-50) to −3 (Int-3) close to the 3’ ss, we developed light gradient boosting machine (LightGBM)-based IntSplice2 models using pathogenic SNVs in the human gene mutation database (HGMD) and ClinVar and common SNVs in dbSNP with 0.01 ≤ minor allelic frequency (MAF) < 0.50. The LightGBM models were generated using features representing splicing cis-elements. The average recall/sensitivity and specificity of IntSplice2 by fivefold cross-validation (CV) of the training dataset were 0.764 and 0.884, respectively. The recall/sensitivity of IntSplice2 was lower than the average recall/sensitivity of 0.800 of IntSplice that we previously made with support vector machine (SVM) modeling for the same intronic positions. In contrast, the specificity of IntSplice2 was higher than the average specificity of 0.849 of IntSplice. For benchmarking (BM) of IntSplice2 with IntSplice, we made a test dataset that was not used to train IntSplice. After excluding the test dataset from the training dataset, we generated IntSplice2-BM and compared it with IntSplice using the test dataset. IntSplice2-BM was superior to IntSplice in all of the seven statistical measures of accuracy, precision, recall/sensitivity, specificity, F1 score, negative predictive value (NPV), and matthews correlation coefficient (MCC). We made the IntSplice2 web service at https://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice2.

Identification of novel single-nucleotide variants altering RNA splicing of PKD1 and PKD2

The development of sequencing techniques identified numerous genetic variants, and accurate evaluation of the clinical significance of these variants facilitates the diagnosis of Mendelian diseases. In the present study, 549 rare single- nucleotide variants of uncertain significance were extracted from the ADPKD and ClinVar databases. MaxEntScan scoresplice is an in silico splicing prediction tool that was used to analyze rare PKD1 and PKD2 variants of unknown significance. An in vitro minigene splicing assay was used to verify 37 splicing-altering candidates that were located within seven residues of the splice donor sequence excluding canonical GT dinucleotides or within 21 residues of the acceptor sequence excluding canonical AG dinucleotides of PKD1 and PKD2. We demonstrated that eight PKD1 variants alter RNA splicing and were predicted to be pathogenic.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Antisense oligonucleotide-based drug development for Cystic Fibrosis patients carrying the 3849+10 kb C-to-T splicing mutation

Background:

Antisense oligonucleotide (ASO)-based drugs for splicing modulation were recently approved for various genetic diseases with unmet need. Here we aimed to develop an ASO-based splicing modulation therapy for Cystic Fibrosis (CF) patients carrying the 3849 + 10 kb C-to-T splicing mutation in the CFTR gene.

Methods:

We have screened, in FRT cells expressing the 3849 + 10 kb C-to-T splicing mutation, ~30 2ʹ-O-Methyl-modified phosphorothioate ASOs, targeted to prevent the recognition and inclusion of a cryptic exon generated due to the mutation. The effect of highly potent ASO candidates on the splicing pattern, protein maturation and CFTR function was further analyzed in well differentiated primary human nasal and bronchial epithelial cells, derived from patients carrying at least one 3849 + 10 kb C-to-T allele.

Results:

A highly potent lead ASO, efficiently delivered by free uptake, was able to significantly increase the level of correctly spliced mRNA and completely restore the CFTR function to wild type levels in cells from a homozygote patient. This ASO led to CFTR function with an average of 43% of wild type levels in cells from various heterozygote patients. Optimized efficiency of the lead ASO was further obtained with 2ʹ-Methoxy Ethyl modification (2ʹMOE).

Conclusion:

The highly efficient splicing modulation and functional correction, achieved by free uptake of the selected lead ASO in various patients, demonstrate the ASO therapeutic potential benefit for CF patients carrying splicing mutations and is aimed to serve as the basis for our current clinical development.

A world-wide analysis of reduced sensitivity to DMI fungicides in the banana pathogen Pseudocercospora fijiensis

BACKGROUND

Pseudocercospora fijiensis is the causal agent of the black leaf streak disease (BLSD) of banana. Bananas are important global export commodities and a major staple food. Their susceptibility to BLSD pushes disease management towards excessive fungicide use, largely relying on multisite inhibitors and sterol demethylation inhibitors (DMIs). These fungicides are ubiquitous in plant disease control, targeting the CYP51 enzyme. We examined sensitivity to DMIs in P. fijiensis field isolates collected from various major banana production zones in Colombia, Costa Rica, Dominican Republic, Ecuador, the Philippines, Guadalupe, Martinique and Cameroon and determined the underlying genetic reasons for the observed phenotypes.

RESULTS

We observed a continuous range of sensitivity towards the DMI fungicides difenoconazole, epoxiconazole and propiconazole with clear cross-sensitivity. Sequence analyses of PfCYP51 in 266 isolates showed 28 independent amino acid substitutions, nine of which correlated with reduced sensitivity to DMIs. In addition to the mutations, we observed up to six insertions in the Pfcyp51 promoter. Such promoter insertions contain repeated elements with a palindromic core and correlate with the enhanced expression of Pfcyp51 and hence with reduced DMI sensitivity. Wild-type isolates from unsprayed bananas fields did not contain any promoter insertions.

CONCLUSION

The presented data significantly contribute to understanding of the evolution and global distribution of DMI resistance mechanisms in P. fijiensis field populations and facilitate the prediction of different DMI efficacy. The overall reduced DMI sensitivity calls for the deployment of a wider range of solutions for sustainable control of this major banana disease. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.

COL1A1, COL4A3, TIMP2 and TGFB1 polymorphisms in cervical insufficiency

OBJECTIVES

To investigate the association between selected single nucleotide polymorphisms (SNPs) with cervical insufficiency and its relationship with obstetric history.

METHODS

Twenty-eight women with cervical insufficiency (case group) and 29 non-pregnant women (control group) were included. The SNPs sequenced included rs2586490 in collagen type I alpha 1 chain (COL1A1), rs1882435 in collagen type IV alpha 3 chain (COL4A3), rs2277698 in metallopeptidase inhibitor 2 (TIMP2), and rs1800468 in transforming growth factor beta 1 (TGFB1).

RESULTS

We found a higher frequency of the normal allele in the control group (65.5%) and the homozygous mutated genotype in the case group (64.3%) for rs2586490 in COL1A1 (p=0.023). An unplanned finding in the cervical insufficiency group was a higher gestational age of delivery (median≥38 weeks) in the mutated allele than in the wild-type genotype (median of 28.2 weeks) for rs2857396, which is also in the COL1A1 gene (p=0.011).

CONCLUSIONS

The findings of the present study corroborate the hypothesis that cervical insufficiency has a genetic component and probably involves genes encoding proteins in the extracellular matrix, in addition to inflammatory processes.

Diversity in STK4 Deficiency and Review of the Literature

BACKGROUND

Serine-threonine kinase-4 (STK4) deficiency is an autosomal recessive (AR) combined immunodeficiency (CID).

OBJECTIVE

We aimed to define characteristic clinical and laboratory features to aid the differential diagnosis and determine the most suitable therapy.

METHODS

In addition to nine patients diagnosed, we reviewed 15 patients from medical literature. We compared B lymphocyte subgroups of our cohort with age-matched healthy controls.

RESULTS

In our cohort, the median age at symptom onset and age of diagnosis are 6years-8months (mo)(6-248mo) and 7years-5mo (6-260mo), respectively. The main clinical findings were infections (9/9), autoimmune/inflammatory diseases (7/9), and atopy (4/9). CD4 lymphopenia (9/9), lymphopenia (7/9), intermittent eosinophilia (4/9), transient neutropenia (3/9), low immunoglobulin (Ig) M (4/9), and high IgE (4/9) were common. Decreased recent thymic emigrants, naive and central memory T cells, albeit increased effector memory T cells were present. The increase in plasmablasts (p=0.003) and the decrease in switched memory B cells (p=0.022) were significant. Out of a total of 24 patients, cutaneous viral infections (n=20), recurrent pneumonia (n=18), Epstein Barr Virus (EBV)-associated lymphoproliferation (n=11), atopic dermatitis (n=10), autoimmune cytopenia (n=7), and lymphoma (n=6) were frequently seen. Lymphopenia, CD4 lymphopenia, high Ig G, A, and E were the most common laboratory characteristics.

CONCLUSION

The differential diagnosis with AR-hyperimmunoglobulin E syndrome is crucial as atopy and CD4 lymphopenia are prominent in both diseases. Immunoglobulins and antibacterial/antiviral prophylaxis are the mainstays of treatment. Clinicians may use immunomodulatory therapies during inflammatory/autoimmune complications. However, more data is needed to recommend hematopoietic stem cell transplantation (HSCT) as a safe therapy.

Enhancing Diagnosis Through RNA Sequencing

The diagnostic rate of comprehensive genomic sequencing remains only 25% to 30% due to the difficulty in interpreting variants of uncertain significance and noncoding mutations and in elucidating downstream effects of these and other genetic changes. Unlike DNA sequencing, RNA sequencing (RNAseq) reveals the functional consequence of genetic variation through the detection of abnormal gene expression levels, differences in gene splicing, and allele-specific expression. RNAseq can provide nearly 40% improvement in diagnostic rates depending on disease and tissue source. In this burgeoning era of precision medicine, RNAseq offers a powerful tool to improve diagnostic rates and understand disease mechanisms.

Clinical, Laboratory Features and Clinical Courses of Patients with Wiskott Aldrich Syndrome and X-linked Thrombocytopenia-A single center study

Objective: Wiskott Aldrich Syndrome is an X-linked primary immunodeficiency disorder characterized by microthrombocytopenia, severe immunodeficiency, and eczema. To define clinical-laboratory features, genetic defects (known/novel) of 23 patients of Wiskott Aldrich Syndrome/X-linked Thrombocytopenia (WAS/XLT) cohort, establish relationships between molecular defects and clinical features if present, evaluate patients who underwent hematopoietic stem cell transplantation (HSCT) and did not.Methods: Qualitative analysis from patients' hospital files and Sanger sequencing for molecular diagnosis was performed. Twenty-two WAS patients and one XLT patient were included in the study.Results: The median age of diagnosis was 15 months (2.5-172 months). The most common symptom was otitis media and all patients had microthrombocytopenia. Autoimmune findings were detected in 34.7% (8 patients) of the patients; three patients (13%) had positive anti-nuclear antibody (ANA), three patients (13%) hemolytic anemia, one patient autoimmune neutropenia, two patients vasculitis, and one patient demyelinating polyneuropathy. Nine of the 23 (39,1%) patients had HSCT with nearly 90% success. We identified 13 different mutations in our cohort; seven were novel.Conclusions: HSCT is the only curative treatment for WAS. The study confirms that early diagnosis is very important for the success of therapy, so we must increase awareness in society and physicians to keep an eye out for clues. Our study cohort and follow-up period are not sufficient to establish phenotype-genotype correlation, so a larger cohort from various centers with longer follow-up will be more decisive.

A variant in TMPRSS2 is associated with decreased disease severity in COVID-19

Background

Mortality due to COVID-19 caused by SARS-CoV-2 infection varies among populations. Functional relevance of genetic variations in Angiotensin-converting enzyme 2 (ACE2) and Transmembrane serine protease 2 (TMPRSS2), two crucial host factors for viral entry, might explain some of this variation.

Methods

In this comparative study in Indian subjects, we recruited 510 COVID-19 patients and retrieved DNA from 520 controls from a repository. Associations between variants in ACE2 and TMPRSS2 with disease severity were identified by whole exome sequencing (WES, n = 20) and targeted genotyping (n = 1010). Molecular dynamic simulations (MDS) were performed to explore functional relevance of the variants. Cleavage of spike glycoprotein by wild and variant TMPRSS2 was determined in HEK293T cells. Potential effects of confounders on the association between genotype and disease severity were tested (Mantel-Haenszel test).

Results

WES identified deleterious variant in TMPRSS2 (rs12329760, G > A, p. V160M). The minor allele frequency (MAF) was 0·27 in controls, 0·31 in asymptomatic, 0·21 in mild-to-moderately affected and 0·19 in severely affected COVID-19 patients. Risk of severity increased with decreasing MAF: Asymptomatic: Odds ratio-0·69 (95% CI–0·52–0·93; p = 0·01); mild-to-moderate: Odds ratio-1·89 (95% CI–1·22–2.92;p = 0·004) and severe: Odds ratio-1·79 (95% CI–1·11–2.88;p = 0·01). No confounding effect of diabetes and hypertension were observed on the risk of developing severe COVID-19 disease with respect to genotype. MDS revealed decreased stability of TMPRSS2 with 160 M variant. Spike glycoprotein cleavage by TMPRSS2 reduced ~2·4-fold in cells expressing 160 M variant.

Conclusion

We demonstrate association of TMPRSS2 variant rs12329760 with decreased disease severity in COVID-19 patients from India.

Genetic evaluation of the variants using MassARRAY in non-small cell lung cancer among North Indians

Lung cancer is genetically diverse and a major health burden. Non-small cell lung cancer (NSCLC) accounts for 80% of total lung cancer cases and 20% cases are Small cell lung cancer (SCLC). The present case–control association study focused on the cost effective high throughput genotyping using Agena MassARRAY matrix-assisted laser desorption/ionization-time of flight, mass spectrometry (MALDI-TOF) platform to analyze the genetic association of candidate genetic variants. We performed multiplex PCR and genotyped twelve single nucleotide polymorphisms (SNPs) in 723 samples (162 NSCLC cases and 592 healthy controls). These genetic variants were selected from literature for their association with various cancers worldwide and this is the first study from the region to examine these critically important genetic variants. With prospective case–control association study design, twelve variants from ten genes were evaluated. Amongst these six variants, TCF21 (rs12190287), ERCC1 (rs2298881, 11615), ERCC5 (rs751402), ARNTL (rs4757151), BRIP1 (rs4986764) showed significant association with NSCLC risk (p ≤ 0.003) in Jammu and Kashmir population. In-silico findings of these genetic variants showed remarkable functional roles that needs in-vitro validations. It is further anticipated that such case control studies will help us in understanding the missing heritability of non-small cell lung cancer.

MLH1 intronic variants mapping to + 5 position of splice donor sites lead to deleterious effects on RNA splicing

Germline pathogenic variants in the DNA mismatch repair genes (MMR): MLH1, MSH2, MSH6, and PMS2, are causative of Lynch syndrome (LS). However, many of the variants mapping outside the invariant splice site positions (IVS ± 1, IVS ± 2) are classified as variants of unknown significance (VUS). Three such variants (MLH1 c.588+5G>C, c.588+5G>T and c.677+5G>A) were identified in 8 unrelated LS families from Argentina, Brazil and Chile. Herein, we collected clinical information on these families and performed segregation analysis and RNA splicing studies to assess the implication of these VUS in LS etiology. Pedigrees showed a clear pattern of variant co-segregation with colorectal cancer and/or other LS-associated malignancies. Tumors presented deficient expression of MLH1-PMS2 proteins in 7/7 of the LS families, and MSI-high status in 3/3 cases. Moreover, RNA analyses revealed that c.588+5G>C and c.588+5G>T induce skipping of exon 7 whereas c.677+5G>A causes skipping of exon 8. In sum, we report that the combined clinical findings in the families and the molecular studies provided the evidences needed to demonstrate that the three MLH1 variants are causative of LS and to classify c.588+5G>C and c.677+5G>A as class 5 (pathogenic), and c.588+5G>T as class 4 (likely-pathogenic). Our findings underline the importance of performing clinical and family analyses, as well as RNA splicing assays in order to determine the clinical significance of intronic variants, and contribute to the genetic counseling and clinical management of patients and their relatives.

SRSF1 promotes the inclusion of exon 3 of SRA1 and the invasion of hepatocellular carcinoma cells by interacting with exon 3 of SRA1pre-mRNA

Steroid receptor RNA activator 1 (SRA1) has been described as a novel transcriptional co-activator that affects the migration of cancer cells. Through RT-PCR, we identified that skipping exon 3 of SRA1 produces two isoforms, including the truncated short isoform, SRA1-S, and the long isoform, SRA1-L. However, the effect of these two isomers on the migration of HCC cells, as well as the specific mechanism of exon 3 skipping remain unclear. In this study, we found up regulated expression of SRSF1 and SRA1-L in highly metastatic HCCLM3, as well as in HCCs with SRSF1 demonstrating the strongest correlation with SRA1-L. In contrast, we observed a constitutively low expression of SRA1-S and SRSF1 in lowly metastatic HepG2 cells. Overexpression of SRSF1 or SRA1-L promoted migration and invasion by increasing the expression of CD44, while SRA1-S reversed the effect of SRSF1 and SRA1-L in vitro. In addition, lung metastasis in mice revealed that, knockdown of SRSF1 or SRA1-L inhibited the migration of HCC cells, while SRA1-L overexpression abolished the effect of SRSF1 knockout and instead promoted HCC cells migration in vivo. More importantly, RNA immunoprecipitation and Cross-link immunoprecipitation analyses showed that SRSF1 interacts with exon 3 of SRA1 to up regulate the expression of SRA1-L in HCC cells. RNA pull-down results indicated that SRSF1 could also bind to exon 3 of SRA1 in vitro. Finally, minigene -MS2 mutation experiments showed that mutation of the SRA1 exon 3 binding site for SRSF1 prevented the binding of SRA1 pre-mRNA. In summary, our results provide experimental evidence that SRA1 exon 3 inclusion is up regulated by SRSF1 to promote tumor invasion and metastasis in hepatocellular carcinoma.

Rapidly Growing Protein-Centric Technologies to Extensively Identify Protein-RNA Interactions: Application to the Analysis of Co-Transcriptional RNA Processing

During mRNA transcription, diverse RNA-binding proteins (RBPs) are recruited to RNA polymerase II (RNAP II) transcription machinery. These RBPs bind to distinct sites of nascent RNA to co-transcriptionally operate mRNA processing. Recent studies have revealed a close relationship between transcription and co-transcriptional RNA processing, where one affects the other’s activity, indicating an essential role of protein–RNA interactions for the fine-tuning of mRNA production. Owing to their limited amount in cells, the detection of protein–RNA interactions specifically assembled on the transcribing RNAP II machinery still remains challenging. Currently, cross-linking and immunoprecipitation (CLIP) has become a standard method to detect in vivo protein–RNA interactions, although it requires a large amount of input materials. Several improved methods, such as infrared-CLIP (irCLIP), enhanced CLIP (eCLIP), and target RNA immunoprecipitation (tRIP), have shown remarkable enhancements in the detection efficiency. Furthermore, the utilization of an RNA editing mechanism or proximity labeling strategy has achieved the detection of faint protein–RNA interactions in cells without depending on crosslinking. This review aims to explore various methods being developed to detect endogenous protein–RNA interaction sites and discusses how they may be applied to the analysis of co-transcriptional RNA processing.

Hotspot exons are common targets of splicing perturbations

High-throughput splicing assays have demonstrated that many exonic variants can disrupt splicing; however, splice-disrupting variants distribute non-uniformly across genes. We propose the existence of exons that are particularly susceptible to splice-disrupting variants, which we refer to as hotspot exons. Hotspot exons are also more susceptible to splicing perturbation through drug treatment and knock-down of RNA-binding proteins. We develop a classifier for exonic splice-disrupting variants and use it to infer hotspot exons. We estimate that 1400 exons in the human genome are hotspots. Using panels of splicing reporters, we demonstrate how the ability of an exon to tolerate a mutation is inversely proportional to the strength of its neighboring splice sites.

Splicing-disrupting mutations are linked to diseases. By employing a machine learning approach, the authors show that certain exons, termed hotspot exons, are enriched for splicing-disruption variants and susceptible to splicing perturbations.

Mutations in cis that affect mRNA synthesis, processing and translation

Genetic mutations that cause hereditary diseases usually affect the composition of the transcribed mRNA and its encoded protein, leading to instability of the mRNA and/or the protein. Sometimes, however, such mutations affect the synthesis, the processing or the translation of the mRNA, with similar disastrous effects. We here present an overview of mRNA synthesis, its posttranscriptional modification and its translation into protein. We then indicate which elements in these processes are known to be affected by pathogenic mutations, but we restrict our review to mutations in cis, in the DNA of the gene that encodes the affected protein. These mutations can be in enhancer or promoter regions of the gene, which act as binding sites for transcription factors involved in pre-mRNA synthesis. We also describe mutations in polyadenylation sequences and in splice site regions, exonic and intronic, involved in intron removal. Finally, we include mutations in the Kozak sequence in mRNA, which is involved in protein synthesis. We provide examples of genetic diseases caused by mutations in these DNA regions and refer to databases to help identify these regions. The over-all knowledge of mRNA synthesis, processing and translation is essential for improvement of the diagnosis of patients with genetic diseases.