Analysis of the canine brain transcriptome with an emphasis on the hypothalamus and cerebral cortex.

The diversity of dog breeds make the domestic dog a valuable model for identifying genes responsible for many phenotypic and behavioral traits. The brain, in particular, is a region of interest for the analysis of molecular changes that are involved in dog-specific behavioral phenotypes. However, such studies are handicapped due to incomplete annotation of the dog genome. We present a high-coverage transcriptome of the dog brain using RNA-Seq. Two areas of the brain, hypothalamus and cerebral cortex, were selected for their roles in cognition, emotion, and neuroendocrine functions. We detected many novel features of the dog transcriptome, including 13,799 novel exons, 51,357 exons with unique 5' or 3' modifications, and many novel alternative splicing events. We provide some examples of novel features in genes that are related to domestication, including ADCY8, SMOC2, and PRNP. We also found 247 novel protein-coding genes and 328 noncoding RNAs, including 57 long noncoding RNAs that represent the first empirical evidence for a large fraction of noncoding RNAs in the dog. In addition, we analyze both gene expression and alternative splicing differences between the hypothalamus and cerebral cortex and find that there is very little overlap between genes that are differentially alternatively spliced and genes that are differentially expressed. We thereby suggest that researchers who want to pinpoint the genetic causes for dog breed-specific traits and diseases should not confine their studies to gene expression alone, but should consider other factors such as alternative splicing and changes in untranslated regions.

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.

The effect of intron length on exon creation ratios during the evolution of mammalian genomes.

Recent studies report that alternatively spliced exons tend to occur in longer introns, which is attributed to the length constraints for splice site pairing for the two major splicing mechanisms, intron definition versus exon definition. Using genome-wide studies of EST and microarray data from human and mouse, we have analyzed the distribution of various subsets of alternatively spliced exons, based on their inclusion level and evolutionary history, versus increasing intron length. Alternative exons may be included in either a major or minor fraction of all transcripts (known as major-form and minor-form exons, respectively). We find that major-form exons are seven- to eightfold more likely to be contained in short introns (<400 nt) than minor-form exons, which occur preferentially in longer introns. Since minor-form exons are more likely to be novel (approximately 75%), this implied that novel exons arise more frequently in longer introns. To test this hypothesis, we used whole genome alignments to classify exons according to their phylogenetic age. We find that older exons, i.e., exons that are conserved in all mammals, predominate at shorter intron lengths, for both major- and minor-form exons. In contrast, exons that arose recently during primate evolution are more prevalent at longer intron lengths (>1000 nt). This suggests that the observed correlation of longer intron lengths with alternatively spliced exons may be at least partly due to biases in the probability of exon creation, which is higher in long introns.

Vascular abnormalities in mice deficient for the G protein-coupled receptor GPR4 that functions as a pH sensor.

GPR4 is a G protein-coupled receptor expressed in the vasculature, lung, kidney, and other tissues. In vitro ectopic overexpression studies implicated GPR4 in sensing extracellular pH changes leading to cyclic AMP (cAMP) production. To investigate its biological roles in vivo, we generated GPR4-deficient mice by homologous recombination. Whereas GPR4-null adult mice appeared phenotypically normal, neonates showed a higher frequency of perinatal mortality. The average litter size from GPR4(-/-) intercrosses was approximately 30% smaller than that from GPR4(+/+) intercrosses on N3 and N5 C57BL/6 genetic backgrounds. A fraction of knockout embryos and neonates had spontaneous hemorrhages, dilated and tortuous subcutaneous blood vessels, and defective vascular smooth muscle cell coverage. Mesangial cells in kidney glomeruli were also significantly reduced in GPR4-null neonates. Some neonates exhibited respiratory distress with airway lining cell metaplasia. To examine whether GPR4 is functionally involved in vascular pH sensing, an ex vivo aortic ring assay was used under defined pH conditions. Compared to wild-type aortas, microvessel outgrowth from GPR4-null aortas was less inhibited by acidic extracellular pH. Treatment with an analog of cAMP, a downstream effector of GPR4, abolished microvessel outgrowth bypassing the GPR4-knockout phenotype. These results suggest that GPR4 deficiency leads to partially penetrant vascular abnormalities during development and that this receptor functions in blood vessel pH sensing.

The ASAP II database: analysis and comparative genomics of alternative splicing in 15 animal species.

We have greatly expanded the Alternative Splicing Annotation Project (ASAP) database: (i) its human alternative splicing data are expanded approximately 3-fold over the previous ASAP database, to nearly 90,000 distinct alternative splicing events; (ii) it now provides genome-wide alternative splicing analyses for 15 vertebrate, insect and other animal species; (iii) it provides comprehensive comparative genomics information for comparing alternative splicing and splice site conservation across 17 aligned genomes, based on UCSC multigenome alignments; (iv) it provides an approximately 2- to 3-fold expansion in detection of tissue-specific alternative splicing events, and of cancer versus normal specific alternative splicing events. We have also constructed a novel database linking orthologous exons and orthologous introns between genomes, based on multigenome alignment of 17 animal species. It can be a valuable resource for studies of gene structure evolution. ASAP II provides a new web interface enabling more detailed exploration of the data, and integrating comparative genomics information with alternative splicing data. We provide a set of tools for advanced data-mining of ASAP II with Pygr (the Python Graph Database Framework for Bioinformatics) including powerful features such as graph query, multigenome alignment query, etc. ASAP II is available at http://www.bioinformatics.ucla.edu/ASAP2.

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.

Evidence that public database records for many cancer-associated genes reflect a splice form found in tumors and lack normal splice forms.

Alternative splicing is widespread in the human genome, and it appears that many genes display different splice forms in cancerous tissue than in normal human tissues. However, since cDNAs for many cancer-associated genes were originally cloned from tumor samples, it is important to ask whether this repertoire of cDNAs provides a complete or representative picture of the transcript isoforms found in normal tissues. To answer this, we used bioinformatics and RT-PCR to identify novel splice forms, focusing on in-frame exonskips, for a panel of 50 cancer-associated genes in normal tissue samples. These data show that in nearly two-thirds of the genes, normal tissues expressed previously unknown splice forms, of which 40% were normally a dominant splice form. Surprisingly, the tumor-associated splice forms were twice as likely to be represented in GenBank than their normal tissue-associated splice forms, most likely because 70% of the mRNAs in GenBank for these genes were cloned from tumor samples. As an example, we describe a novel normal splice form of IKBbeta, an important regulator of the NFkappaB pathway. Our data suggest that systematic re-evaluation of cancer genes' splice forms in normal tissue will yield insights into their distinct functions in normal tissues and in cancer. Our database contains 1308 novel normal splice forms, including many known cancer genes.

Detecting tissue-specific regulation of alternative splicing as a qualitative change in microarray data.

Alternative splicing has recently emerged as a major mechanism of regulation in the human genome, occurring in perhaps 40-60% of human genes. Thus, microarray studies of functional regulation could, in principle, be extended to detect not only the changes in the overall expression of a gene, but also changes in its splicing pattern between different tissues. However, since changes in the total expression of a gene and changes in its alternative splicing can be mixed in complex ways among a set of samples, separating these effects can be difficult, and is essential for their accurate assessment. We present a simple and general approach for distinguishing changes in alternative splicing from changes in expression, based on detecting systematic anti-correlation between the log-ratios of two different samples versus a pool containing both samples. We have tested this analysis method on microarray data for five human tissues, generated using a standard microarray platform and experimental protocols shown previously to be sensitive to alternative splicing. Our automatic analysis was able to detect a wide variety of tissue-specific alternative splicing events, such as exon skipping,mutually exclusive exons, alternative 3' and alternative 5' splicing, alternative initiation and alternative termination, all of which were validated by independent reverse-transcriptase PCR experiments, with validation rates of 70-85%. Our analysis method also enables hierarchical clustering of genes and samples by the level of similarity to their alternative splicing patterns, revealing patterns of tissue-specific regulation that are distinct from those obtained by hierarchical clustering of gene expression from the same microarray data. Our data and analysis source code are available from http://www.bioinformatics.ucla.edu/ASAP.

Analysis of alternative splicing with microarrays: successes and challenges.

Recently, DNA microarrays have emerged as potentially powerful tools for analyzing alternative splicing. We briefly review the latest results in this field and highlight the current challenges that they have revealed.

Cytokine control of memory B cell homing machinery.

The germinal center (GC) is a pivotal site for the development of B cell memory. Whereas GC B cells do not chemotax to most chemokines and do not express the adhesion receptors L-selectin, alpha(4)beta(7), and cutaneous lymphocyte Ag (CLA), memory B cells respond to various chemotactic signals and express adhesion receptors. In this study, we show that CD40 ligand, IL-2, and IL-10 together drive this transition of GC B cells to memory phenotype in vitro, up-regulating memory B cell markers, chemotactic responses to CXC ligand (CXCL)12, CXCL13, and CCL19, and expression of adhesion receptors L-selectin, alpha(4)beta(7), and CLA. Moreover, addition of IL-4 modulates this transition, preventing chemotactic responses to CXCL12 and CXCL13 (but not to CCL19), and inhibiting the re-expression of L-selectin, but not of CLA or alpha(4)beta(7). CCR7 expression, responsiveness to CCL19, and L-selectin/alpha(4)beta(7) phenotype are coordinately regulated. Thus, IL-2/IL-10 and IL-4 play important and distinctive roles in developing the migratory capacities of memory B cells.