ASD: a bioinformatics resource on alternative splicing

RegRNA 3.0: expanding regulatory RNA analysis with new features for motif, interaction, and annotation

Functional RNA molecules are crucial for biological processes from gene regulation to protein synthesis, and analyzing functional motifs and elements is essential for understanding RNA regulation. Building on RegRNA 1.0 and 2.0, we present RegRNA 3.0, a sophisticated meta-workflow that integrates 26 computational tools and 28 databases for annotation, enabling one-step and customizable RNA motif predictions. RegRNA streamlines multi-step analysis and enhances result interpretation with interactive visualizations and comprehensive reporting tools. When provided with an RNA sequence, RegRNA 3.0 generates predictions for RNA functional motifs, RNA interaction motifs, and comprehensive RNA annotations. Specifically, RNA functional motifs include core promoter elements, RNA decay, G-quadruplex, and 14 previous types. RNA interaction motifs include newly added RNA-ligand interactions and RNA-binding protein predictions, along with three previous types. RNA annotation includes RNA family classification, blood exosomes RNA, subcellular localizations, A-to-I editing events, modifications, and 3D structures, along with four previously supported features. RegRNA 3.0 accelerates gene regulation and RNA biology discoveries by offering a user-friendly platform for identifying and analyzing RNA motifs and interactions. The web interface has been improved for intuitive visualizations of predicted motifs and structures, with flexible download options in multiple formats. It is available at http://awi.cuhk.edu.cn/∼RegRNA/.

RNA-As-Graphs Motif Atlas—Dual Graph Library of RNA Modules and Viral Frameshifting-Element Applications

RNA motif classification is important for understanding structure/function connections and building phylogenetic relationships. Using our coarse-grained RNA-As-Graphs (RAG) representations, we identify recurrent dual graph motifs in experimentally solved RNA structures based on an improved search algorithm that finds and ranks independent RNA substructures. Our expanded list of 183 existing dual graph motifs reveals five common motifs found in transfer RNA, riboswitch, and ribosomal 5S RNA components. Moreover, we identify three motifs for available viral frameshifting RNA elements, suggesting a correlation between viral structural complexity and frameshifting efficiency. We further partition the RNA substructures into 1844 distinct submotifs, with pseudoknots and junctions retained intact. Common modules are internal loops and three-way junctions, and three submotifs are associated with riboswitches that bind nucleotides, ions, and signaling molecules. Together, our library of existing RNA motifs and submotifs adds to the growing universe of RNA modules, and provides a resource of structures and substructures for novel RNA design.

Splicing of exon 9a in FMR1 transcripts results in a truncated FMRP with altered subcellular distribution

FMRP is an RNA-binding protein, loss of which causes fragile X syndrome (FXS). FMRP has several isoforms resulted from alternative splicing (AS) of fragile X mental retardation 1 (FMR1) gene, but their biological functions are still poorly understood. In the analysis of alternatively spliced FMR1 transcripts in the blood cells from a patient with FXS-like phenotypes (normal CGG repeats and no mutation in coding sequence of FMR1), we identified three novel FMR1 transcripts that include a previously unidentified microexon (46 bp), terming the exon 9a. This microexon exists widely in unaffected individuals, inclusion of which introduces an in-frame termination codon. To address whether these exon 9a-containing transcripts could produce protein by evading nonsense-mediated decay (NMD), Western blot was used to analysis blood cell lysate from unaffected individuals and a 34 kDa protein that consistent in size with the molecular weight of the predicted truncated protein produced from mRNA with this microexon was found. Meanwhile, treatment of peripheral blood mononuclear cells with an inhibitor of NMD (Cycloheximide) did not result in significant increase in exon 9a-containing transcripts. Using confocal immunofluorescence, we found the truncated protein displayed both nuclear and cytoplasmic localization in HEK293T and HeLa cells due to lacking C-terminal domains including KH2, NES, and RGG, while the full-length FMRP protein mainly localized in the cytoplasm. Therefore, we hypothesize that the inclusion of this microexon to generate exon 9a-containing transcripts may regulate the normal functionality of FMRP, and the dysregulation of normal FMRP due to increased exon 9a-containing alternatively spliced transcripts in that patient may be associated with the manifestation of FXS phenotype.

T2DiACoD: A Gene Atlas of Type 2 Diabetes Mellitus Associated Complex Disorders

We performed integrative analysis of genes associated with type 2 Diabetes Mellitus (T2DM) associated complications by automated text mining with manual curation and also gene expression analysis from Gene Expression Omnibus. They were analysed for pathogenic or protective role, trends, interaction with risk factors, Gene Ontology enrichment and tissue wise differential expression. The database T2DiACoD houses 650 genes, and 34 microRNAs associated with T2DM complications. Seven genes AGER, TNFRSF11B, CRK, PON1, ADIPOQ, CRP and NOS3 are associated with all 5 complications. Several genes are studied in multiple years in all complications with high proportion in cardiovascular (75.8%) and atherosclerosis (51.3%). T2DM Patients' skeletal muscle tissues showed high fold change in differentially expressed genes. Among the differentially expressed genes, VEGFA is associated with several complications of T2DM. A few genes ACE2, ADCYAP1, HDAC4, NCF1, NFE2L2, OSM, SMAD1, TGFB1, BDNF, SYVN1, TXNIP, CD36, CYP2J2, NLRP3 with details of protective role are catalogued. Obesity is clearly a dominant risk factor interacting with the genes of T2DM complications followed by inflammation, diet and stress to variable extents. This information emerging from the integrative approach used in this work could benefit further therapeutic approaches. The T2DiACoD is available at www.http://t2diacod.igib.res.in/ .

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Computational Identification of Tissue-Specific Splicing Regulatory Elements in Human Genes from RNA-Seq Data

Alternative splicing is a vital process for regulating gene expression and promoting proteomic diversity. It plays a key role in tissue-specific expressed genes. This specificity is mainly regulated by splicing factors that bind to specific sequences called splicing regulatory elements (SREs). Here, we report a genome-wide analysis to study alternative splicing on multiple tissues, including brain, heart, liver, and muscle. We propose a pipeline to identify differential exons across tissues and hence tissue-specific SREs. In our pipeline, we utilize the DEXSeq package along with our previously reported algorithms. Utilizing the publicly available RNA-Seq data set from the Human BodyMap project, we identified 28,100 differentially used exons across the four tissues. We identified tissue-specific exonic splicing enhancers that overlap with various previously published experimental and computational databases. A complicated exonic enhancer regulatory network was revealed, where multiple exonic enhancers were found across multiple tissues while some were found only in specific tissues. Putative combinatorial exonic enhancers and silencers were discovered as well, which may be responsible for exon inclusion or exclusion across tissues. Some of the exonic enhancers are found to be co-occurring with multiple exonic silencers and vice versa, which demonstrates a complicated relationship between tissue-specific exonic enhancers and silencers.

MiasDB: A Database of Molecular Interactions Associated with Alternative Splicing of Human Pre-mRNAs

Alternative splicing (AS) is pervasive in human multi-exon genes and is a major contributor to expansion of the transcriptome and proteome diversity. The accurate recognition of alternative splice sites is regulated by information contained in networks of protein-protein and protein-RNA interactions. However, the mechanisms leading to splice site selection are not fully understood. Although numerous databases have been built to describe AS, molecular interaction databases associated with AS have only recently emerged. In this study, we present a new database, MiasDB, that provides a description of molecular interactions associated with human AS events. This database covers 938 interactions between human splicing factors, RNA elements, transcription factors, kinases and modified histones for 173 human AS events. Every entry includes the interaction partners, interaction type, experimental methods, AS type, tissue specificity or disease-relevant information, a simple description of the functionally tested interaction in the AS event and references. The database can be queried easily using a web server (http://47.88.84.236/Miasdb). We display some interaction figures for several genes. With this database, users can view the regulation network describing AS events for 12 given genes.

A Reverse-Genetics Mutational Analysis of the Barley HvDWARF Gene Results in Identification of a Series of Alleles and Mutants with Short Stature of Various Degree and Disturbance in BR Biosynthesis Allowing a New Insight into the Process

Brassinosteroids (BRs) are plant steroid hormones, regulating a broad range of physiological processes. The largest amount of data related with BR biosynthesis has been gathered in Arabidopsis thaliana, however understanding of this process is far less elucidated in monocot crops. Up to now, only four barley genes implicated in BR biosynthesis have been identified. Two of them, HvDWARF and HvBRD, encode BR-6-oxidases catalyzing biosynthesis of castasterone, but their relation is not yet understood. In the present study, the identification of the HvDWARF genomic sequence, its mutational and functional analysis and characterization of new mutants are reported. Various types of mutations located in different positions within functional domains were identified and characterized. Analysis of their impact on phenotype of the mutants was performed. The identified homozygous mutants show reduced height of various degree and disrupted skotomorphogenesis. Mutational analysis of the HvDWARF gene with the “reverse genetics” approach allowed for its detailed functional analysis at the level of protein functional domains. The HvDWARF gene function and mutants’ phenotypes were also validated by measurement of endogenous BR concentration. These results allowed a new insight into the BR biosynthesis in barley.

ATtRACT-a database of RNA-binding proteins and associated motifs

RNA-binding proteins (RBPs) play a crucial role in key cellular processes, including RNA transport, splicing, polyadenylation and stability. Understanding the interaction between RBPs and RNA is key to improve our knowledge of RNA processing, localization and regulation in a global manner. Despite advances in recent years, a unified non-redundant resource that includes information on experimentally validated motifs, RBPs and integrated tools to exploit this information is lacking. Here, we developed a database named ATtRACT (available athttp://attract.cnic.es) that compiles information on 370 RBPs and 1583 RBP consensus binding motifs, 192 of which are not present in any other database. To populate ATtRACT we (i) extracted and hand-curated experimentally validated data from CISBP-RNA, SpliceAid-F, RBPDB databases, (ii) integrated and updated the unavailable ASD database and (iii) extracted information from Protein-RNA complexes present in Protein Data Bank database through computational analyses. ATtRACT provides also efficient algorithms to search a specific motif and scan one or more RNA sequences at a time. It also allows discoveringde novomotifs enriched in a set of related sequences and compare them with the motifs included in the database.Database URL:http:// attract. cnic. es.

The Functional Human C-Terminome

All translated proteins end with a carboxylic acid commonly called the C-terminus. Many short functional sequences (minimotifs) are located on or immediately proximal to the C-terminus. However, information about the function of protein C-termini has not been consolidated into a single source. Here, we built a new "C-terminome" database and web system focused on human proteins. Approximately 3,600 C-termini in the human proteome have a minimotif with an established molecular function. To help evaluate the function of the remaining C-termini in the human proteome, we inferred minimotifs identified by experimentation in rodent cells, predicted minimotifs based upon consensus sequence matches, and predicted novel highly repetitive sequences in C-termini. Predictions can be ranked by enrichment scores or Gene Evolutionary Rate Profiling (GERP) scores, a measurement of evolutionary constraint. By searching for new anchored sequences on the last 10 amino acids of proteins in the human proteome with lengths between 3-10 residues and up to 5 degenerate positions in the consensus sequences, we have identified new consensus sequences that predict instances in the majority of human genes. All of this information is consolidated into a database that can be accessed through a C-terminome web system with search and browse functions for minimotifs and human proteins. A known consensus sequence-based predicted function is assigned to nearly half the proteins in the human proteome. Weblink: http://cterminome.bio-toolkit.com.

Data in support of a functional analysis of splicing mutations in the IDS gene and the use of antisense oligonucleotides to exploit an alternative therapy for MPS II

This data article contains insights into the methodology used for the analysis of three exonic mutations altering the splicing of the IDS gene: c.241C>T, c.257C>T and c.1122C>T.

We have performed splicing assays for the wild-type and mutant minigenes corresponding to these substitutions. In addition, bioinformatic predictions of splicing regulatory sequence elements as well as RNA interference and overexpression experiments were conducted.

The interpretation of these data and further extensive experiments into the analysis of these three mutations and also into the methodology applied to correct one of them can be found in “Functional analysis of splicing mutations in the IDS gene and the use of antisense oligonucleotides to exploit an alternative therapy for MPS II” Matos et al. (2015) [1].

Functional analysis of splicing mutations in the IDS gene and the use of antisense oligonucleotides to exploit an alternative therapy for MPS II

Mucopolysaccharidosis II is a lysosomal storage disorder caused by mutations in the IDS gene, including exonic alterations associated with aberrant splicing. In the present work, cell-based splicing assays were performed to study the effects of two splicing mutations in exon 3 of IDS, i.e., c.241C>T and c.257C>T, whose presence activates a cryptic splice site in exon 3 and one in exon 8, i.e., c.1122C>T that despite being a synonymous mutation is responsible for the creation of a new splice site in exon 8 leading to a transcript shorter than usual. Mutant minigene analysis and overexpression assays revealed that SRSF2 and hnRNP E1 might be involved in the use and repression of the constitutive 3' splice site of exon 3 respectively. For the c.1122C>T the use of antisense therapy to correct the splicing defect was explored, but transfection of patient fibroblasts with antisense morpholino oligonucleotides (n=3) and a locked nucleic acid failed to abolish the abnormal transcript; indeed, it resulted in the appearance of yet another aberrant splicing product. Interestingly, the oligonucleotides transfection in control fibroblasts led to the appearance of the aberrant transcript observed in patients' cells after treatment, which shows that the oligonucleotides are masking an important cis-acting element for 5' splice site regulation of exon 8. These results highlight the importance of functional studies for understanding the pathogenic consequences of mis-splicing and highlight the difficulty in developing antisense therapies involving gene regions under complex splicing regulation.

CoSREM: a graph mining algorithm for the discovery of combinatorial splicing regulatory elements

BACKGROUND

Alternative splicing (AS) is a post-transcriptional regulatory mechanism for gene expression regulation. Splicing decisions are affected by the combinatorial behavior of different splicing factors that bind to multiple binding sites in exons and introns. These binding sites are called splicing regulatory elements (SREs). Here we develop CoSREM (Combinatorial SRE Miner), a graph mining algorithm to discover combinatorial SREs in human exons. Our model does not assume a fixed length of SREs and incorporates experimental evidence as well to increase accuracy. CoSREM is able to identify sets of SREs and is not limited to SRE pairs as are current approaches.

RESULTS

We identified 37 SRE sets that include both enhancer and silencer elements. We show that our results intersect with previous results, including some that are experimental. We also show that the SRE set GGGAGG and GAGGAC identified by CoSREM may play a role in exon skipping events in several tumor samples. We applied CoSREM to RNA-Seq data for multiple tissues to identify combinatorial SREs which may be responsible for exon inclusion or exclusion across tissues.

CONCLUSION

The new algorithm can identify different combinations of splicing enhancers and silencers without assuming a predefined size or limiting the algorithm to find only pairs of SREs. Our approach opens new directions to study SREs and the roles that AS may play in diseases and tissue specificity.

Survey of Programs Used to Detect Alternative Splicing Isoforms from Deep Sequencing Data In Silico

Next-generation sequencing techniques have been rapidly emerging. However, the massive sequencing reads hide a great deal of unknown important information. Advances have enabled researchers to discover alternative splicing (AS) sites and isoforms using computational approaches instead of molecular experiments. Given the importance of AS for gene expression and protein diversity in eukaryotes, detecting alternative splicing and isoforms represents a hot topic in systems biology and epigenetics research. The computational methods applied to AS prediction have improved since the emergence of next-generation sequencing. In this study, we introduce state-of-the-art research on AS and then compare the research methods and software tools available for AS based on next-generation sequencing reads. Finally, we discuss the prospects of computational methods related to AS.

Identifying splicing regulatory elements with de Bruijn graphs

Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.

ASPicDB: a database web tool for alternative splicing analysis

Alternative splicing (AS) is a basic molecular phenomenon that increases the functional complexity of higher eukaryotic transcriptomes. Indeed, through AS individual gene loci can generate multiple RNAs from the same pre-mRNA. AS has been investigated in a variety of clinical and pathological studies, such as the transcriptome regulation in cancer. In human, recent works based on massive RNA sequencing indicate that >95 % of pre-mRNAs are processed to yield multiple transcripts. Given the biological relevance of AS, several computational efforts have been done leading to the implementation of novel algorithms and specific specialized databases. Here we describe the web application ASPicDB that allows the recovery of detailed biological information about the splicing mechanism. ASPicDB provides powerful querying systems to interrogate AS events at gene, transcript, and protein levels. Finally, ASPicDB includes web visualization instruments to browse and export results for further off-line analyses.

Molecular and computational analyses of genes involved in mannose 6-phosphate independent trafficking

The newly-synthesized lysosomal enzymes travel to the trans-Golgi network (TGN) and are then driven to the acidic organelle. While the best-known pathway for TGN-to-endosome transport is the delivery of soluble hydrolases by the M6P receptors (MPRs), additional pathways do exist, as showed by the identification of two alternative receptors: LIMP-2, implicated in the delivery of β-glucocerebrosidase; and sortilin, involved in the transport of the sphingolipid activator proteins prosaposin and GM2AP, acid sphingomyelinase and cathepsins D and H. Disruption of the intracellular transport and delivery pathways to the lysosomes may result in lysosomal dysfunction, predictably leading to a range of clinical manifestations of lysosomal storage diseases. However, for a great percentage of patients presenting such manifestations, no condition is successfully diagnosed. To analyse if, in this group, phenotypes could be determined by impairments in the known M6P-independent receptors, we screened the genes that encode for LIMP-2 and sortilin. No pathogenic mutations were identified. Other approaches will be needed to clarify whether sortilin dysfunction may cause disease.

Alternative splicing of TAF6: downstream transcriptome impacts and upstream RNA splice control elements

The TAF6δ pathway of apoptosis can dictate life versus death decisions independently of the status of p53 tumor suppressor. TAF6δ is an inducible pro-apoptotic subunit of the general RNA polymerase II (Pol II) transcription factor TFIID. Alternative splice site choice of TAF6δ has been shown to be a pivotal event in triggering death via the TAF6δ pathway, yet nothing is currently known about the mechanisms that promote TAF6δ splicing. Furthermore the transcriptome impact of the gain of function of TAF6δ versus the loss of function of the major TAF6α splice form remains undefined. Here we employ comparative microarray analysis to show that TAF6δ drives a transcriptome profile distinct from that resulting from depletion of TAF6α. To define the cis-acting RNA elements responsible for TAF6δ alternative splicing we performed a mutational analysis of a TAF6 minigene system. The data point to several new RNA elements that can modulate TAF6δ and also reveal a role for RNA secondary structure in the selection of TAF6δ.

SeeSite: Characterizing Relationships between Splice Junctions and Splicing Enhancers

RNA splicing is a cellular process driven by the interaction between numerous regulatory sequences and binding sites, however, such interactions have been primarily explored by laboratory methods since computational tools largely ignore the relationship between different splicing elements. Current computational methods identify either splice sites or other regulatory sequences, such as enhancers and silencers. We present a novel approach for characterizing co-occurring relationships between splice site motifs and splicing enhancers. Our approach relies on an efficient algorithm for approximately solving Consensus Sequence with Outliers , an NP-complete string clustering problem. In particular, we give an algorithm for this problem that outputs near-optimal solutions in polynomial time. To our knowledge, this is the first formulation and computational attempt for detecting co-occurring sequence elements in RNA sequence data. Further, we demonstrate that SeeSite is capable of showing that certain ESEs are preferentially associated with weaker splice sites, and that there exists a co-occurrence relationship with splice site motifs.

Sequence-specific flexibility organization of splicing flanking sequence and prediction of splice sites in the human genome

More and more reported results of nucleosome positioning and histone modifications showed that DNA structure play a well-established role in splicing. In this study, a set of DNA geometric flexibility parameters originated from molecular dynamics (MD) simulations were introduced to discuss the structure organization around splice sites at the DNA level. The obtained profiles of specific flexibility/stiffness around splice sites indicated that the DNA physical-geometry deformation could be used as an alternative way to describe the splicing junction region. In combination with structural flexibility as discriminatory parameter, we developed a hybrid computational model for predicting potential splicing sites. And the better prediction performance was achieved when the benchmark dataset evaluated. Our results showed that the mechanical deformability character of a splice junction is closely correlated with both the splice site strength and structural information in its flanking sequences.

Review: Alternative Splicing (AS) of Genes As An Approach for Generating Protein Complexity

Prior to the completion of the human genome project, the human genome was thought to have a greater number of genes as it seemed structurally and functionally more complex than other simpler organisms. This along with the belief of “one gene, one protein”, were demonstrated to be incorrect. The inequality in the ratio of gene to protein formation gave rise to the theory of alternative splicing (AS). AS is a mechanism by which one gene gives rise to multiple protein products. Numerous databases and online bioinformatic tools are available for the detection and analysis of AS. Bioinformatics provides an important approach to study mRNA and protein diversity by various tools such as expressed sequence tag (EST) sequences obtained from completely processed mRNA. Microarrays and deep sequencing approaches also aid in the detection of splicing events. Initially it was postulated that AS occurred only in about 5% of all genes but was later found to be more abundant. Using bioinformatic approaches, the level of AS in human genes was found to be fairly high with 35-59% of genes having at least one AS form. Our ability to determine and predict AS is important as disorders in splicing patterns may lead to abnormal splice variants resulting in genetic diseases. In addition, the diversity of proteins produced by AS poses a challenge for successful drug discovery and therefore a greater understanding of AS would be beneficial.

Novel GUCA1A mutations suggesting possible mechanisms of pathogenesis in cone, cone-rod, and macular dystrophy patients

Here, we report two novel GUCA1A (the gene for guanylate cyclase activating protein 1) mutations identified in unrelated Spanish families affected by autosomal dominant retinal degeneration (adRD) with cone and rod involvement. All patients from a three-generation adRD pedigree underwent detailed ophthalmic evaluation. Total genome scan using single-nucleotide polymorphisms and then the linkage analysis were undertaken on the pedigree. Haplotype analysis revealed a 55.37 Mb genomic interval cosegregating with the disease phenotype on chromosome 6p21.31-q15. Mutation screening of positional candidate genes found a heterozygous transition c.250C>T in exon 4 of GUCA1A, corresponding to a novel mutation p.L84F. A second missense mutation, c.320T>C (p.I107T), was detected by screening of the gene in a Spanish patients cohort. Using bioinformatics approach, we predicted that either haploinsufficiency or dominant-negative effect accompanied by creation of a novel function for the mutant protein is a possible mechanism of the disease due to c.250C>T and c.320T>C. Although additional functional studies are required, our data in relation to the c.250C>T mutation open the possibility that transacting factors binding to de novo created recognition site resulting in formation of aberrant splicing variant is a disease model which may be more widespread than previously recognized as a mechanism causing inherited RD.

In vitro correction of a pseudoexon-generating deep intronic mutation in LGMD2A by antisense oligonucleotides and modified small nuclear RNAs

Limb-girdle muscular dystrophy type 2A (LGMD2A) is the most frequent autosomal recessive muscular dystrophy. It is caused by mutations in the calpain-3 (CAPN3) gene. The majority of the mutations described to date are located in the coding sequence of the gene. However, it is estimated that 25% of the mutations are present at exon-intron boundaries and modify the pre-mRNA splicing of the CAPN3 transcript. We have previously described the first deep intronic mutation in the CAPN3 gene: c.1782+1072G>C mutation. This mutation causes the pseudoexonization of an intronic sequence of the CAPN3 gene in the mature mRNA. In the present work, we show that the point mutation generates the inclusion of the pseudoexon in the mRNA using a minigene assay. In search of a treatment that restores normal splicing, splicing modulation was induced by RNA-based strategies, which included antisense oligonucleotides and modified small-nuclear RNAs. The best effect was observed with antisense sequences, which induced pseudoexon skipping in both HeLa cells cotransfected with mutant minigene and in fibroblasts from patients. Finally, transfection of antisense sequences and siRNA downregulation of serine/arginine-rich splicing factor 1 (SRSF1) indicate that binding of this factor to splicing enhancer sequences is involved in pseudoexon activation.

Mechanisms of the androgen receptor splicing in prostate cancer cells

Prostate tumors develop resistance to androgen deprivation therapy (ADT) by multiple mechanisms, one of which is to express constitutively active androgen receptor (AR) splice variants lacking the ligand-binding domain. AR splice variant 7 (AR-V7, also termed AR3) is the most abundantly expressed variant that drives prostate tumor progression under ADT conditions. However, the molecular mechanism by which AR-V7 is generated remains unclear. In this manuscript, we demonstrated that RNA splicing of AR-V7 in response to ADT was closely associated with AR gene transcription initiation and elongation rates. Enhanced AR gene transcription by ADT provides a prerequisite condition that further increases the interactions between AR pre-mRNA and splicing factors. Under ADT conditions, recruitment of several RNA splicing factors to the 3' splicing site for AR-V7 was increased. We identified two RNA splicing enhancers and their binding proteins (U2AF65 and ASF/SF2) that had critical roles in splicing AR pre-mRNA into AR-V7. These data indicate that ADT-induced AR gene transcription rate and splicing factor recruitment to AR pre-mRNA contribute to the enhanced AR-V7 levels in prostate cancer cells.

Cardiac glycosides correct aberrant splicing of IKBKAP-encoded mRNA in familial dysautonomia derived cells by suppressing expression of SRSF3

The ability to modulate the production of the wild-type transcript in cells bearing the splice-altering familial dysautonomia (FD) causing mutation in the IKBKAP gene prompted a study of the impact of a panel of pharmaceuticals on the splicing of this transcript, which revealed the ability of the cardiac glycoside digoxin to increase the production of the wild-type, exon-20-containing, IKBKAP-encoded transcript and the full-length IκB-kinase-complex-associated protein in FD-derived cells. Characterization of the cis elements and trans factors involved in the digoxin-mediated effect on splicing reveals that this response is dependent on an SRSF3 binding site(s) located in the intron 5' of the alternatively spliced exon and that digoxin mediates its effect by suppressing the level of the SRSF3 protein. Characterization of the digoxin-mediated effect on the RNA splicing process was facilitated by the identification of several RNA splicing events in which digoxin treatment mediates the enhanced inclusion of exonic sequence. Moreover, we demonstrate the ability of digoxin to impact the splicing process in neuronal cells, a cell type profoundly impacted by FD. This study represents the first demonstration that digoxin possesses splice-altering capabilities that are capable of reversing the impact of the FD-causing mutation. These findings support the clinical evaluation of the impact of digoxin on the FD patient population.

Hypoxia-dependent mRNA expression pattern of splicing factor YT521 and its impact on oncological important target gene expression

The ubiquitously expressed splicing factor YT521 (YTHDC1) is characterized by alternatively spliced isoforms with regulatory impact on cancer-associated gene expression. Our recent findings account for the prognostic significance of YT521 in endometrial cancer. In this study, we investigated the hypoxia-dependency of YT521 expression as well as its differential isoform activities on oncological important target genes. YT521's potential regulatory influence on splicing was investigated by a minigene assay for the specific target gene CD44. Functional splicing analysis was performed by YT521 knock-down or overexpression, respectively. In addition, YT521 expression was determined under hypoxia. The two protein-generating YT521 mRNA isoforms 1 and 2 caused a comparable, specific induction of CD44v alternative splicing (P < 0.01). In a number of oncological target genes, YT521 upregulation significantly altered BRCA2 expression pattern, while YT521 knock-down created a significant regulatory impact on PGR expression, respectively. Hypoxia induced a specific switch towards the processing of two non-protein-coding mRNA variants, of which one is described for the first time in this study. The presented study underlines the comparable regulatory potential of both YT521 isoforms 1 and 2, on the investigated target genes in vivo and in vitro. Hypoxia induces a specific switch in YT521 expression pattern towards the two non-protein coding mRNA variants, the already characterized isoform 3 and the newly discovered exon 8-skipping isoform. The altered YT521 alternative splicing is functionally coupled with nonsense-mediated decay and can be interpreted as regulated unproductive splicing and transcription with consecutive impact on the processing of specific cancer-associated genes, such as BRCA2 and PGR.

A broadly applicable high-throughput screening strategy identifies new regulators of Dlg4 (Psd-95) alternative splicing

Most mammalian genes produce multiple mRNA isoforms derived from alternative pre-mRNA splicing, with each alternative exon controlled by a complex network of regulatory factors. The identification of these regulators can be laborious and is usually carried out one factor at a time. We have developed a broadly applicable high-throughput screening method that simultaneously identifies multiple positive and negative regulators of a particular exon. Two minigene reporters were constructed: One produces green fluorescent protein (GFP) from the mRNA including an exon, and red fluorescent protein (RFP) from the mRNA lacking the exon; the other switches these fluorescent products of exon inclusion and exclusion. Combining results from these two reporters eliminates many false positives and greatly enriches for true splicing regulators. After extensive optimization of this method, we performed a gain-of-function screen of 15,779 cDNA clones and identified 40 genes affecting exon 18 of Discs large homolog 4 (Dlg4; also known as post-synaptic density protein 95 [Psd-95]). We confirmed that 28 of the 34 recoverable clones alter reporter splicing in RT-PCR assays. Remarkably, 18 of the identified genes encode splicing factors or RNA binding proteins, including PTBP1, a previously identified regulator of this exon. Loss-of-function experiments examining endogenous Dlg4 transcripts validated the effects of five of eight genes tested in independent cell lines, and two genes were further confirmed to regulate Dlg4 splicing in primary neurons. These results identify multiple new regulators of Dlg4 splicing, and validate an approach to isolating splicing regulators for almost any cassette exon from libraries of cDNAs, shRNAs, or small molecules.

Alternative splicing for diseases, cancers, drugs, and databases

Alternative splicing is a major diversification mechanism in the human transcriptome and proteome. Several diseases, including cancers, have been associated with dysregulation of alternative splicing. Thus, correcting alternative splicing may restore normal cell physiology in patients with these diseases. This paper summarizes several alternative splicing-related diseases, including cancers and their target genes. Since new cancer drugs often target spliceosomes, several clinical drugs and natural products or their synthesized derivatives were analyzed to determine their effects on alternative splicing. Other agents known to have modulating effects on alternative splicing during therapeutic treatment of cancer are also discussed. Several commonly used bioinformatics resources are also summarized.

Prediction of mutant mRNA splice isoforms by information theory-based exon definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon-skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R(i), which measures binding-site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R(i),total) is the sum of the R(i) values of its acceptor and donor splice sites, adjusted for the self-information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR(i,total)) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis (http://splice.uwo.ca) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.

An enhanced computational platform for investigating the roles of regulatory RNA and for identifying functional RNA motifs

Background

Functional RNA molecules participate in numerous biological processes, ranging from gene regulation to protein synthesis. Analysis of functional RNA motifs and elements in RNA sequences can obtain useful information for deciphering RNA regulatory mechanisms. Our previous work, RegRNA, is widely used in the identification of regulatory motifs, and this work extends it by incorporating more comprehensive and updated data sources and analytical approaches into a new platform.

Methods and results

An integrated web-based system, RegRNA 2.0, has been developed for comprehensively identifying the functional RNA motifs and sites in an input RNA sequence. Numerous data sources and analytical approaches are integrated, and several types of functional RNA motifs and sites can be identified by RegRNA 2.0: (i) splicing donor/acceptor sites; (ii) splicing regulatory motifs; (iii) polyadenylation sites; (iv) ribosome binding sites; (v) rho-independent terminator; (vi) motifs in mRNA 5'-untranslated region (5'UTR) and 3'UTR; (vii) AU-rich elements; (viii) C-to-U editing sites; (ix) riboswitches; (x) RNA cis-regulatory elements; (xi) transcriptional regulatory motifs; (xii) user-defined motifs; (xiii) similar functional RNA sequences; (xiv) microRNA target sites; (xv) non-coding RNA hybridization sites; (xvi) long stems; (xvii) open reading frames; (xviii) related information of an RNA sequence. User can submit an RNA sequence and obtain the predictive results through RegRNA 2.0 web page.

Conclusions

RegRNA 2.0 is an easy to use web server for identifying regulatory RNA motifs and functional sites. Through its integrated user-friendly interface, user is capable of using various analytical approaches and observing results with graphical visualization conveniently. RegRNA 2.0 is now available at http://regrna2.mbc.nctu.edu.tw.

A prevalent and three novel mutations in CYP11B1 gene identified in Chinese patients with 11-beta hydroxylase deficiency

11β-Hydroxylase deficiency (11β-OHD), caused by CYP11B1 mutations, is characterized by hyporeninemic, hypokalemic hypertension and hyperandrogenism. We identified a prevalent and three novel mutations of CYP11B1 gene in nine patients with classic 11β-OHD.

SUBJECTS AND METHODS

Nine patients with 11β-OHD from unrelated families were recruited. The complications of 11β-OHD occurred in three patients who never received glucocorticoid treatment. CYP11B1 gene was sequenced and 11β-hydroxylase enzymatic activities were assessed in vitro. A haplotype analysis was performed to determine a common ancestor for those subjects who carried the same p.R454C mutation.

RESULTS

CYP11B1 gene mutations were identified in all patients, with a prevalent (p.R454C) and three novel mutations (p.V148G, IVS7-9C>A, c.1359_1360insG). The p.R141X, p.V148G, c.1359_1360insG and p.R454C mutations retained 4.9%, 3.9%, 3.7%, 4.5% of residual enzymatic activity, respectively. Five of nine patients carried p.R454C mutation, which was only reported in Chinese 11OHD patients. Haplotype analysis showed that this mutation might be inherited from a common ancestor.

CONCLUSION

The enzymatic activities for p.R141X, p.V148G, c.1359_1360insG and p.R454C mutants were almost completely abolished, which corresponds to classic form of 11β-OHD. The observations of a prevalent mutation and three novel mutations might have potential clinical utility for genetic counseling and prenatal diagnosis in Chinese 11β-OHD patients.

A novel splice site mutation of CDHR1 in a consanguineous Israeli Christian Arab family segregating autosomal recessive cone-rod dystrophy

PURPOSE

To investigate the genetic basis for autosomal recessive cone-rod dystrophy in a consanguineous Israeli Christian Arab family.

METHODS

Patients underwent a detailed ophthalmic examination, including funduscopy, electroretinography (ERG), visual field testing, and optical coherence tomography. Genome-wide homozygosity mapping using a single nucleotide polymorphism array was performed to identify homozygous regions shared between the two affected individuals. Mutation screening of the underlying gene was performed with direct sequencing. In silico analysis was used to predict the effect of the mutation on splicing.

RESULTS

The family included two affected individuals. Clinical findings included progressive deterioration of visual acuity, photophobia, defective color vision, loss of central visual fields, pigmentary deposits localized mainly in the peripheral retina, a thinned and atrophic macular region, retinal vessel attenuation, absent ERG cone responses, and reduced ERG rod responses. Homozygosity mapping revealed several homozygous intervals shared among the affected individuals. One, a 12Mb interval on chromosome 10, included the CDHR1 gene. Direct sequencing revealed a single base transversion, c.1485+2T>G, located in the conserved donor splice site of Intron 13. This mutation cosegregated with the disease in the family, and was not detected in 208 Israeli Christian Arab control chromosomes. In silico analysis predicted that this mutation eliminates the Intron 13 donor splice site.

CONCLUSIONS

Only three distinct pathogenic mutations of CDHR1 have been reported to date in patients with autosomal recessive retinal degeneration. Here we report a novel splice site mutation of CDHR1, c.1485+2T>G, underlying autosomal recessive cone-rod dystrophy in a consanguineous Israeli Christian Arab family. This report expands the spectrum of pathogenic mutations of the CDHR1 gene.

Environmental enrichment upregulates micro-RNA-183 and alters acetylcholinesterase splice variants to reduce anxiety-like behavior in the little Indian field mouse (Mus booduga)

Environmental enrichment (EE) has an influential role in reducing behavioral reactivity to stress. We previously observed that EE reduces the anxiety-like behavior in the field mouse Mus booduga accompanied by a reduction in the expression of molecules involved in the stress pathway. In this study, we demonstrate the effect of different housing condition on regulation of micro-RNA-183-SC35-mediated splicing of acetylcholinesterase (AChE). Adult male M. booduga were captured from an agricultural field and housed under nonenriched standard conditions (SC) for 7 days and considered as directly from the wild (DW). On day 8, individuals were randomly assigned to three groups; DW, SC, and EE. The DW group's anxiety-like behavior was assessed in the elevated plus maze (EPM) and open field test (OFT). The SC and EE groups were transferred to their respective conditions and housed for another 30 days. The mice housed in EE showed less anxiety-like behavior on EPM and in OFT compared with DW and SC mice. Interestingly, miR-183 expression was increased following exposure to EPM in EE mice but not in SC mice. Subsequently, the upregulated miR-183 expression suppresses the SC35 expression and shifting of splicing from AChE-S (synaptic) to AChE-R (read-through) form, whereas standard housing condition downregulate miR-183 and induces the splicing of AChE. The upregulated AChE-R form possibly terminates ACh transmission, which is reflected in the level of anxiety-like behavior. Overall, the present study suggests that EE effectively regulates the miR-183 pathway to reduce anxiety-like behavior.

SpliceAid-F: a database of human splicing factors and their RNA-binding sites

A comprehensive knowledge of all the factors involved in splicing, both proteins and RNAs, and of their interaction network is crucial for reaching a better understanding of this process and its functions. A large part of relevant information is buried in the literature or collected in various different databases. By hand-curated screenings of literature and databases, we retrieved experimentally validated data on 71 human RNA-binding splicing regulatory proteins and organized them into a database called ‘SpliceAid-F’ (http://www.caspur.it/SpliceAidF/). For each splicing factor (SF), the database reports its functional domains, its protein and chemical interactors and its expression data. Furthermore, we collected experimentally validated RNA–SF interactions, including relevant information on the RNA-binding sites, such as the genes where these sites lie, their genomic coordinates, the splicing effects, the experimental procedures used, as well as the corresponding bibliographic references. We also collected information from experiments showing no RNA–SF binding, at least in the assayed conditions.

In total, SpliceAid-F contains 4227 interactions, 2590 RNA-binding sites and 1141 ‘no-binding’ sites, including information on cellular contexts and conditions where binding was tested.

The data collected in SpliceAid-F can provide significant information to explain an observed splicing pattern as well as the effect of mutations in functional regulatory elements.

RNA sequencing reveals a diverse and dynamic repertoire of the Xenopus tropicalis transcriptome over development

The Xenopus embryo has provided key insights into fate specification, the cell cycle, and other fundamental developmental and cellular processes, yet a comprehensive understanding of its transcriptome is lacking. Here, we used paired end RNA sequencing (RNA-seq) to explore the transcriptome of Xenopus tropicalis in 23 distinct developmental stages. We determined expression levels of all genes annotated in RefSeq and Ensembl and showed for the first time on a genome-wide scale that, despite a general state of transcriptional silence in the earliest stages of development, approximately 150 genes are transcribed prior to the midblastula transition. In addition, our splicing analysis uncovered more than 10,000 novel splice junctions at each stage and revealed that many known genes have additional unannotated isoforms. Furthermore, we used Cufflinks to reconstruct transcripts from our RNA-seq data and found that ∼13.5% of the final contigs are derived from novel transcribed regions, both within introns and in intergenic regions. We then developed a filtering pipeline to separate protein-coding transcripts from noncoding RNAs and identified a confident set of 6686 noncoding transcripts in 3859 genomic loci. Since the current reference genome, XenTro3, consists of hundreds of scaffolds instead of full chromosomes, we also performed de novo reconstruction of the transcriptome using Trinity and uncovered hundreds of transcripts that are missing from the genome. Collectively, our data will not only aid in completing the assembly of the Xenopus tropicalis genome but will also serve as a valuable resource for gene discovery and for unraveling the fundamental mechanisms of vertebrate embryogenesis.

A novel MCPH1 isoform complements the defective chromosome condensation of human MCPH1-deficient cells

Biallelic mutations in MCPH1 cause primary microcephaly (MCPH) with the cellular phenotype of defective chromosome condensation. MCPH1 encodes a multifunctional protein that notably is involved in brain development, regulation of chromosome condensation, and DNA damage response. In the present studies, we detected that MCPH1 encodes several distinct transcripts, including two major forms: full-length MCPH1 (MCPH1-FL) and a second transcript lacking the six 3′ exons (MCPH1Δe9–14). Both variants show comparable tissue-specific expression patterns, demonstrate nuclear localization that is mediated independently via separate NLS motifs, and are more abundant in certain fetal than adult organs. In addition, the expression of either isoform complements the chromosome condensation defect found in genetically MCPH1-deficient or MCPH1 siRNA-depleted cells, demonstrating a redundancy of both MCPH1 isoforms for the regulation of chromosome condensation. Strikingly however, both transcripts are regulated antagonistically during cell-cycle progression and there are functional differences between the isoforms with regard to the DNA damage response; MCPH1-FL localizes to phosphorylated H2AX repair foci following ionizing irradiation, while MCPH1Δe9–14 was evenly distributed in the nucleus. In summary, our results demonstrate here that MCPH1 encodes different isoforms that are differentially regulated at the transcript level and have different functions at the protein level.

The role of SH3BP2 in the pathophysiology of cherubism

Cherubism is a rare bone dysplasia that is characterized by symmetrical bone resorption limited to the jaws. Bone lesions are filled with soft fibrous giant cell-rich tissue that can expand and cause severe facial deformity. The disorder typically begins in children at ages of 2-5 years and the bone resorption and facial swelling continues until puberty; in most cases the lesions regress spontaneously thereafter. Most patients with cherubism have germline mutations in the gene encoding SH3BP2, an adapter protein involved in adaptive and innate immune response signaling. A mouse model carrying a Pro416Arg mutation in SH3BP2 develops osteopenia and expansile lytic lesions in bone and some soft tissue organs. In this review we discuss the genetics of cherubism, the biological functions of SH3BP2 and the analysis of the mouse model. The data suggest that the underlying cause for cherubism is a systemic autoinflammatory response to physiologic challenges despite the localized appearance of bone resorption and fibrous expansion to the jaws in humans.

Statistical and Computational Methods for High-Throughput Sequencing Data Analysis of Alternative Splicing

The burgeoning field of high-throughput sequencing significantly improves our ability to understand the complexity of transcriptomes. Alternative splicing, as one of the most important driving forces for transcriptome diversity, can now be studied at an unprecedent resolution. Efficient and powerful computational and statistical methods are in urgent need to facilitate the characterization and quantification of alternative splicing events. Here we discuss methods in splice junction read mapping, and methods in exon-centric or isoform-centric quantification of alternative splicing. In addition, we discuss HITS-CLIP and splicing QTL analyses which are novel high-throughput sequencing based approaches in the dissection of splicing regulation.

Detection, annotation and visualization of alternative splicing from RNA-Seq data with SplicingViewer

Alternative splicing is a crucial mechanism by which diverse gene products can be generated from a limited number of genes, and is thought to be involved in complex orchestration of eukaryotic gene expression. Next-generation sequencing technologies, with reduced time and cost, provide unprecedented opportunities for deep interrogation of alternative splicing at the genome-wide scale. In this study, an integrated software SplicingViewer has been developed for unambiguous detection, annotation and visualization of splice junctions and alternative splicing events from RNA-Seq data. Specifically, it allows easy identification and characterization of splice junctions, and holds a versatile computational pipeline for in-depth annotation and classification of alternative splicing with different patterns. Moreover, it provides a user-friendly environment in which an alternative splicing landscape can be displayed in a straightforward and flexible manner. In conclusion, SplicingViewer can be widely used for studying alternative splicing easily and efficiently. SplicingViewer can be freely accessed at http://bioinformatics.zj.cn/splicingviewer.

Detection of splicing events and multiread locations from RNA-seq data based on a geometric-tail (GT) distribution of intron length

Background

RNA sequencing (RNA-seq) measures gene expression levels and permits splicing analysis. Many existing aligners are capable of mapping millions of sequencing reads onto a reference genome. For reads that can be mapped to multiple positions along the reference genome (multireads), these aligners may either randomly assign them to a location, or discard them altogether. Either way could bias downstream analyses. Meanwhile, challenges remain in the alignment of reads spanning across splice junctions. Existing splicing-aware aligners that rely on the read-count method in identifying junction sites are inevitably affected by sequencing depths.

Results

The distance between aligned positions of paired-end (PE) reads or two parts of a spliced read is dependent on the experiment protocol and gene structures. We here proposed a new method that employs an empirical geometric-tail (GT) distribution of intron lengths to make a rational choice in multireads selection and splice-sites detection, according to the aligned distances from PE and sliced reads.

Conclusions

GT models that combine sequence similarity from alignment, and together with the probability of length distribution, could accurately determine the location of both multireads and spliced reads.

UG repeats/TDP-43 interactions near 5' splice sites exert unpredictable effects on splicing modulation

TDP-43 is a nuclear protein implicated in the pathogenesis of several neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal lobar degeneration, with broad involvement in numerous stages of RNA processing ranging from transcription to translation. In diseased neurons, TDP-43 mostly aggregates in the cytoplasm, suggesting that a loss of protein function in the nucleus may play an important role in neurodegeneration. A better understanding of TDP-43 general nuclear functions is therefore an essential step to evaluate this possibility. Presently, the TDP-43 best-characterized functional property is its ability to modulate pre-mRNA splicing when binding in proximity of 3'SS acceptor sequences. In this work, using a variety of artificial and natural splicing substrates, we have investigated the effects of TDP-43 binding to UG repeats in the vicinity of 5'SS donor sequences. In general, our results show that UG repeats are not powerful splicing regulatory elements when located near to exonic 5'SS sequences. However, in cases like the BRCA1, ETF1, and RXRG genes, TDP-43 binding to natural UG-repeated sequences can act as either an activator or a suppressor of 5'SS recognition, depending on splice site strength and on the presence of additional splicing regulatory sequences. The results of this analysis suggest that a role of UG repeats/TDP-43 in 5'SS recognition may exists and may become critical in the presence of mutations that weaken the 5'SS. The general rule that can be drawn at the moment is that the importance of UG repeats near 5' splice sites should always be experimentally validated on a case-by-case basis.

Gene expression analysis of whole blood, peripheral blood mononuclear cells, and lymphoblastoid cell lines from the Framingham Heart Study

Despite a growing number of reports of gene expression analysis from blood-derived RNA sources, there have been few systematic comparisons of various RNA sources in transcriptomic analysis or for biomarker discovery in the context of cardiovascular disease (CVD). As a pilot study of the Systems Approach to Biomarker Research (SABRe) in CVD Initiative, this investigation used Affymetrix Exon arrays to characterize gene expression of three blood-derived RNA sources: lymphoblastoid cell lines (LCL), whole blood using PAXgene tubes (PAX), and peripheral blood mononuclear cells (PBMC). Their performance was compared in relation to identifying transcript associations with sex and CVD risk factors, such as age, high-density lipoprotein, and smoking status, and the differential blood cell count. We also identified a set of exons that vary substantially between participants, but consistently in each RNA source. Such exons are thus stable phenotypes of the participant and may potentially become useful fingerprinting biomarkers. In agreement with previous studies, we found that each of the RNA sources is distinct. Unlike PAX and PBMC, LCL gene expression showed little association with the differential blood count. LCL, however, was able to detect two genes related to smoking status. PAX and PBMC identified Y-chromosome probe sets similarly and slightly better than LCL.

TILLING: a shortcut in functional genomics

Recent advances in large-scale genome sequencing projects have opened up new possibilities for the application of conventional mutation techniques in not only forward but also reverse genetics strategies. TILLING (Targeting Induced Local Lesions IN Genomes) was developed a decade ago as an alternative to insertional mutagenesis. It takes advantage of classical mutagenesis, sequence availability and high-throughput screening for nucleotide polymorphisms in a targeted sequence. The main advantage of TILLING as a reverse genetics strategy is that it can be applied to any species, regardless of its genome size and ploidy level. The TILLING protocol provides a high frequency of point mutations distributed randomly in the genome. The great mutagenic potential of chemical agents to generate a high rate of nucleotide substitutions has been proven by the high density of mutations reported for TILLING populations in various plant species. For most of them, the analysis of several genes revealed 1 mutation/200-500 kb screened and much higher densities were observed for polyploid species, such as wheat. High-throughput TILLING permits the rapid and low-cost discovery of new alleles that are induced in plants. Several research centres have established a TILLING public service for various plant species. The recent trends in TILLING procedures rely on the diversification of bioinformatic tools, new methods of mutation detection, including mismatch-specific and sensitive endonucleases, but also various alternatives for LI-COR screening and single nucleotide polymorphism (SNP) discovery using next-generation sequencing technologies. The TILLING strategy has found numerous applications in functional genomics. Additionally, wide applications of this throughput method in basic and applied research have already been implemented through modifications of the original TILLING strategy, such as Ecotilling or Deletion TILLING.

The polypyrimidine tract binding protein regulates desaturase alternative splicing and PUFA composition

The Δ6 desaturase, encoded by FADS2, plays a crucial role in omega-3 and omega-6 fatty acid synthesis. These fatty acids are essential components of the central nervous system, and they act as precursors for eicosanoid signaling molecules and as direct modulators of gene expression. The polypyrimidine tract binding protein (PTB or hnRNP I) is a splicing factor that regulates alternative pre-mRNA splicing. Here, PTB is shown to bind an exonic splicing silencer element and repress alternative splicing of FADS2 into FADS2 AT1. PTB and FADS2AT1 were inversely correlated in neonatal baboon tissues, implicating PTB as a major regulator of tissue-specific FADS2 splicing. In HepG2 cells, PTB knockdown modulated alternative splicing of FADS2, as well as FADS3, a putative desaturase of unknown function. Omega-3 fatty acids decreased by nearly one half relative to omega-6 fatty acids in PTB knockdown cells compared with controls, with a particularly strong decrease in eicosapentaenoic acid (EPA) concentration and its ratio to arachidonic acid (ARA). This is a rare demonstration of a mechanism specifically altering the cellular omega-3 to omega-6 fatty acid ratio without any change in diet/media. These findings reveal a novel role for PTB, regulating availability of membrane components and eicosanoid precursors for cell signaling.

A conserved tissue-specific homeodomain-less isoform of MEIS1 is downregulated in colorectal cancer

Colorectal cancer is one of the most common cancers in developed nations and is the result of both environmental and genetic factors. Many of the genetic lesions observed in colorectal cancer alter expression of homeobox genes, which encode homeodomain transcription factors. The MEIS1 homeobox gene is known to be involved in several hematological malignancies and solid tumors and recent evidence suggests that expression of the MEIS1 transcript is altered in colorectal cancer. Despite this potential connection, little is known about the role of the gene in the intestines. We probed murine gastrointestinal tissue samples with an N-terminal Meis1 antibody, revealing expression of two previously described isoforms, as well as two novel Meis1 products. A 32 kD Meis1 product was expressed in the nuclei of non-epithelial cells in the stomach and colon, while a 27 kD product was expressed in the cytoplasm of epithelial cells in the proximal colon. Our data suggest that the 27 kD and 32 kD Meis1 proteins are both forms of the Meis1d protein, a homeodomain-less isoform whose transcript was previously identified in cDNA screens. Both the MEIS1D transcript and protein were expressed in human colon mucosa. Expression of the MEIS1D protein was downregulated in 83% (10/12) of primary colorectal cancer samples compared to matched normal mucosa, indicating that MEIS1D is a biomarker of colorectal tumorigenesis. The decreased expression of MEIS1D in colon tumors also suggests that this conserved homeodomain-less isoform may act as a tumor suppressor in human colorectal cancer.