Gene structure prediction and alternative splicing analysis using genomically aligned ESTs

Spliceosomal complex components are critical for adjusting the C:N balance during high-light acclimation

Plant acclimation to an ever-changing environment is decisive for growth, reproduction, and survival. Light availability limits biomass production on both ends of the intensity spectrum. Therefore, the adjustment of plant metabolism is central to high-light (HL) acclimation, and the accumulation of photoprotective anthocyanins is commonly observed. However, mechanisms and factors regulating the HL acclimation response are less clear. Two Arabidopsis mutants of spliceosome components exhibiting a pronounced anthocyanin overaccumulation in HL were isolated from a forward genetic screen for new factors crucial for plant acclimation. Time-resolved physiological, transcriptome, and metabolome analysis revealed a vital function of the spliceosome components for rapidly adjusting gene expression and metabolism. Deficiency of INCREASED LEVEL OF POLYPLOIDY1 (ILP1), NTC-RELATED PROTEIN1 (NTR1), and PLEIOTROPIC REGULATORY LOCUS1 (PRL1) resulted in a marked overaccumulation of carbohydrates and strongly diminished amino acid biosynthesis in HL. While not generally limited in N-assimilation, ilp1, ntr1, and prl1 showed higher glutamate levels and reduced amino acid biosynthesis in HL. The comprehensive analysis reveals a function of the spliceosome components in the conditional regulation of the carbon:nitrogen balance and the accumulation of anthocyanins during HL acclimation. The importance of gene expression, metabolic regulation, and re-direction of carbon towards anthocyanin biosynthesis for HL acclimation are discussed.

Small but mighty: microexons in glucose homeostasis

Many molecular mechanisms underlying blood glucose homeostasis remain elusive. Juan-Mateu et al. find that pancreatic islet cells utilize a regulatory program, originally identified in neurons, that involves alternative splicing of microexons in genes important for insulin secretion or diabetes risk.

Long non-coding RNAs (lncRNAs) in spermatogenesis and male infertility

BACKGROUND

Long non-coding RNAs (lncRNAs) have a size of more than 200 bp and are known to regulate a host of crucial cellular processes like proliferation, differentiation and apoptosis by regulating gene expression. While small noncoding RNAs (ncRNAs) such as miRNAs, siRNAs, Piwi-interacting RNAs have been extensively studied in male germ cell development, the role of lncRNAs in spermatogenesis remains largely unknown.

OBJECTIVE

In this article, we have reviewed the biology and role of lncRNAs in spermatogenesis along with the tools available for data analysis.

RESULTS AND CONCLUSIONS

Till date, three microarray and four RNA-seq studies have been undertaken to identify lncRNAs in mouse testes or germ cells. These studies were done on pre-natal, post-natal, adult testis, and different germ cells to identify lncRNAs regulating spermatogenesis. In case of humans, five RNA-seq studies on different germ cell populations, including two on sperm, were undertaken. We compared three studies on human germ cells to identify common lncRNAs and found 15 lncRNAs (LINC00635, LINC00521, LINC00174, LINC00654, LINC00710, LINC00226, LINC00326, LINC00494, LINC00535, LINC00616, LINC00662, LINC00668, LINC00467, LINC00608, and LINC00658) to show consistent differential expression across these studies. Some of the targets of these lncRNAs included CENPB, FAM98B, GOLGA6 family, RPGR, TPM2, GNB5, KCNQ10T1, TAZ, LIN28A, CDKN2B, CDKN2A, CDKN1A, CDKN1B, CDKN1C, EZH2, SUZ12, VEGFA genes. A lone study on human male infertility identified 9879 differentially expressed lncRNAs with three (lnc32058, lnc09522, and lnc98497) of them showing specific and high expression in immotile sperm in comparison to normal motile sperm. A few lncRNAs (Mrhl, Drm, Spga-lncRNAs, NLC1-C, HongrES2, Tsx, LncRNA-tcam1, Tug1, Tesra, AK015322, Gm2044, and LncRNA033862) have been functionally validated for their roles in spermatogenesis. Apart from rodents and humans, studies on sheep and bull have also identified lncRNAs potentially important for spermatogenesis. A number of these non-coding RNAs are strong candidates for further research on their roles in spermatogenesis.

A Modifying Autoantigen in Graves' Disease

The TSH receptor (TSHR) is the major autoantigen in Graves' disease (GD). Bioinformatic analyses predict the existence of several human TSHR isoforms from alternative splicing, which can lead to the coexpression of multiple receptor forms. The most abundant of these is TSHRv1.3. In silico modeling of TSHRv1.3 demonstrated the structural integrity of this truncated receptor isoform and its potential binding of TSH. Tissue profiling revealed wide expression of TSHRv1.3, with a predominant presence in thyroid, bone marrow, thymus, and adipose tissue. To gain insight into the role of this v1.3 receptor isoform in thyroid pathophysiology, we cloned the entire open reading frame into a mammalian expression vector. Immunoprecipitation studies demonstrated that both TSHR-stimulating antibody and human TSH could bind v1.3. Furthermore, TSHRv1.3 inhibited the stimulatory effect of TSH and TSHR-Ab MS-1 antibody on TSHR-induced cAMP generation in a dose-dependent manner. To confirm the antigenicity of v1.3, we used a peptide ELISA against two different epitopes. Of 13 GD samples, 11 (84.6%) were positive for a carboxy terminal peptide and 10 (76.9%) were positive with a junction region peptide. To demonstrate that intracellular v1.3 could serve as an autoantigen and modulate disease, we used double-transfected Chinese hamster ovary cells that expressed both green fluorescent protein (GFP)-tagged TSHRv1.3 and full-length TSHR. We then induced cell stress and apoptosis using a TSHR monoclonal antibody and observed the culture supernatant contained v1.3-GFP protein, demonstrating the release of the intracellular receptor variant by this mechanism.

SplicedFamAlign: CDS-to-gene spliced alignment and identification of transcript orthology groups

BACKGROUND

The inference of splicing orthology relationships between gene transcripts is a basic step for the prediction of transcripts and the annotation of gene structures in genomes. The splicing structure of a sequence refers to the exon extremity information in a CDS or the exon-intron extremity information in a gene sequence. Splicing orthologous CDS are pairs of CDS with similar sequences and conserved splicing structures from orthologous genes. Spliced alignment that consists in aligning a spliced cDNA sequence against an unspliced genomic sequence, constitutes a promising, yet unexplored approach for the identification of splicing orthology relationships. Existing spliced alignment algorithms do not exploit the information on the splicing structure of the input sequences, namely the exon structure of the cDNA sequence and the exon-intron structure of the genomic sequences. Yet, this information is often available for coding DNA sequences (CDS) and gene sequences annotated in databases, and it can help improve the accuracy of the computed spliced alignments. To address this issue, we introduce a new spliced alignment problem and a method called SplicedFamAlign (SFA) for computing the alignment of a spliced CDS against a gene sequence while accounting for the splicing structures of the input sequences, and then the inference of transcript splicing orthology groups in a gene family based on spliced alignments.

RESULTS

The experimental results show that SFA outperforms existing spliced alignment methods in terms of accuracy and execution time for CDS-to-gene alignment. We also show that the performance of SFA remains high for various levels of sequence similarity between input sequences, thanks to accounting for the splicing structure of the input sequences. It is important to notice that unlike all current spliced alignment methods that are meant for cDNA-to-genome alignments and can be used for CDS-to-gene alignments, SFA is the first method specifically designed for CDS-to-gene alignments.

CONCLUSION

We show the usefulness of SFA for the comparison of genes and transcripts within a gene family for the purpose of analyzing splicing orthologies. It can also be used for gene structure annotation and alternative splicing analyses. SplicedFamAlign was implemented in Python. Source code is freely available at https://github.com/UdeS-CoBIUS/SpliceFamAlign .

Long Non-Coding RNAs in the Regulation of Gene Expression: Physiology and Disease

The identification of RNAs that are not translated into proteins was an important breakthrough, defining the diversity of molecules involved in eukaryotic regulation of gene expression. These non-coding RNAs can be divided into two main classes according to their length: short non-coding RNAs, such as microRNAs (miRNAs), and long non-coding RNAs (lncRNAs). The lncRNAs in association with other molecules can coordinate several physiological processes and their dysfunction may impact in several pathologies, including cancer and infectious diseases. They can control the flux of genetic information, such as chromosome structure modulation, transcription, splicing, messenger RNA (mRNA) stability, mRNA availability, and post-translational modifications. Long non-coding RNAs present interaction domains for DNA, mRNAs, miRNAs, and proteins, depending on both sequence and secondary structure. The advent of new generation sequencing has provided evidences of putative lncRNAs existence; however, the analysis of transcriptomes for their functional characterization remains a challenge. Here, we review some important aspects of lncRNA biology, focusing on their role as regulatory elements in gene expression modulation during physiological and disease processes, with implications in host and pathogens physiology, and their role in immune response modulation.

ChemGenome2.1: An Ab Initio Gene Prediction Software

Gene prediction, also known as gene identification, gene finding, gene recognition, or gene discovery, is among one of the important problems of molecular biology and is receiving increasing attention due to the advent of large-scale genome sequencing projects. We designed an ab initio model (called ChemGenome) for gene prediction in prokaryotic genomes based on physicochemical characteristics of codons. In this chapter, we present the methodology of the latest version of this model ChemGenome2.1 (CG2.1). The first module of the protocol builds a three-dimensional vector from three calculated quantities for each codon-the double-helical trinucleotide base pairing energy, the base pair stacking energy, and an index of the propensity of a codon for protein-nucleic acid interactions. As this three-dimensional vector moves along any genome, the net orientation of the resultant vector should differ significantly for gene and non-genic regions to make a distinction feasible. The predicted putative protein-coding genes from above parameters are passed through a second module of the protocol which reduces the number of false positives by utilizing a filter based on stereochemical properties of protein sequences. The chemical properties of amino acid side chains taken into consideration are the presence of sp3 hybridized γ carbon atom, hydrogen bond donor ability, short/absence of δ carbon and linearity of the side chains/non-occurrence of bi-dentate forks with terminal hydrogen atoms in the side chain. The final prediction of the potential protein-coding genes is based on the frequency of occurrence of amino acids in the predicted protein sequences and their deviation from the frequency values of Swissprot protein sequences, both at monomer and tripeptide levels. The final screening is based on Z-score. Though CG2.1 is a gene finding tool for prokaryotes, considering the underlying similarity in the chemical and physical properties of DNA among prokaryotes and eukaryotes, we attempted to evaluate its applicability for gene finding in the lower eukaryotes. The results give a hope that the concept of gene finding based on physicochemical model of codons is a viable idea for eukaryotes as well, though, undoubtedly, improvements are needed.

A genome wide analysis of alternative splicing events during the osteogenic differentiation of human cartilage endplate-derived stem cells

Low back pain is a prevalent disease, which leads to suffering and disabilities in a vast number of individuals. Degenerative disc diseases are usually the underlying causes of low back pain. However, the pathogenesis of degenerative disc diseases is highly complex and difficult to determine. Current therapies for degenerative disc diseases are various. In particular, cell-based therapies have proven to be effective and promising. Our research group has previously isolated and identified the cartilage endplate‑derived stem cells. In addition, alternative splicing is a sophisticated regulatory mechanism, which greatly increases cellular complexity and phenotypic diversity of eukaryotic organisms. The present study continued to investigate alternative splicing events in osteogenic differentiation of cartilage endplate‑derived stem cells. An Affymetrix Human Transcriptome Array 2.0 was used to detect splicing changes between the control and differentiated samples. Additionally, molecular function and pathway analysis were also performed. Following rigorous bioinformatics analysis of the data, 3,802 alternatively spliced genes were identified, and 10 of these were selected for validation by reverse transcription‑polymerase chain reaction. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway analysis also revealed numerous enriched GO terms and signaling pathways. To the best of our knowledge, the present study is the first to investigate alternative splicing mechanisms in osteogenic differentiation of stem cells on a genome‑wide scale. The illumination of molecular mechanisms of stem cell osteogenic differentiation may assist the development novel bioengineered methods to treat degenerative disc diseases.

Global Gene Expression Profiling and Alternative Splicing Events during the Chondrogenic Differentiation of Human Cartilage Endplate-Derived Stem Cells

Low back pain (LBP) is a very prevalent disease and degenerative disc diseases (DDDs) usually account for the LBP. However, the pathogenesis of DDDs is complicated and difficult to elucidate. Alternative splicing is a sophisticated regulatory process which greatly increases cellular complexity and phenotypic diversity of eukaryotic organisms. In addition, the cartilage endplate-derived stem cells have been discovered and identified by our research group. In this paper, we continue to investigate gene expression profiling and alternative splicing events during chondrogenic differentiation of cartilage endplate-derived stem cells. We adopted Affymetrix Human Transcriptome Array 2.0 (HTA 2.0) to compare the transcriptional and splicing changes between the control and differentiated samples. RT-PCR and quantitative PCR are used to validate the microarray results. The GO and KEGG pathway analysis was also performed. After bioinformatics analysis of the data, we detected 1953 differentially expressed genes. In terms of alternative splicing, the Splicing Index algorithm was used to select alternatively spliced genes. We detected 4411 alternatively spliced genes. GO and KEGG pathway analysis also revealed several functionally involved biological processes and signaling pathways. To our knowledge, this is the first study to investigate the alternative splicing mechanisms in chondrogenic differentiation of stem cells on a genome-wide scale.

A New Panel of SNP Markers for the Individual Identification of North American Pumas

Pumas Puma concolor are one of the most studied terrestrial carnivores because of their widespread distribution, substantial ecological impacts, and conflicts with humans. Over the past decade, managing pumas has involved extensive efforts including the use of genetic methods. Microsatellites have been the most commonly used genetic markers; however, technical artifacts and little overlap of frequently used loci render large-scale comparison of puma genetic data across studies challenging. Therefore, a panel of genetic markers that can produce consistent genotypes across studies without the need for extensive calibrations is essential for range-wide genetic management of puma populations. Here, we describe the development of PumaPlex, a high-throughput assay to genotype 25 single nucleotide polymorphisms in pumas. We validated PumaPlex in 748 North American pumas Puma concolor couguar, and demonstrated its ability to generate reproducible genotypes and accurately identify individuals. Furthermore, in a test using fecal deoxyribonucleic acid (DNA) samples, we found that PumaPlex produced significantly more genotypes with fewer errors than 12 microsatellite loci, 8 of which are commonly used. Our results demonstrate that PumaPlex is a valuable tool for the genetic monitoring and management of North American puma populations. Given the analytical simplicity, reproducibility, and high-throughput capability of single nucleotide polymorphisms, PumaPlex provides a standard panel of markers that promotes the comparison of genotypes across studies and independent of the genotyping technology used.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Physiology and pathophysiology of the renal Na-K-2Cl cotransporter (NKCC2)

The Na-K-2Cl cotransporter (NKCC2; BSC1) is located in the apical membrane of the epithelial cells of the thick ascending limb of the loop of Henle (TAL). NKCC2 facilitates ∼20–25% of the reuptake of the total filtered NaCl load. NKCC2 is therefore one of the transport proteins with the highest overall reabsorptive capacity in the kidney. Consequently, even subtle changes in NKCC2 transport activity considerably alter the renal reabsorptive capacity for NaCl and eventually lead to perturbations of the salt and water homoeostasis. In addition to facilitating the bulk reabsorption of NaCl in the TAL, NKCC2 transport activity in the macula densa cells of the TAL constitutes the initial step of the tubular-vascular communication within the juxtaglomerular apparatus (JGA); this communications allows the TAL to modulate the preglomerular resistance of the afferent arteriole and the renin secretion from the granular cells of the JGA. This review provides an overview of our current knowledge with respect to the general functions of NKCC2, the modulation of its transport activity by different regulatory mechanisms, and new developments in the pathophysiology of NKCC2-dependent renal NaCl transport.

RNA binding proteins: a common denominator of neuronal function and dysfunction

In eukaryotic cells, gene activity is not directly reflected by protein levels because mRNA processing, transport, stability, and translation are co- and post-transcriptionally regulated. These processes, collectively known as the ribonome, are tightly controlled and carried out by a plethora of trans-acting RNA-binding proteins (RBPs) that bind to specific cis elements throughout the RNA sequence. Within the nervous system, the role of RBPs in brain function turns out to be essential due to the architectural complexity of neurons exemplified by a relatively small somal size and an extensive network of projections and connections. Thus far, RBPs have been shown to be indispensable for several aspects of neurogenesis, neurite outgrowth, synapse formation, and plasticity. Consequently, perturbation of their function is central in the etiology of an ever-growing spectrum of neurological diseases, including fragile X syndrome and the neurodegenerative disorders frontotemporal lobar degeneration and amyotrophic lateral sclerosis.

SpliceVista, a tool for splice variant identification and visualization in shotgun proteomics data

Alternative splicing is a pervasive process in eukaryotic organisms. More than 90% of human genes have alternatively spliced products, and aberrant splicing has been shown to be associated with many diseases. Current methods employed in the detection of splice variants include prediction by clustering of expressed sequence tags, exon microarray, and mRNA sequencing, all methods focusing on RNA-level information. There is a lack of tools for analyzing splice variants at the protein level. Here, we present SpliceVista, a tool for splice variant identification and visualization based on mass spectrometry proteomics data. SpliceVista retrieves gene structure and translated sequences from alternative splicing databases and maps MS-identified peptides to splice variants. The visualization module plots the exon composition of each splice variant and aligns identified peptides with transcript positions. If quantitative mass spectrometry data are used, SpliceVista plots the quantitative patterns for each peptide and provides users with the option to cluster peptides based on their quantitative patterns. SpliceVista can identify splice-variant-specific peptides, providing the possibility for variant-specific analysis. The tool was tested on two experimental datasets (PXD000065 and PXD000134). In A431 cells treated with gefitinib, 2983 splice-variant-specific peptides corresponding to 939 splice variants were identified. Through comparison of splice-variant-centric, protein-centric, and gene-centric quantification, several genes (e.g. EIF4H) were found to have differentially regulated splice variants after gefitinib treatment. The same discrepancy between protein-centric and splice-centric quantification was detected in the other dataset, in which induced pluripotent stem cells were compared with parental fibroblast and human embryotic stem cells. In addition, SpliceVista can be used to visualize novel splice variants inferred from peptide-level evidence. In summary, SpliceVista enables visualization, detection, and differential quantification of protein splice variants that are often missed in current proteomics pipelines.

Dietary salt intake modulates differential splicing of the Na-K-2Cl cotransporter NKCC2

Both sodium reabsorption in the thick ascending limb of the loop of Henle (TAL) and macula densa salt sensing crucially depend on the function of the Na/K/2Cl cotransporter NKCC2. The NKCC2 gene gives rise to at least three different full-length NKCC2 isoforms derived from differential splicing. In the present study, we addressed the influence of dietary salt intake on the differential splicing of NKCC2. Mice were subjected to diets with low-salt, standard salt, and high-salt content for 7 days, and NKCC2 isoform mRNA abundance was determined. With decreasing salt intake, we found a reduced abundance of the low-affinity isoform NKCC2A and an increase in the high-affinity isoform NKCC2B in the renal cortex and the outer stripe of the outer medulla. This shift from NKCC2A to NKCC2B during a low-salt diet could be mimicked by furosemide in vivo and in cultured kidney slices. Furthermore, the changes in NKCC2 isoform abundance during a salt-restricted diet were partly mediated by the actions of angiotensin II on AT1 receptors, as determined using chronic angiotensin II infusion. In contrast to changes in oral salt intake, water restriction (48 h) and water loading (8% sucrose solution) increased and suppressed the expression of all NKCC2 isoforms, without changing the distribution pattern of the single isoforms. In summary, the differential splicing of NKCC2 pre-mRNA is modulated by dietary salt intake, which may be mediated by changes in intracellular ion composition. Differential splicing of NKCC2 appears to contribute to the adaptive capacity of the kidney to cope with changes in reabsorptive needs.

Generation and analysis of the expressed sequence tags from the mycelium of Ganoderma lucidum

Ganoderma lucidum (G. lucidum) is a medicinal mushroom renowned in East Asia for its potential biological effects. To enable a systematic exploration of the genes associated with the various phenotypes of the fungus, the genome consortium of G. lucidum has carried out an expressed sequence tag (EST) sequencing project. Using a Sanger sequencing based approach, 47,285 ESTs were obtained from in vitro cultures of G. lucidum mycelium of various durations. These ESTs were further clustered and merged into 7,774 non-redundant expressed loci. The features of these expressed contigs were explored in terms of over-representation, alternative splicing, and natural antisense transcripts. Our results provide an invaluable information resource for exploring the G. lucidum transcriptome and its regulation. Many cases of the genes over-represented in fast-growing dikaryotic mycelium are closely related to growth, such as cell wall and bioactive compound synthesis. In addition, the EST-genome alignments containing putative cassette exons and retained introns were manually curated and then used to make inferences about the predominating splice-site recognition mechanism of G. lucidum. Moreover, a number of putative antisense transcripts have been pinpointed, from which we noticed that two cases are likely to reveal hitherto undiscovered biological pathways. To allow users to access the data and the initial analysis of the results of this project, a dedicated web site has been created at http://csb2.ym.edu.tw/est/.

Exon array analysis of alternative splicing of genes in SOD1G93A transgenic mice

Alternative splicing is a common strategy for creating functional diversities of proteins. While conventional identification of splice variants generally targets individual genes in amyotrophic lateral sclerosis, we present a novel exon-centric array that allows genome-wide identification of splice variants and concurrently provides analysis of gene expression. Compare 1 was asymptomatic SOD1G93A transgenic mice with nontransgenic littermates; compare 2 was symptomatic with asymptomatic transgenic mice. RT-PCR was performed to validate. Pathway and GO analysis were performed on abnormal genes. These findings could guide us to demonstrated the potential influence of mutant human CuZn-SOD1 and of splicing regulation in pathological processes.

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

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

High-throughput functional microRNA profiling using recombinant AAV-based microRNA sensor arrays

There is a lack of methods for high-throughput functional microRNA (miRNA) profiling. In this chapter, we describe a recombinant adeno-associated virus-based miRNA sensor array (miRNA Asensor array), which is able to profile functional miRNAs in cultured cells. The preparation of an miRNA Asensor array and its usage are discussed.

In silico prediction of mRNA poly(A) sites in Chlamydomonas reinhardtii

Accurately predicting polyadenylation [poly(A)] sites is important for defining the end of genes and understanding gene regulation mechanisms. Alternative polyadenylation (APA) has been demonstrated to play an important role in transcriptome diversity and regulating gene expression. To accurately predict poly(A) and APA sites in Chlamydomonas reinhardtii, a green alga that can be used to produce renewable energy, we proposed a novel model that integrated five methods for representing the features of these sites with a combined classifier. We presented a new grouping method based on pattern assembly to classify the poly(A) sites into four groups. We used five methods, involving the predicted RNA secondary structure, the term frequency-inverse document frequency weight, first-order Markov chain, pentamer ratio and a position weight matrix, to generate the feature space. We then developed a heuristic method to form the combined classifier by weighting multiple classifiers to predict poly(A) sites in each group. The high specificity and sensitivity of this model were demonstrated by testing the four groups of poly(A) sites and the intronic APA sites. The average prediction performance was approximately 8 % higher than the performance of a previous prediction model. For the group without any conserved patterns, the prediction accuracy was 9 % higher than for the accuracy with the previous technique. However, the prediction efficiency of this group was still significantly lower than that of the other groups, indicating the importance of identifying additional signal patterns for poly(A) site prediction. We also predicted the alternative poly(A) sites in introns with good accuracy. This prediction model was designed to be easily expanded with new classifiers or new features. Therefore, this model is applicable to new data or other species. Our model will be useful both in genome annotation because it predicts the end of a mature transcript and in genetic engineering because it enables researchers to eliminate undesirable poly(A) sites.

Identification of novel splice variants and exons of human endothelial cell-specific chemotaxic regulator (ECSCR) by bioinformatics analysis

Recent discovery of biological function of endothelial cell-specific chemotaxic regulator (ECSCR), previously known as endothelial cell-specific molecule 2 (ECSM2), in modulating endothelial cell migration, apoptosis, and angiogenesis, has made it an attractive molecule in vascular research. Thus, identification of splice variants of ECSCR could provide new strategies for better understanding its roles in health and disease. In this study, we performed a series of blast searches on the human EST database with known ECSCR cDNA sequence (Variant 1), and identified additional three splice variants (Variants 2-4). When examining the ECSCR gene in the human genome assemblies, we found a large unknown region between Exons 9 and 11. By PCR amplification and sequencing, we partially mapped Exon 10 within this previously unknown region of the ECSCR gene. Taken together, in addition to previously reported human ECSCR, we identified three novel full-length splice variants potentially encoding different protein isoforms. We further defined a total of twelve exons and nearly all exon-intron boundaries of the gene, of which only eight are annotated in current public databases. Our work provides new information on gene structure and alternative splicing of the human ECSCR, which may imply its functional complexity. This undoubtedly opens new opportunities for future investigation of the biological and pathological significance of these ECSCR splice variants.

Revisit on the evolutionary relationship between alternative splicing and gene duplication

Gene duplications and alternative splicing (AS) isoforms are two widespread types of genetic variations that can facilitate diversification of protein function. A number of studies claimed that after gene duplication, two AS isoforms with differential functions can be 'fixed', respectively, in each of the duplicate copies. This simple 'functional-sharing' hypothesis was recently challenged by Roux and Robinson-Rechavi (2011). Instead, they proposed a more sophisticated hypothesis, invoking that less alternative splicing genes tend to be duplicated more frequently, and single-copy genes are younger than duplicate genes, or the 'duplicability-age' hypothesis for short. In this letter, we show that all these genome-wide analyses of AS isoforms actually did not provide clear-cut evidence to nullify the basic idea of functional-sharing hypothesis. After updating our understanding of genome-wide alternative splicing, duplicability and CNV (copy number variation), we argue that the foundation of the duplicability-age hypothesis remains to be justified carefully. Finally, we suggest that a better approach to resolving this controversy is the correspondence analysis of indels (insertions and deletions) between duplicate genes to the genomic exon-intron structure, which can be used to experimentally test the effect of functional-sharing hypothesis.

mRNA transcript diversity creates new opportunities for pharmacological intervention

Most protein coding genes generate multiple RNA transcripts through alternative splicing, variable 3' and 5'UTRs, and RNA editing. Although drug design typically targets the main transcript, alternative transcripts can have profound physiological effects, encoding proteins with distinct functions or regulatory properties. Formation of these alternative transcripts is tissue-selective and context-dependent, creating opportunities for more effective and targeted therapies with reduced adverse effects. Moreover, genetic variation can tilt the balance of alternative versus constitutive transcripts or generate aberrant transcripts that contribute to disease risk. In addition, environmental factors and drugs modulate RNA splicing, affording new opportunities for the treatment of splicing disorders. For example, therapies targeting specific mRNA transcripts with splice-site-directed oligonucleotides that correct aberrant splicing are already in clinical trials for genetic disorders such as Duchenne muscular dystrophy. High-throughput sequencing technologies facilitate discovery of novel RNA transcripts and protein isoforms, applications ranging from neuromuscular disorders to cancer. Consideration of a gene's transcript diversity should become an integral part of drug design, development, and therapy.

PIntron: a fast method for detecting the gene structure due to alternative splicing via maximal pairings of a pattern and a text

BACKGROUND

A challenging issue in designing computational methods for predicting the gene structure into exons and introns from a cluster of transcript (EST, mRNA) sequences, is guaranteeing accuracy as well as efficiency in time and space, when large clusters of more than 20,000 ESTs and genes longer than 1 Mb are processed. Traditionally, the problem has been faced by combining different tools, not specifically designed for this task.

RESULTS

We propose a fast method based on ad hoc procedures for solving the problem. Our method combines two ideas: a novel algorithm of proved small time complexity for computing spliced alignments of a transcript against a genome, and an efficient algorithm that exploits the inherent redundancy of information in a cluster of transcripts to select, among all possible factorizations of EST sequences, those allowing to infer splice site junctions that are largely confirmed by the input data. The EST alignment procedure is based on the construction of maximal embeddings, that are sequences obtained from paths of a graph structure, called embedding graph, whose vertices are the maximal pairings of a genomic sequence T and an EST P. The procedure runs in time linear in the length of P and T and in the size of the output.The method was implemented into the PIntron package. PIntron requires as input a genomic sequence or region and a set of EST and/or mRNA sequences. Besides the prediction of the full-length transcript isoforms potentially expressed by the gene, the PIntron package includes a module for the CDS annotation of the predicted transcripts.

CONCLUSIONS

PIntron, the software tool implementing our methodology, is available at http://www.algolab.eu/PIntron under GNU AGPL. PIntron has been shown to outperform state-of-the-art methods, and to quickly process some critical genes. At the same time, PIntron exhibits high accuracy (sensitivity and specificity) when benchmarked with ENCODE annotations.

BEAT: Bioinformatics Exon Array Tool to store, analyze and visualize Affymetrix GeneChip Human Exon Array data from disease experiments

BACKGROUND

It is known from recent studies that more than 90% of human multi-exon genes are subject to Alternative Splicing (AS), a key molecular mechanism in which multiple transcripts may be generated from a single gene. It is widely recognized that a breakdown in AS mechanisms plays an important role in cellular differentiation and pathologies. Polymerase Chain Reactions, microarrays and sequencing technologies have been applied to the study of transcript diversity arising from alternative expression. Last generation Affymetrix GeneChip Human Exon 1.0 ST Arrays offer a more detailed view of the gene expression profile providing information on the AS patterns. The exon array technology, with more than five million data points, can detect approximately one million exons, and it allows performing analyses at both gene and exon level. In this paper we describe BEAT, an integrated user-friendly bioinformatics framework to store, analyze and visualize exon arrays datasets. It combines a data warehouse approach with some rigorous statistical methods for assessing the AS of genes involved in diseases. Meta statistics are proposed as a novel approach to explore the analysis results. BEAT is available at http://beat.ba.itb.cnr.it.

RESULTS

BEAT is a web tool which allows uploading and analyzing exon array datasets using standard statistical methods and an easy-to-use graphical web front-end. BEAT has been tested on a dataset with 173 samples and tuned using new datasets of exon array experiments from 28 colorectal cancer and 26 renal cell cancer samples produced at the Medical Genetics Unit of IRCCS Casa Sollievo della Sofferenza.To highlight all possible AS events, alternative names, accession Ids, Gene Ontology terms and biochemical pathways annotations are integrated with exon and gene level expression plots. The user can customize the results choosing custom thresholds for the statistical parameters and exploiting the available clinical data of the samples for a multivariate AS analysis.

CONCLUSIONS

Despite exon array chips being widely used for transcriptomics studies, there is a lack of analysis tools offering advanced statistical features and requiring no programming knowledge. BEAT provides a user-friendly platform for a comprehensive study of AS events in human diseases, displaying the analysis results with easily interpretable and interactive tables and graphics.

Heterogeneous nuclear ribonucleoprotein hrp36 acts as an alternative splicing repressor in Litopenaeus vannamei Dscam

Heterogeneous nuclear ribonucleoproteins (hnRNPs) are highly conserved from nematode to mammal because they play an important role in several aspects of pre-mRNA maturation, including RNA packaging and alternative splicing. In Drosophila, the hnRNP A1 homolog hrp36 regulates alternative splicing in several genes, including the Down syndrome cell adhesion molecule (Dscam), which produces tens of thousands of isoforms from one locus. In this study, the first hrp36 gene was identified and characterized from Litopenaeus vannamei (Lvhrp36). Its open reading frame (ORF) contains 1101 bp encoding 366 amino acids. The deduced Lvhrp36 protein includes two copies of the RNA recognition motif (RRM), a C-terminal glycine-rich domain (GRD), the highly degenerate RNP consensus sequences RNP-1 and RNP-2, and two RGG boxes. Tissue tropism analysis indicated that Lvhrp36 is expressed ubiquitously and at high levels in most tissues. dsRNA silencing of shrimp Lvhrp36 in vivo induced abnormal exon inclusions in LvDscam, especially in the Ig3 variable region. In the Ig3 region, a total of 14 different combinations were arranged in three different types of abnormal inclusion pattern. This compares to a single combination (one abnormal pattern) in Ig2 and three different combinations (one abnormal pattern) in Ig7. This is the first evidence to suggest that hrp36 may be involved in the regulation of the Ig7 variable region, and it is noteworthy because, at least in Drosophila, there are no hrp36 binding sequences in the Ig7 exon cluster. The above aberrant events were not observed in all of the Lvhrp36-silenced shrimp, and even when they occurred, the normal patterns of inclusion were far more common. We hypothesize that this continued prevalence of normal inclusions was probably due to other unsilenced proteins that were able to rescue Lvhrp36's functionality. Taken together, our results suggest that Lvhrp36 acts as a splicing repressor that regulates alternative splicing events in the Ig2, Ig3 and Ig7 variable regions of shrimp L. vannamei Dscam.

Alternative splicing in self-renewal of embryonic stem cells

Much of embryonic stem cell biology has focused on transcriptional expression and regulation of genes that could mediate its unique potential in self-renewal or pluripotency. In alignment with our present understanding on the genetic, protein, and epigenetic factors that may direct cell fate, we present a short overview of the often overlooked contribution of alternative splice variants to regulatory diversity. Progressing beyond the limitations of a fixed genomic sequence, alternative splicing offers an additional layer of complexity to produce protein variants that may differ in function and localization that can direct embryonic stem cells to specific differentiation pathways. In light of the number of variants that can be produced at key ES cell genes alone, it is challenging to consider how much more multifaceted transcriptional regulation truly is, and if this can be captured more fully in future works.

Differential expression of novel adiponectin receptor-1 transcripts in skeletal muscle of subjects with normal glucose tolerance and type 2 diabetes

OBJECTIVE

Adiponectin receptor-1 (AdipoR1) expression in skeletal muscle has been suggested to play an important role in insulin resistance and diabetes. We aimed at evaluating the presence of novel AdiopR1 splice variants in human muscle and their regulation under physiological and pathophysiological states.

RESEARCH DESIGN AND METHODS

AdipoR1 5'UTR mRNA transcripts, predicted from bioinformatics data, were evaluated in fetal and adult human tissues. Expression and function of the identified transcripts were assessed in cultured human skeletal muscle cells and in muscle biopsies obtained from individuals with normal glucose tolerance (NGT) and type 2 diabetes (n = 49).

RESULTS

Screening of potential AdipoR1 5'UTR splice variants revealed a novel highly abundant muscle transcript (R1T3) in addition to the previously described transcript (R1T1). Unlike R1T1, R1T3 expression was significantly increased during fetal development and myogenesis, paralleled with increased AdipoR1 protein expression. The 5'UTR of R1T3 was found to contain upstream open reading frames that repress translation of downstream coding sequences. Conversely, AdipoR1 3'UTR was associated with enhanced translation efficiency during myoblast-myotube differentiation. A marked reduction in muscle expression of R1T3, R1T1, and R1T3-to-R1T1 ratio was observed in individuals with type 2 diabetes compared with expression levels of NGT subjects, paralleled with decreased expression of the differentiation marker myogenin. Among NGT subjects, R1T3 expression was positively correlated with insulin sensitivity.

CONCLUSIONS

These results indicate that AdipoR1 receptor expression in human skeletal muscle is subjected to posttranscriptional regulation, including alternative splicing and translational control. These mechanisms play an important role during myogenesis and may be important for whole-body insulin sensitivity.

Novel RBPJ transcripts identified in human amniotic fluid cells

The NOTCH signaling pathway plays important roles in stem cell maintenance, cell-fate determination and differentiation during development. Following ligand binding, the cleaved NOTCH intracellular domain (NICD) interacts directly with the recombinant signal binding protein for immunoglobulin kappa J region (RBPJ) transcription factor and the resulting complex targets gene expression in the nucleus. To date, four human RBPJ isoforms have been described in Entrez Gene, varying in the first 5'coding exons. Using an improved protocol, we were able to further identify all four known and five novel RBPJ transcript variants in human amniotic fluid (AF) cells, a cell type known for its stem cell characteristics. In addition, we used human embryonal carcinoma (EC) NTera2/D1 (NT2) cells and NT2-derived neuron and astrocytes to compare the expression pattern of RBPJ transcripts. Further examination of RBPJ transcripts showed that the novel splice variants contain open reading frames in-frame with the known isoforms, suggesting that they can putatively generate similar function proteins. All known and novel RBPJ transcripts contain the putative nuclear localization signal (NLS), an important component of RBPJ-mediated gene regulation.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background

Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression.

Results

The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system.

Conclusions

In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

Alternative splicing of the antitrypsin gene in the silkworm, Bombyx mori

Alternative splicing plays an important role in expanding protein diversity. In the present study, different splice variants of the antitrypsin gene (sw-AT) in the silkworm were identified by bioinformatics analyses using expressed sequence tags and genomic information. Four splice variants were obtained by RT-PCR with suitably designed primers, confirmed by sequencing, and designated as sw-AT-1, sw-AT-2, sw-AT-3, and sw-AT-4. The sw-AT gene contains 10 exons and nine introns. The splice variants differ in exon 9, with sw-AT-1, sw-AT-2, and sw-AT-3 using different versions of the exon, namely exon 9a, 9b, and 9c, respectively. In sw-AT-4, exon 9 consists of the combination of exons 9b and 9c. The expression patterns of the four isoforms in different tissues, at different developmental stages, and under different stress conditions (temperature, starvation, and mycotic infection) were characterized and quantified. The sw-AT isoforms showed tissue-specific expression patterns, with sw-AT-1 present in almost all tissues and sw-AT-4 found in only a few tissues. The four isoforms were predominantly expressed in the fat body, body wall, and testes of larvae, and exhibited similar expression profiles during development of the fat body. Among the stress treatments, low temperature had the greatest effect on isoform expression, and expression was also upregulated with mycotic infection.

Identification of two APOBEC3F splice variants displaying HIV-1 antiviral activity and contrasting sensitivity to Vif

Approximately half of all human genes undergo alternative mRNA splicing. This process often yields homologous gene products exhibiting diverse functions. Alternative splicing of APOBEC3G (A3G) and APOBEC3F (A3F), the major host resistance factors targeted by the HIV-1 protein Vif, has not been explored. We investigated the effects of alternative splicing on A3G/A3F gene expression and antiviral activity. Three alternatively spliced A3G mRNAs and two alternatively spliced A3F mRNAs were detected in peripheral blood mononuclear cells in each of 10 uninfected, healthy donors. Expression of these splice variants was altered in different cell subsets and in response to cellular stimulation. Alternatively spliced A3G variants were insensitive to degradation by Vif but displayed no antiviral activity against HIV-1. Conversely, alternative splicing of A3F produced a 37-kDa variant lacking exon 2 (A3FΔ2) that was prominently expressed in macrophages and monocytes and was resistant to Vif-mediated degradation. Alternative splicing also produced a 24-kDa variant of A3F lacking exons 2-4 (A3FΔ2-4) that was highly sensitive to Vif. Both A3FΔ2 and A3FΔ2-4 displayed reduced cytidine deaminase activity and moderate antiviral activity. These alternatively spliced A3F gene products, particularly A3FΔ2, were incorporated into HIV virions, albeit at levels less than wild-type A3F. Thus, alternative splicing of A3F mRNA generates truncated antiviral proteins that differ sharply in their sensitivity to Vif.

Bioinformatic analysis of TE-spliced new exons within human, mouse and zebrafish genomes

Recent studies indicate major roles for transposable elements (TEs) in alternative splicing. In this study, we conducted genome-wide alternative splicing analyses focusing on new internal exon birth derived from TEs in human, mouse, and zebrafish genomes. We identified two different exon sets, TE-spliced exons and non-TE-spliced exons. The proportion of TE-spliced exons was nearly twice as high as the proportion of non-TE-spliced exons in the coding sequence (CDS) region. Detailed analysis of various families of TEs in three different species of TE-spliced exons revealed a different pattern in zebrafish. In our analysis, we could identify the functional role of TE insertions in the vertebrate genome affecting mRNA splicing machinery. Their effects can be directly linked to the shift from constitutive to alternative splicing during primate evolution. Our results indicate that TEs have a significant effect on shaping new internal exons in human, mouse, and zebrafish transcriptomes.

Global regulation of alternative splicing during myogenic differentiation

Recent genome-wide analyses have elucidated the extent of alternative splicing (AS) in mammals, often focusing on comparisons of splice isoforms between differentiated tissues. However, regulated splicing changes are likely to be important in biological transitions such as cellular differentiation, or response to environmental stimuli. To assess the extent and significance of AS in myogenesis, we used splicing-sensitive microarray analysis of differentiating C2C12 myoblasts. We identified 95 AS events that undergo robust splicing transitions during C2C12 differentiation. More than half of the splicing transitions are conserved during differentiation of avian myoblasts, suggesting the products and timing of transitions are functionally significant. The majority of splicing transitions during C2C12 differentiation fall into four temporal patterns and were dependent on the myogenic program, suggesting that they are integral components of myogenic differentiation. Computational analyses revealed enrichment of many sequence motifs within the upstream and downstream intronic regions near the alternatively spliced regions corresponding to binding sites of splicing regulators. Western analyses demonstrated that several splicing regulators undergo dynamic changes in nuclear abundance during differentiation. These findings show that within a developmental context, AS is a highly regulated and conserved process, suggesting a major role for AS regulation in myogenic differentiation.

Splicing variants impact in thyroid normal physiology and pathological conditions

RNA splicing is an essential, precisely regulated process that occurs after gene transcription and before mRNA translation, in which introns may be removed and exons, retained. Variability in splicing patterns is a major source of protein diversity from the genome and function to generate a tremendously diverse proteome from a relatively small number of genes. Changes in splice site choice can determine different effects on the encoded protein. Small changes in peptide sequence can alter ligand binding, enzymatic activity, allosteric regulation, or protein localization. Errors in splicing regulation have been implicated in a number of different disease states. This study reviewed the mechanisms of splicing and their repercussion in endocrinology, emphasizing its importance in some thyroid physiological and pathological conditions.

Proteogenomics to discover the full coding content of genomes: a computational perspective

Proteogenomics has emerged as a field at the junction of genomics and proteomics. It is a loose collection of technologies that allow the search of tandem mass spectra against genomic databases to identify and characterize protein-coding genes. Proteogenomic peptides provide invaluable information for gene annotation, which is difficult or impossible to ascertain using standard annotation methods. Examples include confirmation of translation, reading-frame determination, identification of gene and exon boundaries, evidence for post-translational processing, identification of splice-forms including alternative splicing, and also, prediction of completely novel genes. For proteogenomics to deliver on its promise, however, it must overcome a number of technological hurdles, including speed and accuracy of peptide identification, construction and search of specialized databases, correction of sampling bias, and others. This article reviews the state of the art of the field, focusing on the current successes, and the role of computation in overcoming these challenges. We describe how technological and algorithmic advances have already enabled large-scale proteogenomic studies in many model organisms, including arabidopsis, yeast, fly, and human. We also provide a preview of the field going forward, describing early efforts in tackling the problems of complex gene structures, searching against genomes of related species, and immunoglobulin gene reconstruction.

Localization and functional characterization of the human NKCC2 isoforms

AIM

Salt reabsorption across the apical membrane of cells in the thick ascending limb (TAL) of Henle is primarily mediated by the bumetanide-sensitive Na(+)/K(+)/2Cl(-) cotransporter NKCC2. Three full-length splice variants of NKCC2 (NKCC2B, NKCC2A and NKCC2F) have been described. The NKCC2 isoforms have specific localizations and transport characteristics, as assessed for rabbit, rat and mouse. In the present study, we aimed to address the localization and transport characteristics of the human NKCC2 isoforms.

METHODS

RT-PCR, in situ hybridization and uptake studies in Xenopus oocytes were performed to characterize human NKCC2 isoforms.

RESULTS

All three classical NKCC2 isoforms were detected in the human kidney; in addition, we found splice variants with tandem duplicates of the variable exon 4. Contrary to rodents, in which NKCC2F is the most abundant NKCC2 isoform, NKCC2A was the dominant isoform in humans; similarly, isoform-specific in situ hybridization showed high expression levels of human NKCC2A along the TAL. Compared to NKCC2B and NKCC2F, human NKCC2A had the lowest Cl(-) affinity as determined by (86)Rb(+) uptake studies in oocytes. All NKCC2 isoforms were more efficiently inhibited by bumetanide than by furosemide. A sequence analysis of the amino acids encoded by exon 4 variants revealed high similarities between human and rodent NKCC2 isoforms, suggesting that differences in ion transport characteristics between species may be related to sequence variations outside the highly conserved sequence encoded by exon 4.

CONCLUSION

The human NKCC2 is an example of how differential splicing forms the basis for a diversification of transporter protein function.

A study of alternative splicing in the pig

Background

Since at least half of the genes in mammalian genomes are subjected to alternative splicing, alternative pre-mRNA splicing plays an important contribution to the complexity of the mammalian proteome. Expressed sequence tags (ESTs) provide evidence of a great number of possible alternative isoforms. With the EST resource for the domestic pig now containing more than one million porcine ESTs, it is possible to identify alternative splice forms of the individual transcripts in this species from the EST data with some confidence.

Results

The pig EST data generated by the Sino-Danish Pig Genome project has been assembled with publicly available ESTs and made available in the PigEST database. Using the Distiller package 2,515 EST clusters with candidate alternative isoforms were identified in the EST data with high confidence. In agreement with general observations in human and mouse, we find putative splice variants in about 30% of the contigs with more than 50 ESTs. Based on the criteria that a minimum of two EST sequences confirmed each splice event, a list of 100 genes with the most distinct tissue-specific alternative splice events was generated from the list of candidates. To confirm the tissue specificity of the splice events, 10 genes with functional annotation were randomly selected from which 16 individual splice events were chosen for experimental verification by quantitative PCR (qPCR). Six genes were shown to have tissue specific alternatively spliced transcripts with expression patterns matching those of the EST data. The remaining four genes had tissue-restricted expression of alternative spliced transcripts. Five out of the 16 splice events that were experimentally verified were found to be putative pig specific.

Conclusions

In accordance with human and rodent studies we estimate that approximately 30% of the porcine genes undergo alternative splicing. We found a good correlation between EST predicted tissue-specificity and experimentally validated splice events in different porcine tissue. This study indicates that a cluster size of around 50 ESTs is optimal for in silico detection of alternative splicing. Although based on a limited number of splice events, the study supports the notion that alternative splicing could have an important impact on species differentiation since 31% of the splice events studied appears to be species specific.

Genetically modified mice-successes and failures of a widely used technology

Genetically modified mice, created by random integration of a transgene into the genome or by targeted mutation of a specific gene, have proven to be extremely powerful tools for studying gene function in vivo. In this article, we give (1) a short overview of the traditional methods in mouse transgenesis and (2) a discussion of the problems with these methods, (3) more recent methods that were developed to overcome these problems, and (4) an outlook on future directions in gene targeting.

Mining expressed sequence tag (EST) libraries for cancer-associated genes

Originally established in the beginning of the 1990s as a direct route to gene finding, expressed sequence tags (ESTs) still lend themselves as a means to analyze gene expression in almost all human tissues. The type of questions that can be addressed using public EST libraries ranges from tissue-specific gene profiling to the comparison between tissues in diseased and healthy states. Thanks to a multitude of web-based online bioinformatics resources, mining in EST libraries is not restricted to experts in the field of data analysis, but can readily be performed by the medical or life scientist. In this chapter, a couple of cases studies are presented that guide the scientist to the most useful online resources so that they can conduct their own research.

Massive expansions of Dscam splicing diversity via staggered homologous recombination during arthropod evolution

The arthropod Down syndrome cell adhesion molecule (Dscam) gene can generate tens of thousands of protein isoforms via combinatorial splicing of numerous alternative exons encoding immunoglobulin variable domains organized into three clusters referred to as the exon 4, 6, and 9 clusters. Dscam protein diversity is important for nervous system development and immune functions. We have performed extensive phylogenetic analyses of Dscam from 20 arthropods (each containing between 46 and 96 alternative exons) to reconstruct the detailed history of exon duplication and loss events that built this remarkable system over 450 million years of evolution. Whereas the structure of the exon 4 cluster is ancient, the exon 6 and 9 clusters have undergone massive, independent expansions in each insect lineage. An analysis of nearly 2000 duplicated exons enabled detailed reconstruction of the timing, location, and boundaries of these duplication events. These data clearly show that new Dscam exons have arisen continuously throughout arthropod evolution and that this process is still occurring in the exon 6 and 9 clusters. Recently duplicated regions display boundaries corresponding to a single exon and the adjacent intron. The boundaries, homology, location, clustering, and relative frequencies of these duplication events strongly suggest that staggered homologous recombination is the major mechanism by which new Dscam exons evolve. These data provide a remarkably detailed picture of how complex gene structure evolves and reveal the molecular mechanism behind this process.