SplicerAV: a tool for mining microarray expression data for changes in RNA processing

Background  

Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets.     

ExpressionPlot: a web-based framework for analysis of RNA-Seq and microarray gene expression data

RNA-Seq and microarray platforms have emerged as important tools for detecting changes in gene expression and RNA processing in biological samples. We present ExpressionPlot, a software package consisting of a default back end, which prepares raw sequencing or Affymetrix microarray data, and a web-based front end, which offers a biologically centered interface to browse, visualize, and compare different data sets. Download and installation instructions, a user's manual, discussion group, and a prototype are available at .     

Interactively optimizing signal-to-noise ratios in expression profiling: project-specific algorithm selection and detection p-value weighting in Affymetrix microarrays

MOTIVATION: The most commonly utilized microarrays for mRNA profiling (Affymetrix) include 'probe sets' of a series of perfect match and mismatch probes (typically 22 oligonucleotides per probe set). There are an increasing number of reported 'probe set algorithms' that differ in their interpretation of a probe set to derive a single normalized 'signal' representative of expression of each mRNA. These algorithms are known to differ in accuracy and sensitivity, and optimization has been done using a small set of standardized control microarray data. We hypothesized that different mRNA profiling projects have varying sources and degrees of confounding noise, and that these should alter the choice of a specific probe set algorithm. Also, we hypothesized that use of the Microarray Suite (MAS) 5.0 probe set detection p-value as a weighting function would improve the performance of all probe set algorithms. RESULTS: We built an interactive visual analysis software tool (HCE2W) to test and define parameters in Affymetrix analyses that optimize the ratio of signal (desired biological variable) versus noise (confounding uncontrolled variables). Five probe set algorithms were studied with and without statistical weighting of probe sets using the MAS 5.0 probe set detection p-values. The signal-to-noise ratio optimization method was tested in two large novel microarray datasets with different levels of confounding noise, a 105 sample U133A human muscle biopsy dataset (11 groups: mutation-defined, extensive noise), and a 40 sample U74A inbred mouse lung dataset (8 groups: little noise). Performance was measured by the ability of the specific probe set algorithm, with and without detection p-value weighting, to cluster samples into the appropriate biological groups (unsupervised agglomerative clustering with F-measure values). Of the total random sampling analyses, 50% showed a highly statistically significant difference between probe set algorithms by ANOVA [F(4,10) > 14, p < 0.0001], with weighting by MAS 5.0 detection p-value showing significance in the mouse data by ANOVA [F(1,10) > 9, p < 0.013] and paired t-test [t(9) = -3.675, p = 0.005]. Probe set detection p-value weighting had the greatest positive effect on performance of dChip difference model, ProbeProfiler and RMA algorithms. Importantly, probe set algorithms did indeed perform differently depending on the specific project, most probably due to the degree of confounding noise. Our data indicate that significantly improved data analysis of mRNA profile projects can be achieved by optimizing the choice of probe set algorithm with the noise levels intrinsic to a project, with dChip difference model with MAS 5.0 detection p-value continuous weighting showing the best overall performance in both projects. Furthermore, both existing and newly developed probe set algorithms should incorporate a detection p-value weighting to improve performance. AVAILABILITY: The Hierarchical Clustering Explorer 2.0 is available at http://www.cs.umd.edu/hcil/hce/ Murine arrays (40 samples) are publicly available at the PEPR resource (http://microarray.cnmcresearch.org/pgadatatable.asp http://pepr.cnmcresearch.org Chen et al., 2004).     

mRNA isoform diversity can obscure detection of miRNA-mediated control of translation

Reporter-based studies support inhibition of translation at the level of initiation as a substantial component of the miRNA mechanism, yet recent global analyses have suggested that they predominantly act through decreasing target mRNA stability. Cells commonly coexpress several processing isoforms of an mRNA, which may also differ in their regulatory untranslated regions (UTR). In particular, cancer cells are known to express high levels of short 3' UTR isoforms that evade miRNA-mediated regulation, whereas longer 3' UTRs predominate in nontransformed cells. To test whether mRNA isoform diversity can obscure detection of miRNA-mediated control at the level of translation, we assayed the responses of 11 endogenous let-7 targets to inactivation of this miRNA in HeLa cells, an intensively studied model system. We show that translational regulation in many cases appears to be modest when measuring the composite polysome profile of all extant isoforms of a given mRNA by density ultracentrifugation. In contrast, we saw clear effects at the level of translation initiation for multiple examples when selectively profiling mRNA isoforms carrying the 5' or 3' untranslated regions that were actually permissive to let-7 action, or when let-7 and a second targeting miRNA were jointly manipulated. Altogether, these results highlight a caveat to the mechanistic interpretation of data from global miRNA target analyses in transformed cells. Importantly, they reaffirm the importance of translational control as part of the miRNA mechanism in animal cells.     

Generation of troponin T isoforms by alternative RNA splicing in avian skeletal muscle. Conserved and divergent features in birds and mammals

We describe the isolation and sequence analysis of quail muscle cDNA clones encoding two closely related isoforms of the striated muscle contractile protein, troponin T. The cDNAs represent two troponin T mRNAs that exhibit an unusual sequence relationship. The two mRNAs have identical sequences over hundreds of nucleotides including 3' untranslated regions, but they differ dramatically in a discrete, internally located block of 38 nucleotides. The two alternative sequences of this 38-nucleotide block encode two different but related versions of amino acid residues 230-242, near the C terminus of the protein. These results are consistent with a novel mechanism of troponin T isoform generation by alternative mRNA splicing pathways from a single gene containing two different exons corresponding to amino acids 229-242, as recently proposed by Medford et al. (Medford, R. M., Nguyen, H. T., Destree, A. T., Summers, E., and Nadal-Ginard, B. (1984) Cell 38, 409-421). This proposal was based on analysis of a rat troponin T genomic DNA clone and a cDNA clone corresponding to one of the two alternatively spliced mRNAs. Our analysis of quail troponin T cDNA clones, apparently corresponding to two alternatively spliced mRNA species, provides important new evidence for this novel mechanism of troponin T isoform generation and reveals the differential splicing mechanism to be of great antiquity, antedating the bird-mammal divergence. One of the quail alternative isoform sequences clearly corresponds to one of the rat sequences, but the other quail alternative sequence does not correspond to either of the rat sequences. This result suggests a greater complexity of troponin T gene structure or a greater diversity of troponin T isoform genes than is currently known, and also has implications for the functional significance of the troponin T protein isoform heterogeneity. Comparison of quail and mammal alternative isoform sequences also reveals strongly conserved features which suggest that all the isoform alternative amino acid sequences are variations on a common structural theme.     

Genome-Wide Identification of Alternative Splice Forms Down-Regulated by Nonsense-Mediated mRNA Decay in Drosophila

Alternative mRNA splicing adds a layer of regulation to the expression of thousands of genes in Drosophila melanogaster. Not all alternative splicing results in functional protein; it can also yield mRNA isoforms with premature stop codons that are degraded by the nonsense-mediated mRNA decay (NMD) pathway. This coupling of alternative splicing and NMD provides a mechanism for gene regulation that is highly conserved in mammals. NMD is also active in Drosophila, but its effect on the repertoire of alternative splice forms has been unknown, as has the mechanism by which it recognizes targets. Here, we have employed a custom splicing-sensitive microarray to globally measure the effect of alternative mRNA processing and NMD on Drosophila gene expression. We have developed a new algorithm to infer the expression change of each mRNA isoform of a gene based on the microarray measurements. This method is of general utility for interpreting splicing-sensitive microarrays and high-throughput sequence data. Using this approach, we have identified a high-confidence set of 45 genes where NMD has a differential effect on distinct alternative isoforms, including numerous RNA–binding and ribosomal proteins. Coupled alternative splicing and NMD decrease expression of these genes, which may in turn have a downstream effect on expression of other genes. The NMD–affected genes are enriched for roles in translation and mitosis, perhaps underlying the previously observed role of NMD factors in cell cycle progression. Our results have general implications for understanding the NMD mechanism in fly. Most notably, we found that the NMD–target mRNAs had significantly longer 3′ untranslated regions (UTRs) than the nontarget isoforms of the same genes, supporting a role for 3′ UTR length in the recognition of NMD targets in fly.     

Alternatively spliced products of the human kinesin light chain 1 (KNS2) gene

Conventional kinesin is a microtubule-based molecular motor involved in the transport of membranous and non-membranous cargoes. The kinesin holoenzyme exists as a heterotetramer, consisting of two heavy chain and two light chain subunits. It is thought that one function of the light chains is to interact with the cargo. Alternative splicing of kinesin light chain pre-mRNA has been observed in lower organisms, although evidence for alternative splicing of the human gene has not been reported. We have identified 19 variants of the human KNS2 gene ( KLC1 ) that are generated by alternative splicing of downstream exons, but calculate that KNS2 has the potential to produce 285 919 spliceforms. Corresponding spliceforms of the mouse KLC1 gene were also identified. The alternative exons are all located 3' of exon 12 and the novel spliceforms produce both alternative carboxy termini and alternative 3' untranslated regions. The observation of multiple light chain isoforms is consistent with their proposed role in specific cargo attachment.