Estimation of alternative splicing variability in human populations

Differences in alternative splicing events in the adaptive strategies of two Daphnia galeata genotypes induced by fish kairomones

BACKGROUND

Daphnia galeata is a suitable model organism for investigating predator-induced defense. Genes and pathways exhibiting differential expression between fish kairomone-treated and untreated groups in D. galeata have been identified. However, understanding of the significance of alternative splicing, a crucial process of the regulation of gene expression in eukaryotes, to this mechanism remains limited. This study measured life-history traits and conducted short-read RNA sequencing and long-read isoform sequencing of two Korean D. galeata genotypes (KB1 and KE2) to uncover the genetic mechanism underlying their phenotypic plasticity under predation stress.

RESULTS

KB1 exhibited strategies to enhance fertility and decrease body length when exposed to fish kairomones, while KE2 deployed an adaptive strategy to increase body length. Full-length transcriptomes from KB1 and KE2 yielded 65,736 and 57,437 transcripts, respectively, of which 32 differentially expressed transcripts (DETs) were shared under predation stress across both genotypes. Prominent DETs common to both genotypes were related to energy metabolism and the immune system. Additionally, differential alternative splicing (DAS) events were detected in both genotypes in response to fish kairomones. DAS genes shared between both genotypes may indicate their significant role in the post-transcriptional stress response to fish predation. Calpain-3, involved in digestion and nutrient absorption, was identified as a DAS gene in both genotypes when exposed to fish kairomones. In addition, the gene encoding thymosin beta, which is related to growth, was found to be a statistically significant DAS only in KB1, while that encoding ultraspiracle protein, also associated with growth, was only identified in KE2. Moreover, transcripts encoding proteins such as EGF-like domain-containing protein, vitellogenin fused with superoxide dismutase, and others were identified overlapping between DAS events and DETs and potentially elucidating their association with the observed phenotypic variation in each genotype.

CONCLUSIONS

Our findings highlight the crucial role of alternative splicing in modulating transcriptome landscape under predation stress in D. galeata, emphasizing the requirement for integrating gene expression and splicing analyses in evolutionary adaptation studies.

A fast non-parametric test of association for multiple traits

The increasing availability of multidimensional phenotypic data in large cohorts of genotyped individuals requires efficient methods to identify genetic effects on multiple traits. Permutational multivariate analysis of variance (PERMANOVA) offers a powerful non-parametric approach. However, it relies on permutations to assess significance, which hinders the analysis of large datasets. Here, we derive the limiting null distribution of the PERMANOVA test statistic, providing a framework for the fast computation of asymptotic p values. Our asymptotic test presents controlled type I error and high power, often outperforming parametric approaches. We illustrate its applicability in the context of QTL mapping and GWAS.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Serum lysyl oxidase concentration increases in long-standing systemic sclerosis: Can lysyl oxidase change over time?

Objectives

This study aims to investigate the association of serum lysyl oxidase (LOX) levels with systemic sclerosis (SSc), to examine the relationship between LOX and disease onset, and to evaluate the probable effects of hyperlipidemia on the circulating levels of LOX among patients with SSc.

Patients and methods

Between May 2017 and November 2018, a total of 39 patients with SSc (2 males, 37 females; mean age: 46.6±12.3 years; range, 18 to 65 years) and 35 healthy controls (4 males, 31 females; mean age: 43.1±14.1 years; range, 18 to 65 years) were included. Serum LOX concentration was measured using the enzyme-linked immunoassay in triplicate.

Results

We found higher levels of serum LOX in patients with SSc compared to healthy controls. There was a significant relationship between serum LOX levels and disease onset. Patients with long-standing disease demonstrated increased levels of LOX in the blood compared to the recent-onset group. Hyperlipidemia did not have a significant effect on circulating levels of LOX. There was a significant negative correlation between LOX levels and modified Rodnan Skin Score in the subgroup of patients with skin involvement only and in patients without gastrointestinal involvement.

Conclusion

Our study findings show an increased level of LOX protein level in the blood of patients diagnosed with SSc. Hyperlipidemia seems not to affect the concentrations of LOX in the peripheral blood of patients with SSc.

Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome

Alternative splicing (AS) is a fundamental step in eukaryotic mRNA biogenesis. Here, we develop an efficient and reproducible pipeline for the discovery of genetic variants that affect AS (splicing QTLs, sQTLs). We use it to analyze the GTEx dataset, generating a comprehensive catalog of sQTLs in the human genome. Downstream analysis of this catalog provides insight into the mechanisms underlying splicing regulation. We report that a core set of sQTLs is shared across multiple tissues. sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many also affect the expression of the same or other genes, uncovering regulatory loci that act through different mechanisms. sQTLs tend to be located in post-transcriptionally spliced introns, which would function as hotspots for splicing regulation. While many variants affect splicing patterns by altering the sequence of splice sites, many more modify the binding sites of RNA-binding proteins. Genetic variants affecting splicing can have a stronger phenotypic impact than those affecting gene expression.

RNA Sequencing of Human Peripheral Blood Cells Indicates Upregulation of Immune-Related Genes in Huntington's Disease

Huntington's disease (HD) is an autosomal dominantly inherited neurodegenerative disorder caused by a trinucleotide repeat expansion in the Huntingtin gene. As disease-modifying therapies for HD are being developed, peripheral blood cells may be used to indicate disease progression and to monitor treatment response. In order to investigate whether gene expression changes can be found in the blood of individuals with HD that distinguish them from healthy controls, we performed transcriptome analysis by next-generation sequencing (RNA-seq). We detected a gene expression signature consistent with dysregulation of immune-related functions and inflammatory response in peripheral blood from HD cases vs. controls, including induction of the interferon response genes, IFITM3, IFI6 and IRF7. Our results suggest that it is possible to detect gene expression changes in blood samples from individuals with HD, which may reflect the immune pathology associated with the disease.

A limited set of transcriptional programs define major cell types

We have produced RNA sequencing data for 53 primary cells from different locations in the human body. The clustering of these primary cells reveals that most cells in the human body share a few broad transcriptional programs, which define five major cell types: epithelial, endothelial, mesenchymal, neural, and blood cells. These act as basic components of many tissues and organs. Based on gene expression, these cell types redefine the basic histological types by which tissues have been traditionally classified. We identified genes whose expression is specific to these cell types, and from these genes, we estimated the contribution of the major cell types to the composition of human tissues. We found this cellular composition to be a characteristic signature of tissues and to reflect tissue morphological heterogeneity and histology. We identified changes in cellular composition in different tissues associated with age and sex, and found that departures from the normal cellular composition correlate with histological phenotypes associated with disease.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure

BACKGROUND

Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses.

RESULTS

By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing.

CONCLUSIONS

Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Enteric infection induces Lark-mediated intron retention at the 5' end of Drosophila genes

BACKGROUND

RNA splicing is a key post-transcriptional mechanism that generates protein diversity and contributes to the fine-tuning of gene expression, which may facilitate adaptation to environmental challenges. Here, we employ a systems approach to study alternative splicing changes upon enteric infection in females from classical Drosophila melanogaster strains as well as 38 inbred lines.

RESULTS

We find that infection leads to extensive differences in isoform ratios, which results in a more diverse transcriptome with longer 5' untranslated regions (5'UTRs). We establish a role for genetic variation in mediating inter-individual splicing differences, with local splicing quantitative trait loci (local-sQTLs) being preferentially located at the 5' end of transcripts and directly upstream of splice donor sites. Moreover, local-sQTLs are more numerous in the infected state, indicating that acute stress unmasks a substantial number of silent genetic variants. We observe a general increase in intron retention concentrated at the 5' end of transcripts across multiple strains, whose prevalence scales with the degree of pathogen virulence. The length, GC content, and RNA polymerase II occupancy of these introns with increased retention suggest that they have exon-like characteristics. We further uncover that retained intron sequences are enriched for the Lark/RBM4 RNA binding motif. Interestingly, we find that lark is induced by infection in wild-type flies, its overexpression and knockdown alter survival, and tissue-specific overexpression mimics infection-induced intron retention.

CONCLUSION

Our collective findings point to pervasive and consistent RNA splicing changes, partly mediated by Lark/RBM4, as being an important aspect of the gut response to infection.

Transcriptome variation in human populations and its potential application in forensics

This review presents the state-of-the-art in the forensic application of genetic methods driven by the research in population transcriptomics. In the first part of the review, the constraints of using classical genomic markers are shortly reviewed. In the second part, the developments in the field of inter-population diversity at the transcriptomic level are presented. Subsequently, a potential of population-specific transcriptomic markers in forensic science applications, including ascertaining population affiliation of human samples and cell mixtures separation, are presented.

Gene expression variability: the other dimension in transcriptome analysis

Transcriptome sequencing is a powerful technique to study molecular changes that underlie the differences in physiological conditions and disease progression. A typical question that is posed in such studies is finding genes with significant changes between sample groups. In this respect expression variability is regarded as a nuisance factor that is primarily of technical origin and complicates the data analysis. However, it is becoming apparent that the biological variation in gene expression might be an important molecular phenotype that can affect physiological parameters. In this review we explore the recent literature on technical and biological variability in gene expression, sources of expression variability, (epi-)genetic hallmarks, and evolutionary constraints in genes with robust and variable gene expression. We provide an overview of recent findings on effects of external cues, such as diet and aging, on expression variability and on other biological phenomena that can be linked to it. We discuss metrics and tools that were developed for quantification of expression variability and highlight the importance of future studies in this direction. To assist the adoption of expression variability analysis, we also provide a detailed description and computer code, which can easily be utilized by other researchers. We also provide a reanalysis of recently published data to highlight the value of the analysis method.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

Towards a complete map of the human long non-coding RNA transcriptome

Gene maps, or annotations, enable us to navigate the functional landscape of our genome. They are a resource upon which virtually all studies depend, from single-gene to genome-wide scales and from basic molecular biology to medical genetics. Yet present-day annotations suffer from trade-offs between quality and size, with serious but often unappreciated consequences for downstream studies. This is particularly true for long non-coding RNAs (lncRNAs), which are poorly characterized compared to protein-coding genes. Long-read sequencing technologies promise to improve current annotations, paving the way towards a complete annotation of lncRNAs expressed throughout a human lifetime.

Alternatively spliced variants in Atlantic cod (Gadus morhua) support response to variable salinity environment

Analysis of gill transcriptome of the Atlantic cod from the Baltic Sea demonstrated that alternatively spliced (AS) variants may be actively involved in the process of adaptation to altered salinity. Some AS variants of different genes, like phospholipase A2 group IVC (PLA2G4C), appeared only in fish exposed to altered salinity, while other isoforms of the same genes were present in all experimental groups. Novel sequence arrangements represent 89% of all AS in the Baltic cod compared to the Atlantic population. Profiles of modified pathways suggest that regulation by AS can afford specific changes of genes expressed in response to the environment. The AS variants appear to be involved in the response to stress by modifications of signalling in apoptosis pathways, an innate immunological response and pro-inflammatory process. Present results support the hypothesis that developing new AS variants could support genome complexity and reinforce the ability to fast adapt to local environments.

Changes in Alternative Splicing as Pharmacodynamic Markers for Sudemycin D6

OBJECTIVE

The aim of the study was to define pharmacodynamic markers for sudemycin D6, an experimental cancer drug that changes alternative splicing in human blood.

METHODS

Blood samples from 12 donors were incubated with sudemycin D6 for up to 24 hours, and at several time points total RNA from lymphocytes was prepared and the pre-messenger RNA (mRNA) splicing patterns were analyzed with reverse transcription-polymerase chain reaction.

RESULTS

Similar to immortalized cells, blood lymphocytes change alternative splicing due to sudemycin D6 treatment. However, lymphocytes in blood respond slower than immortalized cultured cells.

CONCLUSIONS

Exon skipping in the DUSP11 and SRRM1 pre-mRNAs are pharmacodynamic markers for sudemycin D6 treatment and show effects beginning at 9 hours after treatment.

snoRNA and piRNA expression levels modified by tobacco use in women with lung adenocarcinoma

Lung cancer is one of the most frequent types of cancer worldwide. Most patients are diagnosed at advanced stage and thus have poor prognosis. Smoking is a risk factor for lung cancer, however most smokers do not develop lung cancer while 20% of women with lung adenocarcinoma are non-smokers. Therefore, it is possible that these two groups present differences besides the smoking status, including differences in their gene expression signature. The altered expression patterns of non-coding RNAs in complex diseases make them potential biomarkers for diagnosis and treatment. We analyzed data from differentially and constitutively expressed PIWI-interacting RNAs and small nucleolar RNAs from publicly available small RNA high-throughput sequencing data in search of an expression pattern of non-coding RNA that could differentiate these two groups. Here, we report two sets of differentially expressed small non-coding RNAs identified in normal and tumoral tissues of women with lung adenocarcinoma, that discriminate between smokers and non-smokers. Our findings may offer new insights on metabolic alterations caused by tobacco and may be used for early diagnosis of lung cancer.

Evidence of selection on splicing-associated loci in human populations and relevance to disease loci mapping

We performed a whole-genome scan of genetic variants in splicing regulatory elements (SREs) and evaluated the extent to which natural selection has shaped extant patterns of variation in SREs. We investigated the degree of differentiation of single nucleotide polymorphisms (SNPs) in SREs among human populations and applied long-range haplotype- and multilocus allelic differentiation-based methods to detect selection signatures. We describe an approach, sampling a large number of loci across the genome from functional classes and using the consensus from multiple tests, for identifying candidates for selection signals. SRE SNPs in various SNP functional classes show different patterns of population differentiation compared with their non-SRE counterparts. Intronic regions display a greater enrichment for extreme population differentiation among the potentially tissue-dependent transcript ratio quantitative trait loci (trQTLs) than SRE SNPs in general and includ outlier trQTLs for cross-population composite likelihood ratio, suggesting that incorporation of context annotation for regulatory variation may lead to improved detection of signature of selection on these loci. The proportion of extremely rare SNPs disrupting SREs is significantly higher in European than in African samples. The approach developed here will be broadly useful for studies of function and disease-associated variation in the human genome.

Translational contributions to tissue specificity in rhythmic and constitutive gene expression

Background

The daily gene expression oscillations that underlie mammalian circadian rhythms show striking differences between tissues and involve post-transcriptional regulation. Both aspects remain poorly understood. We have used ribosome profiling to explore the contribution of translation efficiency to temporal gene expression in kidney and contrasted our findings with liver data available from the same mice.

Results

Rhythmic translation of constantly abundant messenger RNAs (mRNAs) affects largely non-overlapping transcript sets with distinct phase clustering in the two organs. Moreover, tissue differences in translation efficiency modulate the timing and amount of protein biosynthesis from rhythmic mRNAs, consistent with organ specificity in clock output gene repertoires and rhythmicity parameters. Our comprehensive datasets provided insights into translational control beyond temporal regulation. Between tissues, many transcripts show differences in translation efficiency, which are, however, of markedly smaller scale than mRNA abundance differences. Tissue-specific changes in translation efficiency are associated with specific transcript features and, intriguingly, globally counteracted and compensated transcript abundance variations, leading to higher similarity at the level of protein biosynthesis between both tissues.

Conclusions

We show that tissue specificity in rhythmic gene expression extends to the translatome and contributes to define the identities, the phases and the expression levels of rhythmic protein biosynthesis. Moreover, translational compensation of transcript abundance divergence leads to overall higher similarity at the level of protein production across organs. The unique resources provided through our study will serve to address fundamental questions of post-transcriptional control and differential gene expression in vivo.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1222-2) contains supplementary material, which is available to authorized users.

Clustering of mRNA-Seq data based on alternative splicing patterns

Sequencing of messenger RNA (mRNA) can provide estimates of the levels of individual isoforms within the cell. It remains to adapt many standard statistical methods commonly used for analyzing gene expression levels to take advantage of this additional information. One novel question is whether we can find clusters of samples that are distinguished not by their gene expression but by their isoform usage. We propose a novel approach for clustering mRNA-Seq data that identifies such clusters. We show via simulation that our methods are more sensitive to finding clusters based on isoform usage than standard clustering techniques. We demonstrate its performance by finding a technical artifact that resulted in different batches having different isoform usage patterns, and illustrate its usage on several The Cancer Genome Atlas datasets.

Dissecting the Effect of Genetic Variation on the Hepatic Expression of Drug Disposition Genes across the Collaborative Cross Mouse Strains

A central challenge in pharmaceutical research is to investigate genetic variation in response to drugs. The Collaborative Cross (CC) mouse reference population is a promising model for pharmacogenomic studies because of its large amount of genetic variation, genetic reproducibility, and dense recombination sites. While the CC lines are phenotypically diverse, their genetic diversity in drug disposition processes, such as detoxification reactions, is still largely uncharacterized. Here we systematically measured RNA-sequencing expression profiles from livers of 29 CC lines under baseline conditions. We then leveraged a reference collection of metabolic biotransformation pathways to map potential relations between drugs and their underlying expression quantitative trait loci (eQTLs). By applying this approach on proximal eQTLs, including eQTLs acting on the overall expression of genes and on the expression of particular transcript isoforms, we were able to construct the organization of hepatic eQTL-drug connectivity across the CC population. The analysis revealed a substantial impact of genetic variation acting on drug biotransformation, allowed mapping of potential joint genetic effects in the context of individual drugs, and demonstrated crosstalk between drug metabolism and lipid metabolism. Our findings provide a resource for investigating drug disposition in the CC strains, and offer a new paradigm for integrating biotransformation reactions to corresponding variations in DNA sequences.

Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans

Background

Long non-coding RNAs (lncRNAs) are increasingly implicated as gene regulators and may ultimately be more numerous than protein-coding genes in the human genome. Despite large numbers of reported lncRNAs, reference annotations are likely incomplete due to their lower and tighter tissue-specific expression compared to mRNAs. An unexplored factor potentially confounding lncRNA identification is inter-individual expression variability. Here, we characterize lncRNA natural expression variability in human primary granulocytes.

Results

We annotate granulocyte lncRNAs and mRNAs in RNA-seq data from 10 healthy individuals, identifying multiple lncRNAs absent from reference annotations, and use this to investigate three known features (higher tissue-specificity, lower expression, and reduced splicing efficiency) of lncRNAs relative to mRNAs. Expression variability was examined in seven individuals sampled three times at 1- or more than 1-month intervals. We show that lncRNAs display significantly more inter-individual expression variability compared to mRNAs. We confirm this finding in two independent human datasets by analyzing multiple tissues from the GTEx project and lymphoblastoid cell lines from the GEUVADIS project. Using the latter dataset we also show that including more human donors into the transcriptome annotation pipeline allows identification of an increasing number of lncRNAs, but minimally affects mRNA gene number.

Conclusions

A comprehensive annotation of lncRNAs is known to require an approach that is sensitive to low and tight tissue-specific expression. Here we show that increased inter-individual expression variability is an additional general lncRNA feature to consider when creating a comprehensive annotation of human lncRNAs or proposing their use as prognostic or disease markers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0873-8) contains supplementary material, which is available to authorized users.

Long non-coding RNA expression in primary human monocytes

Long non-coding RNAs (lncRNAs) have been shown to possess a wide range of functions in both cellular and developmental processes including cancers. Although some of the lncRNAs have been implicated in the regulation of the immune response, the exact function of the large majority of lncRNAs still remains unknown. In this study, we characterized the lncRNAs in human primary monocytes, an essential component of the innate immune system. We performed RNA sequencing of monocytes from four individuals and combined our data with eleven other publicly available datasets. Our analysis led to identification of ~8000 lncRNAs of which >1000 have not been previously reported in monocytes. PCR-based validation of a subset of the identified novel long intergenic noncoding RNAs (lincRNAs) revealed distinct expression patterns. Our study provides a landscape of lncRNAs in monocytes, which could facilitate future experimental studies to characterize the functions of these molecules in the innate immune system.

Tissue Restricted Splice Junctions Originate Not Only from Tissue-Specific Gene Loci, but Gene Loci with a Broad Pattern of Expression

Cellular mechanisms that achieve protein diversity in eukaryotes are multifaceted, including transcriptional components such as RNA splicing. Through alternative splicing, a single protein-coding gene can generate multiple mRNA transcripts and protein isoforms, some of which are tissue-specific. We have conducted qualitative and quantitative analyses of the Bodymap 2.0 messenger RNA-sequencing data from 16 human tissue samples and identified 209,363 splice junctions. Of these, 22,231 (10.6%) were not previously annotated and 21,650 (10.3%) were expressed in a tissue-restricted pattern. Tissue-restricted alternative splicing was found to be widespread, with approximately 65% of expressed multi-exon genes containing at least one tissue-specific splice junction. Interestingly, we observed many tissue-specific splice junctions not only in genes expressed in one or a few tissues, but also from gene loci with a broad pattern of expression.

Regulatory Divergence of Transcript Isoforms in a Mammalian Model System

Phenotypic differences between species are driven by changes in gene expression and, by extension, by modifications in the regulation of the transcriptome. Investigation of mammalian transcriptome divergence has been restricted to analysis of bulk gene expression levels and gene-internal splicing. Using allele-specific expression analysis in inter-strain hybrids of Mus musculus, we determined the contribution of multiple cellular regulatory systems to transcriptome divergence, including: alternative promoter usage, transcription start site selection, cassette exon usage, alternative last exon usage, and alternative polyadenylation site choice. Between mouse strains, a fifth of genes have variations in isoform usage that contribute to transcriptomic changes, half of which alter encoded amino acid sequence. Virtually all divergence in isoform usage altered the post-transcriptional regulatory instructions in gene UTRs. Furthermore, most genes with isoform differences between strains contain changes originating from multiple regulatory systems. This result indicates widespread cross-talk and coordination exists among different regulatory systems. Overall, isoform usage diverges in parallel with and independently to gene expression evolution, and the cis and trans regulatory contribution to each differs significantly.

Human genomics. The human transcriptome across tissues and individuals

Transcriptional regulation and posttranscriptional processing underlie many cellular and organismal phenotypes. We used RNA sequence data generated by Genotype-Tissue Expression (GTEx) project to investigate the patterns of transcriptome variation across individuals and tissues. Tissues exhibit characteristic transcriptional signatures that show stability in postmortem samples. These signatures are dominated by a relatively small number of genes—which is most clearly seen in blood—though few are exclusive to a particular tissue and vary more across tissues than individuals. Genes exhibiting high interindividual expression variation include disease candidates associated with sex, ethnicity, and age. Primary transcription is the major driver of cellular specificity, with splicing playing mostly a complementary role; except for the brain, which exhibits a more divergent splicing program. Variation in splicing, despite its stochasticity, may play in contrast a comparatively greater role in defining individual phenotypes.

RNA-Seq: Improving Our Understanding of Retinal Biology and Disease

Over the past several years, rapid technological advances have allowed for a dramatic increase in our knowledge and understanding of the transcriptional landscape, because of the ability to study gene expression in greater depth and with more detail than previously possible. To this end, RNA-Seq has quickly become one of the most widely used methods for studying transcriptomes of tissues and individual cells. Unlike previously favored analysis methods, RNA-Seq is extremely high-throughput, and is not dependent on an annotated transcriptome, laying the foundation for novel genetic discovery. Additionally, RNA-Seq derived transcriptomes provide a basis for widening the scope of research to identify potential targets in the treatment of retinal disease.

A generalized dSpliceType framework to detect differential splicing and differential expression events using RNA-Seq

Transcriptomes are routinely compared in term of a list of differentially expressed genes followed by functional enrichment analysis. Due to the technology limitations of microarray, the molecular mechanisms of differential expression is poorly understood. Using RNA-seq data, we propose a generalized dSpliceType framework to systematically investigate the synergistic and antagonistic effects of differential splicing and differential expression. We applied the method to two public RNA-seq data sets and compared the transcriptomes between treatment and control conditions. The generalized dSpliceType detects and prioritizes a list of genes that are differentially expressed and/or spliced. In particular, the multivariate dSpliceType is among the fist to utilize sequential dependency of normalized base-wise read coverage signals and capture biological variability among replicates using a multivariate statistical model. We compared dSpliceType with two other methods in terms of five most common types of differential splicing events between two conditions using RNA-Seq. dSpliceType is free, available from http://dsplicetype.sourceforge.net/.

Analysis of long non-coding RNAs highlights tissue-specific expression patterns and epigenetic profiles in normal and psoriatic skin

BACKGROUND

Although analysis pipelines have been developed to use RNA-seq to identify long non-coding RNAs (lncRNAs), inference of their biological and pathological relevance remains a challenge. As a result, most transcriptome studies of autoimmune disease have only assessed protein-coding transcripts.

RESULTS

We used RNA-seq data from 99 lesional psoriatic, 27 uninvolved psoriatic, and 90 normal skin biopsies, and applied computational approaches to identify and characterize expressed lncRNAs. We detect 2,942 previously annotated and 1,080 novel lncRNAs which are expected to be skin specific. Notably, over 40% of the novel lncRNAs are differentially expressed and the proportions of differentially expressed transcripts among protein-coding mRNAs and previously-annotated lncRNAs are lower in psoriasis lesions versus uninvolved or normal skin. We find that many lncRNAs, in particular those that are differentially expressed, are co-expressed with genes involved in immune related functions, and that novel lncRNAs are enriched for localization in the epidermal differentiation complex. We also identify distinct tissue-specific expression patterns and epigenetic profiles for novel lncRNAs, some of which are shown to be regulated by cytokine treatment in cultured human keratinocytes.

CONCLUSIONS

Together, our results implicate many lncRNAs in the immunopathogenesis of psoriasis, and our results provide a resource for lncRNA studies in other autoimmune diseases.

Enhanced transcriptome maps from multiple mouse tissues reveal evolutionary constraint in gene expression

Mice have been a long-standing model for human biology and disease. Here we characterize, by RNA sequencing, the transcriptional profiles of a large and heterogeneous collection of mouse tissues, augmenting the mouse transcriptome with thousands of novel transcript candidates. Comparison with transcriptome profiles in human cell lines reveals substantial conservation of transcriptional programmes, and uncovers a distinct class of genes with levels of expression that have been constrained early in vertebrate evolution. This core set of genes captures a substantial fraction of the transcriptional output of mammalian cells, and participates in basic functional and structural housekeeping processes common to all cell types. Perturbation of these constrained genes is associated with significant phenotypes including embryonic lethality and cancer. Evolutionary constraint in gene expression levels is not reflected in the conservation of the genomic sequences, but is associated with conserved epigenetic marking, as well as with characteristic post-transcriptional regulatory programme, in which sub-cellular localization and alternative splicing play comparatively large roles.

The analysis of mammalian transcriptomes could provide new insights into human biology. Here the authors carry out RNA sequencing in a large collection of mouse tissues and compare these data to human transcriptome profiles, identifying a set of constrained genes that carry out basic cellular functions with remarkably constant expression levels across tissues and species.

Variation in genomic landscape of clear cell renal cell carcinoma across Europe

The incidence of renal cell carcinoma (RCC) is increasing worldwide, and its prevalence is particularly high in some parts of Central Europe. Here we undertake whole-genome and transcriptome sequencing of clear cell RCC (ccRCC), the most common form of the disease, in patients from four different European countries with contrasting disease incidence to explore the underlying genomic architecture of RCC. Our findings support previous reports on frequent aberrations in the epigenetic machinery and PI3K/mTOR signalling, and uncover novel pathways and genes affected by recurrent mutations and abnormal transcriptome patterns including focal adhesion, components of extracellular matrix (ECM) and genes encoding FAT cadherins. Furthermore, a large majority of patients from Romania have an unexpected high frequency of A:T>T:A transversions, consistent with exposure to aristolochic acid (AA). These results show that the processes underlying ccRCC tumorigenesis may vary in different populations and suggest that AA may be an important ccRCC carcinogen in Romania, a finding with major public health implications.

RNA-Seq Data: A Complexity Journey

A paragraph from the highlights of “Transcriptomics: Throwing light on dark matter” by L. Flintoft (Nature Reviews Genetics 11, 455, 2010), says: “Reports over the past few years of extensive transcription throughout eukaryotic genomes have led to considerable excitement. However, doubts have been raised about the methods that have detected this pervasive transcription and about how much of it is functional.” Since the appearance of the ENCODE project and due to follow-up work, a shift from the pervasive transcription observed from RNA-Seq data to its functional validation is gradually occurring. However, much less attention has been turned to the problem of deciphering the complexity of transcriptome data, which determines uncertainty with regard to identification, quantification and differential expression of genes and non-coding RNAs. The aim of this mini-review is to emphasize transcriptome-related problems of direct and inverse nature for which novel inference approaches are needed.

Identification of genetic variants associated with alternative splicing using sQTLseekeR

Identification of genetic variants affecting splicing in RNA sequencing population studies is still in its infancy. Splicing phenotype is more complex than gene expression and ought to be treated as a multivariate phenotype to be recapitulated completely. Here we represent the splicing pattern of a gene as the distribution of the relative abundances of a gene’s alternative transcript isoforms. We develop a statistical framework that uses a distance-based approach to compute the variability of splicing ratios across observations, and a non-parametric analogue to multivariate analysis of variance. We implement this approach in the R package sQTLseekeR and use it to analyze RNA-Seq data from the Geuvadis project in 465 individuals. We identify hundreds of single nucleotide polymorphisms (SNPs) as splicing QTLs (sQTLs), including some falling in genome-wide association study SNPs. By developing the appropriate metrics, we show that sQTLseekeR compares favorably with existing methods that rely on univariate approaches, predicting variants that behave as expected from mutations affecting splicing.

RNA sequencing has enabled the global analysis of both gene expression levels and splicing events. Here, the authors develop a multivariate approach that is able to identify SNPs that influence splicing, and investigate the overlap of these with functional domains across the genome, including previously identified GWAS signals.

Transcriptome sequencing from diverse human populations reveals differentiated regulatory architecture

Large-scale sequencing efforts have documented extensive genetic variation within the human genome. However, our understanding of the origins, global distribution, and functional consequences of this variation is far from complete. While regulatory variation influencing gene expression has been studied within a handful of populations, the breadth of transcriptome differences across diverse human populations has not been systematically analyzed. To better understand the spectrum of gene expression variation, alternative splicing, and the population genetics of regulatory variation in humans, we have sequenced the genomes, exomes, and transcriptomes of EBV transformed lymphoblastoid cell lines derived from 45 individuals in the Human Genome Diversity Panel (HGDP). The populations sampled span the geographic breadth of human migration history and include Namibian San, Mbuti Pygmies of the Democratic Republic of Congo, Algerian Mozabites, Pathan of Pakistan, Cambodians of East Asia, Yakut of Siberia, and Mayans of Mexico. We discover that approximately 25.0% of the variation in gene expression found amongst individuals can be attributed to population differences. However, we find few genes that are systematically differentially expressed among populations. Of this population-specific variation, 75.5% is due to expression rather than splicing variability, and we find few genes with strong evidence for differential splicing across populations. Allelic expression analyses indicate that previously mapped common regulatory variants identified in eight populations from the International Haplotype Map Phase 3 project have similar effects in our seven sampled HGDP populations, suggesting that the cellular effects of common variants are shared across diverse populations. Together, these results provide a resource for studies analyzing functional differences across populations by estimating the degree of shared gene expression, alternative splicing, and regulatory genetics across populations from the broadest points of human migration history yet sampled.

Previous gene expression studies have identified factors influencing population-level variation in gene regulation. However, these efforts have been limited to a small set of well-studied populations. By leveraging the high resolution of RNA sequencing and broad population sampling, we survey the landscape of transcriptome variation across a globally distributed set of seven populations that span a breadth of human genetic variation and major dispersal events. We assess differences in gene expression, transcript structure, and regulatory variation. We find only 44 transcripts that show significant differences in expression, likely as a result of the small sample size, but we find that 25% of the variance in gene expression is due to population differences. This is a larger fraction than previously observed, and it is likely due to the greater breadth of human diversity assayed in this study. We also find that population-specific variance is mostly due to transcription variability rather than the configuration of expressed gene products. Additionally, known common regulatory variants have similar effects across populations including those we study here. These data and results serve as a resource cataloging the wide array of gene expression regulation affecting population variation among diverse groups, improving our understanding of transcriptional diversity.

Expression quantitative trait analyses to identify causal genetic variants for type 2 diabetes susceptibility

Type 2 diabetes (T2D) is a common metabolic disorder which is caused by multiple genetic perturbations affecting different biological pathways. Identifying genetic factors modulating the susceptibility of this complex heterogeneous metabolic phenotype in different ethnic and racial groups remains challenging. Despite recent success, the functional role of the T2D susceptibility variants implicated by genome-wide association studies (GWAS) remains largely unknown. Genetic dissection of transcript abundance or expression quantitative trait (eQTL) analysis unravels the genomic architecture of regulatory variants. Availability of eQTL information from tissues relevant for glucose homeostasis in humans opens a new avenue to prioritize GWAS-implicated variants that may be involved in triggering a causal chain of events leading to T2D. In this article, we review the progress made in the field of eQTL research and knowledge gained from those studies in understanding transcription regulatory mechanisms in human subjects. We highlight several novel approaches that can integrate eQTL analysis with multiple layers of biological information to identify ethnic-specific causal variants and gene-environment interactions relevant to T2D pathogenesis. Finally, we discuss how the eQTL analysis mediated search for “missing heritability” may lead us to novel biological and molecular mechanisms involved in susceptibility to T2D.

Transcriptome analyses and differential gene expression in a non-model fish species with alternative mating tactics

Background

Social dominance is important for the reproductive success of males in many species. In the black-faced blenny (Tripterygion delaisi) during the reproductive season, some males change color and invest in nest making and defending a territory, whereas others do not change color and ‘sneak’ reproductions when females lay their eggs. Using RNAseq, we profiled differential gene expression between the brains of territorial males, sneaker males, and females to study the molecular signatures of male dimorphism.

Results

We found that more genes were differentially expressed between the two male phenotypes than between males and females, suggesting that during the reproductive period phenotypic plasticity is a more important factor in differential gene expression than sexual dimorphism. The territorial male overexpresses genes related to synaptic plasticity and the sneaker male overexpresses genes involved in differentiation and development.

Conclusions

Previously suggested candidate genes for social dominance in the context of alternative mating strategies seem to be predominantly species-specific. We present a list of novel genes which are differentially expressed in Tripterygion delaisi. This is the first genome-wide study for a molecular non-model species in the context of alternative mating strategies and provides essential information for further studies investigating the molecular basis of social dominance.

Splicing in the human brain

It has become increasingly clear over the past decade that RNA has important functions in human cells beyond its role as an intermediate translator of DNA to protein. It is now known that RNA plays highly specific roles in pathways involved in regulatory, structural, and catalytic functions. The complexity of RNA production and regulation has become evident with the advent of high-throughput methods to study the transcriptome. Deep sequencing has revealed an enormous diversity of RNA types and transcript isoforms in human cells. The transcriptome of the human brain is particularly interesting as it contains more expressed genes than other tissues and also displays an extreme diversity of transcript isoforms, indicating that highly complex regulatory pathways are present in the brain. Several of these regulatory proteins are now identified, including RNA-binding proteins that are neuron specific. RNA-binding proteins also play important roles in regulating the splicing process and the temporal and spatial isoform production. While significant progress has been made in understanding the human transcriptome, many questions still remain regarding the basic mechanisms of splicing and subcellular localization of RNA. A long-standing question is to what extent the splicing of pre-mRNA is cotranscriptional and posttranscriptional, respectively. Recent data, including studies of the human brain, indicate that splicing is primarily cotranscriptional in human cells. This chapter describes the current understanding of splicing and splicing regulation in the human brain and discusses the recent global sequence-based analyses of transcription and splicing.

Transcriptome characterization by RNA sequencing identifies a major molecular and clinical subdivision in chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) has heterogeneous clinical and biological behavior. Whole-genome and -exome sequencing has contributed to the characterization of the mutational spectrum of the disease, but the underlying transcriptional profile is still poorly understood. We have performed deep RNA sequencing in different subpopulations of normal B-lymphocytes and CLL cells from a cohort of 98 patients, and characterized the CLL transcriptional landscape with unprecedented resolution. We detected thousands of transcriptional elements differentially expressed between the CLL and normal B cells, including protein-coding genes, noncoding RNAs, and pseudogenes. Transposable elements are globally derepressed in CLL cells. In addition, two thousand genes—most of which are not differentially expressed—exhibit CLL-specific splicing patterns. Genes involved in metabolic pathways showed higher expression in CLL, while genes related to spliceosome, proteasome, and ribosome were among the most down-regulated in CLL. Clustering of the CLL samples according to RNA-seq derived gene expression levels unveiled two robust molecular subgroups, C1 and C2. C1/C2 subgroups and the mutational status of the immunoglobulin heavy variable (IGHV) region were the only independent variables in predicting time to treatment in a multivariate analysis with main clinico-biological features. This subdivision was validated in an independent cohort of patients monitored through DNA microarrays. Further analysis shows that B-cell receptor (BCR) activation in the microenvironment of the lymph node may be at the origin of the C1/C2 differences.

Functional transcriptomics in the post-ENCODE era

The last decade has seen tremendous effort committed to the annotation of the human genome sequence, most notably perhaps in the form of the ENCODE project. One of the major findings of ENCODE, and other genome analysis projects, is that the human transcriptome is far larger and more complex than previously thought. This complexity manifests, for example, as alternative splicing within protein-coding genes, as well as in the discovery of thousands of long noncoding RNAs. It is also possible that significant numbers of human transcripts have not yet been described by annotation projects, while existing transcript models are frequently incomplete. The question as to what proportion of this complexity is truly functional remains open, however, and this ambiguity presents a serious challenge to genome scientists. In this article, we will discuss the current state of human transcriptome annotation, drawing on our experience gained in generating the GENCODE gene annotation set. We highlight the gaps in our knowledge of transcript functionality that remain, and consider the potential computational and experimental strategies that can be used to help close them. We propose that an understanding of the true overlap between transcriptional complexity and functionality will not be gained in the short term. However, significant steps toward obtaining this knowledge can now be taken by using an integrated strategy, combining all of the experimental resources at our disposal.

Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals

Understanding the consequences of regulatory variation in the human genome remains a major challenge, with important implications for understanding gene regulation and interpreting the many disease-risk variants that fall outside of protein-coding regions. Here, we provide a direct window into the regulatory consequences of genetic variation by sequencing RNA from 922 genotyped individuals. We present a comprehensive description of the distribution of regulatory variation--by the specific expression phenotypes altered, the properties of affected genes, and the genomic characteristics of regulatory variants. We detect variants influencing expression of over ten thousand genes, and through the enhanced resolution offered by RNA-sequencing, for the first time we identify thousands of variants associated with specific phenotypes including splicing and allelic expression. Evaluating the effects of both long-range intra-chromosomal and trans (cross-chromosomal) regulation, we observe modularity in the regulatory network, with three-dimensional chromosomal configuration playing a particular role in regulatory modules within each chromosome. We also observe a significant depletion of regulatory variants affecting central and critical genes, along with a trend of reduced effect sizes as variant frequency increases, providing evidence that purifying selection and buffering have limited the deleterious impact of regulatory variation on the cell. Further, generalizing beyond observed variants, we have analyzed the genomic properties of variants associated with expression and splicing and developed a Bayesian model to predict regulatory consequences of genetic variants, applicable to the interpretation of individual genomes and disease studies. Together, these results represent a critical step toward characterizing the complete landscape of human regulatory variation.

Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals

Understanding the consequences of regulatory variation in the human genome remains a major challenge, with important implications for understanding gene regulation and interpreting the many disease-risk variants that fall outside of protein-coding regions. Here, we provide a direct window into the regulatory consequences of genetic variation by sequencing RNA from 922 genotyped individuals. We present a comprehensive description of the distribution of regulatory variation—by the specific expression phenotypes altered, the properties of affected genes, and the genomic characteristics of regulatory variants. We detect variants influencing expression of over ten thousand genes, and through the enhanced resolution offered by RNA-sequencing, for the first time we identify thousands of variants associated with specific phenotypes including splicing and allelic expression. Evaluating the effects of both long-range intra-chromosomal and trans (cross-chromosomal) regulation, we observe modularity in the regulatory network, with three-dimensional chromosomal configuration playing a particular role in regulatory modules within each chromosome. We also observe a significant depletion of regulatory variants affecting central and critical genes, along with a trend of reduced effect sizes as variant frequency increases, providing evidence that purifying selection and buffering have limited the deleterious impact of regulatory variation on the cell. Further, generalizing beyond observed variants, we have analyzed the genomic properties of variants associated with expression and splicing and developed a Bayesian model to predict regulatory consequences of genetic variants, applicable to the interpretation of individual genomes and disease studies. Together, these results represent a critical step toward characterizing the complete landscape of human regulatory variation.

Transcriptome and genome sequencing uncovers functional variation in humans

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues

Alternative splicing is widely recognized for its roles in regulating genes and creating gene diversity. However, despite many efforts, the repertoire of gene splicing variation is still incompletely characterized, even in humans. Here we describe a new computational system, ASprofile, and its application to RNA-seq data from Illumina's Human Body Map project (>2.5 billion reads).  Using the system, we identified putative alternative splicing events in 16 different human tissues, which provide a dynamic picture of splicing variation across the tissues. We detected 26,989 potential exon skipping events representing differences in splicing patterns among the tissues. A large proportion of the events (>60%) were novel, involving new exons (~3000), new introns (~16000), or both. When tracing these events across the sixteen tissues, only a small number (4-7%) appeared to be differentially expressed ('switched') between two tissues, while 30-45% showed little variation, and the remaining 50-65% were not present in one or both tissues compared.  Novel exon skipping events appeared to be slightly less variable than known events, but were more tissue-specific. Our study represents the first effort to build a comprehensive catalog of alternative splicing in normal human tissues from RNA-seq data, while providing insights into the role of alternative splicing in shaping tissue transcriptome differences. The catalog of events and the ASprofile software are freely available from the Zenodo repository ( http://zenodo.org/record/7068; doi: 10.5281/zenodo.7068) and from our web site http://ccb.jhu.edu/software/ASprofile.

Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues

Alternative splicing is widely recognized for its roles in regulating genes and creating gene diversity. However, despite many efforts, the repertoire of gene splicing variation is still incompletely characterized, even in humans. Here we describe a new computational system, ASprofile, and its application to RNA-seq data from Illumina’s Human Body Map project (>2.5 billion reads).  Using the system, we identified putative alternative splicing events in 16 different human tissues, which provide a dynamic picture of splicing variation across the tissues. We detected 26,989 potential exon skipping events representing differences in splicing patterns among the tissues. A large proportion of the events (>60%) were novel, involving new exons (~3000), new introns (~16000), or both. When tracing these events across the sixteen tissues, only a small number (4-7%) appeared to be differentially expressed (‘switched’) between two tissues, while 30-45% showed little variation, and the remaining 50-65% were not present in one or both tissues compared.  Novel exon skipping events appeared to be slightly less variable than known events, but were more tissue-specific. Our study represents the first effort to build a comprehensive catalog of alternative splicing in normal human tissues from RNA-seq data, while providing insights into the role of alternative splicing in shaping tissue transcriptome differences. The catalog of events and the ASprofile software are freely available from the Zenodo repository

( http://zenodo.org/record/7068; doi: 10.5281/zenodo.7068) and from our web site http://ccb.jhu.edu/software/ASprofile.

Transcriptome analyses of the human retina identify unprecedented transcript diversity and 3.5 Mb of novel transcribed sequence via significant alternative splicing and novel genes

BACKGROUND

The retina is a complex tissue comprised of multiple cell types that is affected by a diverse set of diseases that are important causes of vision loss. Characterizing the transcripts, both annotated and novel, that are expressed in a given tissue has become vital for understanding the mechanisms underlying the pathology of disease.

RESULTS

We sequenced RNA prepared from three normal human retinas and characterized the retinal transcriptome at an unprecedented level due to the increased depth of sampling provided by the RNA-seq approach. We used a non-redundant reference transcriptome from all of the empirically-determined human reference tracks to identify annotated and novel sequences expressed in the retina. We detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3' and 5' alternate splice sites, and 28,271 exon skipping events. We also identified 116 potential novel genes. These data represent a significant addition to the annotated human transcriptome. For example, the novel exons detected increase the number of identified exons by 3%. Using a high-throughput RNA capture approach to validate 14,696 of these novel transcriptome features we found that 99% of the putative novel events can be reproducibly detected. Further, 15-36% of the novel splicing events maintain an open reading frame, suggesting they produce novel protein products.

CONCLUSIONS

To our knowledge, this is the first application of RNA capture to perform large-scale validation of novel transcriptome features. In total, these analyses provide extensive detail about a previously uncharacterized level of transcript diversity in the human retina.

Transcriptome analysis of human tissues and cell lines reveals one dominant transcript per gene

BACKGROUND

RNA sequencing has opened new avenues for the study of transcriptome composition. Significant evidence has accumulated showing that the human transcriptome contains in excess of a hundred thousand different transcripts. However, it is still not clear to what extent this diversity prevails when considering the relative abundances of different transcripts from the same gene.

RESULTS

Here we show that, in a given condition, most protein coding genes have one major transcript expressed at significantly higher level than others, that in human tissues the major transcripts contribute almost 85 percent to the total mRNA from protein coding loci, and that often the same major transcript is expressed in many tissues. We detect a high degree of overlap between the set of major transcripts and a recently published set of alternatively spliced transcripts that are predicted to be translated utilizing proteomic data. Thus, we hypothesize that although some minor transcripts may play a functional role, the major ones are likely to be the main contributors to the proteome. However, we still detect a non-negligible fraction of protein coding genes for which the major transcript does not code a protein.

CONCLUSIONS

Overall, our findings suggest that the transcriptome from protein coding loci is dominated by one transcript per gene and that not all the transcripts that contribute to transcriptome diversity are equally likely to contribute to protein diversity. This observation can help to prioritize candidate targets in proteomics research and to predict the functional impact of the detected changes in variation studies.

TrueSight: a new algorithm for splice junction detection using RNA-seq

RNA-seq has proven to be a powerful technique for transcriptome profiling based on next-generation sequencing (NGS) technologies. However, due to the short length of NGS reads, it is challenging to accurately map RNA-seq reads to splice junctions (SJs), which is a critically important step in the analysis of alternative splicing (AS) and isoform construction. In this article, we describe a new method, called TrueSight, which for the first time combines RNA-seq read mapping quality and coding potential of genomic sequences into a unified model. The model is further utilized in a machine-learning approach to precisely identify SJs. Both simulations and real data evaluations showed that TrueSight achieved higher sensitivity and specificity than other methods. We applied TrueSight to new high coverage honey bee RNA-seq data to discover novel splice forms. We found that 60.3% of honey bee multi-exon genes are alternatively spliced. By utilizing gene models improved by TrueSight, we characterized AS types in honey bee transcriptome. We believe that TrueSight will be highly useful to comprehensively study the biology of alternative splicing.

The evolutionary landscape of alternative splicing in vertebrate species

How species with similar repertoires of protein-coding genes differ so markedly at the phenotypic level is poorly understood. By comparing organ transcriptomes from vertebrate species spanning ~350 million years of evolution, we observed significant differences in alternative splicing complexity between vertebrate lineages, with the highest complexity in primates. Within 6 million years, the splicing profiles of physiologically equivalent organs diverged such that they are more strongly related to the identity of a species than they are to organ type. Most vertebrate species-specific splicing patterns are cis-directed. However, a subset of pronounced splicing changes are predicted to remodel protein interactions involving trans-acting regulators. These events likely further contributed to the diversification of splicing and other transcriptomic changes that underlie phenotypic differences among vertebrate species.

Reuse of public genome-wide gene expression data

Our understanding of gene expression has changed dramatically over the past decade, largely catalysed by technological developments. High-throughput experiments - microarrays and next-generation sequencing - have generated large amounts of genome-wide gene expression data that are collected in public archives. Added-value databases process, analyse and annotate these data further to make them accessible to every biologist. In this Review, we discuss the utility of the gene expression data that are in the public domain and how researchers are making use of these data. Reuse of public data can be very powerful, but there are many obstacles in data preparation and analysis and in the interpretation of the results. We will discuss these challenges and provide recommendations that we believe can improve the utility of such data.

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.