Discovery of tissue-specific exons using comprehensive human exon microarrays

Melanoma spheroids grown under neural crest cell conditions are highly plastic migratory/invasive tumor cells endowed with immunomodulator function

Background

The aggressiveness of melanoma tumors is likely to rely on their well-recognized heterogeneity and plasticity. Melanoma comprises multi-subpopulations of cancer cells some of which may possess stem cell-like properties. Although useful, the sphere-formation assay to identify stem cell-like or tumor initiating cell subpopulations in melanoma has been challenged, and it is unclear if this model can predict a functional phenotype associated with aggressive tumor cells.

Methodology/Principal Findings

We analyzed the molecular and functional phenotypes of melanoma spheroids formed in neural crest cell medium. Whether from metastatic or advanced primary tumors, spheroid cells expressed melanoma-associated markers. They displayed higher capacity to differentiate along mesenchymal lineages and enhanced expression of SOX2, NANOG, KLF4, and/or OCT4 transcription factors, but not enhanced self-renewal or tumorigenicity when compared to their adherent counterparts. Gene expression profiling attributed a neural crest cell signature to these spheroids and indicated that a migratory/invasive and immune-function modulating program could be associated with these cells. In vitro assays confirmed that spheroids display enhanced migratory/invasive capacities. In immune activation assays, spheroid cells elicited a poorer allogenic response from immune cells and inhibited mitogen-dependent T cells activation and proliferation more efficiently than their adherent counterparts. Our findings reveal a novel immune-modulator function of melanoma spheroids and suggest specific roles for spheroids in invasion and in evasion of antitumor immunity.

Conclusion/Significance

The association of a more plastic, invasive and evasive, thus a more aggressive tumor phenotype with melanoma spheroids reveals a previously unrecognized aspect of tumor cells expanded as spheroid cultures. While of limited efficiency for melanoma initiating cell identification, our melanoma spheroid model predicted aggressive phenotype and suggested that aggressiveness and heterogeneity of melanoma tumors can be supported by subpopulations other than cancer stem cells. Therefore, it could be constructive to investigate melanoma aggressiveness, relevant to patients and clinical transferability.

Identification of Tat-SF1 cellular targets by exon array analysis reveals dual roles in transcription and splicing

Tat specific factor 1 (Tat-SF1) interacts with components of both the transcription and splicing machineries and has been classified as a transcription-splicing factor. Although its function as an HIV-1 dependency factor has been investigated, relatively little is known about the cellular functions of Tat-SF1. To identify target genes of Tat-SF1, we utilized a combination of RNAi and exon-specific microarrays. These arrays, which survey genome-wide changes in transcript and individual exon levels, revealed 450 genes with transcript level changes upon Tat-SF1 depletion. Strikingly, 98% of these target genes were down-regulated upon depletion, indicating that Tat-SF1 generally activates gene expression. We also identified 89 genes that showed differential exon level changes after Tat-SF1 depletion. The 89 genes showed evidence of many different types of alternative exon use consistent with the regulation of transcription initiation sites and RNA processing. Minimal overlap between genes with transcript-level and exon-level changes suggests that Tat-SF1 does not functionally couple transcription and splicing. Biological processes significantly enriched with transcript- and exon-level targets include the cell cycle and nucleic acid metabolism; the insulin signaling pathway was enriched with Tat-SF1 transcript-level targets but not exon-level targets. Additionally, a hexamer, ATGCCG, was over-represented in the promoter region of genes showing changes in transcription initiation upon Tat-SF1 depletion. This may represent a novel motif that Tat-SF1 recognizes during transcription. Together, these findings suggest that Tat-SF1 functions independently in transcription and splicing of cellular genes.

Genome-wide analysis of alternative splicing in medulloblastoma identifies splicing patterns characteristic of normal cerebellar development

Alternative splicing is an important mechanism for the generation of protein diversity at a post-transcriptional level. Modifications in the splicing patterns of several genes have been shown to contribute to the malignant transformation of different tissue types. In this study, we used the Affymetrix Exon arrays to investigate patterns of differential splicing between pediatric medulloblastomas and normal cerebellum on a genome-wide scale. Of the 1,262 genes identified as potentially generating tumor-associated splice forms, we selected 14 examples of differential splicing of known cassette exons and successfully validated 11 of them by reverse transcriptase PCR. The pattern of differential splicing of three validated events was characteristic for the molecular subset of sonic hedgehog (Shh)-driven medulloblastomas, suggesting that their unique gene signature includes the expression of distinctive transcript variants. Generally, we observed that tumor and normal fetal cerebellar samples shared significantly lower exon inclusion rates than normal adult cerebellum. We investigated whether tumor-associated splice forms were expressed in primary cultures of Shh-dependent mouse cerebellar granule cell precursors (GCP) and found that Shh caused a decrease in the cassette exon inclusion rate of five of the seven tested genes. Furthermore, we observed a significant increase in exon inclusion between postnatal days 7 and 14 of mouse cerebellar development, at the time when GCPs mature into postmitotic neurons. We conclude that inappropriate splicing frequently occurs in human medulloblastomas and may be linked to the activation of developmental signaling pathways and a failure of cerebellar precursor cells to differentiate.

Alternative splicing in bone following mechanical loading

It is estimated that more than 90% of human genes express multiple mRNA transcripts due to alternative splicing. Consequently, the proteins produced by different splice variants will likely have different functions and expression levels. Several genes with splice variants are known in bone, with functions that affect osteoblast function and bone formation. The primary goal of this study was to evaluate the extent of alternative splicing in a bone subjected to mechanical loading and subsequent bone formation. We used the rat forelimb loading model, in which the right forelimb was loaded axially for 3 min, while the left forearm served as a non-loaded control. Animals were subjected to loading sessions every day, with 24 h between sessions. Ulnae were sampled at 11 time points, from 4 h to 32days after beginning loading. RNA was isolated and mRNA abundance was measured at each time point using Affymetrix exon arrays (GeneChip® Rat Exon 1.0 ST Arrays). An ANOVA model was used to identify potential alternatively spliced genes across the time course, and five alternatively spliced genes were validated with qPCR: Akap12, Fn1, Pcolce, Sfrp4, and Tpm1. The number of alternatively spliced genes varied with time, ranging from a low of 68 at 12h to a high of 992 at 16d. We identified genes across the time course that encoded proteins with known functions in bone formation, including collagens, matrix proteins, and components of the Wnt/β-catenin and TGF-β signaling pathways. We also identified alternatively spliced genes encoding cytokines, ion channels, muscle-related genes, and solute carriers that do not have a known function in bone formation and represent potentially novel findings. In addition, a functional characterization was performed to categorize the global functions of the alternatively spliced genes in our data set. In conclusion, mechanical loading induces alternative splicing in bone, which may play an important role in the response of bone to mechanical loading.

Comprehensive analysis of alternative splicing and functionality in neuronal differentiation of P19 cells

BACKGROUND

Alternative splicing, which produces multiple mRNAs from a single gene, occurs in most human genes and contributes to protein diversity. Many alternative isoforms are expressed in a spatio-temporal manner, and function in diverse processes, including in the neural system.

METHODOLOGY/PRINCIPAL FINDINGS

The purpose of the present study was to comprehensively investigate neural-splicing using P19 cells. GeneChip Exon Array analysis was performed using total RNAs purified from cells during neuronal cell differentiation. To efficiently and readily extract the alternative exon candidates, 9 filtering conditions were prepared, yielding 262 candidate exons (236 genes). Semiquantitative RT-PCR results in 30 randomly selected candidates suggested that 87% of the candidates were differentially alternatively spliced in neuronal cells compared to undifferentiated cells. Gene ontology and pathway analyses suggested that many of the candidate genes were associated with neural events. Together with 66 genes whose functions in neural cells or organs were reported previously, 47 candidate genes were found to be linked to 189 events in the gene-level profile of neural differentiation. By text-mining for the alternative isoform, distinct functions of the isoforms of 9 candidate genes indicated by the result of Exon Array were confirmed.

CONCLUSIONS/SIGNIFICANCE

Alternative exons were successfully extracted. Results from the informatics analyses suggested that neural events were primarily governed by genes whose expression was increased and whose transcripts were differentially alternatively spliced in the neuronal cells. In addition to known functions in neural cells or organs, the uninvestigated alternative splicing events of 11 genes among 47 candidate genes suggested that cell cycle events are also potentially important. These genes may help researchers to differentiate the roles of alternative splicing in cell differentiation and cell proliferation.

Human transcriptome array for high-throughput clinical studies

A 6.9 million-feature oligonucleotide array of the human transcriptome [Glue Grant human transcriptome (GG-H array)] has been developed for high-throughput and cost-effective analyses in clinical studies. This array allows comprehensive examination of gene expression and genome-wide identification of alternative splicing as well as detection of coding SNPs and noncoding transcripts. The performance of the array was examined and compared with mRNA sequencing (RNA-Seq) results over multiple independent replicates of liver and muscle samples. Compared with RNA-Seq of 46 million uniquely mappable reads per replicate, the GG-H array is highly reproducible in estimating gene and exon abundance. Although both platforms detect similar expression changes at the gene level, the GG-H array is more sensitive at the exon level. Deeper sequencing is required to adequately cover low-abundance transcripts. The array has been implemented in a multicenter clinical program and has generated high-quality, reproducible data. Considering the clinical trial requirements of cost, sample availability, and throughput, the GG-H array has a wide range of applications. An emerging approach for large-scale clinical genomic studies is to first use RNA-Seq to the sufficient depth for the discovery of transcriptome elements relevant to the disease process followed by high-throughput and reliable screening of these elements on thousands of patient samples using custom-designed arrays.

Jumonji domain-containing protein 6 (Jmjd6) is required for angiogenic sprouting and regulates splicing of VEGF-receptor 1

JmjC domain-containing proteins play a crucial role in the control of gene expression by acting as protein hydroxylases or demethylases, thereby controlling histone methylation or splicing. Here, we demonstrate that silencing of Jumonji domain-containing protein 6 (Jmjd6) impairs angiogenic functions of endothelial cells by changing the gene expression and modulating the splicing of the VEGF-receptor 1 (Flt1). Reduction of Jmjd6 expression altered splicing of Flt1 and increased the levels of the soluble form of Flt1, which binds to VEGF and placental growth factor (PlGF) and thereby inhibits angiogenesis. Saturating VEGF or PlGF or neutralizing antibodies directed against soluble Flt1 rescued the angiogenic defects induced by Jmjd6 silencing. Jmjd6 interacts with the splicing factors U2AF65 that binds to Flt1 mRNA. In conclusion, Jmjd6 regulates the splicing of Flt1, thereby controlling angiogenic sprouting.

Differential gene expression in liver and small intestine from lactating rats compared to age-matched virgin controls detects increased mRNA of cholesterol biosynthetic genes

Background

Lactation increases energy demands four- to five-fold, leading to a two- to three-fold increase in food consumption, requiring a proportional adjustment in the ability of the lactating dam to absorb nutrients and to synthesize critical biomolecules, such as cholesterol, to meet the dietary needs of both the offspring and the dam. The size and hydrophobicity of the bile acid pool increases during lactation, implying an increased absorption and disposition of lipids, sterols, nutrients, and xenobiotics. In order to investigate changes at the transcriptomics level, we utilized an exon array and calculated expression levels to investigate changes in gene expression in the liver, duodenum, jejunum, and ileum of lactating dams when compared against age-matched virgin controls.

Results

A two-way mixed models ANOVA was applied to detect differentially expressed genes. Significance calls were defined as a p < 0.05 for the overall physiologic state effect (lactation vs. control), and a within tissue pairwise comparison of p < 0.01. The proportion of false positives, an estimate of the ratio of false positives in the list of differentially expressed genes, was calculated for each tissue. The number of differentially expressed genes was 420 in the liver, 337 in the duodenum, 402 in the jejunum, and 523 in the ileum. The list of differentially expressed genes was in turn analyzed by Ingenuity Pathways Analysis (IPA) to detect biological pathways that were overrepresented. In all tissues, sterol regulatory element binding protein (Srebp)-regulated genes involved in cholesterol synthesis showed increased mRNA expression, with the fewest changes detected in the jejunum. We detected increased Scap mRNA in the liver only, suggesting an explanation for the difference in response to lactation between the liver and small intestine. Expression of Cyp7a1, which catalyzes the rate limiting step in the bile acid biosynthetic pathway, was also significantly increased in liver. In addition, decreased levels of mRNA associated with T-cell signaling were found in the jejunum and ileum. Several members of the Solute Carrier (SLC) and Adenosine Triphosphate Binding Cassette (ABC) superfamilies of membrane transporters were found to be differentially expressed; these genes may play a role in differences in nutrient and xenobiotic absorption and disposition. mRNA expression of SLC39a4_predicted, a zinc transporter, was increased in all tissues, suggesting that it is involved in increased zinc uptake during lactation. Microarray data are available through GEO under GSE19175.

Conclusions

We detected differential expression of mRNA from several pathways in lactating dams, including upregulation of the cholesterol biosynthetic pathway in liver and intestine, consistent with Srebp activation. Differential T-Cell signaling in the two most distal regions of the small intestine (ileum and jejunum) was also noted, as well as differential expression of transporters that likely play a key role in nutrient uptake.

Systems analysis of alternative splicing and its regulation

Alternative splicing (AS) has emerged as a key mechanism that accounts for gene expression diversity in metazoan organisms. Splicing is tightly regulated by a repertoire of RNA and protein factors and RNA sequence elements that function in a cooperative manner. Systems-level experimental and computational approaches have been instrumental in establishing comprehensive profiles of transcript variants generated by AS. In addition, systems biology approaches are starting to define how combinatorial splicing regulation shapes the complex splicing phenotypes observed in different tissue types and developmental stages and under different conditions. Here, we review recent progress in these areas.

A molecular inversion probe assay for detecting alternative splicing

BACKGROUND

A sensitive, high-throughput method for monitoring pre-mRNA splicing on a genomic scale is needed to understand the spectrum of alternatively spliced mRNA in human cells.

RESULTS

We adapted Molecular Inversion Probes (MIPs), a padlock-probe based technology, for the multiplexed capture and quantitation of individual splice events in human tissues. Individual MIP capture probes can be quantified using either DNA microarrays or high-throughput sequencing, which permits independent assessment of each spliced junction. Using our methodology we successfully identified 100% of our positive controls and showed that there is a strong correlation between the data from our alternative splicing MIP (asMIP) assay and quantitative PCR.

CONCLUSION

The asMIP assay provides a sensitive, accurate and multiplexed means for measuring pre-mRNA splicing. Fully optimized, we estimate that the assay could accommodate a throughput of greater than 20,000 splice junctions in a single reaction. This would represent a significant improvement over existing technologies.

Exon array analysis using re-defined probe sets results in reliable identification of alternatively spliced genes in non-small cell lung cancer

Background

Treatment of non-small cell lung cancer with novel targeted therapies is a major unmet clinical need. Alternative splicing is a mechanism which generates diverse protein products and is of functional relevance in cancer.

Results

In this study, a genome-wide analysis of the alteration of splicing patterns between lung cancer and normal lung tissue was performed. We generated an exon array data set derived from matched pairs of lung cancer and normal lung tissue including both the adenocarcinoma and the squamous cell carcinoma subtypes. An enhanced workflow was developed to reliably detect differential splicing in an exon array data set. In total, 330 genes were found to be differentially spliced in non-small cell lung cancer compared to normal lung tissue. Microarray findings were validated with independent laboratory methods for CLSTN1, FN1, KIAA1217, MYO18A, NCOR2, NUMB, SLK, SYNE2, TPM1, (in total, 10 events) and ADD3, which was analysed in depth. We achieved a high validation rate of 69%. Evidence was found that the activity of FOX2, the splicing factor shown to cause cancer-specific splicing patterns in breast and ovarian cancer, is not altered at the transcript level in several cancer types including lung cancer.

Conclusions

This study demonstrates how alternatively spliced genes can reliably be identified in a cancer data set. Our findings underline that key processes of cancer progression in NSCLC are affected by alternative splicing, which can be exploited in the search for novel targeted therapies.

Improvements to previous algorithms to predict gene structure and isoform concentrations using Affymetrix Exon arrays

BACKGROUND

Exon arrays provide a way to measure the expression of different isoforms of genes in an organism. Most of the procedures to deal with these arrays are focused on gene expression or on exon expression. Although the only biological analytes that can be properly assigned a concentration are transcripts, there are very few algorithms that focus on them. The reason is that previously developed summarization methods do not work well if applied to transcripts. In addition, gene structure prediction, i.e., the correspondence between probes and novel isoforms, is a field which is still unexplored.

RESULTS

We have modified and adapted a previous algorithm to take advantage of the special characteristics of the Affymetrix exon arrays. The structure and concentration of transcripts -some of them possibly unknown- in microarray experiments were predicted using this algorithm. Simulations showed that the suggested modifications improved both specificity (SP) and sensitivity (ST) of the predictions. The algorithm was also applied to different real datasets showing its effectiveness and the concordance with PCR validated results.

CONCLUSIONS

The proposed algorithm shows a substantial improvement in the performance over the previous version. This improvement is mainly due to the exploitation of the redundancy of the Affymetrix exon arrays. An R-Package of SPACE with the updated algorithms have been developed and is freely available.

Alternative splicing-related factor YT521: an independent prognostic factor in endometrial cancer

BACKGROUND

YT521 is a splicing factor involved in alternative splicing regulation of several tumor biological important genes. Two messenger RNA (mRNA) isoforms due to YT521 exon6 alternative splicing exist, with so far unknown functional consequences. Further evidence exists for a direct influence of YT521 expression in tumorigenesis because its mRNA level is changed in tumors compared with physiological tissue. We investigated the potential impact of YT521 expression on tumor biological parameters in endometrial cancer (EC).

METHODS

Real-time reverse transcription-polymerase chain reaction specifically detecting YT521 exon6-retention and exon6-skipping mRNA isoforms and immunohistochemistry were performed in a cohort of 130 EC tissue samples.

RESULTS

Whereas YT521 exon6-retention mRNA was detectable in 86 (66.2%), the exon6-skipping isoform mRNA was expressed in only 8 (6.2%) of all EC samples. On the protein level, 104 (80%) of EC samples showed nuclear expression. The mRNA levels of exon6-skipping isoform were not correlated to any of the clinicopathological parameters of EC. In contrast, YT521 exon6-retention mRNA expression was positively correlated to metastasis (R = 0.196, P = 0.026) and inversely correlated to the protein expression levels (R = -0.205, P = 0.019). In univariate analyses, higher levels of YT521 exon6-retention mRNA were correlated to a poorer progression-free survival (P = 0.003), and this is confirmed by multivariate analyses (P = 0.019). The negative YT521 protein expression was correlated to poorer overall and disease-specific survival (P = 0.036 and P = 0.034), respectively, in univariate analyses. They are also confirmed by multivariate analyses (P = 0.021 and P = 0.010, respectively).

CONCLUSIONS

We characterized for the first time in a clinical setting a new but rare exon6-skipping mRNA splicing isoform of YT521. Furthermore, we identified YT521 as a potential new independent prognostic factor for patients with EC: the lack of YT521 protein in tumor cells was highly predictive for a poor overall and disease-specific survival and independent from the histological subtypes.

Detection of clinically relevant exonic copy-number changes by array CGH

Array comparative genomic hybridization (aCGH) is a powerful tool for the molecular elucidation and diagnosis of disorders resulting from genomic copy-number variation (CNV). However, intragenic deletions or duplications--those including genomic intervals of a size smaller than a gene--have remained beyond the detection limit of most clinical aCGH analyses. Increasing array probe number improves genomic resolution, although higher cost may limit implementation, and enhanced detection of benign CNV can confound clinical interpretation. We designed an array with exonic coverage of selected disease and candidate genes and used it clinically to identify losses or gains throughout the genome involving at least one exon and as small as several hundred base pairs in size. In some patients, the detected copy-number change occurs within a gene known to be causative of the observed clinical phenotype, demonstrating the ability of this array to detect clinically relevant CNVs with subkilobase resolution. In summary, we demonstrate the utility of a custom-designed, exon-targeted oligonucleotide array to detect intragenic copy-number changes in patients with various clinical phenotypes.

Cotranscriptional exon skipping in the genotoxic stress response

Pre-mRNA splicing is functionally coupled to transcription, and genotoxic stresses can enhance alternative exon inclusion by affecting elongating RNA polymerase II. We report here that various genotoxic stress inducers, including camptothecin (CPT), inhibit the interaction between Ewing's sarcoma proto-oncoprotein (EWS), an RNA polymerase II-associated factor, and YB-1, a spliceosome-associated factor. This results in the cotranscriptional skipping of several exons of the MDM2 gene, which encodes the main p53 ubiquitin ligase. This reversible exon skipping participates in the regulation of MDM2 expression that may contribute to the accumulation of p53 during stress exposure and its rapid shut-off when stress is removed. Finally, a splicing-sensitive microarray identified numerous exons that are skipped in response to CPT and EWS-YB-1 depletion. These data demonstrate genotoxic stress-induced alteration of the communication between the transcriptional and splicing machineries, which results in widespread exon skipping and plays a central role in the genotoxic stress response.

Role of SFRS13A in low-density lipoprotein receptor splicing

Low-density lipoprotein receptor (LDLR) is a major apolipoprotein E (APOE) receptor and thereby is critical to cholesterol homeostasis and, possibly, Alzheimer disease (AD) development. We previously identified a single nucleotide polymorphism (SNP), rs688:C>T, that modulates LDLR exon 12 splicing and is associated with cholesterol levels in premenopausal women and with Alzheimer disease in men. To gain additional insights into LDLR splicing regulation, we seek to identify splicing factors that modulate LDLR splicing efficiency. By using an in vitro minigene study, we first found that ectopic expression of SFRS3 (SRp20), SFRS13A (SRp38), SFRS13A-2 (SRp38-2), and RBMX (hnRNP G) robustly decreased LDLR splicing efficiency. Although SFRS3 and SFRS13A specifically increased the LDLR transcript lacking exon 11, SFRS13A-2 and RBMX primarily increased the LDLR isoform lacking both exons 11 and 12. When we evaluated the relationship between the expression of these splicing factors and LDLR splicing in human brain and liver specimens, we found that overall SFRS13A expression was significantly associated with LDLR splicing efficiency in vivo. We interpret these results as suggesting that SFRS13A regulates LDLR splicing efficiency and may therefore emerge as a modulator of cholesterol homeostasis.

Implementation of exon arrays: alternative splicing during T-cell proliferation as determined by whole genome analysis

Background

The contribution of alternative splicing and isoform expression to cellular response is emerging as an area of considerable interest, and the newly developed exon arrays allow for systematic study of these processes. We use this pilot study to report on the feasibility of exon array implementation looking to replace the 3' in vitro transcription expression arrays in our laboratory.

One of the most widely studied models of cellular response is T-cell activation from exogenous stimulation. Microarray studies have contributed to our understanding of key pathways activated during T-cell stimulation. We use this system to examine whole genome transcription and alternate exon usage events that are regulated during lymphocyte proliferation in an attempt to evaluate the exon arrays.

Results

Peripheral blood mononuclear cells form healthy donors were activated using phytohemagglutinin, IL2 and ionomycin and harvested at 5 points over a 7 day period. Flow cytometry measured cell cycle events and the Affymetrix exon array platform was used to identify the gene expression and alternate exon usage changes. Gene expression changes were noted in a total of 2105 transcripts, and alternate exon usage identified in 472 transcript clusters. There was an overlap of 263 transcripts which showed both differential expression and alternate exon usage over time. Gene ontology enrichment analysis showed a broader range of biological changes in biological processes for the differentially expressed genes, which include cell cycle, cell division, cell proliferation, chromosome segregation, cell death, component organization and biogenesis and metabolic process ontologies. The alternate exon usage ontological enrichments are in metabolism and component organization and biogenesis. We focus on alternate exon usage changes in the transcripts of the spliceosome complex. The real-time PCR validation rates were 86% for transcript expression and 71% for alternate exon usage.

Conclusions

This study illustrates that the Exon array technology has the potential to provide information on both transcript expression and isoform usage, with very little increase in expense.

Position-dependent alternative splicing activity revealed by global profiling of alternative splicing events regulated by PTB

To gain global insights into the role of the well-known repressive splicing regulator PTB, we analyzed the consequences of PTB knockdown in HeLa cells using high-density oligonucleotide splice-sensitive microarrays. The major class of identified PTB-regulated splicing event was PTB-repressed cassette exons, but there was also a substantial number of PTB-activated splicing events. PTB-repressed and PTB-activated exons showed a distinct arrangement of motifs with pyrimidine-rich motif enrichment within and upstream of repressed exons but downstream of activated exons. The N-terminal half of PTB was sufficient to activate splicing when recruited downstream of a PTB-activated exon. Moreover, insertion of an upstream pyrimidine tract was sufficient to convert a PTB-activated exon to a PTB-repressed exon. Our results show that PTB, an archetypal splicing repressor, has variable splicing activity that predictably depends upon its binding location with respect to target exons.

The CASPR2 cell adhesion molecule functions as a tumor suppressor gene in glioma

Genomic translocations have been implicated in cancer. In this study, we performed a screen for genetic translocations in gliomas based on exon-level expression profiles. We identified a translocation in the contactin-associated protein-like 2 (CASPR2) gene, encoding a cell adhesion molecule. CASPR2 mRNA was fused to an expressed sequence tag that likely is part of the nuclear receptor coactivator 1 gene. Despite high mRNA expression levels, no CASPR2 fusion protein was detected. In a set of 25 glioblastomas and 22 oligodendrogliomas, mutation analysis identified two additional samples with genetic alterations in the CASPR2 gene and all three identified genetic alterations are likely to reduce CASPR2 protein expression levels. Methylation of the CASPR2 gene was also observed in gliomas and glioma cell lines. CASPR2-overexpressing cells showed decreased proliferation rates, likely because of an increase in apoptosis. Moreover, high CASPR2 mRNA expression level is positively correlated with survival and is an independent prognostic factor. These results indicate that CASPR2 acts as a tumor suppressor gene in glioma.

Bioinformatics of the urinary proteome

Proteomics-based biomarker discovery studies usually entail the isolation of peptide fragments from candidate biomarkers of interest. Detection of such peptides from biological or clinical samples and identification of the corresponding full-length protein and the gene encoding that protein provide the means to gather a wealth of information. This information, termed annotation because it is attached to the gene or protein sequence under study, describes relationships to human disease, cytogenetic map position, protein domains, protein-protein and small molecule interactions, tissues or cell types in which the gene is expressed, as well as several other aspects of gene and protein function. Bioinformatics tools are employed and genome databases are mined to retrieve this information. Coupled with extensive gene and protein annotation, detected peptides are better placed in a biological context with respect to the health status of the subject. Examples of the status include cancers (bladder, kidney), metabolic disorders (diabetes and kidney function), and the nutritional state of the subject.

A framework for oligonucleotide microarray preprocessing

MOTIVATION

The availability of flexible open source software for the analysis of gene expression raw level data has greatly facilitated the development of widely used preprocessing methods for these technologies. However, the expansion of microarray applications has exposed the limitation of existing tools.

RESULTS

We developed the oligo package to provide a more general solution that supports a wide range of applications. The package is based on the BioConductor principles of transparency, reproducibility and efficiency of development. It extends the existing tools and leverages existing code for visualization, accessing data and widely used preprocessing routines. The oligo package implements a unified paradigm for preprocessing data and interfaces with other BioConductor tools for downstream analysis. Our infrastructure is general and can be used by other BioConductor packages.

AVAILABILITY

The oligo package is freely available through BioConductor, http://www.bioconductor.org.

Heritability in the efficiency of nonsense-mediated mRNA decay in humans

Background

In eukaryotes mRNA transcripts of protein-coding genes in which an intron has been retained in the coding region normally result in premature stop codons and are therefore degraded through the nonsense-mediated mRNA decay (NMD) pathway. There is evidence in the form of selective pressure for in-frame stop codons in introns and a depletion of length three introns that this is an important and conserved quality-control mechanism. Yet recent reports have revealed that the efficiency of NMD varies across tissues and between individuals, with important clinical consequences.

Principal Findings

Using previously published Affymetrix exon microarray data from cell lines genotyped as part of the International HapMap project, we investigated whether there are heritable, inter-individual differences in the abundance of intron-containing transcripts, potentially reflecting differences in the efficiency of NMD. We identified intronic probesets using EST data and report evidence of heritability in the extent of intron expression in 56 HapMap trios. We also used a genome-wide association approach to identify genetic markers associated with intron expression. Among the top candidates was a SNP in the DCP1A gene, which forms part of the decapping complex, involved in NMD.

Conclusions

While we caution that some of the apparent inter-individual difference in intron expression may be attributable to different handling or treatments of cell lines, we hypothesize that there is significant polymorphism in the process of NMD, resulting in heritable differences in the abundance of intronic mRNA. Part of this phenotype is likely to be due to a polymorphism in a decapping enzyme on human chromosome 3.

Systems biology-based approaches to understand HIV-exposed uninfected women

Worldwide HIV infects women more frequently than men, and it is clear that not all exposed to HIV become infected. Several populations of HIV-exposed uninfected (EU) women have been identified, including discordant couples and sex workers. Understanding what provides natural protection in EU women is critical in vaccine or microbicide development. However, correlates of protection in these women are still unclear. Most studies have used classical methods, examining single genes or cellular factors, a mainstay for traditional immunobiology. This reductionist approach may be limited in the information it can provide. Novel technologies are now available that allow us to take a "systems biology" approach, which allows the study of a complex biological system and identifies factors that may provide protection against HIV infection. Herein we report developments in discovery-based systems biology approaches in EU women and how this broadens our understanding of natural protection against HIV-1.

Exon-level microarray analyses identify alternative splicing programs in breast cancer

Protein isoforms produced by alternative splicing (AS) of many genes have been implicated in several aspects of cancer genesis and progression. These observations motivated a genome-wide assessment of AS in breast cancer. We accomplished this by measuring exon level expression in 31 breast cancer and nonmalignant immortalized cell lines representing luminal, basal, and claudin-low breast cancer subtypes using Affymetrix Human Junction Arrays. We analyzed these data using a computational pipeline specifically designed to detect AS with a low false-positive rate. This identified 181 splice events representing 156 genes as candidates for AS. Reverse transcription-PCR validation of a subset of predicted AS events confirmed 90%. Approximately half of the AS events were associated with basal, luminal, or claudin-low breast cancer subtypes. Exons involved in claudin-low subtype-specific AS were significantly associated with the presence of evolutionarily conserved binding motifs for the tissue-specific Fox2 splicing factor. Small interfering RNA knockdown of Fox2 confirmed the involvement of this splicing factor in subtype-specific AS. The subtype-specific AS detected in this study likely reflects the splicing pattern in the breast cancer progenitor cells in which the tumor arose and suggests the utility of assays for Fox-mediated AS in cancer subtype definition and early detection. These data also suggest the possibility of reducing the toxicity of protein-targeted breast cancer treatments by targeting protein isoforms that are not present in limiting normal tissues.

Single molecule analysis of combinatorial splicing

Alternative splicing creates diverse mRNA isoforms from single genes and thereby enhances complexity of transcript structure and of gene function. We describe a method called spliceotyping, which translates combinatorial mRNA splicing patterns along transcripts into a library of binary strings of nucleic acid tags that encode the exon composition of individual mRNA molecules. The exon inclusion pattern of each analyzed transcript is thus represented as binary data, and the abundance of different splice variants is registered by counts of individual molecules. The technique is illustrated in a model experiment by analyzing the splicing patterns of the adenovirus early 1A gene and the beta actin reference transcript. The method permits many genes to be analyzed in parallel and it will be valuable for elucidating the complex effects of combinatorial splicing.

Genetic variation in microRNA networks: the implications for cancer research

Many studies have highlighted the role that microRNAs have in physiological processes and how their deregulation can lead to cancer. More recently, it has been proposed that the presence of single nucleotide polymorphisms in microRNA genes, their processing machinery and target binding sites affects cancer risk, treatment efficacy and patient prognosis. In reviewing this new field of cancer biology, we describe the methodological approaches of these studies and make recommendations for which strategies will be most informative in the future.

AltAnalyze and DomainGraph: analyzing and visualizing exon expression data

Alternative splicing is an important mechanism for increasing protein diversity. However, its functional effects are largely unknown. Here, we present our new software workflow composed of the open-source application AltAnalyze and the Cytoscape plugin DomainGraph. Both programs provide an intuitive and comprehensive end-to-end solution for the analysis and visualization of alternative splicing data from Affymetrix Exon and Gene Arrays at the level of proteins, domains, microRNA binding sites, molecular interactions and pathways. Our software tools include easy-to-use graphical user interfaces, rigorous statistical methods (FIRMA, MiDAS and DABG filtering) and do not require prior knowledge of exon array analysis or programming. They provide new methods for automatic interpretation and visualization of the effects of alternative exon inclusion on protein domain composition and microRNA binding sites. These data can be visualized together with affected pathways and gene or protein interaction networks, allowing a straightforward identification of potential biological effects due to alternative splicing at different levels of granularity. Our programs are available at http://www.altanalyze.org and http://www.domaingraph.de. These websites also include extensive documentation, tutorials and sample data.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing

Background

There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important.

Results

Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage.

Conclusions

We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Evolution of alternative splicing in primate brain transcriptomes

Alternative splicing is a predominant form of gene regulation in higher eukaryotes. The evolution of alternative splicing provides an important mechanism for the acquisition of novel gene functions. In this work, we carried out a genome-wide phylogenetic survey of lineage-specific splicing patterns in the primate brain, via high-density exon junction array profiling of brain transcriptomes of humans, chimpanzees and rhesus macaques. We identified 509 genes showing splicing differences among these species. RT-PCR analysis of 40 exons confirmed the predicted splicing evolution of 33 exons. Of these 33 exons, outgroup analysis using rhesus macaques confirmed 13 exons with human-specific increase or decrease in transcript inclusion levels after humans diverged from chimpanzees. Some of the human-specific brain splicing patterns disrupt domains critical for protein-protein interactions, and some modulate translational efficiency of their host genes. Strikingly, for exons showing splicing differences across species, we observed a significant increase in the rate of silent substitutions within exons, coupled with accelerated sequence divergence in flanking introns. This indicates that evolution of cis-regulatory signals is a major contributor to the emergence of human-specific splicing patterns. In one gene (MAGOH), using minigene reporter assays, we demonstrated that the combination of two human-specific cis-sequence changes created its human-specific splicing pattern. Together, our data reveal widespread human-specific changes of alternative splicing in the brain and suggest an important role of splicing in the evolution of neuronal gene regulation and functions.

Genome-wide analysis of alternative splicing evolution among Mus subspecies

Alternative splicing, the combination of different exons to produce a variety of transcripts from a single gene, contributes enormously to transcriptome diversity in mammals, and the majority of genes encode alternatively spliced products. Previous research comparing mouse, rat and human has shown that a significant proportion of splice forms are not conserved across species, suggesting that alternative transcripts are an important source of evolutionary novelty. Here, we studied the evolution of alternative splicing in the early stages of species divergence in the house mouse. We sequenced the testis transcriptomes of three Mus musculus subspecies and Mus spretus using Illumina technology. On the basis of a genome-wide analysis of read coverage differences among subspecies, we identified several hundred candidate alternatively spliced regions. We conservatively estimate that 6.5% of testis-expressed genes show alternative splice differences between at least one pair of M. musculus subspecies, a proportion slightly higher than the proportion of genes differentially expressed among subspecies. These results suggest that differences in both the structure and abundance of transcripts contribute to early transcriptome divergence.

Different effects of the probe summarization algorithms PLIER and RMA on high-level analysis of Affymetrix exon arrays

BACKGROUND

Alternative splicing is an important mechanism that increases protein diversity and functionality in higher eukaryotes. Affymetrix exon arrays are a commercialized platform used to detect alternative splicing on a genome-wide scale. Two probe summarization algorithms, PLIER (Probe Logarithmic Intensity Error) and RMA (Robust Multichip Average), are commonly used to compute gene-level and exon-level expression values. However, a systematic comparison of these two algorithms on their effects on high-level analysis of the arrays has not yet been reported.

RESULTS

In this study, we showed that PLIER summarization led to over-estimation of gene-level expression changes, relative to exon-level expression changes, in two-group comparisons. Consequently, it led to detection of substantially more skipped exons on up-regulated genes, as well as substantially more included (i.e., non-skipped) exons on down-regulated genes. In contrast, this bias was not observed for RMA-summarized data. By using a published human tissue dataset, we compared the tissue-specific expression and splicing detected by Affymetrix exon arrays with those detected based on expressed sequence databases. We found the tendency of PLIER was not supported by the expressed sequence data.

CONCLUSION

We showed that the tendency of PLIER in detection of alternative splicing is likely caused by a technical bias in the approach, rather than a biological bias. Moreover, we observed abnormal summarization results when using the PLIER algorithm, indicating that mathematical problems, such as numerical instability, may affect PLIER performance.

Multidimensional gene set analysis of genomic data

Understanding the functional implications of changes in gene expression, mutations, etc., is the aim of most genomic experiments. To achieve this, several functional profiling methods have been proposed. Such methods study the behaviour of different gene modules (e.g. gene ontology terms) in response to one particular variable (e.g. differential gene expression). In spite to the wealth of information provided by functional profiling methods, a common limitation to all of them is their inherent unidimensional nature. In order to overcome this restriction we present a multidimensional logistic model that allows studying the relationship of gene modules with different genome-scale measurements (e.g. differential expression, genotyping association, methylation, copy number alterations, heterozygosity, etc.) simultaneously. Moreover, the relationship of such functional modules with the interactions among the variables can also be studied, which produces novel results impossible to be derived from the conventional unidimensional functional profiling methods. We report sound results of gene sets associations that remained undetected by the conventional one-dimensional gene set analysis in several examples. Our findings demonstrate the potential of the proposed approach for the discovery of new cell functionalities with complex dependences on more than one variable.

Understanding neuronal connectivity through the post-transcriptional toolkit

Post-transcriptional regulatory mechanisms have emerged as a critical component underlying the diversification and spatiotemporal control of the proteome during the establishment of precise neuronal connectivity. These mechanisms have been shown to be important for virtually all stages of assembling a neural network, from neurite guidance, branching, and growth to synapse morphogenesis and function. From the moment a gene is transcribed, it undergoes a series of post-transcriptional regulatory modifications in the nucleus and cytoplasm until its final deployment as a functional protein. Initially, a message is subjected to extensive structural regulation through alternative splicing, which is capable of greatly expanding the protein repertoire by generating, in some cases, thousands of functionally distinct isoforms from a single gene locus. Then, RNA packaging into neuronal transport granules and recognition by RNA-binding proteins and/or microRNAs is capable of restricting protein synthesis to selective locations and under specific input conditions. This ability of the post-transcriptional apparatus to expand the informational content of a cell and control the deployment of proteins in both spatial and temporal dimensions is a feature well adapted for the extreme morphological properties of neural cells. In this review, we describe recent advances in understanding how post-transcriptional regulatory mechanisms refine the proteomic complexity required for the assembly of intricate and specific neural networks.

The comparison of different pre- and post-analysis filters for determination of exon-level alternative splicing events using Affymetrix arrays

Understanding the biologic significance of alternative splicing has been impeded by the difficulty in systematically identifying and validating transcript isoforms. Current exon array workflows suggest several different filtration steps to reduce the number of tests and increase the detection of alternative splicing events. In this study, we examine the effects of the suggested pre-analysis filtration by detection above background P value or signal intensity. This is followed post-analytically by restriction of exon expression to a fivefold change between groups, limiting the analysis to known alternative splicing events, or using the intersection of the results from different algorithms. Combinations of the filters are also examined. We find that none of the filtering methods reduces the number of technical false-positive calls identified by visual inspection. These include edge effects, nonresponsive probe sets, and inclusion of intronic and untranslated region probe sets into transcript annotations. Modules for filtering the exon microarray data on the basis of annotation features are needed. We propose new approaches to data filtration that would reduce the number of technical false-positives and therefore, impact the time spent performing visual inspection of the exon arrays.

The ancient mammalian KRAB zinc finger gene cluster on human chromosome 8q24.3 illustrates principles of C2H2 zinc finger evolution associated with unique expression profiles in human tissues

Background

Expansion of multi-C2H2 domain zinc finger (ZNF) genes, including the Krüppel-associated box (KRAB) subfamily, paralleled the evolution of tetrapodes, particularly in mammalian lineages. Advances in their cataloging and characterization suggest that the functions of the KRAB-ZNF gene family contributed to mammalian speciation.

Results

Here, we characterized the human 8q24.3 ZNF cluster on the genomic, the phylogenetic, the structural and the transcriptome level. Six (ZNF7, ZNF34, ZNF250, ZNF251, ZNF252, ZNF517) of the seven locus members contain exons encoding KRAB domains, one (ZNF16) does not. They form a paralog group in which the encoded KRAB and ZNF protein domains generally share more similarities with each other than with other members of the human ZNF superfamily. The closest relatives with respect to their DNA-binding domain were ZNF7 and ZNF251. The analysis of orthologs in therian mammalian species revealed strong conservation and purifying selection of the KRAB-A and zinc finger domains. These findings underscore structural/functional constraints during evolution. Gene losses in the murine lineage (ZNF16, ZNF34, ZNF252, ZNF517) and potential protein truncations in primates (ZNF252) illustrate ongoing speciation processes. Tissue expression profiling by quantitative real-time PCR showed similar but distinct patterns for all tested ZNF genes with the most prominent expression in fetal brain. Based on accompanying expression signatures in twenty-six other human tissues ZNF34 and ZNF250 revealed the closest expression profiles. Together, the 8q24.3 ZNF genes can be assigned to a cerebellum, a testis or a prostate/thyroid subgroup. These results are consistent with potential functions of the ZNF genes in morphogenesis and differentiation. Promoter regions of the seven 8q24.3 ZNF genes display common characteristics like missing TATA-box, CpG island-association and transcription factor binding site (TFBS) modules. Common TFBS modules partly explain the observed expression pattern similarities.

Conclusions

The ZNF genes at human 8q24.3 form a relatively old mammalian paralog group conserved in eutherian mammals for at least 130 million years. The members persisted after initial duplications by undergoing subfunctionalizations in their expression patterns and target site recognition. KRAB-ZNF mediated repression of transcription might have shaped organogenesis in mammalian ontogeny.

PositionMatcher: A Fast Custom-Annotation Tool for Short DNA Sequences

Microarray probes and reads from massively parallel sequencing technologies are two most widely used genomic tags for a transcriptome study. Names and underlying technologies might differ, but expression technologies share a common objective-to obtain mRNA abundance values at the gene level, with high sensitivity and specificity. However, the initial tag annotation becomes obsolete as more insight is gained into biological references (genome, transcriptome, SNP, etc.). While novel alignment algorithms for short reads are being released every month, solutions for rapid annotation of tags are rare. We have developed a generic matching algorithm that uses genomic positions for rapid custom-annotation of tags with a time complexity O(nlogn). We demonstrate our algorithm on the custom annotation of Illumina massively parallel sequencing reads and Affymetrix microarray probes and identification of alternatively spliced regions.

Sam68 sequestration and partial loss of function are associated with splicing alterations in FXTAS patients

Fragile X-associated Tremor/Ataxia Syndrome (FXTAS) is a neurodegenerative disorder caused by expansion of 55-200 CGG repeats in the 5'-UTR of the FMR1 gene. FXTAS is characterized by action tremor, gait ataxia and impaired executive cognitive functioning. It has been proposed that FXTAS is caused by titration of RNA-binding proteins by the expanded CGG repeats. Sam68 is an RNA-binding protein involved in alternative splicing regulation and its ablation in mouse leads to motor coordination defects. Here, we report that mRNAs containing expanded CGG repeats form large and dynamic intranuclear RNA aggregates that recruit several RNA-binding proteins sequentially, first Sam68, then hnRNP-G and MBNL1. Importantly, Sam68 is sequestered by expanded CGG repeats and thereby loses its splicing-regulatory function. Consequently, Sam68-responsive splicing is altered in FXTAS patients. Finally, we found that regulation of Sam68 tyrosine phosphorylation modulates its localization within CGG aggregates and that tautomycin prevents both Sam68 and CGG RNA aggregate formation. Overall, these data support an RNA gain-of-function mechanism for FXTAS neuropathology, and suggest possible target routes for treatment options.

Nonsense-mediated RNA decay regulation by cellular stress: implications for tumorigenesis

Nonsense-mediated RNA decay (NMD) has long been viewed as an important constitutive mechanism to rapidly eliminate mutated mRNAs. More recently, it has been appreciated that NMD also degrades multiple nonmutated transcripts and that NMD can be regulated by wide variety of cellular stresses. Many of the stresses that inhibit NMD, including cellular hypoxia and amino acid deprivation, are experienced in cells exposed to hostile microenvironments, and several NMD-targeted transcripts promote cellular adaptation in response to these environmental stresses. Because adaptation to the microenvironment is crucial in tumorigenesis, and because NMD targets many mutated tumor suppressor gene transcripts, the regulation of NMD may have particularly important implications in cancer. This review briefly outlines the mechanisms by which transcripts are identified and targeted by NMD and reviews the evidence showing that NMD is a regulated process that can dynamically alter gene expression. Although much of the focus in NMD research has been in identifying the proteins that play a role in NMD and identifying NMD-targeted transcripts, recent data about the potential functional significance of NMD regulation, including the stabilization of alternatively spliced mRNA isoforms, the validation of mRNAs as bona fide NMD targets, and the role of NMD in tumorigenesis, are explored.

Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments

Alternative splicing, polyadenylation of pre-messenger RNA molecules and differential promoter usage can produce a variety of transcript isoforms whose respective expression levels are regulated in time and space, thus contributing specific biological functions. However, the repertoire of mammalian alternative transcripts and their regulation are still poorly understood. Second-generation sequencing is now opening unprecedented routes to address the analysis of entire transcriptomes. Here, we developed methods that allow the prediction and quantification of alternative isoforms derived solely from exon expression levels in RNA-Seq data. These are based on an explicit statistical model and enable the prediction of alternative isoforms within or between conditions using any known gene annotation, as well as the relative quantification of known transcript structures. Applying these methods to a human RNA-Seq dataset, we validated a significant fraction of the predictions by RT-PCR. Data further showed that these predictions correlated well with information originating from junction reads. A direct comparison with exon arrays indicated improved performances of RNA-Seq over microarrays in the prediction of skipped exons. Altogether, the set of methods presented here comprehensively addresses multiple aspects of alternative isoform analysis. The software is available as an open-source R-package called Solas at http://cmb.molgen.mpg.de/2ndGenerationSequencing/Solas/.

The role of HMGCR alternative splicing in statin efficacy

Statins, or 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR) inhibitors, are widely prescribed to lower plasma cholesterol levels and reduce cardiovascular disease risk. Despite the well-documented efficacy of statins, there is large interindividual variation in response. Using a panel of immortalized lymphocyte cell lines incubated with simvastatin, we recently found that the magnitude of expression of an alternatively spliced HMGCR transcript lacking exon 13 was inversely correlated with in vivo reductions of total cholesterol, low-density lipoprotein cholesterol, apoB, and triglycerides after statin treatment of the individuals from whom the cells were derived. This review will discuss the potential significance of alternative splicing as a mechanism contributing to variation in statin efficacy as well as the use of immortalized lymphocyte cell lines for identifying pharmacogenetically relevant polymorphisms and molecular mechanisms.

AGRONOMICS1: a new resource for Arabidopsis transcriptome profiling

Transcriptome profiling has become a routine tool in biology. For Arabidopsis (Arabidopsis thaliana), the Affymetrix ATH1 expression array is most commonly used, but it lacks about one-third of all annotated genes present in the reference strain. An alternative are tiling arrays, but previous designs have not allowed the simultaneous analysis of both strands on a single array. We introduce AGRONOMICS1, a new Affymetrix Arabidopsis microarray that contains the complete paths of both genome strands, with on average one 25mer probe per 35-bp genome sequence window. In addition, the new AGRONOMICS1 array contains all perfect match probes from the original ATH1 array, allowing for seamless integration of the very large existing ATH1 knowledge base. The AGRONOMICS1 array can be used for diverse functional genomics applications such as reliable expression profiling of more than 30,000 genes, detection of alternative splicing, and chromatin immunoprecipitation coupled to microarrays (ChIP-chip). Here, we describe the design of the array and compare its performance with that of the ATH1 array. We find results from both microarrays to be of similar quality, but AGRONOMICS1 arrays yield robust expression information for many more genes, as expected. Analysis of the ATH1 probes on AGRONOMICS1 arrays produces results that closely mirror those of ATH1 arrays. Finally, the AGRONOMICS1 array is shown to be useful for ChIP-chip experiments. We show that heterochromatic H3K9me2 is strongly confined to the gene body of target genes in euchromatic chromosome regions, suggesting that spreading of heterochromatin is limited outside of pericentromeric regions.

Splicing factor and exon profiling across human tissues

It has been shown that alternative splicing is especially prevalent in brain and testis when compared to other tissues. To test whether there is a specific propensity of these tissues to generate splicing variants, we used a single source of high-density microarray data to perform both splicing factor and exon expression profiling across 11 normal human tissues. Paired comparisons between tissues and an original exon-based statistical group analysis demonstrated after extensive RT-PCR validation that the cerebellum, testis, and spleen had the largest proportion of differentially expressed alternative exons. Variations at the exon level correlated with a larger number of splicing factors being expressed at a high level in the cerebellum, testis and spleen than in other tissues. However, this splicing factor expression profile was similar to a more global gene expression pattern as a larger number of genes had a high expression level in the cerebellum, testis and spleen. In addition to providing a unique resource on expression profiling of alternative splicing variants and splicing factors across human tissues, this study demonstrates that the higher prevalence of alternative splicing in a subset of tissues originates from the larger number of genes, including splicing factors, being expressed than in other tissues.

Exon-based clustering of murine breast tumor transcriptomes reveals alternative exons whose expression is associated with metastasis

In the field of bioinformatics, exon profiling is a developing area of disease-associated transcriptome analysis. In this study, we performed a microarray-based transcriptome analysis at the single exon level in mouse 4T1 primary mammary tumors with different metastatic capabilities. A novel bioinformatics platform was developed that identified 679 genes with differentially expressed exons in 4T1 tumors, many of which were involved in cell morphology and movement. Of 152 alternative exons tested by reverse transcription-PCR, 97 were validated as differentially expressed in primary tumors with different metastatic capability. This analysis revealed candidate progression genes, hinting at variations in protein functions by alternate exon usage. In a parallel effort, we developed a novel exon-based clustering analysis and identified alternative exons in tumor transcriptomes that were associated with dissemination of primary tumor cells to sites of pulmonary metastasis. This analysis also revealed that the splicing events identified by comparing primary tumors were not aberrant events. Lastly, we found that a subset of differentially spliced variant transcripts identified in the murine model was associated with poor prognosis in a large clinical cohort of patients with breast cancer. Our findings illustrate the utility of exon profiling to define novel theranostic markers for study in cancer progression and metastasis.

Alternative splicing: global insights

Following the original reports of pre-mRNA splicing in 1977, it was quickly realized that splicing together of different combinations of splice sites--alternative splicing--allows individual genes to generate more than one mRNA isoform. The full extent of alternative splicing only began to be revealed once large-scale genome and transcriptome sequencing projects began, rapidly revealing that alternative splicing is the rule rather than the exception. Recent technical innovations have facilitated the investigation of alternative splicing at a global scale. Splice-sensitive microarray platforms and deep sequencing allow quantitative profiling of very large numbers of alternative splicing events, whereas global analysis of the targets of RNA binding proteins reveals the regulatory networks involved in post-transcriptional gene control. Combined with sophisticated computational analysis, these new approaches are beginning to reveal the so-called 'RNA code' that underlies tissue and developmentally regulated alternative splicing, and that can be disrupted by disease-causing mutations.

Splicing, cis genetic variation and disease

Splicing is a post-transcriptional modification of RNA during which introns are removed and exons are joined. Most of the mammalian genes undergo constitutive and alternative splicing events. In addition to the strong signals of the splice sites, splicing is influenced at a distance by a range of trans factors that interact with cis regulatory elements and influence the spliceosome. The intention of the present mini-review is to give some insights into the complexity of this interaction and to introduce the consequences of some kinds of detrimental genetic variation on alternative splicing and disease.

RNA processing and its regulation: global insights into biological networks

In recent years views of eukaryotic gene expression have been transformed by the finding that enormous diversity can be generated at the RNA level. Advances in technologies for characterizing RNA populations are revealing increasingly complete descriptions of RNA regulation and complexity; for example, through alternative splicing, alternative polyadenylation and RNA editing. New biochemical strategies to map protein-RNA interactions in vivo are yielding transcriptome-wide insights into mechanisms of RNA processing. These advances, combined with bioinformatics and genetic validation, are leading to the generation of functional RNA maps that reveal the rules underlying RNA regulation and networks of biologically coherent transcripts. Together these are providing new insights into molecular cell biology and disease.

Determination of alternate splicing events using the Affymetrix Exon 1.0 ST arrays

Alternative splicing plays an important role in regulation of normal cellular function. Alternative splicing of pre-mRNA leads to the diversity of downstream protein products in the cell. The Affymetrix Exon arrays allow for a high throughput evaluation of the differences in spliced mRNA expressed in a biological system. In this study, we describe a method using this technology to study the generation of alternative mRNA transcripts in breast cancer cells that differ in the levels of a particular integrin, alpha3beta1.

Predicting functional alternative splicing by measuring RNA selection pressure from multigenome alignments

High-throughput methods such as EST sequencing, microarrays and deep sequencing have identified large numbers of alternative splicing (AS) events, but studies have shown that only a subset of these may be functional. Here we report a sensitive bioinformatics approach that identifies exons with evidence of a strong RNA selection pressure ratio (RSPR)--i.e., evolutionary selection against mutations that change only the mRNA sequence while leaving the protein sequence unchanged--measured across an entire evolutionary family, which greatly amplifies its predictive power. Using the UCSC 28 vertebrate genome alignment, this approach correctly predicted half to three-quarters of AS exons that are known binding targets of the NOVA splicing regulatory factor, and predicted 345 strongly selected alternative splicing events in human, and 262 in mouse. These predictions were strongly validated by several experimental criteria of functional AS such as independent detection of the same AS event in other species, reading frame-preservation, and experimental evidence of tissue-specific regulation: 75% (15/20) of a sample of high-RSPR exons displayed tissue specific regulation in a panel of ten tissues, vs. only 20% (4/20) among a sample of low-RSPR exons. These data suggest that RSPR can identify exons with functionally important splicing regulation, and provides biologists with a dataset of over 600 such exons. We present several case studies, including both well-studied examples (GRIN1) and novel examples (EXOC7). These data also show that RSPR strongly outperforms other approaches such as standard sequence conservation (which fails to distinguish amino acid selection pressure from RNA selection pressure), or pairwise genome comparison (which lacks adequate statistical power for predicting individual exons).

The epithelial splicing factors ESRP1 and ESRP2 positively and negatively regulate diverse types of alternative splicing events

Cell-type and tissue-specific alternative splicing events are regulated by combinatorial control involving both abundant RNA binding proteins as well as those with more discrete expression and specialized functions. Epithelial Splicing Regulatory Proteins 1 and 2 (ESRP1 and ESRP2) are recently discovered epithelial-specific RNA binding proteins that promote splicing of the epithelial variant of the FGFR2, ENAH, CD44 and CTNND1 transcripts. To catalogue a larger set of splicing events under the regulation of the ESRPs we profiled splicing changes induced by RNA interference-mediated knockdown of ES RP1 and ES RP2 expression in a human epithelial cell line using the splicing sensitive Affymetrix Exon ST1.0 Arrays. Analysis of the microarray data resulted in the identification of over a hundred candidate ESRP regulated splicing events. We were able to independently validate 38 of these targets by RT-PCR. The ESRP regulated events encompass all known types of alternative splicing events, most prominent being alternative cassette exons and splicing events leading to alternative 3' terminal exons. Importantly, a number of these regulated splicing events occur in gene transcripts that encode proteins with well-described roles in the regulation of actin cytoskeleton organization, cell-cell adhesion, cell polarity and cell migration. In sum, this work reveals a novel list of transcripts differentially spliced in epithelial and mesenchymal cells, implying that coordinated alternative splicing plays a critical role in determination of cell type identity. These results further establish ESRP1 and ESRP2 as global regulators of an epithelial splicing regulatory network.