Single molecule profiling of alternative pre-mRNA splicing

CD44 alternative splicing senses intragenic DNA methylation in tumors via direct and indirect mechanisms

DNA methylation (meDNA) is a modulator of alternative splicing, and splicing perturbations are involved in tumorigenesis nearly as frequently as DNA mutations. However, the impact of meDNA on tumorigenesis via splicing-mediated mechanisms has not been thoroughly explored. Here, we found that HCT116 colon carcinoma cells inactivated for the DNA methylases DNMT1/3b undergo a partial epithelial to mesenchymal transition associated with increased CD44 variant exon skipping. These skipping events are directly mediated by the loss of intragenic meDNA and the chromatin factors MBD1/2/3 and HP1γ and are also linked to phosphorylation changes in elongating RNA polymerase II. The role of meDNA in alternative splicing was confirmed by using the dCas9/DNMT3b tool. We further tested whether the meDNA level could have predictive value in the MCF10A model for breast cancer progression and in patients with acute lymphoblastic leukemia (B ALL). We found that a small number of differentially spliced genes, mostly involved in splicing and signal transduction, are correlated with the local modulation of meDNA. Our observations suggest that, although DNA methylation has multiple avenues to affect alternative splicing, its indirect effect may also be mediated through alternative splicing isoforms of these meDNA sensors.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Highly sensitive and multiplexed quantification of mRNA splice variants by the direct ligation of DNA probes at the exon junction and universal PCR amplification

Alternative messenger RNA (mRNA) splicing is a basic mechanism of gene regulation. In general, reverse transcription and polymerase based primer extension limit the sensitivity and selectivity of the current detection of mRNA splice variants, respectively. Here, we show that, using the ligation of two properly designed probes at the exon junction combined with universal PCR amplification, as little as a single copy of a mRNA splice variant per cell can be accurately determined, and the dynamic range covers six orders of magnitude. Three mRNA splice variants were measured from total RNA samples derived from different cell lines. Moreover, by encoding the ligation probes with different lengths, multiplexed mRNA splice variants can be simultaneously detected in one-tube PCR amplification using electrophoretic separation.

A homozygous PMS2 founder mutation with an attenuated constitutional mismatch repair deficiency phenotype

BACKGROUND

Inherited mutations in DNA mismatch repair genes predispose to different cancer syndromes depending on whether they are mono-allelic or bi-allelic. This supports a causal relationship between expression level in the germline and phenotype variation. As a model to study this relationship, our study aimed to define the pathogenic characteristics of a recurrent homozygous coding variant in PMS2 displaying an attenuated phenotype identified by clinical genetic testing in seven Inuit families from Northern Quebec.

METHODS

Pathogenic characteristics of the PMS2 mutation NM_000535.5:c.2002A>G were studied using genotype-phenotype correlation, single-molecule expression detection and single genome microsatellite instability analysis.

RESULTS

This PMS2 mutation generates a de novo splice site that competes with the authentic site. In homozygotes, expression of the full-length protein is reduced to a level barely detectable by conventional diagnostics. Median age at primary cancer diagnosis is 22 years among 13 NM_000535.5:c.2002A>G homozygotes, versus 8 years in individuals carrying bi-allelic truncating mutations. Residual expression of full-length PMS2 transcript was detected in normal tissues from homozygotes with cancers in their 20s.

CONCLUSIONS

Our genotype-phenotype study of c.2002A>G illustrates that an extremely low level of PMS2 expression likely delays cancer onset, a feature that could be exploited in cancer preventive intervention.

Review: Alternative Splicing (AS) of Genes As An Approach for Generating Protein Complexity

Prior to the completion of the human genome project, the human genome was thought to have a greater number of genes as it seemed structurally and functionally more complex than other simpler organisms. This along with the belief of “one gene, one protein”, were demonstrated to be incorrect. The inequality in the ratio of gene to protein formation gave rise to the theory of alternative splicing (AS). AS is a mechanism by which one gene gives rise to multiple protein products. Numerous databases and online bioinformatic tools are available for the detection and analysis of AS. Bioinformatics provides an important approach to study mRNA and protein diversity by various tools such as expressed sequence tag (EST) sequences obtained from completely processed mRNA. Microarrays and deep sequencing approaches also aid in the detection of splicing events. Initially it was postulated that AS occurred only in about 5% of all genes but was later found to be more abundant. Using bioinformatic approaches, the level of AS in human genes was found to be fairly high with 35-59% of genes having at least one AS form. Our ability to determine and predict AS is important as disorders in splicing patterns may lead to abnormal splice variants resulting in genetic diseases. In addition, the diversity of proteins produced by AS poses a challenge for successful drug discovery and therefore a greater understanding of AS would be beneficial.

Tau alternative splicing in familial and sporadic tauopathies

Six tau isoforms differing in their affinity for microtubules are produced by alternative splicing from the MAPT (microtubule-associated protein tau) gene in adult human brain. Several MAPT mutations causing the familial tauopathy, FTDP-17 (frontotemporal dementia with parkinsonism linked to chromosome 17), affect alternative splicing of exon 10, encoding a microtubule-binding motif. Advanced RNA analysis methods have suggested that levels of exon 10-containing MAPT mRNA are elevated in Alzheimer's disease. Furthermore, the MAPT H1 haplotype, associated with Alzheimer's disease, promotes exon 10 inclusion in MAPT mRNA. Thus an accurate regulation of tau alternative splicing is critical for the maintenance of neuronal viability, and its alteration might be a contributing factor to Alzheimer's disease. Tau alternative splicing could represent a target for therapeutic intervention to delay the progression of pathology in familial as well as sporadic tauopathies.

Human mucin MUC1 RNA undergoes different types of alternative splicing resulting in multiple isoforms

MUC1 is a transmembrane mucin with important functions in normal and transformed cells, carried out by the extracellular domain or the cytoplasmic tail. A characteristic feature of the MUC1 extracellular domain is the variable number of tandem repeats (VNTR) region. Alternative splicing may regulate MUC1 expression and possibly function. We developed an RT-PCR method for efficient isolation of MUC1 mRNA isoforms that allowed us to evaluate the extent of alternative splicing of MUC1 and elucidate some of the rules that govern this process. We cloned and analyzed 21, 24, and 36 isoforms from human tumor cell lines HeLa, MCF7, and Jurkat, respectively, and 16 from normal activated human T cells. Among the 78 MUC1 isoforms we isolated, 76 are new and different cells showed varied MUC1 expression patterns. The VNTR region of exon 2 was recognized as an intron with a fixed 5' splice site but variable 3' splice sites. We also report that the 3506 A/G SNP in exon 2 can regulate 3' splice sites selection in intron 1 and produce different MUC1 short isoform proteins. Furthermore, the SNP A to G mutation was also observed in vivo, during de novo tumor formation in MUC1(+/-)Kras(G12D/+)Pten(loxP/loxP) mice. No specific functions have been associated with previously reported short isoforms. We now report that one new G SNP-associated isoform MUC1/Y-LSP, but not the A SNP-associated isoform MUC1/Y, inhibits tumor growth in immunocompetent but not immunocompromised mice.

Improving oligonucleotide fingerprinting of rRNA genes by implementation of polony microarray technology

Improvements to oligonucleotide fingerprinting of rRNA genes (OFRG) were obtained by implementing polony microarray technology. OFRG is an array-based method for analyzing microbial community composition. Polonies are discrete clusters of DNA, produced by solid-phase PCR in hydrogels, and derived from individual, spatially isolated DNA molecules. The advantages of a polony-based OFRG method include higher throughput and reductions in the PCR-induced errors and compositional skew inherent in all other PCR-based community composition methods, including high-throughput sequencing of rRNA genes. Given the similarities between polony microarrays and certain aspects of sequencing methods such as the Illumina platform, we suggest that if concepts presented in this study were implemented in high-throughput sequencing protocols, a reduction of PCR-induced errors and compositional skew may be realized.

Massively parallel haplotyping on microscopic beads for the high-throughput phase analysis of single molecules

In spite of the many advances in haplotyping methods, it is still very difficult to characterize rare haplotypes in tissues and different environmental samples or to accurately assess the haplotype diversity in large mixtures. This would require a haplotyping method capable of analyzing the phase of single molecules with an unprecedented throughput. Here we describe such a haplotyping method capable of analyzing in parallel hundreds of thousands single molecules in one experiment. In this method, multiple PCR reactions amplify different polymorphic regions of a single DNA molecule on a magnetic bead compartmentalized in an emulsion drop. The allelic states of the amplified polymorphisms are identified with fluorescently labeled probes that are then decoded from images taken of the arrayed beads by a microscope. This method can evaluate the phase of up to 3 polymorphisms separated by up to 5 kilobases in hundreds of thousands single molecules. We tested the sensitivity of the method by measuring the number of mutant haplotypes synthesized by four different commercially available enzymes: Phusion, Platinum Taq, Titanium Taq, and Phire. The digital nature of the method makes it highly sensitive to detecting haplotype ratios of less than 1:10,000. We also accurately quantified chimera formation during the exponential phase of PCR by different DNA polymerases.

Single molecule analysis of c-myb alternative splicing reveals novel classifiers for precursor B-ALL

The c-Myb transcription factor, a key regulator of proliferation and differentiation in hematopoietic and other cell types, has an N-terminal DNA binding domain and a large C-terminal domain responsible for transcriptional activation, negative regulation and determining target gene specificity. Overexpression and rearrangement of the c-myb gene (MYB) has been reported in some patients with leukemias and other types of cancers, implicating activated alleles of c-myb in the development of human tumors. Alternative RNA splicing can produce variants of c-myb with qualitatively distinct transcriptional activities that may be involved in transformation and leukemogenesis. Here, by performing a detailed, single molecule assay we found that c-myb alternative RNA splicing was elevated and much more complex in leukemia samples than in cell lines or CD34+ hematopoietic progenitor cells from normal donors. The results revealed that leukemia samples express more than 60 different c-myb splice variants, most of which have multiple alternative splicing events and were not detectable by conventional microarray or PCR approaches. For example, the single molecule assay detected 21 and 22 splice variants containing the 9B and 9S exons, respectively, most of which encoded unexpected variant forms of c-Myb protein. Furthermore, the detailed analysis identified some splice variants whose expression correlated with poor survival in a small cohort of precursor B-ALL samples. Our findings indicate that single molecule assays can reveal complexities in c-myb alternative splicing that have potential as novel biomarkers and could help explain the role of c-Myb variants in the development of human leukemia.

Key features of the two-intron Saccharomyces cerevisiae gene SUS1 contribute to its alternative splicing

Alternative pre-mRNA splicing allows dramatic expansion of the eukaryotic proteome and facilitates cellular response to changes in environmental conditions. The Saccharomyces cerevisiae gene SUS1, which encodes a protein involved in mRNA export and histone H2B deubiquitination, contains two introns; non-canonical sequences in the first intron contribute to its retention, a common form of alternative splicing in plants and fungi. Here we show that the pattern of SUS1 splicing changes in response to environmental change such as temperature elevation, and the retained intron product is subject to nonsense-mediated decay. The activities of different splicing factors determine the pattern of SUS1 splicing, including intron retention and exon skipping. Unexpectedly, removal of the 3′ intron is affected by splicing of the upstream intron, suggesting that cross-exon interactions influence intron removal. Production of different SUS1 isoforms is important for cellular function, as we find that the temperature sensitivity and histone H2B deubiquitination defects observed in sus1Δ cells are only partially suppressed by SUS1 cDNA, but SUS1 that is able to undergo splicing complements these phenotypes. These data illustrate a role for S. cerevisiae alternative splicing in histone modification and cellular function and reveal important mechanisms for splicing of yeast genes containing multiple introns.

Fine-mapping resolves Eae23 into two QTLs and implicates ZEB1 as a candidate gene regulating experimental neuroinflammation in rat

Background

To elucidate mechanisms involved in multiple sclerosis (MS), we studied genetic regulation of experimental autoimmune encephalomyelitis (EAE) in rats, assuming a conservation of pathogenic pathways. In this study, we focused on Eae23, originally identified to regulate EAE in a (LEW.1AV1xPVG.1AV1)F2 cross. Our aim was to determine whether one or more genes within the 67 Mb region regulate EAE and to define candidate risk genes.

Methodology/Principal Findings

We used high resolution quantitative trait loci (QTL) analysis in the 10th generation (G10) of an advanced intercross line (AIL) to resolve Eae23 into two QTLs that independently regulate EAE, namely Eae23a and Eae23b. We established a congenic strain to validate the effect of this region on disease. PVG alleles in Eae23 resulted in significant protection from EAE and attenuated CNS inflammation/demyelination. Disease amelioration was accompanied with increased levels of Foxp3⁺ cells in the CNS of the congenic strain compared to DA. We then focused on candidate gene investigation in Eae23b, a 9 Mb region linked to all clinical phenotypes. Affymetrix exon arrays were used to study expression of the genes in Eae23b in the parental strains, where none showed differential expression. However, we found lower expression of exon 4 of ZEB1, which is specific for splice-variant Zfhep1. ZEB1 is an interleukin 2 (IL2) repressor involved in T cell development. The splice-specific variance prompted us to next analyze the expression of ZEB1 and its two splice variants, Zfhep1 and Zfhep2, in both lymph node and spleen. We demonstrated that ZEB1 splice-variants are differentially expressed; severity of EAE and higher IL2 levels were associated with down-regulation of Zfhep1 and up-regulation of Zfhep2.

Conclusions/Significance

We speculate that the balance between splice-variants of ZEB1 could influence the regulation of EAE. Further functional studies of ZEB1 and the splice-variants may unravel novel pathways contributing to MS pathogenesis and inflammation in general.

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.

Characterization of DNA polymerase beta splicing variants in gastric cancer: the most frequent exon 2-deleted isoform is a non-coding RNA

DNA repair polymerase beta (Pol beta) gene variants are frequently associated with tumor tissues. In this study a search for Pol beta mutants and splice variants was conducted in matched normal and tumor gastric tissues and blood samples from healthy donors. No tumor associated mutations were found while a variety of alternative Pol beta splicing variants were detected with high frequency in all the specimens analysed. Quantitative PCR of the Pol beta variant lacking exon 2 (Ex2Delta) and the isoforms with exon 11 skipping allowed to clarify that these variants are not tumor- neither tissue-specific and their levels vary greatly among different individuals. The most frequent Ex2Delta variant was further characterized. We clearly demonstrated that this variant does not encode protein, as detected by both western blotting and immunofluorescence analysis of human AGS cells expressing HA-tagged Ex2Delta. The lack of translation was confirmed by comparing the DNA gap-filling capacity and alkylation sensitivity of wild type and Pol beta null murine fibroblasts expressing the human Ex2Delta variant. We showed that the Ex2Delta transcript is polyadenylated and its half-life is significantly longer than that of the wild type mRNA as inferred by treating AGS cells with actinomycin D. Moreover, we found that it localizes to polyribosomes suggesting a role as post-transcriptional regulator. This study identifies a new type of DNA repair variants that do not give rise to functional proteins but to non-coding RNAs that could either modulate target mRNAs or represent unproductive splicing events.

Stochastic noise in splicing machinery

The number of known alternative human isoforms has been increasing steadily with the amount of available transcription data. To date, over 100 000 isoforms have been detected in EST libraries, and at least 75% of human genes have at least one alternative isoform. In this paper, we propose that most alternative splicing events are the result of noise in the splicing process. We show that the number of isoforms and their abundance can be predicted by a simple stochastic noise model that takes into account two factors: the number of introns in a gene and the expression level of a gene. The results strongly support the hypothesis that most alternative splicing is a consequence of stochastic noise in the splicing machinery, and has no functional significance. The results are also consistent with error rates tuned to ensure that an adequate level of functional product is produced and to reduce the toxic effect of accumulation of misfolding proteins. Based on simulation of sampling of virtual cDNA libraries, we estimate that error rates range from 1 to 10% depending on the number of introns and the expression level of a gene.

Hepatitis C virus infection of T cells inhibits proliferation and enhances fas-mediated apoptosis by down-regulating the expression of CD44 splicing variant 6

BACKGROUND

A lymphotropic hepatitis C virus strain (HCV, SB strain, hereafter "SB-HCV") has been shown to infect established T cell lines (Molt-4 and Jurkat) and primary human naive CD4(+) T cells. During T cell development and activation, transient expression of CD44 splicing variant 6 (CD44v6) plays a significant role.

METHODS

SB-HCV was used to infect Molt-4 cells, and their cellular proliferation and CD44 expression was examined.

RESULTS

SB-HCV-infected Molt-4 cells expressed a significantly lower level of the CD44v6 isoform. The infected cells could be divided into 2 carboxyfluorescein succinimidyl ester (CFSE) groups, CFSE-high (indicating low proliferation activity; 34.2% of the cells) and CFSE-low (indicating high proliferation activity; 62.5% of the cells), whereas uninfected cells consisted of only a CFSE-low population. Of the CFSE-high cells, 82.4% were positive for the HCV protein NS5A, whereas only 1.2% of the CFSE-low cells were positive for this protein. Among the HCV proteins, NS5A alone caused the down-regulation of CD44v6 expression. After cells were stimulated with phorbol myristate acetate, the amount of phosphorylated mitogen-activated protein (MAP) kinase was significantly reduced in CFSE-high, SB-HCV-infected Molt-4 cells. After Fas ligand stimulation, SB-HCV-infected Molt-4 cells had increased cleavage of caspase 8 and 3 and enhanced apoptosis, compared with the rates of cleavage and apoptosis in control groups, indicating that SB-HCV infection increased Fas-mediated apoptosis.

CONCLUSION

HCV replication in T cells suppresses cellular proliferation and enhances susceptibility to Fas signaling by inhibiting CD44v6 signaling and expression.

Quantitative analysis of somatic mitochondrial DNA mutations by single-cell single-molecule PCR

Mitochondrial genome integrity is an important issue in somatic mitochondrial genetics. Development of quantitative methods is indispensable to somatic mitochondrial genetics as quantitative studies are required to characterize heteroplasmy and mutation processes, as well as their effects on phenotypic developments. Quantitative studies include the identification and measurement of the load of pathogenic and non-pathogenic clonal mutations, screening mitochondrial genomes for mutations in order to determine the mutation spectra and characterize an ongoing mutation process. Single-molecule PCR (smPCR) has been shown to be an effective method that can be applied to all areas of quantitative studies. It has distinct advantages over conventional vector-based cloning techniques avoiding the well-known PCR-related artifacts such as the introduction of artificial mutations, preferential allelic amplifications, and "jumping" PCR. smPCR is a straightforward and robust method, which can be effectively used for molecule-by-molecule mutational analysis, even when mitochondrial whole genome (mtWG) analysis is involved. This chapter describes the key features of the smPCR method and provides three examples of its applications in single-cell analysis: di-plex smPCR for deletion quantification, smPCR cloning for clonal point mutation quantification, and smPCR cloning for whole genome sequencing (mtWGS).

Novel splice variants derived from the receptor tyrosine kinase superfamily are potential therapeutics for rheumatoid arthritis

INTRODUCTION

Despite the advent of biological therapies for the treatment of rheumatoid arthritis, there is a compelling need to develop alternative therapeutic targets for nonresponders to existing treatments. Soluble receptors occur naturally in vivo, such as the splice variant of the cell surface receptor for vascular endothelial growth factor (VEGF)--a key regulator of angiogenesis in rheumatoid arthritis. Bioinformatics analyses predict that the majority of human genes undergo alternative splicing, generating proteins--many of which may have regulatory functions. The objective of the present study was to identify alternative splice variants (ASV) from cell surface receptor genes, and to determine whether the novel proteins encoded exert therapeutic activity in an in vivo model of arthritis.

METHODS

To identify novel splice variants, we performed RT-PCR using an mRNA pool representing major human tissue types and tumors. Novel ASV were identified by alignment of each cloned sequence to its respective genomic sequence in comparison with full-length transcripts. To test whether these ASV have biologic activity, we characterized a subset of them for ligand binding, and for efficacy in an animal model of arthritis. The in vivo study was accomplished using adenoviruses expressing secreted ASV.

RESULTS

We cloned 60 novel human ASV from 21 genes, encoding cell surface receptors--many of which are known to be important in the regulation of angiogenesis. The ASV were characterized by exon extension, intron retention and alternative exon utilization. Efficient expression and secretion of selected ASV--corresponding to VEGF receptor type 1, VEGF receptor type 2, VEGF receptor type 3, angiopoietin receptor Tie1, Met (receptor for hepatocyte growth factor), colony-stimulating factor 1 receptor, platelet-derived growth factor receptor beta, fibroblast growth factor receptor 1, Kit, and RAGE--was demonstrated, together with binding to their cognate ligands. Importantly, ASV derived from VEGF receptor type 1 and Tie1, and to a lesser extent from VEGF receptor type 2 and fibroblast growth factor receptor 1, reduced clinical signs of arthritis in vivo. The reduction was paralleled by decreased joint inflammation and destruction.

CONCLUSION

The present study shows that unique ASV derived from receptors that play key roles in angiogenesis--namely, VEGF receptor type 1 and, for the first time, Tie1--can markedly reduce arthritis severity. More broadly, our results demonstrate that ASV are a source of novel proteins with therapeutic potential in diseases in which angiogenesis and cellular hyperplasia play a central role, such as rheumatoid arthritis.

Functional importance of different patterns of correlation between adjacent cassette exons in human and mouse

Background

Alternative splicing expands transcriptome diversity and plays an important role in regulation of gene expression. Previous studies focus on the regulation of a single cassette exon, but recent experiments indicate that multiple cassette exons within a gene may interact with each other. This interaction can increase the potential to generate various transcripts and adds an extra layer of complexity to gene regulation. Several cases of exon interaction have been discovered. However, the extent to which the cassette exons coordinate with each other remains unknown.

Results

Based on EST data, we employed a metric of correlation coefficients to describe the interaction between two adjacent cassette exons and then categorized these exon pairs into three different groups by their interaction (correlation) patterns. Sequence analysis demonstrates that strongly-correlated groups are more conserved and contain a higher proportion of pairs with reading frame preservation in a combinatorial manner. Multiple genome comparison further indicates that different groups of correlated pairs have different evolutionary courses: (1) The vast majority of positively-correlated pairs are old, (2) most of the weakly-correlated pairs are relatively young, and (3) negatively-correlated pairs are a mixture of old and young events.

Conclusion

We performed a large-scale analysis of interactions between adjacent cassette exons. Compared with weakly-correlated pairs, the strongly-correlated pairs, including both the positively and negatively correlated ones, show more evidence that they are under delicate splicing control and tend to be functionally important. Additionally, the positively-correlated pairs bear strong resemblance to constitutive exons, which suggests that they may evolve from ancient constitutive exons, while negatively and weakly correlated pairs are more likely to contain newly emerging exons.

Polony DNA sequencing

Polony DNA sequencing provides an inexpensive, accurate, high-throughput way to resequence genomes of interest by comparison to a reference genome. Mate-paired in vitro shotgun genomic libraries are produced and clonally amplified on microbeads by emulsion PCR. These serve as templates for sequencing by fluorescent nonamer ligation reactions on a microscope slide. Each sequencing run results in millions of 26-bp reads that can be aligned to the reference genome, allowing the identification of differences between sequences.

Implementation of alternating excitation schemes in a biochip-reader for quasi-simultaneous multi-color single-molecule detection

We report here the development of a device for single-molecule biochip readout using fast alternating excitation. The technology extends standard imaging cytometry by offering additional color channels in excitation. To enable the study of mobile objects, e.g. actively transported vesicles in living cells or freely diffusing lipids in a lipid bilayer, the frequency of the illumination pulses was chosen high enough to virtually freeze the motion of the biomolecules, as they are shifted through the illuminated area. The synchronization of sample illumination, scanning and line-camera readout yield two quasi-simultaneously recorded images covering the same sample region. Diffraction-limited resolution and high localization precision for point-light sources down to approximately 10 nm was shown by scanning immobilized 100 nm fluorescence latex beads. Ultra-sensitivity was demonstrated by imaging single fluorescent streptavidin molecules diffusing in a fluid lipid bilayer. Two-color streptavidin labeled with Cy3 and Cy5 could be easily identified in the two respective excitation channels; high accordance in the dye positions confirms the applicability for colocalization studies of moving objects. Finally, scans of antibody-receptor interactions in large populations of live cells illustrate the feasibility of this method for biochip application.

Polony analysis of gene expression in ES cells and blastocysts

Expression profiling of stem cells is challenging due to their small numbers and heterogeneity. The PCR colony (polony) approach has theoretical advantages as an assay for stem cells but has not been applied to small numbers of cells. An assay has been developed that is sensitive enough to detect mRNAs from small numbers of ES cells and from fractions of a single mouse blastocyst. Genes assayed include Oct3, Rex1, Nanog, Cdx2 and GLUT-1. The assay is highly sensitive so that multiple mRNAs from a single blastocyst were easily detected in the same assay. In its present version, the assay is an attractive alternative to conventional RT–PCR for profiling small populations of stem cells. The assay is also amenable to improvements that will increase its sensitivity and ability to analyze many cDNAs simultaneously.

Multisite and bidirectional exonic splicing enhancer in CD44 alternative exon v3

The human CD44 gene encodes multiple isoforms of a transmembrane protein that differ in their extracellular domains as a result of alternative splicing of its variable exons. Expression of CD44 is tightly regulated according to the type and physiological status of a cell, with expression of high molecular weight isoforms by inclusion of variable exons and low molecular weight isoforms containing few or no variable exons. Human CD44 variable exon 3 (v3) can follow a specific alternative splicing route different from that affecting other variable exons. Here we map and functionally describe the splicing enhancer element within CD44 exon v3 which regulates its inclusion in the final mRNA. The v3 splicing enhancer is a multisite bipartite element consisting of a tandem nonamer, the XX motif, and an heptamer, the Y motif, located centrally in the exon. Each of the three sites of this multisite enhancer partially retains its splicing enhancing capacity independently from each other in CD44 and shows full enhancing function in gene contexts different from CD44. We further demonstrate that these motifs act cooperatively as at least two motifs are needed to maintain exon inclusion. Their action is differential with respect to the splice-site target abutting v3. The first X motif acts on the 3' splice site, the second X motif acts on both splice sites (as a bidirectional exonic splicing enhancer), and the Y motif acts on the 5' splice site. We also show that the multisite v3 splicing enhancer is functional irrespective of flanking intron length and spatial organization within v3.

Genome-wide analysis of alternative pre-mRNA splicing

Alternative splicing of mRNA precursors allows the synthesis of multiple mRNAs from a single primary transcript, significantly expanding the information content and regulatory possibilities of higher eukaryotic genomes. High-throughput enabling technologies, particularly large-scale sequencing and splicing-sensitive microarrays, are providing unprecedented opportunities to address key questions in this field. The picture emerging from these pioneering studies is that alternative splicing affects most human genes and a significant fraction of the genes in other multicellular organisms, with the potential to greatly influence the evolution of complex genomes. A combinatorial code of regulatory signals and factors can deploy physiologically coherent programs of alternative splicing that are distinct from those regulated at other steps of gene expression. Pre-mRNA splicing and its regulation play important roles in human pathologies, and genome-wide analyses in this area are paving the way for improved diagnostic tools and for the identification of novel and more specific pharmaceutical targets.

Identification of differentially regulated splice variants and novel exons in glial brain tumors using exon expression arrays

Aberrant splice variants are involved in the initiation and/or progression of glial brain tumors. We therefore set out to identify splice variants that are differentially expressed between histologic subgroups of gliomas. Splice variants were identified using a novel platform that profiles the expression of virtually all known and predicted exons present in the human genome. Exon-level expression profiling was done on 26 glioblastomas, 22 oligodendrogliomas, and 6 control brain samples. Our results show that Human Exon arrays can identify subgroups of gliomas based on their histologic appearance and genetic aberrations. We next used our expression data to identify differentially expressed splice variants. In two independent approaches, we identified 49 and up to 459 exons that are differentially spliced between glioblastomas and oligodendrogliomas, a subset of which (47% and 33%) were confirmed by reverse transcription-PCR (RT-PCR). In addition, exon level expression profiling also identified >700 novel exons. Expression of approximately 67% of these candidate novel exons was confirmed by RT-PCR. Our results indicate that exon level expression profiling can be used to molecularly classify brain tumor subgroups, can identify differentially regulated splice variants, and can identify novel exons. The splice variants identified by exon level expression profiling may help to detect the genetic changes that cause or maintain gliomas and may serve as novel treatment targets.

Single molecule profiling of tau gene expression in Alzheimer's disease

Tau is a microtubule-associated protein that is important for establishing and maintaining neuronal morphology. In addition to its role in normal cells, tau protein is involved in many neurodegenerative diseases, e.g. Alzheimer's disease (AD) and frontotemporal dementia, as the main component of intraneuronal aggregates. Alternative splicing of tau gene in the brain can give rise to at least six protein variants. A causative role of skewed tau exon 10 inclusion has been defined in frontotemporal dementia; however, no link was established between the aberrant splicing of tau and AD. Here, we applied a single-molecule-based technology, polymerase colony or polony, to simultaneously monitor tau splicing variant and haplotype profile in sporadic AD and normal brains. We found that the coordinated expression of tau exons 2 and 10 is altered in AD. Additional investigations of cis and trans mechanisms of this observation revealed a decreased protein expression of a known tau splicing factor, htra2-beta-1 in AD, thereby implicating a trans mechanism. Our results demonstrate that dysregulation of combinatorial splicing might serve as a signature for aging-related diseases, and the polony assay could be widely adapted for the study of other tauopathies. Furthermore, splicing-based therapeutics is an emerging area of drug development, and a well-defined and quantitative assay for monitoring single-gene transcriptome will be relevant for such development.

Quantitative monitoring by polymerase colony assay of known mutations resistant to ABL kinase inhibitors

Resistance to molecularly targeted chemotherapy, and the development of novel agents that are active against resistant forms of target proteins create the need for a sensitive and quantitative assay to monitor drug-resistant mutations in patients to guide treatment and assess response. Here, we describe an application of the polymerase colony (polony) method to identify and quantify known point mutations in the BCR-ABL oncogene in patients with chronic myelogenous leukemia who evolve resistance to ABL kinase inhibitors. The assay can detect mutations with a sensitivity of 10(-4), quantify the burden of drug-resistant cells, and simultaneously monitor the dynamics of several coexisting mutations. As a proof of concept, we analysed blood samples from three patients undergoing therapy with ABL kinase inhibitors and found that the patients' response to therapy correlated with our molecular monitoring. We were also able to detect mutations emerging in patients long before clinical relapse. Therefore, the polony assay could be applied to a larger patient sample to assess the utility of early mutation detection in patient-specific treatment decisions. Finally, this methodology could be a valuable research tool to shed light on the natural behavior of mutations pre-existing kinase inhibitors therapy and either disappearing over time or slowly taking over.

Alternative splicing in cancer: Noise, functional, or systematic?

Pre-messenger RNA splicing is a fine-tuned process that generates multiple functional variants from individual genes. Various cell types and developmental stages regulate alternative splicing patterns differently in their generation of specific gene functions. In cancers, splicing is significantly altered, and understanding the underlying mechanisms and patterns in cancer will shed new light onto cancer biology. Cancer-specific transcript variants are promising biomarkers and targets for diagnostic, prognostic, and treatment purposes. In this review, we explore how alternative splicing cannot simply be considered as noise or an innocent bystander, but is actively regulated or deregulated in cancers. A special focus will be on aspects of cell biology and biochemistry of alternative splicing in cancer cells, addressing differences in splicing mechanisms between normal and malignant cells. The systems biology of splicing is only now applied to the field of cancer research. We explore functional annotations for some of the most intensely spliced gene classes, and provide a literature mining and clustering that reflects the most intensely investigated genes. A few well-established cancer-specific splice events, such as the CD44 antigen, are used to illustrate the potential behind the exploration of the mechanisms of their regulation. Accordingly, we describe the functional connection between the regulatory machinery (i.e., the spliceosome and its accessory proteins) and their global impact on qualitative transcript variation that are only now emerging from the use of genomic technologies such as microarrays. These studies are expected to open an entirely new level of genetic information that is currently still poorly understood.

RIP-Chip: the isolation and identification of mRNAs, microRNAs and protein components of ribonucleoprotein complexes from cell extracts

RNA targets of multitargeted RNA-binding proteins (RBPs) can be studied by various methods including mobility shift assays, iterative in vitro selection techniques and computational approaches. These techniques, however, cannot be used to identify the cellular context within which mRNAs associate, nor can they be used to elucidate the dynamic composition of RNAs in ribonucleoprotein (RNP) complexes in response to physiological stimuli. But by combining biochemical and genomics procedures to isolate and identify RNAs associated with RNA-binding proteins, information regarding RNA-protein and RNA-RNA interactions can be examined more directly within a cellular context. Several protocols--including the yeast three-hybrid system and immunoprecipitations that use physical or chemical cross-linking--have been developed to address this issue. Cross-linking procedures in general, however, are limited by inefficiency and sequence biases. The approach outlined here, termed RNP immunoprecipitation-microarray (RIP-Chip), allows the identification of discrete subsets of RNAs associated with multi-targeted RNA-binding proteins and provides information regarding changes in the intracellular composition of mRNPs in response to physical, chemical or developmental inducements of living systems. Thus, RIP-Chip can be used to identify subsets of RNAs that have related functions and are potentially co-regulated, as well as proteins that are associated with them in RNP complexes. Using RIP-Chip, the identification and/or quantification of RNAs in RNP complexes can be accomplished within a few hours or days depending on the RNA detection method used.

Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long disordered regions

Identifying relationships between function, amino acid sequence, and protein structure represents a major challenge. In this study, we propose a bioinformatics approach that identifies functional keywords in the Swiss-Prot database that correlate with intrinsic disorder. A statistical evaluation is employed to rank the significance of these correlations. Protein sequence data redundancy and the relationship between protein length and protein structure were taken into consideration to ensure the quality of the statistical inferences. Over 200,000 proteins from the Swiss-Prot database were analyzed using this approach. The predictions of intrinsic disorder were carried out using PONDR VL3E predictor of long disordered regions that achieves an accuracy of above 86%. Overall, out of the 710 Swiss-Prot functional keywords that were each associated with at least 20 proteins, 238 were found to be strongly positively correlated with predicted long intrinsically disordered regions, whereas 302 were strongly negatively correlated with such regions. The remaining 170 keywords were ambiguous without strong positive or negative correlation with the disorder predictions. These functions cover a large variety of biological activities and imply that disordered regions are characterized by a wide functional repertoire. Our results agree well with literature findings, as we were able to find at least one illustrative example of functional disorder or order shown experimentally for the vast majority of keywords showing the strongest positive or negative correlation with intrinsic disorder. This work opens a series of three papers, which enriches the current view of protein structure-function relationships, especially with regards to functionalities of intrinsically disordered proteins, and provides researchers with a novel tool that could be used to improve the understanding of the relationships between protein structure and function. The first paper of the series describes our statistical approach, outlines the major findings, and provides illustrative examples of biological processes and functions positively and negatively correlated with intrinsic disorder.

Multivariate analysis and visualization of splicing correlations in single-gene transcriptomes

BACKGROUND

RNA metabolism, through 'combinatorial splicing', can generate enormous structural diversity in the proteome. Alternative domains may interact, however, with unpredictable phenotypic consequences, necessitating integrated RNA-level regulation of molecular composition. Splicing correlations within transcripts of single genes provide valuable clues to functional relationships among molecular domains as well as genomic targets for higher-order splicing regulation.

RESULTS

We present tools to visualize complex splicing patterns in full-length cDNA libraries. Developmental changes in pair-wise correlations are presented vectorially in 'clock plots' and linkage grids. Higher-order correlations are assessed statistically through Monte Carlo analysis of a log-linear model with an empirical-Bayes estimate of the true probabilities of observed and unobserved splice forms. Log-linear coefficients are visualized in a 'spliceprint,' a signature of splice correlations in the transcriptome. We present two novel metrics: the linkage change index, which measures the directional change in pair-wise correlation with tissue differentiation, and the accuracy index, a very simple goodness-of-fit metric that is more sensitive than the integrated squared error when applied to sparsely populated tables, and unlike chi-square, does not diverge at low variance. Considerable attention is given to sparse contingency tables, which are inherent to single-gene libraries.

CONCLUSION

Patterns of splicing correlations are revealed, which span a broad range of interaction order and change in development. The methods have a broad scope of applicability, beyond the single gene--including, for example, multiple gene interactions in the complete transcriptome.

Detection of minor drug-resistant populations by parallel allele-specific sequencing

We developed a highly sensitive parallel allele-specific sequencing (PASS) assay to simultaneously analyze a large number of viral genomes and detect minor drug-resistant populations at approximately 0.1-0.01% levels. Using this assay on samples from individuals infected with human immunodeficiency viruses (HIV), we successfully detected and quantified minor populations of drug-resistant viruses and performed linkage analysis of multiple-drug resistance mutations. This assay may serve as a useful tool to study drug resistance in HIV and other infectious agents.

Probe selection and expression index computation of Affymetrix Exon Arrays

Background

There is great current interest in developing microarray platforms for measuring mRNA abundance at both gene level and exon level. The Affymetrix Exon Array is a new high-density gene expression microarray platform, with over six million probes targeting all annotated and predicted exons in a genome. An important question for the analysis of exon array data is how to compute overall gene expression indexes. Because of the complexity of the design of exon array probes, this problem is different in nature from summarizing gene-level expression from traditional 3′ expression arrays.

Methodology/Principal Findings

In this manuscript, we use exon array data from 11 human tissues to study methods for computing gene-level expression. We showed that for most genes there is a subset of exon array probes having highly correlated intensities across multiple samples. We suggest that these probes could be used as reliable indicators of overall gene expression levels. We developed a probe selection algorithm to select such a subset of highly correlated probes for each gene, and computed gene expression indexes using the selected probes.

Conclusions/Significance

Our results demonstrate that probe selection improves gene expression estimates from exon arrays. The selected probes can be used in future analyses of other exon array datasets to compute gene expression indexes.

SNPSplicer: systematic analysis of SNP-dependent splicing in genotyped cDNAs

Functional annotation of SNPs (as generated by HapMap (http://www.hapmap.org) for instance) is a major challenge. SNPs that lead to single amino acid substitutions, stop codons, or frameshift mutations can be readily interpreted, but these represent only a fraction of known SNPs. Many SNPs are located in sequences of splicing relevance-the canonical splice site consensus sequences, exonic and intronic splice enhancers or silencers (exonic splice enhancer [ESE], intronic splice enhancer [ISE], exonic splicing silencer [ESS], and intronic splicing silencer [ISS]), and others. We propose using sets of matching DNA and complementary DNA (cDNA) as a screening method to investigate the potential splice effects of SNPs in RT-PCR experiments with tissue material from genotyped sources. We have developed a software solution (SNPSplicer; http://www.ikmb.uni-kiel.de/snpsplicer) that aids in the rapid interpretation of such screening experiments. The utility of the approach is illustrated for SNPs affecting the donor splice sites (rs2076530:A>G, rs3816989:G>A) leading to the use of a cryptic splice site and exon skipping, respectively, and an exonic splice enhancer SNP (rs2274987:C/T), leading to inclusion of a new exon. We anticipate that this methodology may help in the functional annotation of SNPs in a more high-throughput fashion.

Profiling the array of Ca(v)3.1 variants from the human T-type calcium channel gene CACNA1G: alternative structures, developmental expression, and biophysical variations

We describe the regulated transcriptome of CACNA1G, a human gene for T-type Ca(v)3.1 calcium channels that is subject to extensive alternative RNA splicing. Fifteen sites of transcript variation include 2 alternative 5'-UTR promoter sites, 2 alternative 3'-UTR polyadenylation sites, and 11 sites of alternative splicing within the open reading frame. A survey of 1580 fetal and adult human brain full-length complementary DNAs reveals a family of 30 distinct transcripts, including multiple functional forms that vary in expression with development. Statistical analyses of fetal and adult transcript populations reveal patterns of linkages among intramolecular splice site configurations that change dramatically with development. A shift from nearly independent, biased splicing in fetal transcripts to strongly concerted splicing in adult transcripts suggests progressive activation of multiple "programs" of splicing regulation that reorganize molecular structures in differentiating cells. Patch-clamp studies of nine selected variants help relate splicing regulation to permutations of the gating parameters most likely to modify T-channel physiology in expressing neurons. Gating behavior reflects combinatorial interactions between variable domains so that molecular phenotype depends on ensembles of coselected domains, consistent with the observed emergence of concerted splicing during development. We conclude that the structural gene and networks of splicing regulatory factors define an integrated system for the phenotypic variation of Ca(v)3.1 biophysics during nervous system development.

An expectation-maximization algorithm for probabilistic reconstructions of full-length isoforms from splice graphs

Reconstructing full-length transcript isoforms from sequence fragments (such as ESTs) is a major interest and challenge for bioinformatic analysis of pre-mRNA alternative splicing. This problem has been formulated as finding traversals across the splice graph, which is a directed acyclic graph (DAG) representation of gene structure and alternative splicing. In this manuscript we introduce a probabilistic formulation of the isoform reconstruction problem, and provide an expectation-maximization (EM) algorithm for its maximum likelihood solution. Using a series of simulated data and expressed sequences from real human genes, we demonstrate that our EM algorithm can correctly handle various situations of fragmentation and coupling in the input data. Our work establishes a general probabilistic framework for splice graph-based reconstructions of full-length isoforms.

Regulation of CD44 alternative splicing by SRm160 and its potential role in tumor cell invasion

The multiple isoforms of the transmembrane glycoprotein CD44 are produced by alternative RNA splicing. Expression of CD44 isoforms containing variable 5 exon (v5) correlates with enhanced malignancy and invasiveness of some tumors. Here we demonstrate that SRm160, a splicing coactivator, regulates CD44 alternative splicing in a Ras-dependent manner. Overexpression of SRm160 stimulates inclusion of CD44 v5 when Ras is activated. Conversely, small interfering RNA (siRNA)-mediated silencing of SRm160 significantly reduces v5 inclusion. Immunoprecipitation shows association of SRm160 with Sam68, a protein that also stimulates v5 inclusion in a Ras-dependent manner, suggesting that these two proteins interact to regulate CD44 splicing. Importantly, siRNA-mediated depletion of CD44 v5 decreases tumor cell invasion. Reduction of SRm160 by siRNA transfection downregulates the endogenous levels of CD44 isoforms, including v5, and correlates with a decrease in tumor cell invasiveness.

Long-range polony haplotyping of individual human chromosome molecules

We report a method for multilocus long-range haplotyping on human chromosome molecules in vitro based on the DNA polymerase colony (polony) technology. By immobilizing thousands of intact chromosome molecules within a polyacrylamide gel on a microscope slide and performing multiple amplifications from single molecules, we determined long-range haplotypes spanning a 153-Mb region of human chromosome 7 and found evidence of rare mitotic recombination events in human lymphocytes. Furthermore, the parallel nature of DNA polony technology allows efficient haplotyping on pooled DNAs from a population on one slide, with a throughput three orders of magnitudes higher than current molecular haplotyping methods. Linkage disequilibrium statistics established by our pooled DNA haplotyping method are more accurate than statistically inferred haplotypes. This haplotyping method is well suited for candidate gene-based association studies as well as for investigating the pattern of recombination in mammalian cells.

Misregulation of tau alternative splicing in neurodegeneration and dementia

Tau is a microtubule-associated protein that fulfills several functions critical for neuronal formation and health. Tau discharges its functions by producing multiple isoforms via intricately regulated alternative splicing. These isoforms modulate tau function in normal brain by altering the domains of the protein, thereby influencing its conformation and post-translational modifications and hence its affinity for microtubules and other ligands. Disturbances in tau expression result in disruption of the neuronal cytoskeleton and formation of pathological tau structures (neurofibrillary tangles) found in brains of dementia sufferers. More specifically, aberrations in tau splicing regulation directly cause several neurodegenerative diseases that lead to dementia. This review briefly presents our cumulative knowledge of tau splicing regulation in connection with the alterations in tau splicing seen in neurodegeneration.

Detection of alternatively spliced transcripts in leukemia cell lines by minisequencing on microarrays

BACKGROUND

Recent genome-wide expression studies suggest that approximately 80% of the 25,000 human genes undergo alternative splicing. Alternative splicing may be associated with human diseases, particularly with cancer, but the molecular disease mechanisms are poorly understood. Convenient, novel methods for multiplexed detection of alternatively spliced transcripts are needed.

METHODS

We devised a new approach for detecting splice variants based on a tag-microarray minisequencing system, originally developed for genotyping single-nucleotide polymorphisms. We established the system for multiplexed detection of 61 alternatively spliced transcripts in a panel of 19 cancer-related genes and used it to dissect the splicing patterns in cancer and endothelial cells.

RESULTS

Our microarray system detected 82% of the splice variants screened for, including both simple and complex splice variants, in at least 1 of the leukemia cell types analyzed. The intraassay CV values for our method ranged from 0.01 to 0.34 (mean, 0.13) for 5 replicate measurements. Our system allowed semiquantitative comparison of the splicing patterns between the cell lines. Similar, but not identical, patterns of alternative splicing were observed among the leukemia cell lines. Size analysis of the PCR products subjected to the tag-array minisequencing system and real-time PCR with exon-junction probes verified the results from the microarray system.

CONCLUSIONS

The microarray-based method is a robust and easily accessible tool for parallel detection of alternatively spliced transcripts of multiple genes. It can be used for studying alternative splicing in cancer progression and for following up drug treatment, and it may be a useful tool in clinical diagnostics for cancer and other disorders.

A polar mechanism coordinates different regions of alternative splicing within a single gene

Alternative splicing plays a key role in generating protein diversity. Transfections with minigenes revealed coordination between two distant, alternatively spliced exons in the same gene. Mutations that either inhibit or stimulate inclusion of the upstream alternative exon deeply affect inclusion of the downstream one. However, similar mutations at the downstream alternative exon have little effect on the upstream one. This polar effect is promoter specific and is enhanced by inhibition of transcriptional elongation. Consistently, cells from mutant mice with either constitutive or null inclusion of a fibronectin alternative exon revealed coordination with a second alternative splicing region, located far downstream. Using allele-specific RT-PCR, we demonstrate that this coordination occurs in cis and is also affected by transcriptional elongation rates. Bioinformatics supports the generality of these findings, indicating that 25% of human genes contain multiple alternative splicing regions and identifying several genes with nonrandom distribution of mRNA isoforms at two alternative regions.

Expressible molecular colonies

Carrying out polymerase chain reaction in a gel layer generates a 2-D pattern of DNA colonies comprising pure genetic clones. Here we demonstrate that transcription, translation and protein folding can be performed in the same gel. The resulting nucleoprotein colonies mimic living cells by serving as compartments in which the synthesized RNAs and proteins co-localize with their templates. Yet, due to the absence of penetration barriers, such a molecular colony display allows cloned genes to be directly tested for the encoded functions. Now, the results imply that virtually any manipulations with genes and their expression products can be accomplished in vitro.