Alternative splicing: a potential source of functional innovation in the eukaryotic genome

Systematic approach of polyploidy as an evolutionary genetic and genomic phenomenon in horticultural crops

Polyploidy is thought to be an evolutionary and systematic mechanism for gene flow and phenotypic advancement in flowering plants. It is a natural phenomenon that promotes diversity by creating new permutations enhancing the prime potentials as compared to progenitors. Two different pathways have been recognized in studying polyploidy in nature; mitotic or somatic chromosome doubling and cytogenetics variation. Secondly, the vital influence of being polyploid is its heritable property (unreduced reproductive cells) formed during first and second-division restitution (FDR & SDR). Different approaches either chemical (Colchicine, Oryzalin, Caffeine, Trifuralin, or phosphoric amides) or gaseous i.e. Nitrous oxide have been deliberated as strong polyploidy causing agents. A wide range of cytogenetic practices like chromosomes study, ploidy, genome analysis, and plant morphology and anatomy have been studied in different plant species. Flow cytometry for ploidy and chromosome analysis through fluorescence and genomic in situ hybridization (FISH & GISH) are the basic methods to evaluate heredity substances sampled from leaves and roots. Many horticultural crops have been developed successfully and released commercially for consumption. Moreover, some deep detailed studies are needed to check the strong relationship between unique morphological features and genetic makeup concerning genes and hormonal expression in a strong approach.

Alternative splicing of lncRNA LAIR fine-tunes the regulation of neighboring yield-related gene LRK1 expression in rice

The diversity in alternative splicing of long noncoding RNAs (lncRNAs) poses a challenge for functional annotation of lncRNAs. Moreover, little is known on the effects of alternatively spliced lncRNAs on crop yield. In this study, we cloned nine isoforms resulting from the alternative splicing of the lncRNA LAIR in rice. The LAIR isoforms are generated via alternative 5'/3' splice sites and different combinations of specific introns. All LAIR isoforms activate the expression of the neighboring LRK1 gene and enhance yield-related rice traits. In addition, there are slight differences in the binding ability of LAIR isoforms to the epigenetic modification-related proteins OsMOF and OsWDR5, which affect the enrichment of H4K16ac and H3K4me3 at the LRK1 locus, and consequently fine-tune the regulation of LRK1 expression and yield-related traits. These differences in binding may be caused by polymorphic changes to the RNA secondary structure resulting from alternative splicing. It was also observed that the composition of LAIR isoforms was sensitive to abiotic stress. These findings suggest that the alternative splicing of LAIR leads to the formation of a functional transcript population that precisely regulates yield-related gene expression, which may be relevant for phenotypic polymorphism-based crop breeding under changing environmental conditions.

Identification and functional regulation of three alternative splicing isoforms of the fthl27 gene in miiuy croaker, Miichthys miiuy

Alternative splicing is an important basic mechanism for eukaryotes to control gene expression. Different forms of alternative splicing may lead to the production of protein subtypes with different functions, leading to the expansion of protein diversity in organisms, affecting cell production and metabolism, and is even related to the occurrence of many diseases. Many studies have shown that ferritin is usually associated with inflammation, vascular proliferation, and tumors, which is the focus of immunological research. It not only plays a role in iron metabolism and storage in the body, but also plays an important regulatory role in pathways related to immune and inflammatory regulation. However, there are few studies on alternative splicing events of the ferritin gene nowadays. Therefore, this study identified three different splicing isoforms in its ferritin gene fthl27 of Miichthys miiuy through Sanger sequencing, qRT-PCR, and other experimental techniques, and we found that three different splicing isoforms of the ferritin gene fthl27 in M. Miiuy cells showed an upregulation trend after being stimulated by Lipopolysaccharide (LPS) and poly (I: C). The experiment also found that the three isoforms may have different regulatory effects on the expression of inflammatory factors and antiviral immune factors, playing an important role in the innate immune response of fish.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

Whole-genome sequencing reveals the genetic mechanisms of domestication in classical inbred mice

Background

The laboratory mouse was domesticated from the wild house mouse. Understanding the genetics underlying domestication in laboratory mice, especially in the widely used classical inbred mice, is vital for studies using mouse models. However, the genetic mechanism of laboratory mouse domestication remains unknown due to lack of adequate genomic sequences of wild mice.

Results

We analyze the genetic relationships by whole-genome resequencing of 36 wild mice and 36 inbred strains. All classical inbred mice cluster together distinctly from wild and wild-derived inbred mice. Using nucleotide diversity analysis, Fst, and XP-CLR, we identify 339 positively selected genes that are closely associated with nervous system function. Approximately one third of these positively selected genes are highly expressed in brain tissues, and genetic mouse models of 125 genes in the positively selected genes exhibit abnormal behavioral or nervous system phenotypes. These positively selected genes show a higher ratio of differential expression between wild and classical inbred mice compared with all genes, especially in the hippocampus and frontal lobe. Using a mutant mouse model, we find that the SNP rs27900929 (T>C) in gene Astn2 significantly reduces the tameness of mice and modifies the ratio of the two Astn2 (a/b) isoforms.

Conclusion

Our study indicates that classical inbred mice experienced high selection pressure during domestication under laboratory conditions. The analysis shows the positively selected genes are closely associated with behavior and the nervous system in mice. Tameness may be related to the Astn2 mutation and regulated by the ratio of the two Astn2 (a/b) isoforms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02772-1.

A survey of transcriptome complexity using full-length isoform sequencing in the tea plant Camellia sinensis

Tea is one of the most popular beverages and its leaves are rich in catechins, contributing to the diverse flavor as well as beneficial for human health. However, the study of the post-transcriptional regulatory mechanism affecting the synthesis of catechins remains insufficient. Here, we sequenced the transcriptome using PacBio sequencing technology and obtained 63,111 full-length high-quality isoforms, including 1302 potential novel genes and 583 highly reliable fusion transcripts. We also identified 1204 lncRNAs with high quality, containing 188 known and 1016 novel lncRNAs. In addition, 311 mis-annotated genes were corrected based on the high-quality Isoseq reads. A large number of alternative splicing (AS) events (3784) and alternative polyadenylation (APA) genes (18,714) were analyzed, accounting for 8.84% and 43.7% of the total annotated genes, respectively. We also found that 2884 genes containing AS and APA features exhibited higher expression levels than other genes. These genes are mainly involved in amino acid biosynthesis, carbon fixation in photosynthetic organisms, phenylalanine, tyrosine, tryptophan biosynthesis, and pyruvate metabolism, suggesting that they play an essential role in the catechins content of tea polyphenols. Our results further improved the level of genome annotation and indicated that post-transcriptional regulation plays a crucial part in synthesizing catechins.

RNA sequencing role and application in clinical diagnostic

Although whole-exome sequencing and whole-genome sequencing has tremendously improved our understanding of the genetic etiology of human disorders, about half of the patients still do not receive a molecular diagnosis. The high fraction of variants with uncertain significance and the challenges of interpretation of noncoding variants have urged scientists to implement RNA sequencing (RNA-seq) in the diagnostic approach as a high throughput assay to complement genomic data with functional evidence. RNA-seq data can be used to identify aberrantly spliced genes, detect allele-specific expression, and identify gene expression outliers. Amongst eight studies utilizing RNA-seq, a mean diagnostic uplift of 15% has been reported. Here, we provide an overview of how RNA-seq has been implemented to aid in identifying the causal variants of Mendelian disorders.

Principles and correction of 5'-splice site selection

In Eukarya, immature mRNA transcripts (pre-mRNA) often contain coding sequences, or exons, interleaved by non-coding sequences, or introns. Introns are removed upon splicing, and further regulation of the retained exons leads to alternatively spliced mRNA. The splicing reaction requires the stepwise assembly of the spliceosome, a macromolecular machine composed of small nuclear ribonucleoproteins (snRNPs). This review focuses on the early stage of spliceosome assembly, when U1 snRNP defines each intron 5'-splice site (5'ss) in the pre-mRNA. We first introduce the splicing reaction and the impact of alternative splicing on gene expression regulation. Thereafter, we extensively discuss splicing descriptors that influence the 5'ss selection by U1 snRNP, such as sequence determinants, and interactions mediated by U1-specific proteins or U1 small nuclear RNA (U1 snRNA). We also include examples of diseases that affect the 5'ss selection by U1 snRNP, and discuss recent therapeutic advances that manipulate U1 snRNP 5'ss selectivity with antisense oligonucleotides and small-molecule splicing switches.

Genome-wide Analysis of Alternative Gene Splicing Associated with Virulence in the Brown Planthopper Nilaparvata lugens (Hemiptera: Delphacidae)

Alternative splicing of protein coding genes plays a profound role in phenotypic variation for many eukaryotic organisms. The development of high-throughput sequencing and bioinformatics algorithms provides the possibility of genome-wide identification of alternative splicing events in eukaryotes. However, for the brown planthopper Nilaparvata lugens, a destructive pest of rice crops, whole-genome distribution of alternative splicing events and the role of alternative splicing in the phenotypic plasticity of virulence have not previously been estimated. Here, we developed an analysis pipeline to identify alternative splicing events in the genome of N. lugens. Differential expression analysis and functional annotation were performed on datasets related to different virulence phenotypes. In total, 27,880 alternative splicing events corresponding to 9,787 multi-exon genes were detected in N. lugens. Among them, specifically expressed alternative splicing transcripts in the virulent Mudgo population were enriched in metabolic process categories, while transcripts in the avirulent TN1 population were enriched in regulator activity categories. In addition, genes encoding odorant receptor, secreted saliva protein and xenobiotic metabolic P450 monooxygenase showed different splicing patterns between Mudgo population and TN1 population. Host change experiment also revealed that an isoform of a P450 gene could be specially induced by the stimulation of resistant rice variety Mudgo. This research pioneered a genome-wide study of alternative gene splicing in the rice brown planthopper. Differences in alternative splicing between virulent and avirulent populations indicated that alternative splicing might play an important role in the formation of virulence phenotypes in N. lugens.

The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis

Safflower (Carthamus tinctorius L.), a member of the Asteraceae, is a popular crop due to its high linoleic acid (LA) and flavonoid (such as hydroxysafflor yellow A) contents. Here, we report the first high-quality genome assembly (contig N50 of 21.23 Mb) for the 12 pseudochromosomes of safflower using single-molecule real-time sequencing, Hi-C mapping technologies and a genetic linkage map. Phyloge nomic analysis showed that safflower diverged from artichoke (Cynara cardunculus) and sunflower (Helianthus annuus) approximately 30.7 and 60.5 million years ago, respectively. Comparative genomic analyses revealed that uniquely expanded gene families in safflower were enriched for those predicted to be involved in lipid metabolism and transport and abscisic acid signalling. Notably, the fatty acid desaturase 2 (FAD2) and chalcone synthase (CHS) families, which function in the LA and flavonoid biosynthesis pathways, respectively, were expanded via tandem duplications in safflower. CarFAD2-12 was specifically expressed in seeds and was vital for high-LA content in seeds, while tandemly duplicated CarFAD2 genes were up-regulated in ovaries compared to CarFAD2-12, which indicates regulatory divergence of FAD2 in seeds and ovaries. CarCHS1, CarCHS4 and tandem-duplicated CarCHS5˜CarCHS6, which were up-regulated compared to other CarCHS members at early stages, contribute to the accumulation of major flavonoids in flowers. In addition, our data reveal multiple alternative splicing events in gene families related to fatty acid and flavonoid biosynthesis. Together, these results provide a high-quality reference genome and evolutionary insights into the molecular basis of fatty acid and flavonoid biosynthesis in safflower.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families

Sclerotinia sclerotiorum is a notorious generalist plant pathogen that threatens more than 600 host plants, including wild and cultivated species. The molecular bases underlying the broad compatibility of S. sclerotiorum with its hosts is not fully elucidated. In contrast to higher plants and animals, alternative splicing (AS) is not well studied in plant-pathogenic fungi. AS is a common regulated cellular process that increases cell protein and RNA diversity. In this study, we annotated spliceosome genes in the genome of S. sclerotiorum and characterized their expression in vitro and during the colonization of six host species. Several spliceosome genes were differentially expressed in planta, suggesting that AS was altered during infection. Using stringent parameters, we identified 1,487 S. sclerotiorum genes differentially expressed in planta and exhibiting alternative transcripts. The most common AS events during the colonization of all plants were retained introns and the alternative 3' receiver site. We identified S. sclerotiorum genes expressed in planta for which (a) the relative accumulation of alternative transcripts varies according to the host being colonized and (b) alternative transcripts harbour distinct protein domains. This notably included 42 genes encoding predicted secreted proteins showing high-confidence AS events. This study indicates that AS events are taking place in the plant pathogenic fungus S. sclerotiorum during the colonization of host plants and could generate functional diversity in the repertoire of proteins secreted by S. sclerotiorum during infection.

Comparative Transcriptome Profiling of Skeletal Muscle from Black Muscovy Duck at Different Growth Stages Using RNA-seq

In China, the production for duck meat is second only to that of chicken, and the demand for duck meat is also increasing. However, there is still unclear on the internal mechanism of regulating skeletal muscle growth and development in duck. This study aimed to identity candidate genes related to growth of duck skeletal muscle and explore the potential regulatory mechanism. RNA-seq technology was used to compare the transcriptome of skeletal muscles in black Muscovy ducks at different developmental stages (day 17, 21, 27, 31, and 34 of embryos and postnatal 6-month-olds). The SNPs and InDels of black Muscovy ducks at different growth stages were mainly in “INTRON”, “SYNONYMOUS_CODING”, “UTR_3_PRIME”, and “DOWNSTREAM”. The average number of AS in each sample was 37,267, mainly concentrated in TSS and TTS. Besides, a total of 19 to 5377 DEGs were detected in each pairwise comparison. Functional analysis showed that the DEGs were mainly involved in the processes of cell growth, muscle development, and cellular activities (junction, migration, assembly, differentiation, and proliferation). Many of DEGs were well known to be related to growth of skeletal muscle in black Muscovy duck, such as MyoG, FBXO1, MEF2A, and FoxN2. KEGG pathway analysis identified that the DEGs were significantly enriched in the pathways related to the focal adhesion, MAPK signaling pathway and regulation of the actin cytoskeleton. Some DEGs assigned to these pathways were potential candidate genes inducing the difference in muscle growth among the developmental stages, such as FAF1, RGS8, GRB10, SMYD3, and TNNI2. Our study identified several genes and pathways that may participate in the regulation of skeletal muscle growth in black Muscovy duck. These results should serve as an important resource revealing the molecular basis of muscle growth and development in duck.

The role of alternative splicing in cancer: From oncogenesis to drug resistance

Alternative splicing is a tightly regulated process whereby non-coding sequences of pre-mRNA are removed and protein-coding segments are assembled in diverse combinations, ultimately giving rise to proteins with distinct or even opposing functions. In the past decade, whole genome/transcriptome sequencing studies revealed the high complexity of splicing regulation, which occurs co-transcriptionally and is influenced by chromatin status and mRNA modifications. Consequently, splicing profiles of both healthy and malignant cells display high diversity and alternative splicing was shown to be widely deregulated in multiple cancer types. In particular, mutations in pre-mRNA regulatory sequences, splicing regulators and chromatin modifiers, as well as differential expression of splicing factors are important contributors to cancer pathogenesis. It has become clear that these aberrations contribute to many facets of cancer, including oncogenic transformation, cancer progression, response to anticancer drug treatment as well as resistance to therapy. In this respect, alternative splicing was shown to perturb the expression a broad spectrum of relevant genes involved in drug uptake/metabolism (i.e. SLC29A1, dCK, FPGS, and TP), activation of nuclear receptor pathways (i.e. GR, AR), regulation of apoptosis (i.e. MCL1, BCL-X, and FAS) and modulation of response to immunotherapy (CD19). Furthermore, aberrant splicing constitutes an important source of novel cancer biomarkers and the spliceosome machinery represents an attractive target for a novel and rapidly expanding class of therapeutic agents. Small molecule inhibitors targeting SF3B1 or splice factor kinases were highly cytotoxic against a wide range of cancer models, including drug-resistant cells. Importantly, these effects are enhanced in specific cancer subsets, such as splicing factor-mutated and c-MYC-driven tumors. Furthermore, pre-clinical studies report synergistic effects of spliceosome modulators in combination with conventional antitumor agents. These strategies based on the use of low dose splicing modulators could shift the therapeutic window towards decreased toxicity in healthy tissues. Here we provide an extensive overview of the latest findings in the field of regulation of splicing in cancer, including molecular mechanisms by which cancer cells harness alternative splicing to drive oncogenesis and evade anticancer drug treatment as well as splicing-based vulnerabilities that can provide novel treatment opportunities. Furthermore, we discuss current challenges arising from genome-wide detection and prediction methods of aberrant splicing, as well as unravelling functional relevance of the plethora of cancer-related splicing alterations.

Properties of two spliceoforms of major ampullate spidroin 1 reveal unique functions of N-linker region

Spiders produce diverse silk fibers with distinct properties for daily survival. Among these silk fibers, dragline silk spun by major ampullate gland is used for bridgelines and web radii, exhibiting both outstanding tensile strength and extensibility. Although more and more full-length major ampullate spidroin gene sequences have been reported, the research regarding alternative splicing events of spidroins are rare. Here we describe two spliceoforms of major ampullate spidroin 1 (MaSp1) from Araneus ventricosus, and both of them are lack of central repetitive region. The minor isoform only has terminal regions. For the major isoform, however, the N-linker and terminal regions are all retained. Furthermore, we investigated the functions of N-linker structure of A. ventricosus MaSp1, based on the properties of the two spliceoforms. The dimer level of major isoform (MaSp1-2) is higher than that of the minor isoform (MaSp1-1). Moreover, the MaSp1-2 protein display higher melting temperature (T_m) than MaSp1-1, and the MaSp1-2 fibers exhibit higher tensile strength than MaSp1-1 fibers. These studies demonstrate that the N-linker region promotes the formation of intermolecular disulphide bond, suggesting a strategy to enhance the thermostabilization and mechanical properties of spidroins.

Trans-Disciplinary Responses to Climate Change: Lessons from Rice-Based Systems in Asia

Climate change will continue to have a largely detrimental impact on the agricultural sector worldwide because of predicted rising temperatures, variable rainfall, and an increase in extreme weather events. Reduced crop yields will lead to higher food prices and increased hardship for low income populations, especially in urban areas. Action on climate change is one of the Sustainable Development Goals (SDG 13) and is linked to the Paris Climate Agreement. The research challenge posed by climate change is so complex that a trans-disciplinary response is required, one that brings together researchers, practitioners, and policy-makers in networks where the lines between “research” and “development” become deliberately blurred. Fostering such networks will require researchers, throughout the world, not only to work across disciplines but also to pursue new South–North and South–South partnerships incorporating policy-makers and practitioners. We use our diverse research experiences to describe the emergence of such networks, such as the Direct Seeded Rice Consortium (DSRC) in South and Southeast Asia, and to identify lessons on how to facilitate and strengthen the development of trans-disciplinary responses to climate change.

Full-length transcriptome sequencing from multiple tissues of duck, Anas platyrhynchos

Duck (Anas platyrhynchos), one of the most economically important waterfowl, is an ideal model for studying the immune protection mechanism of birds. An incomplete duck reference genome and very limited availability of full-length cDNAs has hindered the identification of alternatively spliced transcripts and slowed down many basic studies in ducks. We applied PacBio Iso-Seq technologies to multiple tissues from duck for use in transcriptome sequencing. We obtained 199,993 full-length transcripts and comprehensively annotated these transcripts. 23,755 lncRNAs were predicted from all identified transcripts and 35,031 alternative splicing events, which divided into 5 models, were accurately predicted from 3,346 genes. Our data constitute a large increase in the known number of both lncRNA, and alternatively spliced transcripts of duck and plays an important role in improving current genome annotation. In addition, the data will be extremely useful for functional studies in other birds.

RNA-Seq Perspectives to Improve Clinical Diagnosis

In recent years, high-throughput next-generation sequencing technology has allowed a rapid increase in diagnostic capacity and precision through different bioinformatics processing algorithms, tools, and pipelines. The identification, annotation, and classification of sequence variants within different target regions are now considered a gold standard in clinical genetic diagnosis. However, this procedure lacks the ability to link regulatory events such as differential splicing to diseases. RNA-seq is necessary in clinical routine in order to interpret and detect among others splicing events and splicing variants, as it would increase the diagnostic rate by up to 10-35%. The transcriptome has a very dynamic nature, varying according to tissue type, cellular conditions, and environmental factors that may affect regulatory events such as splicing and the expression of genes or their isoforms. RNA-seq offers a robust technical analysis of this complexity, but it requires a profound knowledge of computational/statistical tools that may need to be adjusted depending on the disease under study. In this article we will cover RNA-seq analyses best practices applied to clinical routine, bioinformatics procedures, and present challenges of this approach.

Single-molecule long-read sequencing facilitates shrimp transcriptome research

Although shrimp are of great economic importance, few full-length shrimp transcriptomes are available. Here, we used Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology to generate transcripts from the Pacific white shrimp (Litopenaeus vannamei). We obtained 322,600 full-length non-chimeric reads, from which we generated 51,367 high-quality unique full-length transcripts. We corrected errors in the SMRT sequences by comparison with Illumina-produced short reads. We successfully annotated 81.72% of all unique SMRT transcripts against the NCBI non-redundant database, 58.63% against Swiss-Prot, 45.38% against Gene Ontology, 32.57% against Clusters of Orthologous Groups of proteins (COG), and 47.83% against Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Across all transcripts, we identified 3,958 long non-coding RNAs (lncRNAs) and 80,650 simple sequence repeats (SSRs). Our study provides a rich set of full-length cDNA sequences for L. vannamei, which will greatly facilitate shrimp transcriptome research.

Transcription-Associated Mutation Promotes RNA Complexity in Highly Expressed Genes-A Major New Source of Selectable Variation

Alternatively spliced transcript isoforms are thought to play a critical role for functional diversity. However, the mechanism generating the enormous diversity of spliced transcript isoforms remains unknown, and its biological significance remains unclear. We analyzed transcriptomes in saker falcons, chickens, and mice to show that alternative splicing occurs more frequently, yielding more isoforms, in highly expressed genes. We focused on hemoglobin in the falcon, the most abundantly expressed genes in blood, finding that alternative splicing produces 10-fold more isoforms than expected from the number of splice junctions in the genome. These isoforms were produced mainly by alternative use of de novo splice sites generated by transcription-associated mutation (TAM), not by the RNA editing mechanism normally invoked. We found that high expression of globin genes increases mutation frequencies during transcription, especially on nontranscribed DNA strands. After DNA replication, transcribed strands inherit these somatic mutations, creating de novo splice sites, and generating multiple distinct isoforms in the cell clone. Bisulfate sequencing revealed that DNA methylation may counteract this process by suppressing TAM, suggesting DNA methylation can spatially regulate RNA complexity. RNA profiling showed that falcons living on the high Qinghai-Tibetan Plateau possess greater global gene expression levels and higher diversity of mean to high abundance isoforms (reads per kilobases per million mapped reads ≥18) than their low-altitude counterparts, and we speculate that this may enhance their oxygen transport capacity under low-oxygen environments. Thus, TAM-induced RNA diversity may be physiologically significant, providing an alternative strategy in lifestyle evolution.

Insertion of 275-bp SINE into first intron of PDIA4 gene is associated with litter size in Xiang pigs

The aim of the study was to investigate the SINE polymorphism in Xiang, Kele, Qianbei black, Jiangkouluobo, Large White, and Duroc pig breeds. The PCR based detection method was conducted to identify the short interspersed nuclear element (SINE) polymorphism in the PDIA4 gene. There were greater frequencies of the SINE-/- genotypes in Xiang pigs (55.9%) as compared with other pig breed groups. There was an association between this 275 bp SINE polymorphism and litter size (P = 0.003). The homozygous SINE+/+ genotype of the PDIA4 gene had a 1.45-piglets larger litter sizes compared to those with the homozygous SINE-/- genotype. Furthermore, there were assessments of mRNA and protein abundances as a result of PDIA4 gene expression in the ovaries of Xiang pigs for the three different SINE genotypes, and the results indicated that relative abundances of PDIA4 mRNA and protein was greater for the SINE-/- and SINE-/+ genotypes compared with the SINE+/+ genotype (P < 0.05). These findings suggested that the 275 bp SINE polymorphism might change the expression of the PDIA4 gene and could be used as a candidate DNA marker for the selection for litter size in Xiang pigs.

Genome-wide transcriptome analysis identifies alternative splicing regulatory network and key splicing factors in mouse and human psoriasis

Psoriasis is a chronic inflammatory disease that affects the skin, nails, and joints. For understanding the mechanism of psoriasis, though, alternative splicing analysis has received relatively little attention in the field. Here, we developed and applied several computational analysis methods to study psoriasis. Using psoriasis mouse and human datasets, our differential alternative splicing analyses detected hundreds of differential alternative splicing changes. Our analysis of conservation revealed many exon-skipping events conserved between mice and humans. In addition, our splicing signature comparison analysis using the psoriasis datasets and our curated splicing factor perturbation RNA-Seq database, SFMetaDB, identified nine candidate splicing factors that may be important in regulating splicing in the psoriasis mouse model dataset. Three of the nine splicing factors were confirmed upon analyzing the human data. Our computational methods have generated predictions for the potential role of splicing in psoriasis. Future experiments on the novel candidates predicted by our computational analysis are expected to provide a better understanding of the molecular mechanism of psoriasis and to pave the way for new therapeutic treatments.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Evidence for a Strong Correlation Between Transcription Factor Protein Disorder and Organismic Complexity

Studies of diverse phylogenetic lineages reveal that protein disorder increases in concert with organismic complexity but that differences nevertheless exist among lineages. To gain insight into this phenomenology, we analyzed all of the transcription factor (TF) families for which sequences are known for 17 species spanning bacteria, yeast, algae, land plants, and animals and for which the number of different cell types has been reported in the primary literature. Although the fraction of disordered residues in TF sequences is often moderately or poorly correlated with organismic complexity as gauged by cell-type number (r² < 0.5), an unbiased and phylogenetically broad analysis shows that organismic complexity is positively and strongly correlated with the total number of TFs, the number of their spliced variants and their total disordered residues content (r² > 0.8). Furthermore, the correlation between the fraction of disordered residues and cell-type number becomes stronger when confined to the TF families participating in cell cycle, cell size, cell division, cell differentiation, or cell proliferation, and other important developmental processes. The data also indicate that evolutionarily simpler organisms allow for the detection of subtle differences in the conserved IDRs of TFs as well as changes in variable IDRs, which can influence the DNA recognition and multifunctionality of TFs through direct or indirect mechanisms. Although strong correlations cannot be taken as evidence for cause-and-effect relationships, we interpret our data to indicate that increasing TF disorder likely was an important factor contributing to the evolution of organismic complexity and not merely a concurrent unrelated effect of increasing organismic complexity.

Identification of genome-wide non-canonical spliced regions and analysis of biological functions for spliced sequences using Read-Split-Fly

Background

It is generally thought that most canonical or non-canonical splicing events involving U2- and U12 spliceosomes occur within nuclear pre-mRNAs. However, the question of whether at least some U12-type splicing occurs in the cytoplasm is still unclear. In recent years next-generation sequencing technologies have revolutionized the field. The “Read-Split-Walk” (RSW) and “Read-Split-Run” (RSR) methods were developed to identify genome-wide non-canonical spliced regions including special events occurring in cytoplasm. As the significant amount of genome/transcriptome data such as, Encyclopedia of DNA Elements (ENCODE) project, have been generated, we have advanced a newer more memory-efficient version of the algorithm, “Read-Split-Fly” (RSF), which can detect non-canonical spliced regions with higher sensitivity and improved speed. The RSF algorithm also outputs the spliced sequences for further downstream biological function analysis.

Results

We used open access ENCODE project RNA-Seq data to search spliced intron sequences against the U12-type spliced intron sequence database to examine whether some events could occur as potential signatures of U12-type splicing. The check was performed by searching spliced sequences against 5’ss and 3’ss sequences from the well-known orthologous U12-type spliceosomal intron database U12DB. Preliminary results of searching 70 ENCODE samples indicated that the presence of 5’ss with U12-type signature is more frequent than U2-type and prevalent in non-canonical junctions reported by RSF. The selected spliced sequences have also been further studied using miRBase to elucidate their functionality. Preliminary results from 70 samples of ENCODE datasets show that several miRNAs are prevalent in studied ENCODE samples. Two of these are associated with many diseases as suggested in the literature. Specifically, hsa-miR-1273 and hsa-miR-548 are associated with many diseases and cancers.

Conclusions

Our RSF pipeline is able to detect many possible junctions (especially those with a high RPKM) with very high overall accuracy and relative high accuracy for novel junctions. We have incorporated useful parameter features into the pipeline such as, handling variable-length read data, and searching spliced sequences for splicing signatures and miRNA events. We suggest RSF, a tool for identifying novel splicing events, is applicable to study a range of diseases across biological systems under different experimental conditions.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1801-y) contains supplementary material, which is available to authorized users.

Comparison of Alternative Splicing Junction Detection Tools Using RNA-Seq Data

BACKGROUND

Alternative splicing (AS) is a posttranscriptional process that produces differ-ent transcripts from the same gene and is important to produce diverse protein products in response to environmental stimuli. AS occurs at specific sites on the mRNA sequence, some of which have been de-fined. Multiple bioinformatics tools have been developed to detect AS from experimental data.

OBJECTIVES

The goal of this review is to help researchers use specific tools to aid their research and to develop new AS detection tools based on these previously established tools.

METHOD

We selected 15 AS detection tools that were recently published; we classified and delineated them on several aspects. Also, a performance comparison of these tools with the same starting input was conducted.

RESULT

We reviewed the following categorized features of the tools: Publication information, working principles, generic and distinct workflows, running platform, input data requirement, sequencing depth dependency, reads mapped to multiple locations, isoform annotation basis, precise detected AS types, and performance benchmarks.

CONCLUSION

Through comparisons of these tools, we provide a panorama of the advantages and short-comings of each tool and their scopes of application.

A Bioinformatics-Based Alternative mRNA Splicing Code that May Explain Some Disease Mutations Is Conserved in Animals

Deep sequencing of cDNAs made from spliced mRNAs indicates that most coding genes in many animals and plants have pre-mRNA transcripts that are alternatively spliced. In pre-mRNAs, in addition to invariant exons that are present in almost all mature mRNA products, there are at least 6 additional types of exons, such as exons from alternative promoters or with alternative polyA sites, mutually exclusive exons, skipped exons, or exons with alternative 5′ or 3′ splice sites. Our bioinformatics-based hypothesis is that, in analogy to the genetic code, there is an “alternative-splicing code” in introns and flanking exon sequences, analogous to the genetic code, that directs alternative splicing of many of the 36 types of introns. In humans, we identified 42 different consensus sequences that are each present in at least 100 human introns. 37 of the 42 top consensus sequences are significantly enriched or depleted in at least one of the 36 types of introns. We further supported our hypothesis by showing that 96 out of 96 analyzed human disease mutations that affect RNA splicing, and change alternative splicing from one class to another, can be partially explained by a mutation altering a consensus sequence from one type of intron to that of another type of intron. Some of the alternative splicing consensus sequences, and presumably their small-RNA or protein targets, are evolutionarily conserved from 50 plant to animal species. We also noticed the set of introns within a gene usually share the same splicing codes, thus arguing that one sub-type of splicesosome might process all (or most) of the introns in a given gene. Our work sheds new light on a possible mechanism for generating the tremendous diversity in protein structure by alternative splicing of pre-mRNAs.

Perspectives on the history of evo-devo and the contemporary research landscape in the genomics era

A fundamental question in biology is how the extraordinary range of living organisms arose. In this theme issue, we celebrate how evolutionary studies on the origins of morphological diversity have changed over the past 350 years since the first publication of the Philosophical Transactions of The Royal Society Current understanding of this topic is enriched by many disciplines, including anatomy, palaeontology, developmental biology, genetics and genomics. Development is central because it is the means by which genetic information of an organism is translated into morphology. The discovery of the genetic basis of development has revealed how changes in form can be inherited, leading to the emergence of the field known as evolutionary developmental biology (evo-devo). Recent approaches include imaging, quantitative morphometrics and, in particular, genomics, which brings a new dimension. Articles in this issue illustrate the contemporary evo-devo field by considering general principles emerging from genomics and how this and other approaches are applied to specific questions about the evolution of major transitions and innovations in morphology, diversification and modification of structures, intraspecific morphological variation and developmental plasticity. Current approaches enable a much broader range of organisms to be studied, thus building a better appreciation of the origins of morphological diversity.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Genome-Wide Analysis of Alternative Splicing Provides Insights into Stress Adaptation of the Pacific Oyster

Alternative splicing (AS) is thought to enhance transcriptome diversity dramatically and play an important role in stress adaptation. While well studied in vertebrates, AS remains poorly understood in invertebrates. Here, we used high-throughput RNA-sequencing data to perform a genome-wide survey of AS in the Pacific oyster (Crassostrea gigas), an economically important mollusk that is cultivated worldwide. This analysis identified 8223 AS events corresponding to 4480 genes in the Pacific oyster, suggesting that about 16 % of oyster multiexonic genes undergo AS. We observed that a majority of the identified AS events were related to skipped exons (37.8 %). Then Gene Ontology analysis was conducted to analyze the function of the genes that undergo AS and the genes that produce more than five AS isoforms. After that, the expression of AS isoforms facing temperature, salinity, and air exposure challenge were examined. To validate our bioinformatic-predicted results and examine whether AS affects stress adaptation, we selected heat-shock protein 60 (HSP60) and HSP90 genes, both of which experience AS, for reverse transcription PCR (RT-PCR). We also performed quantitative real-time PCR (qRT-PCR) to determine the relative expression of each AS isoform among different stress adapted populations. Our study indicates that AS events are likely complex in the Pacific oyster and may be related to stress adaptation. These results will complement the predicted gene database of C. gigas and provide an invaluable resource for future functional genomic studies on molluscs.

Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing

Zea mays is an important genetic model for elucidating transcriptional networks. Uncertainties about the complete structure of mRNA transcripts limit the progress of research in this system. Here, using single-molecule sequencing technology, we produce 111,151 transcripts from 6 tissues capturing ∼70% of the genes annotated in maize RefGen_v3 genome. A large proportion of transcripts (57%) represent novel, sometimes tissue-specific, isoforms of known genes and 3% correspond to novel gene loci. In other cases, the identified transcripts have improved existing gene models. Averaging across all six tissues, 90% of the splice junctions are supported by short reads from matched tissues. In addition, we identified a large number of novel long non-coding RNAs and fusion transcripts and found that DNA methylation plays an important role in generating various isoforms. Our results show that characterization of the maize B73 transcriptome is far from complete, and that maize gene expression is more complex than previously thought.

Abnormalities in Alternative Splicing of Apoptotic Genes and Cardiovascular Diseases

Apoptosis is required for normal heart development in the embryo, but has also been shown to be an important factor in the occurrence of heart disease. Alternative splicing of apoptotic genes is currently emerging as a diagnostic and therapeutic target for heart disease. This review addresses the involvement of abnormalities in alternative splicing of apoptotic genes in cardiac disorders including cardiomyopathy, myocardial ischemia and heart failure. Many pro-apoptotic members of the Bcl-2 family have alternatively spliced isoforms that lack important active domains. These isoforms can play a negative regulatory role by binding to and inhibiting the pro-apoptotic forms. Alternative splicing is observed to be increased in various cardiovascular diseases with the level of alternate transcripts increasing elevated in diseased hearts compared to healthy subjects. In many cases these isoforms appear to be the underlying cause of the disease, while in others they may be induced in response to cardiovascular pathologies. Regardless of this, the detection of alternate splicing events in the heart can serve as useful diagnostic or prognostic tools, while those splicing events that seem to play a causative role in cardiovascular disease make attractive future drug targets.

Alternative Splice in Alternative Lice

Genomic and transcriptomics analyses have revealed human head and body lice to be almost genetically identical; although con-specific, they nevertheless occupy distinct ecological niches and have differing feeding patterns. Most importantly, while head lice are not known to be vector competent, body lice can transmit three serious bacterial diseases; epidemictyphus, trench fever, and relapsing fever. In order to gain insights into the molecular bases for these differences, we analyzed alternative splicing (AS) using next-generation sequencing data for one strain of head lice and one strain of body lice. We identified a total of 3,598 AS events which were head or body lice specific. Exon skipping AS events were overrepresented among both head and body lice, whereas intron retention events were underrepresented in both. However, both the enrichment of exon skipping and the underrepresentation of intron retention are significantly stronger in body lice compared with head lice. Genes containing body louse-specific AS events were found to be significantly enriched for functions associated with development of the nervous system, salivary gland, trachea, and ovarian follicle cells, as well as regulation of transcription. In contrast, no functional categories were overrepresented among genes with head louse-specific AS events. Together, our results constitute the first evidence for transcript pool differences in head and body lice, providing insights into molecular adaptations that enabled human lice to adapt to clothing, and representing a powerful illustration of the pivotal role AS can play in functional adaptation.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1914-5) contains supplementary material, which is available to authorized users.

Exploration of alternative splicing events in ten different grapevine cultivars

BACKGROUND

The complex dynamics of gene regulation in plants are still far from being fully understood. Among many factors involved, alternative splicing (AS) in particular is one of the least well documented. For many years, AS has been considered of less relevant in plants, especially when compared to animals, however, since the introduction of next generation sequencing techniques the number of plant genes believed to be alternatively spliced has increased exponentially.

RESULTS

Here, we performed a comprehensive high-throughput transcript sequencing of ten different grapevine cultivars, which resulted in the first high coverage atlas of the grape berry transcriptome. We also developed findAS, a software tool for the analysis of alternatively spliced junctions. We demonstrate that at least 44% of multi-exonic genes undergo AS and a large number of low abundance splice variants is present within the 131.622 splice junctions we have annotated from Pinot noir.

CONCLUSIONS

Our analysis shows that ~70% of AS events have relatively low expression levels, furthermore alternative splice sites seem to be enriched near the constitutive ones in some extent showing the noise of the splicing mechanisms. However, AS seems to be extensively conserved among the 10 cultivars.

Transcriptome analysis provides insights into the regulatory function of alternative splicing in antiviral immunity in grass carp (Ctenopharyngodon idella)

Characterization of the transcriptomic response to infection is an effective approach to understanding the immune mechanisms. Herein we challenged grass carp (Ctenopharyngodon idella) with grass carp reovirus (GCRV) and sequenced four cDNA libraries obtained from head-kidney and spleen by using Illumina Miseq. As a result, we gained a total of 21.52 Gb clean data with 107.96 million reads, and de novo assembled 55,199 unigenes with an average length of 1,470 bp. Comparative transcriptome analysis reveals that 217 unigenes are differentially expressed (fold-change of at least 4) between resistant and susceptible fish in both head-kidney and spleen, and of which 36 unigenes were validated by RT-qPCR experiment. The expression profile of immune-related genes demonstrates that the immune response of spleen is more intense than that of head-kidney. Remarkably, 11,811 unigenes contain multiple transcripts, of which 322 unigenes possess notably differentially expressed transcripts between the four transcriptomic datasets. Furthermore, the splicing transcripts of IL-12p40 and IL-1R1 are firstly found to play diverse roles in the antiviral response of fishes. This study provides a complete transcriptome dataset of C. idella, which is valuable for the studies of immune complexity and, moreover, throws light on the regulatory role of AS in antiviral immunity.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.

Lineage-specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage-specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage-specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage-specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.

Are all of the human exons alternatively spliced?

Alternative mRNA splicing (AS) is a major mechanism for increasing regulatory complexity. A key concept in AS is the distinction between alternatively and constitutively spliced exons (ASEs and CSEs, respectively). ASEs and CSEs have been reported to be differentially regulated, and to have distinct biological properties. However, the recent flood of RNA-sequencing data has obscured the boundary between ASEs and CSEs. Researchers are beginning to question whether ‘authentic CSEs’ do exist, and whether the ASE/CSE distinction is biologically invalid. Here, I examine the influences of increasing transcriptome data on the human ASE/CSE classification and our past understanding of the properties of these two types of exons. Interestingly, although the percentage of human ASEs has increased dramatically in recent years, the overall distinction between ASEs and CSEs remain valid. For example, CSEs are longer, evolve more slowly, and less frequently correspond to intrinsically disordered protein regions than ASEs. In addition, only a relatively small number of human genes have their transcripts composed entirely of ASEs despite the large amount of high-throughput transcriptome information. Therefore, the ‘backbone’ concept of AS, in which CSEs constitute the invariant part and ASEs the flexible part of the transcript, appears to be generally true despite the increasing percentage of ASEs in the human exome.

Potential impact of gene regulatory mechanisms on the evolution of multicellularity in the volvocine algae

A fundamental question in biology is how multicellular organisms can arise from their single-celled precursors. The evolution of multicellularity requires the adoption of new traits in unicellular ancestors that allows the generation of form by, for example, increasing the size and developing new cell types. But what are the genetic, cellular and biochemical bases underlying the evolution of multicellularity? Recent advances in evolutionary developmental biology suggest that the regulation of gene expression by cis-regulatory factors, gene duplication and alternative splicing contribute to phenotypic evolution. These mechanisms enable different degrees of phenotypic divergence and complexity with variation in traits from genomes with similar gene contents. In addition, signaling pathways specific to cell types are developed to guarantee the modulation of cellular and developmental processes matched to the cell types as well as the maintenance of multicellularity.

Perinatal methylmercury exposure perturbs the expression of Plp1 and Cnp splice variants in cerebellum of rat pups

Early life exposure to environmental chemicals can interfere with myelin formation in the developing brain, leading to neurological disorders. The Proteolipid Protein 1 (Plp1), Myelin Basic Protein (Mbp) and 2',3'-Cyclic Nucleotide 3'Phosphodiesterase (Cnp) genes expressed in oligodendrocytes and involved in myelination processes can be useful biomarkers of potential developmental neurotoxicity. In an earlier study, we concluded that the reduction in the expression levels of Mbp splice variants in juvenile rat cerebellum following perinatal methylmercury (MeHg) exposure were compatible with an overall reduction of mature oligodendrocytes population. This observation prompted us to analyze the expression of Plp1 and Cnp in developing rat cerebellum to further confirm and investigate the toxic effects of MeHg on vulnerable oligodendrocytes. Splice variants of Plp1 in human and of Cnp in mouse are curated in NCBI RefSeq database, but not for rat. Lack of annotation of splice variants can pose significant challenge for the reliable quantification of gene expression levels in toxicological studies. Therefore, we applied a "comparative sequence analysis" approach, relying on annotated splice variants in human/mouse and on evolutionary conservation of intron-exon structures, to identify additional splice variants of Plp1 and Cnp in rat. Then, we confirmed their identity by nucleotide sequencing and characterized their temporal expression patterns during brain development by RT-PCR. The measurement of total transcripts and individual splice variants of Plp1 and Cnp in the cerebellum of MeHg-exposed rat pups revealed a relatively similar level of reduction in their expression levels. This study further confirms that perinatal exposure to MeHg can impact oligodendrocytes in pups. Based on these observations, we conclude that monitoring the expression of these oligodendrocyte-enriched genes can be useful to identify toxic chemicals affecting myelination.

Alternative Splicing in the Obligate Biotrophic Oomycete Pathogen Pseudoperonospora cubensis

Pseudoperonospora cubensis is an obligate pathogen and causative agent of cucurbit downy mildew. To help advance our understanding of the pathogenicity of P. cubensis, we used RNA-Seq to improve the quality of its reference genome sequence. We also characterized the RNA-Seq dataset to inventory transcript isoforms and infer alternative splicing during different stages of its development. Almost half of the original gene annotations were improved and nearly 4,000 previously unannotated genes were identified. We also demonstrated that approximately 24% of the expressed genome and nearly 55% of the intron-containing genes from P. cubensis had evidence for alternative splicing. Our analyses revealed that intron retention is the predominant alternative splicing type in P. cubensis, with alternative 5'- and alternative 3'-splice sites occurring at lower frequencies. Representatives of the newly identified genes and predicted alternatively spliced transcripts were experimentally validated. The results presented herein highlight the utility of RNA-Seq for improving draft genome annotations and, through this approach, we demonstrate that alternative splicing occurs more frequently than previously predicted. In total, the current study provides evidence that alternative splicing plays a key role in transcriptome regulation and proteome diversification in plant-pathogenic oomycetes.

Functional Implications of RNA Splicing for Human Long Intergenic Noncoding RNAs

Long intergenic noncoding RNAs (lincRNAs) have been suggested as playing important roles in human gene regulation. The majority of annotated human lincRNAs include multiple exons and are alternatively spliced. However, the connections between alternative RNA splicing (AS) and the functions/regulations of lincRNAs have remained elusive. In this study, we compared the sequence evolution and biological features between single-exonic lincRNAs and multi-exonic lincRNAs (SELs and MELs, respectively) that were present only in the hominoids (hominoid-specific) or conserved in primates (primate-conserved). The MEL exons were further classified into alternatively spliced exons (ASEs) and constitutively spliced exons (CSEs) for evolutionary analyses. Our results indicate that SELs and MELs differed significantly from each other. Firstly, in hominoid-specific lincRNAs, MELs (both CSEs and ASEs) evolved slightly more rapidly than SELs, which evolved approximately at the neutral rate. In primate-conserved lincRNAs, SELs and ASEs evolved slightly more slowly than CSEs and neutral sequences. The evolutionary path of hominid-specific lincRNAs thus seemed to have diverged from that of their more ancestral counterparts. Secondly, both of the exons and transcripts of SELs were significantly longer than those of MELs, and this was probably because SEL transcripts were more resistant to RNA splicing than MELs. Thirdly, SELs were physically closer to coding genes than MELs. Fourthly, SELs were more widely expressed in human tissues than MELs. These results suggested that SELs and MELs represented two biologically distinct groups of genes. In addition, the SEL-MEL and ASE-CSE differences implied that splicing might be important for the functionality or regulations of lincRNAs in primates.

Role of Pnn in alternative splicing of a specific subset of lncRNAs of the corneal epithelium

PURPOSE

GG-H whole transcriptome array analysis suggested involvement of PININ (PNN) in the alternative splicing of multiple long non-coding RNAs (lncRNAs). To further investigate PNN's role in regulating the alternative splicing of lncRNAs in a corneal epithelial context, we performed detailed analyses for detecting and identifying alternatively spliced lncRNAs.

METHODS

Total RNA was isolated from PNN knockdown human corneal epithelial (HCET) cells or Pnn-deficient mouse corneas, and subjected to real-time-PCR (RT-PCR) assays, and the alternatively spliced lncRNAs were counted. Alternatively spliced lncRNAs were detected with in situ hybridization with variant-specific RNA probes on human cornea sections.

RESULTS

Our analysis uncovered PNN's impact on the transcript levels of several lncRNAs including Linc00085 and HAS2-AS1. Interestingly, a mouse ortholog of HAS2-AS1, Has2as, clearly exhibited a differential splicing pattern among three major splice variants in the Pnn-deficient mouse cornea. The sequence analyses and quantification of splice variants of candidate lncRNAs, including RP11-295B20.2, RP11-18I14.1, and RP11-322M19.1, demonstrated complex configuration of their splicing changes, with a significant impact of PNN on the process. Knockdown of PNN in HCET cells led to specific changes in the inclusion of multiple cassette exons as well as in the use of alternative splice sites in RP11-322M19.1 and RP11-18I14.1, resulting in considerable net changes in the ratio between the splice variants. Finally, in situ hybridization analyses revealed the presence of RP11-295G20.2 in the nuclei of corneal epithelial cells, but not in the stromal cells of the human cornea, while RP11-322M19.1 was present in epithelial and non-epithelial cells.

CONCLUSIONS

The data suggest PNN's role in the alternative splicing of a specific subset of lncRNAs might have a significant impact on the corneal epithelium.

Polyploidy and novelty: Gottlieb's legacy

Nearly four decades ago, Roose & Gottlieb (Roose & Gottlieb 1976 Evolution 30, 818-830. (doi:10.2307/2407821)) showed that the recently derived allotetraploids Tragopogon mirus and T. miscellus combined the allozyme profiles of their diploid parents (T. dubius and T. porrifolius, and T. dubius and T. pratensis, respectively). This classic paper addressed the link between genotype and biochemical phenotype and documented enzyme additivity in allopolyploids. Perhaps more important than their model of additivity, however, was their demonstration of novelty at the biochemical level. Enzyme multiplicity-the production of novel enzyme forms in the allopolyploids-can provide an extensive array of polymorphism for a polyploid individual and may explain, for example, the expanded ranges of polyploids relative to their diploid progenitors. In this paper, we extend the concept of evolutionary novelty in allopolyploids to a range of genetic and ecological features. We observe that the dynamic nature of polyploid genomes-with alterations in gene content, gene number, gene arrangement, gene expression and transposon activity-may generate sufficient novelty that every individual in a polyploid population or species may be unique. Whereas certain combinations of these features will undoubtedly be maladaptive, some unique combinations of newly generated variation may provide tremendous evolutionary potential and adaptive capabilities.

Association of eight EST-derived SNPs with carcass and meat quality traits in pigs

The identification of genetic markers associated with important economic traits is fundamental to improving the productivity and quality of livestock. In this investigation, we searched for 177 expressed sequence tags (ESTs) putatively involved in meat quality from the available pig EST database, and detected eight single nucleotide polymorphisms (SNPs) in eight ESTs. We investigated the associations of these SNPs with 18 carcass and meat quality traits in a Landrace × Lantang F2 resource population (n = 257). Association analysis revealed that seven SNPs (except E42) were associated with some of the carcass- and meat quality-related traits. Particularly, significant associations of three SNPs (E53, E82, and E36) with backfat thickness traits were observed. Further, the genetic effects of E53 on four live backfat thickness traits were validated in an independent population (n = 221). More investigations about E53 sequence characteristics were performed, i.e., radiation hybrid (RH) mapping, 3'-RACE, and screening analysis of the positive BAC clones. Our research identified the genetic effects of eight EST-derived SNPs on carcass and meat quality traits, and suggested that E53 may be a useful marker for live backfat thickness traits in pig breeding programs.

NFAT1 and NFAT3 cooperate with HDAC4 during regulation of alternative splicing of PMCA isoforms in PC12 cells

Background

The bulk of human genes undergo alternative splicing (AS) upon response to physiological stimuli. AS is a great source of protein diversity and biological processes and is associated with the development of many diseases. Pheochromocytoma is a neuroendocrine tumor, characterized by an excessive Ca²⁺-dependent secretion of catecholamines. This underlines the importance of balanced control of calcium transport via regulation of gene expression pattern, including different calcium transport systems, such as plasma membrane Ca²⁺-ATPases (PMCAs), abundantly expressed in pheochromocytoma chromaffin cells (PC12 cells). PMCAs are encoded by four genes (Atp2b1, Atp2b2, Atp2b3, Atp2b4), whose transcript products undergo alternative splicing giving almost 30 variants.

Results

In this scientific report, we propose a novel mechanism of regulation of PMCA alternative splicing in PC12 cells through cooperation of the nuclear factor of activated T-cells (NFAT) and histone deacetylases (HDACs). Luciferase assays showed increased activity of NFAT in PC12 cells, which was associated with altered expression of PMCA. RT-PCR experiments suggested that inhibition of the transcriptional activity of NFAT might result in the rearrangement of PMCA splicing variants in PC12 cells. NFAT inhibition led to dominant expression of 2x/c, 3x/a and 4x/a PMCA variants, while in untreated cells the 2w,z/b, 3z,x/b,c,e,f, and 4x/b variants were found as well. Furthermore, chromatin immunoprecipitation experiments showed that NFAT1-HDAC4 or NFAT3-HDAC4 complexes might be involved in regulation of PMCA2x splicing variant generation.

Conclusions

We suggest that the influence of NFAT/HDAC on PMCA isoform composition might be important for altered dopamine secretion by PC12 cells.