Abstract
Understanding the functional consequences of genetic variation, and how it affects complex human disease and quantitative traits, remains a critical challenge for biomedicine. We present an analysis of RNA sequencing data from 1641 samples across 43 tissues from 175 individuals, generated as part of the pilot phase of the Genotype-Tissue Expression (GTEx) project. We describe the landscape of gene expression across tissues, catalog thousands of tissue-specific and shared regulatory expression quantitative trait loci (eQTL) variants, describe complex network relationships, and identify signals from genome-wide association studies explained by eQTLs. These findings provide a systematic understanding of the cellular and biological consequences of human genetic variation and of the heterogeneity of such effects among a diverse set of human tissues.

Abstract
Chimeric RNAs originating from two or more different genes are known to exist not only in cancer, but also in normal tissues, where they can play a role in human evolution. However, the exact mechanism of their formation is unknown. Here, we use RNA sequencing data from 462 healthy individuals representing 5 human populations to systematically identify and in depth characterize 81 RNA tandem chimeric transcripts, 13 of which are novel. We observe that 6 out of these 81 chimeras have been regarded as cancer-specific. Moreover, we show that a prevalence of long introns at the fusion breakpoint is associated with the chimeric transcripts formation. We also find that tandem RNA chimeras have lower abundances as compared to their partner genes. Finally, by combining our results with genomic data from the same individuals we uncover intronic genetic variants associated with the chimeric RNA formation. Taken together our findings provide an important insight into the chimeric transcripts formation and open new avenues of research into the role of intronic genetic variants in post-transcriptional processing events. © 2014 Greger et al.

Abstract
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â(tm) 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. © 2014 Macmillan Publishers Limited.

Abstract
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. © 2014 Hansen et al.

Abstract
Accurate prediction of the functional effect of genetic variation is critical for clinical genome interpretation.We systematically characterized the transcriptome effects of protein-truncating variants, a class of variants expected to have profound effects on gene function, using data from the Genotype-Tissue Expression (GTEx) and Geuvadis projects. We quantitated tissue-specific and positional effects on nonsense-mediated transcript decay and present an improved predictive model for this decay. We directly measured the effect of variants both proximal and distal to splice junctions. Furthermore, we found that robustness to heterozygous gene inactivation is not due to dosage compensation. Our results illustrate the value of transcriptome data in the functional interpretation of genetic variants.

Abstract
Type 2 diabetes involves defective insulin secretion with islet inflammation governed in part by IL-1ß. Prolonged exposure of islets to high concentrations of IL-1ß (>24 h, 20 ng/ml) impairs beta cell function and survival. Conversely, exposure to lower concentrations of IL-1ß for >24 h improves these same parameters. The impact on insulin secretion of shorter exposure times to IL-1ßand the underlying molecular mechanisms are poorly understood and were the focus of this study. Treatment of rat primary beta cells, as well as rat or human whole islets, with 0.1 ng/ml IL-1ß for 2 h increased glucose-stimulated (but not basal) insulin secretion, whereas 20 ng/ml was without effect. Similar differential effects of IL-1ß depending on concentration were observed after 15 min of KCl stimulation but were prevented by diazoxide. Studies on sorted rat beta cells indicated that the enhancement of stimulated secretion by 0.1 ng/ml IL-1ß was mediated by the NF-ßB pathway and c-JUN/JNK pathway acting in parallel to elicit focal adhesion remodeling and the phosphorylation of paxillin independently of upstream regulation by focal adhesion kinase. Because the beneficial effect of IL-1ß was dependent in part upon transcription, gene expression was analyzed by RNAseq. There were 18 genes regulated uniquely by 0.1 but not 20 ng/ml IL-1ß, which are mostly involved in transcription and apoptosis. These results indicate that 2h of exposure of beta cells to a low but not a high concentration of IL-1ß enhances glucose-stimulated insulin secretion through focal adhesion and actin remodeling, as well as modulation of gene expression. © 2015 by The American Society for Biochemistry and Molecular Biology, Inc.