Transcript Isoforms Research Articles

An impaired splicing program underlies differentiation defects in hSOD1G93A neural progenitor cells.

Amyotrophic lateral sclerosis (ALS) is an adult devastating neurodegenerative disease characterized by the loss of upper and lower motor neurons (MNs), resulting in progressive paralysis and death. Genetic animal models of ALS have highlighted dysregulation of synaptic structure and function as a pathogenic feature of ALS-onset and progression. Alternative pre-mRNA splicing (AS), which allows expansion of the coding power of genomes by generating multiple transcript isoforms from each gene, is widely associated with synapse formation and functional specification. Deciphering the link between aberrant splicing regulation and pathogenic features of ALS could pave the ground for novel therapeutic opportunities. Herein, we found that neural progenitor cells (NPCs) derived from the hSOD1G93A mouse model of ALS displayed increased proliferation and propensity to differentiate into neurons. In parallel, hSOD1G93A NPCs showed impaired splicing patterns in synaptic genes, which could contribute to the observed phenotype. Remarkably, master splicing regulators of the switch from stemness to neural differentiation are de-regulated in hSOD1G93A NPCs, thus impacting the differentiation program. Our data indicate that hSOD1G93A mutation impacts on neurogenesis by increasing the NPC pool in the developing mouse cortex and affecting their intrinsic properties, through the establishment of a specific splicing program.

A novel antisense lncRNA, ARBAG harboring an RNA destabilizing GWAS variant for C-peptide dictates the transcript isoforms of GABRA6 in cerebellum.

Human disease-associated genetic variations often map to long non-coding RNA (lncRNA) genes; however, elucidation of their functional impact is challenging. We previously identified a new genetic variant rs4454083 (A/G) residing in exon of an uncharacterized lncRNA ARBAG that strongly associates with plasma levels of C-peptide, a hormone that regulates insulin bioavailability. On the opposite strand, rs4454083 also corresponds to an intron of a cerebellum-specific GABA receptor subunit gene GABRA6 that mediates strengthening of inhibitory synapses by insulin. Here, we show that alleles of rs4454083 modulate transcript levels of the antisense gene, ARBAG, which then controls the expression of the sense gene, GABRA6. Predisposing to low C-peptide, GG (a minor allele genotype across ethnicities) stabilizes ARBAG lncRNA causing higher transcript levels in cerebellum. ARBAG lncRNA abundance leads to cleavage of GABRA6 mRNA at the complementary region, resulting in a dysfunctional GABRA6 protein that would not be recruited for synapse strengthening. Together, our findings in human cerebellar cell-line and induced Pluripotent Stem Cells (iPSCs) demonstrate biological role of a novel lncRNA in determining the ratio of mRNA isoforms of a protein-coding gene and the ability of an embedded variant in modulating lncRNA stability leading to inter-individual differences in protein expression.

Full-length transcript alterations in human bronchial epithelial cells with U2AF1 S34F mutations.

U2AF1 is one of the most recurrently mutated splicing factors in lung adenocarcinoma and has been shown to cause transcriptome-wide pre-mRNA splicing alterations; however, the full-length altered mRNA isoforms associated with the mutation are largely unknown. To better understand the impact U2AF1 has on full-length isoform fate and function, we conducted high-throughput long-read cDNA sequencing from isogenic human bronchial epithelial cells with and without a U2AF1 S34F mutation. We identified 49,366 multi-exon transcript isoforms, more than half of which did not match GENCODE or short-read-assembled isoforms. We found 198 transcript isoforms with significant expression and usage changes relative to WT, only 68% of which were assembled by short reads. Expression of isoforms from immune-related genes is largely down-regulated in mutant cells and without observed splicing changes. Finally, we reveal that isoforms likely targeted by nonsense-mediated decay are down-regulated in U2AF1 S34F cells, suggesting that isoform changes may alter the translational output of those affected genes. Altogether, our work provides a resource of full-length isoforms associated with U2AF1 S34F in lung cells.

Supervised Relevance-Redundancy assessments for feature selection in omics-based classification scenarios

Background and objective:Many classification tasks in translational bioinformatics and genomics are characterized by the high dimensionality of potential features and unbalanced sample distribution among classes. This can affect classifier robustness and increase the risk of overfitting, curse of dimensionality and generalization leaks; furthermore and most importantly, this can prevent obtaining adequate patient stratification required for precision medicine in facing complex diseases, like cancer. Setting up a feature selection strategy able to extract only proper predictive features by removing irrelevant, redundant, and noisy ones is crucial to achieving valuable results on the desired task. Methods:We propose a new feature selection approach, called ReRa, based on supervised Relevance-Redundancy assessments. ReRa consists of a customized step of relevance-based filtering, to identify a reduced subset of meaningful features, followed by a supervised similarity-based procedure to minimize redundancy. This latter step innovatively uses a combination of global and class-specific similarity assessments to remove redundant features while preserving those differentiated across classes, even when these classes are strongly unbalanced. Results:We compared ReRa with several existing feature selection methods to obtain feature spaces on which performing breast cancer patient subtyping using several classifiers: we considered two use cases based on gene or transcript isoform expression. In the vast majority of the assessed scenarios, when using ReRa-selected feature spaces, the performances were significantly increased compared to simple feature filtering, LASSO regularization, or even MRmr — another Relevance-Redundancy method. The two use cases represent an insightful example of translational application, taking advantage of ReRa capabilities to investigate and enhance a clinically-relevant patient stratification task, which could be easily applied also to other cancer types and diseases. Conclusions:ReRa approach has the potential to improve the performance of machine learning models used in an unbalanced classification scenario. Compared to another Relevance-Redundancy approach like MRmr, ReRa does not require tuning the number of preserved features, ensures efficiency and scalability over huge initial dimensionalities and allows re-evaluation of all previously selected features at each iteration of the redundancy assessment, to ultimately preserve only the most relevant and class-differentiated features.

Human intestine and placenta exhibit tissue-specific expression of RAGE isoforms

The receptor for advanced glycation end products (RAGE) is encoded by AGER, a gene that is subjected to tissue-specific alternative splicing. Splice variants of RAGE in intestine and placenta are unknown and contradictory data concerning RAGE protein expression in these tissues have been published. As a basis for future functional studies, we examined RAGE expression in small intestine, colon and placentas. PCR cloning revealed that full-length RAGE is the only RAGE transcript isoform expressed in placenta. In the small intestine, the major transcript isoform detected was RAGE_v1 encoding the C-terminally truncated soluble receptor. In the colon, both full-length RAGE as well as several splice variants were identified. Four antibodies were used to study protein expression by immunoblotting and were carefully validated. Appropriate controls were essential to avoid misinterpretation of bands caused by non-specific reactivity of antibodies. Only one of four antibodies tested detected full-length RAGE in placenta, whereas no RAGE-specific band was detected in intestinal tissues despite loading >30-fold more intestinal tissue than the positive control, human lung. RAGE expression levels in the placenta were 100-fold lower compared with human lung when analyzed by ELISA, and no significant differences in RAGE expression were detected between healthy placentas and placentas from women with preeclampsia, gestational diabetes mellitus, or fetal growth restriction. We conclude that healthy placental chorionic tissue expresses low levels of full-length RAGE, whereas expression of the tissue-specific intestinal isoforms is below the limit of detection. Low RAGE expression levels in combination with a lack of antibody validation may explain the conflicting published results on RAGE protein expression in intestine and placenta.

HiCognition: a visual exploration and hypothesis testing tool for 3D genomics

Genome browsers facilitate integrated analysis of multiple genomics datasets yet visualize only a few regions at a time and lack statistical functions for extracting meaningful information. We present HiCognition, a visual exploration and machine-learning tool based on a new genomic region set concept, enabling detection of patterns and associations between 3D chromosome conformation and collections of 1D genomics profiles of any type. By revealing how transcription and cohesion subunit isoforms contribute to chromosome conformation, we showcase how the flexible user interface and machine learning tools of HiCognition help to understand the relationship between the structure and function of the genome.

Identification of experimentally-supported poly(A) sites in single-cell RNA-seq data with SCINPAS.

Alternative polyadenylation is a main driver of transcriptome diversity in mammals, generating transcript isoforms with different 3' ends via cleavage and polyadenylation at distinct polyadenylation (poly(A)) sites. The regulation of cell type-specific poly(A) site choice is not completely resolved, and requires quantitative poly(A) site usage data across cell types. 3' end-based single-cell RNA-seq can now be broadly used to obtain such data, enabling the identification and quantification of poly(A) sites with direct experimental support. We propose SCINPAS, a computational method to identify poly(A) sites from scRNA-seq datasets. SCINPAS modifies the read deduplication step to favor the selection of distal reads and extract those with non-templated poly(A) tails. This approach improves the resolution of poly(A) site recovery relative to standard software. SCINPAS identifies poly(A) sites in genic and non-genic regions, providing complementary information relative to other tools. The workflow is modular, and the key read deduplication step is general, enabling the use of SCINPAS in other typical analyses of single cell gene expression. Taken together, we show that SCINPAS is able to identify experimentally-supported, known and novel poly(A) sites from 3' end-based single-cell RNA sequencing data.

Hybrid sequencing discloses unique aspects of the transcriptomic architecture in equid alphaherpesvirus 1

This study employed both short-read sequencing (SRS, Illumina) and long-read sequencing (LRS Oxford Nanopore Technologies) platforms to conduct a comprehensive analysis of the equid alphaherpesvirus 1 (EHV-1) transcriptome. The study involved the annotation of canonical mRNAs and their transcript variants, encompassing transcription start site (TSS) and transcription end site (TES) isoforms, in addition to alternative splicing forms. Furthermore, the study revealed the presence of numerous non-coding RNA (ncRNA) molecules, including intergenic and antisense transcripts, produced by EHV-1. An intriguing finding was the abundant production of chimeric transcripts, some of which potentially encode fusion polypeptides. Moreover, EHV-1 exhibited a greater incidence of transcriptional overlaps and splicing compared to related viruses. It is noteworthy that many genes have their unique TESs along with the co-terminal transcription ends, a characteristic scarcely seen in other alphaherpesviruses. The study also identified transcripts that overlap the replication origins of the virus. Moreover, a novel ncRNA, referred to as NOIR, was found to intersect with the 5′-ends of longer transcript isoform specified by the major transactivator genes ORF64 and ORF65, surrounding the OriL. These findings together imply the existence of a key regulatory mechanism that governs both transcription and replication through, among others, a process that involves interference between the DNA and RNA synthesis machineries.

The clinical phenotype of germline RUNX1 mutations in relation to the accompanying somatic variants and RUNX1 isoform expression.

Germline RUNX1 mutations lead to familial platelet disorder with associated myeloid malignancy (FPDMM), characterized by thrombocytopenia, abnormal bleeding, and an elevated risk of developing myelodysplastic neoplasia (MDS) and acute myeloid leukemia (AML) at young age. However, it is not known why or how germline carriers of RUNX1 mutations have a particular propensity to develop myeloid hematologic malignancies, but the acquisition and composition of somatic mutations are believed to initiate and determine disease progression. We present a novel family pedigree that shares a common germline RUNX1R204* variant and exhibits a spectrum of somatic mutations and related myeloid malignancies (MM). RUNX1 mutations are associated with inferior clinical outcome; however, the proband of this family developed MDS with ring sideroblasts (MDS-RS), classified as a low-risk MDS subgroup. His relatively indolent clinical course is likely due to a specific somatic mutation in the SF3B1 gene. While the three main RUNX1isoforms have been ascribed various roles in normal hematopoiesis, they are now being increasingly recognized as involved in myeloid disease. We investigated the RUNX1 transcript isoform patterns in the proband and his sister, who carries the same germline RUNX1R204* variant, and has FPDMM but no MM. We demonstrate a RUNX1a increase in MDS-RS, as previously reported in MM. Interestingly, we identify a striking unbalance of RUNX1b and -c in FPDMM. In conclusion, this report reinforces the relevance of somatic variants on the clinical phenotypic heterogeneity in families with germline RUNX1 deficiency and investigates a potential new role for RUNX1 isoform disequilibrium as a mechanism for development of MM.

High-throughput RNA isoform sequencing using programmed cDNA concatenation.

Full-length RNA-sequencing methods using long-read technologies can capture complete transcript isoforms, but their throughput is limited. We introduce multiplexed arrays isoform sequencing (MAS-ISO-seq), a technique for programmably concatenating complementary DNAs (cDNAs) into molecules optimal for long-read sequencing, increasing the throughput >15-fold to nearly 40 million cDNA reads per run on the Sequel IIe sequencer. When applied to single-cell RNA sequencing of tumor-infiltrating T cells, MAS-ISO-seq demonstrated a 12- to 32-fold increase in the discovery of differentially spliced genes.

The Versatility of Plectin in Cancer: A Pan-Cancer Analysis on Potential Diagnostic and Prognostic Impacts of Plectin Isoforms

Plectin, encoded by PLEC, is a cytoskeletal and scaffold protein with a number of unique isoforms that act on various cellular functions such as cell adhesion, signal transduction, cancer cell invasion, and migration. While plectin has been shown to display high expression and mislocalization in tumor cells, our knowledge of the biological significance of plectin and its isoforms in tumorigenesis remain limited. In this study, we first performed pathway enrichment analysis to identify cancer hallmark proteins associated with plectin. Then, a pan-cancer analysis was performed using RNA-seq data collected from the Cancer Genome Atlas (TCGA) to detect the mRNA expression levels of PLEC and its transcript isoforms, and the prognostic as well as diagnostic significance of the transcript isoforms was evaluated considering cancer stages. We show here that several tissue specific PLEC isoforms are dysregulated in different cancer types and stages but not the expression of PLEC. Among them, PLEC 1d and PLEC 1f are potential biomarker candidates and call for further translational and personalized medicine research. This study makes a contribution as a stride to unravel the molecular mechanisms underpinning plectin isoforms in cancer development and progression by revealing the potent plectin isoforms in different stages of cancer as potential early cancer detection biomarkers. Importantly, uncovering how plectin isoforms guide malignancy and particular cancer types by comprehensive functional studies might open new avenues toward novel cancer therapeutics.

Full-length transcriptome reveals rapid genetic changes in triploid hybrid grass carp derived from female grass carp × male topmouth culter

Distant hybridization refers to the combination of different genomes from different species and further leads to genetic structure changes in hybrids. For next-generation sequencing limitations, the genetic changes with respect to hybrids or polyploids cannot be fully explained. Here, we obtained 117,526 full-length isoforms from a triploid hybrid grass carp (3nGT, 3n = 72) by the crossing of female grass carp (GC, 2n = 48) and male topmouth culter (TC, 2n = 48) using single-molecule long reads sequencing technology. Further, the gene fusion, alternative splicing (AS), transcription of homologous genes, and novel genes were analyzed. The results showed that 1563 fusion genes were detected in 3nGT including 872 heterogenous fusion genes derived from GC and TC genomes. Meanwhile, 175,644 AS events were identified and divided into five models, of which IntronR was the dominant model. In addition, 8964 novel genes were discovered, of which 2802 novel genes were annotated. Finally, scans results of Hox gene family revealed that both GC and TC Hox genes were transcribed in 3nGT, but presented asymmetric transcription patterns with a bias towards the maternal parent GC genome. This study not only provides a rich resource of transcript isoforms for the superior triploid grass carp, but also provides an important insight into rapid genetic changes in hybrid fish.

Characterization of the genomic and transcriptional structure of chicken NRG4 gene.

Neuregulin 4 (NRG4) is an important adipocytokine, which plays crucial roles in maintaining energy balance, regulating glucose and lipid metabolism, and preventing non-alcoholic fatty liver disease in mammals. At present, the genomic organization, transcript and protein isoforms of human NRG4 gene have been fully explored. Previous studies in our laboratory have shown that the NRG4 gene is expressed in chicken adipose tissue, but the chicken NRG4 (cNRG4) genomic structure, transcript and protein isoforms are still unknown. To this end, in this study, the genomic and transcriptional structure of the cNRG4 gene were systematically investigated using rapid amplification of cDNA ends (RACE) and reverse transcription-polymerase chain reaction (RT-PCR). The results showed that the coding region (CDS) of the cNRG4 gene was small, but it had a very complex transcriptional structure characterized by multiple transcription start sites, alternative splicing, intron retention, cryptic exons, and alternative polyadenylation, thus leading to production of four 5?UTR isoforms (cNRG4 A, cNRG4 B, cNRG4 C, and cNRG4 D) and six 3?UTR isoforms (cNRG4 a, cNRG4 b, cNRG4 c, cNRG4 d, cNRG4 e, and cNRG4 f) of the cNRG4 gene. The cNRG4 gene spanned 21,969 bp of genomic DNA (Chr.10:3,490,314~3,512,282) and consisted of 11 exons and 10 introns. Compared with the cNRG4 gene mRNA sequence (NM_001030544.4), two novel exons and one cryptic exon of the cNRG4 gene were identified in this study. Bioinformatics analysis, RT-PCR, cloning and sequencing analysis showed that the cNRG4 gene could encode three protein isoforms (cNRG4-1, cNRG4-2 and cNRG4-3). This study lays a foundation for further research on the function and regulation of the cNRG4 gene.

Adipose tissue plasticity in pheochromocytoma patients suggests a role of the splicing machinery in human adipose browning

Adipose tissue from pheochromocytoma patients acquires brown fat features, making it a valuable model for studying the mechanisms that control thermogenic adipose plasticity in humans. Transcriptomic analyses revealed a massive downregulation of splicing machinery components and splicing regulatory factors in browned adipose tissue from patients, with upregulation of a few genes encoding RNA-binding proteins potentially involved in splicing regulation. These changes were also observed in cell culture models of human brown adipocyte differentiation, confirming a potential involvement of splicing in the cell-autonomous control of adipose browning. The coordinated changes in splicing are associated with a profound modification in the expression levels of splicing-driven transcript isoforms for genes involved in the specialized metabolism of brown adipocytes and those encoding master transcriptional regulators of adipose browning. Splicing control appears to be a relevant component of the coordinated gene expression changes that allow human adipose tissue to acquire a brown phenotype.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe2+ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Human neutrophils are a key target of glucocorticoids

Abstract Neutrophils play central roles in immunity and inflammation. Glucocorticoids (GCs) are the cornerstone of anti-inflammatory and immunosuppressive therapy. While GCs trigger an acute rise in circulating neutrophils, little is known about the mechanisms by which GCs modulate neutrophil biology. We studied the response of 9 primary human cells to in vitro treatment with GCs. Surprisingly, neutrophils were transcriptionally the most responsive cell type, and they express proportionally the most glucocorticoid receptor (GR) transcripts. Long-read RNA-seq revealed 42 full-length GR transcript isoforms in neutrophils, twice as many as had been reported previously in all other cell types. At the protein level, neutrophils also have a unique and diverse set of GR isoforms. In contrast to the widely accepted model of GR signaling, in which GR binds GC molecules in the cytoplasm and then translocates to the nucleus, GR localization in neutrophils is almost exclusively nuclear, even in the absence of GC treatment. To validate our findings in vivo, we then performed a study of healthy volunteers, with blood collected before and serially after GC infusion. We generated a pathway-level map of the human neutrophil transcriptome at baseline and in response to GCs. The in vivo transcriptional response of human neutrophils to GCs was as strong as that observed in vitro, with a predominance of upregulation. Among the most strongly induced genes are several that are either not expressed, or expressed at very low levels, in human neutrophils at baseline. Some of the proteins encoded by these genes may have putative effects on neutrophil function, including SAMSN1 and ZBTB16. These findings indicate that human neutrophils are a key target for GC action.

Nano-DMS-MaP allows isoform-specific RNA structure determination

Genome-wide measurements of RNA structure can be obtained using reagents that react with unpaired bases, leading to adducts that can be identified by mutational profiling on next-generation sequencing machines. One drawback of these experiments is that short sequencing reads can rarely be mapped to specific transcript isoforms. Consequently, information is acquired as a population average in regions that are shared between transcripts, thus blurring the underlying structural landscape. Here, we present nanopore dimethylsulfate mutational profiling (Nano-DMS-MaP)—a method that exploits long-read sequencing to provide isoform-resolved structural information of highly similar RNA molecules. We demonstrate the value of Nano-DMS-MaP by resolving the complex structural landscape of human immunodeficiency virus-1 transcripts in infected cells. We show that unspliced and spliced transcripts have distinct structures at the packaging site within the common 5′ untranslated region, likely explaining why spliced viral RNAs are excluded from viral particles. Thus, Nano-DMS-MaP is a straightforward method to resolve biologically important transcript-specific RNA structures that were previously hidden in short-read ensemble analyses.

Exploring Toxin Genes of Myanmar Russell's Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics.

The Russell's viper (Daboia siamensis) is a medically important venomous snake in Myanmar. Next-generation sequencing (NGS) shows potential to investigate the venom complexity, giving deeper insights into snakebite pathogenesis and possible drug discoveries. mRNA from venom gland tissue was extracted and sequenced on the Illumina HiSeq platform and de novo assembled by Trinity. The candidate toxin genes were identified via the Venomix pipeline. Protein sequences of identified toxin candidates were compared with the previously described venom proteins using Clustal Omega to assess the positional homology among candidates. Candidate venom transcripts were classified into 23 toxin gene families including 53 unique full-length transcripts. C-type lectins (CTLs) were the most highly expressed, followed by Kunitz-type serine protease inhibitors, disintegrins and Bradykinin potentiating peptide/C-type natriuretic peptide (BPP-CNP) precursors. Phospholipase A2, snake venom serine proteases, metalloproteinases, vascular endothelial growth factors, L-amino acid oxidases and cysteine-rich secretory proteins were under-represented within the transcriptomes. Several isoforms of transcripts which had not been previously reported in this species were discovered and described. Myanmar Russell's viper venom glands displayed unique sex-specific transcriptome profiles which were correlated with clinical manifestation of envenoming. Our results show that NGS is a useful tool to comprehensively examine understudied venomous snakes.

Molecular Gene Expression Analysis of Medium Spiny Neurons in Huntington’s Disease (P1-11.015)

<h3>Objective:</h3> To identify molecular pathways underlying caudate, medium spiny neuron (MSN) degeneration in post-mortem Huntington’s Disease (HD) subjects. <h3>Background:</h3> Huntington’s Disease is an autosomal dominant, neurodegenerative disorder marked by aberrant expression of the expanded Huntingtin gene (HTT). Expression of mutant <i>HTT</i> results in selective atrophy of MSNs in the striatum. Previous studies have analyzed the caudate transcriptome in HD with a lack of consistency and an absence of corresponding proteome characterization. <h3>Design/Methods:</h3> Our study includes 20 HD and 20 neurotypical control post-mortem human subjects that are 55% female with a mean age of 57.6 at the time of death. We employed laser capture microdissection (LCM) to isolate the whole cell bodies of MSNs (~5000 neuron enrichments) from each of the 40 samples, facilitating single-cell analysis of the cytosol in addition to MSN nuclei. We performed RNA sequencing and labeled quantitative tandem mass spectrometry (TMT-MS/MS) on the same MSN lysate. Differential molecular analyses reveal RNA and protein isoform changes and phospho-site signaling changes in HD caudate. <h3>Results:</h3> We present MSN transcriptome, proteome, and phospho-proteome differences underlying HD at both the transcript and protein feature levels. Our findings suggest that there are <i>HTT</i> protein-protein interactions specific to MSN degeneration in HD. <h3>Conclusions:</h3> Applying this methodology and analytical paradigm, we expect to find (i) differentially expressed pathways (transcript and protein isoforms) underlying HD in MSN populations (ii) MSN-specific gene and protein interaction networks that have previously not been defined in controls or HD subjects, and (iii) phosphorylation-driven signaling changes in specific protein pathways that are closer to functional gene expression underlying HD. <b>Disclosure:</b> Mr. El-Halawani has nothing to disclose. Miss Martins has nothing to disclose. Dr. A Bharadwaj has nothing to disclose. Dr. Tao has nothing to disclose. The institution of Dr. Chen has received research support from National Science Foundation. Dr. Kleinman has received personal compensation in the range of $500-$4,999 for serving on a Scientific Advisory or Data Safety Monitoring board for Merck. Dr. Kleinman has received personal compensation in the range of $500-$4,999 for serving as an Expert Witness for Austin Sidley law Firm. Daniel Weinberger has nothing to disclose. Dr. Hyde has received personal compensation for serving as an employee of Lieber Institute for Brain Development. The institution of Dr. Hyde has received research support from National Institutes of Health.

Analysis of 200 unrelated individuals with a constitutional NF1 deep intronic pathogenic variant reveals that variants flanking the alternatively spliced NF1 exon 31 [23a] cause a classical neurofibromatosis type 1 phenotype while altering predominantly NF1 isoform type II

Neurofibromatosis type 1 results from loss-of-function NF1 pathogenic variants (PVs). Up to 30% of all NF1 PVs disrupt mRNA splicing, including deep intronic variants. Here, we retrospectively investigated the spectrum of NF1 deep intronic PVs in a cohort of 8,090 unrelated individuals from the University of Alabama at Birmingham (UAB) dataset with a molecularly confirmed neurofibromatosis type 1. All variants were identified through their effect on the NF1 transcript, followed by variant characterization at the DNA-level. A total of 68 distinct variants, which were ≥ 20 nucleotides away from the closest exon–intron junction, were identified in 2.5% unrelated individuals with NF1 (200/8,090). Nine different pathogenic splice variants, identified in 20 probands, led to exonization of different parts of intron 30 [23.2] or 31 [23a]. The two major NF1 transcript isoforms, distinguished by the absence (type I) or presence (type II) of the alternatively spliced cassette exon 31 [23a], are equally expressed in blood in control individuals without NF1 or NF1-affected individuals carrying their PV not in the introns flanking exon 31 [23a]. By fragment and cloning analysis we demonstrated that the exonization of intron 31 [23a] sequences due to deep intronic PV predominantly affects the NF1 isoform II. Seven additional (likely) pathogenic NF1 deep intronic variants not observed in the UAB dataset were found by classification of 36 variants identified by a literature search. Hence, the unique list of these 75 deep intronic (likely) PVs should be included in any comprehensive NF1 testing strategy.