A graph theoretic approach to testing associations between disparate sources of functional genomic data

The last few years have seen the advent of high-throughput technologies to analyze various properties of the transcriptome and proteome of several organisms. The congruency of these different data sources, or lack thereof, can shed light on the mechanisms that govern cellular function. A central challenge for bioinformatics research is to develop a unified framework for combining the multiple sources of functional genomics information and testing associations between them, thus obtaining a robust and integrated view of the underlying biology. We present a graph theoretic approach to test the significance of the association between multiple disparate sources of functional genomics data by proposing two statistical tests, namely edge permutation and node label permutation tests. We demonstrate the use of the proposed tests by finding significant association between a Gene Ontology-derived "predictome" and data obtained from mRNA expression and phenotypic experiments for Saccharomyces cerevisiae. Moreover, we employ the graph theoretic framework to recast a surprising discrepancy presented in Giaever et al. (2002) between gene expression and knockout phenotype, using expression data from a different set of experiments.

Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

Common target genes of palatal and gingival fibroblasts for EMD: the microarray approach.

BACKGROUND AND OBJECTIVE Connective tissue grafts are frequently applied, together with Emdogain(®) , for root coverage. However, it is unknown whether fibroblasts from the gingiva and from the palate respond similarly to Emdogain. The aim of this study was therefore to evaluate the effect of Emdogain(®) on fibroblasts from palatal and gingival connective tissue using a genome-wide microarray approach. MATERIAL AND METHODS Human palatal and gingival fibroblasts were exposed to Emdogain(®) and RNA was subjected to microarray analysis followed by gene ontology screening with Database for Annotation, Visualization and Integrated Discovery functional annotation clustering, Kyoto Encyclopedia of Genes and Genomes pathway analysis and the Search Tool for the Retrieval of Interacting Genes/Proteins functional protein association network. Microarray results were confirmed by quantitative RT-PCR analysis. RESULTS The transcription levels of 106 genes were up-/down-regulated by at least five-fold in both gingival and palatal fibroblasts upon exposure to Emdogain(®) . Gene ontology screening assigned the respective genes into 118 biological processes, six cellular components, eight molecular functions and five pathways. Among the striking patterns observed were the changing expression of ligands targeting the transforming growth factor-beta and gp130 receptor family as well as the transition of mesenchymal epithelial cells. Moreover, Emdogain(®) caused changes in expression of receptors for chemokines, lipids and hormones, and for transcription factors such as SMAD3, peroxisome proliferator-activated receptor gamma and those of the ETS family. CONCLUSION The present data suggest that Emdogain(®) causes substantial alterations in gene expression, with similar patterns observed in palatal and gingival fibroblasts.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.

Doxorubicin Affects Expression of Proteins of Neuronal Pathways in MCF-7 Breast Cancer Cells.

In the present article, we report on the semi-quantitative proteome analysis and related changes in protein expression of the MCF-7 breast cancer cell line following treatment with doxorubicin, using the precursor acquisition independent from ion count (PAcIFIC) mass spectrometry method. PAcIFIC represents a cost-effective and easy-to-use proteomics approach, enabling for deep proteome sequencing with minimal sample handling. The acquired proteomic data sets were searched for regulated Reactome pathways and Gene Ontology annotation terms using a new algorithm (SetRank). Using this approach, we identified pathways with significant changes (≤0.05), such as chromatin organization, DNA binding, embryo development, condensed chromosome, sequence-specific DNA binding, response to oxidative stress and response to toxin, as well as others. These sets of pathways are already well-described as being susceptible to chemotherapeutic drugs. Additionally, we found pathways related to neuron development, such as central nervous system neuron differentiation, neuron projection membrane and SNAP receptor activity. These later pathways might indicate biological mechanisms on the molecular level causing the known side-effect of doxorubicin chemotherapy, characterized as cognitive impairment, also called 'chemo brain'. Mass spectrometry data are available via ProteomeXchange with identifier PXD002998.

Experiment on marine fungus Microascus brevicaulis strain

The marine fungus Microascus brevicaulis strain LF580 is a non-model secondary metabolite producer with high yields of the two secondary metabolites scopularide A and B, which exhibit distinct activities against tumour cell lines. A mutant strain was obtained using UV mutagenesis, showing besides higher production levels faster growth and differences in pellet formation. Comparative proteomics were applied to gain deeper understanding of the regulation of production and of the physiology of this fungus and its mutant. For this purpose, an optimised protein extraction protocol was established. Here, we show the first proteome study of a marine fungus. In total, 4759 proteins were identified. The central metabolic pathway of LF580 could be mapped by using KEGG pathway analysis and GO annotation. Using iTRAQ labelling, 318 proteins were shown to be significantly regulated in the mutant strain: 189 were down- and 129 upregulated. Proteomics are a powerful tool for the understanding of regulatory aspects: The differences on proteome level could be attributed to a limited nutrient availability in wild type strain due to a strong pellet formation. This information can be applied to optimisation on strain and process level. The linkage between nutrient limitation and pellet formation in the non-model fungus M. brevicaulis is in consensus with the knowledge on model organisms like Aspergillus niger and Penicillium chrysogenum.

Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to CO2-driven acidification during the initiation of calcification

The impact of ocean acidification (OA) on coral calcification, a subject of intense current interest, is poorly understood in part because of the presence of symbionts in adult corals. Early life history stages of Acropora spp. provide an opportunity to study the effects of elevated CO(2) on coral calcification without the complication of symbiont metabolism. Therefore, we used the Illumina RNAseq approach to study the effects of acute exposure to elevated CO(2) on gene expression in primary polyps of Acropora millepora, using as reference a novel comprehensive transcriptome assembly developed for this study. Gene ontology analysis of this whole transcriptome data set indicated that CO(2) -driven acidification strongly suppressed metabolism but enhanced extracellular organic matrix synthesis, whereas targeted analyses revealed complex effects on genes implicated in calcification. Unexpectedly, expression of most ion transport proteins was unaffected, while many membrane-associated or secreted carbonic anhydrases were expressed at lower levels. The most dramatic effect of CO(2) -driven acidification, however, was on genes encoding candidate and known components of the skeletal organic matrix that controls CaCO(3) deposition. The skeletal organic matrix effects included elevated expression of adult-type galaxins and some secreted acidic proteins, but down-regulation of other galaxins, secreted acidic proteins, SCRiPs and other coral-specific genes, suggesting specialized roles for the members of these protein families and complex impacts of OA on mineral deposition. This study is the first exhaustive exploration of the transcriptomic response of a scleractinian coral to acidification and provides an unbiased perspective on its effects during the early stages of calcification.

Transcriptional changes and oxidative stress associated with the synergistic interaction between potato virus X and potato virus Y and their relationship with symptom expression.

Many virus diseases of economic importance to agriculture result from mixtures of different pathogens invading the host at a given time. This contrasts with the relatively scarce studies available on the molecular events associated with virus---host interactions in mixed infections. Compared with single infections, co-infection of Nicotiana benthamiana with Potato virus X (PVX) and Potato virus Y (PVY) resulted in increased systemic symptoms (synergism) that led to necrosis of the newly emerging leaves and death of the plant. A comparative transcriptional analysis was undertaken to identify quantitative and qualitative differences in gene expression during this synergistic infection and correlate these changes with the severe symptoms it caused. Global transcription profiles of doubly infected leaves were compared with those from singly infected leaves using gene ontology enrichment analysis and metabolic pathway annotator software. Functional gene categories altered by the double infection comprise suites of genes regulated coordinately, which are associated with chloroplast functions (downregulated), protein synthesis and degradation (upregulated), carbohydrate metabolism (upregulated), and response to biotic stimulus and stress (upregulated). The expressions of reactive oxygen species?generating enzymes as well as several mitogen-activated protein kinases were also significantly induced. Accordingly, synergistic infection induced a severe oxidative stress in N. benthamiana leaves, as judged by increases in lipid peroxidation and by the generation of superoxide radicals in chloroplasts, which correlated with the misregulation of antioxidative genes in microarray data. Interestingly, expression of genes encoding oxylipin biosynthesis was uniquely upregulated by the synergistic infection. Virus-induced gene silencing of ?-dioxygenase1 delayed cell death during PVX?PVY infection.

Investigação de mecanismos genéticos e epigenéticos em distúrbios do crescimento humano

Parcela considerável de pacientes com distúrbios de crescimento não têm a causa de seus quadros clínicos estabelecida, incluindo aproximadamente 50% dos pacientes com diagnóstico clínico de síndrome de Silver−Russell (SRS) e 10-20% dos pacientes com síndrome de Beckwith-Wiedemann (BWS). O objetivo deste estudo foi investigar as causas genéticas e epigenéticas de distúrbios de crescimento, de etiologia desconhecida, numa contribuição para o entendimento de mecanismos que regulam o crescimento. O estudo compreendeu: (1) a investigação de microdesequilíbrios cromossômicos, por aCGH; (2) a análise do perfil de expressão alelo-específica de genes sujeitos a imprinting (IG), por pirossequenciamento (PSQ) ou sequenciamento de Sanger; (3) a investigação do padrão de metilação global em pacientes com restrição de crescimento, utilizando microarray de metilação. A casuística constituiu-se de 41 pacientes não aparentados, com distúrbios de crescimento, de etiologia desconhecida: (1) 25, com hipótese diagnóstica de SRS; (2) seis, com restrição de crescimento intrauterino e peso ao nascimento abaixo do 10º percentil, associados a outros sinais clínicos; (3) sete, com hipótese diagnóstica de BWS; e (4) três, com macrossomia pré-natal ou pós-natal, associada a outros sinais. A investigação de microdesequilíbrios cromossômicos foi realizada em 40 pacientes. Foram detectadas 58 variantes raras em 30/40 pacientes (75%): 40 foram consideradas provavelmente benignas (18 pacientes, 45%), 12, com efeito patogênico desconhecido (11 pacientes, 27,5%), duas, provavelmente patogênicas (um paciente, 2,5%) e quatro, patogênicas (três pacientes, 7,5%). Essas frequências são comparáveis àquelas descritas em estudos que investigaram CNV em grupos de pacientes com distúrbios de crescimento e outras alterações congênitas, incluindo SRS, e mostram a importância da investigação de microdesequilíbrios cromossômicos nesses pacientes. A diversidade dos microdesequilíbrios cromossômicos identificados é reflexo da heterogeneidade clínica das casuísticas. Neste estudo, muitos dos pacientes com hipótese diagnóstica de SRS e BWS apresentavam sinais clínicos atípicos, explicando a ausência neles das alterações (epi)genéticas que causam essas síndromes. A identificação de CNV características de outras síndromes reflete a sobreposição de sinais clínicos com BWS e SRS. A análise do perfil de expressão alelo-específica de IG foi realizada em um subgrupo de 18 pacientes com restrição de crescimento. Trinta IG com função em proliferação celular, crescimento fetal ou neurodesenvolvimento foram inicialmente selecionados. Após seleção de SNP transcritos com alta frequência na população, genotipagem de pacientes, genitores e indivíduos controle, determinação da expressão dos IG em sangue periférico e seu padrão de expressão (mono ou bialélico), 13 IG, expressos no sangue, tiveram a expressão alelo-específica avaliada, sete deles por PSQ e seis por sequenciamento de Sanger. Alterações no perfil de expressão de dois genes, de expressão normalmente paterna, foram detectadas em 4/18 pacientes (22%). Este estudo é o primeiro a utilizar pirossequenciamento e sequenciamento de Sanger na avaliação do perfil de expressão alelo-específica de IG, em pacientes com restrição de crescimento. Apesar de terem limitações, ambas as técnicas mostraram-se robustas e revelaram alterações de expressão alélica interessantes; entretanto, a relação dessas alterações com o quadro clínico dos pacientes permanece por esclarecer. A investigação da metilação global do DNA foi realizada em subgrupo de 21 pacientes com restrição de crescimento e em 24 indivíduos controle. Dois tipos de análise foram realizados: (1) análise diferencial de grupo e (2) análise diferencial individual. Na primeira análise, em que foi comparado o padrão de metilação do grupo de pacientes com quadro clínico sugestivo de SRS (n=16) com o do grupo controle (n=24), não houve indicação de hipo ou hipermetilação global no grupo SRS. Na segunda análise, foi comparado o padrão de metilação de cada um dos 21 pacientes com restrição de crescimento e dos 24 indivíduos controle, com o padrão de metilação do grupo controle. O número médio de CpG hipermetilados e de segmentos diferencialmente metilados (SDM) foi significativamente maior nos pacientes. Foram identificados 82 SDM hipermetilados, estando 57 associados a gene(s) (69,5%), em 16 pacientes, e 51 SDM hipometilados, 41 deles associados a gene(s) (80,4%), em 10 pacientes. A análise de ontologia genética dos 61 genes associados aos SDM hipo ou hipermetilados nos pacientes destacou genes que atuam no desenvolvimento e na morfogênese do sistema esquelético e de órgãos fetais, e na regulação da transcrição gênica e de processos metabólicos. Alterações de metilação em genes que atuam em processos de proliferação e diferenciação celulares e crescimento foram identificadas em 9/20 dos pacientes (45%), sugerindo implicação clínica. Não foi detectada alteração epigenética comum aos pacientes com diagnóstico clínico de SRS, explicável provavelmente pela heterogeneidade clínica. A investigação de metilação global, utilizando microarray, produziu novos dados que podem contribuir para a compreensão de mecanismos moleculares que influenciam o crescimento pré- e pós-natal. Na translocação aparentemente equilibrada - t(5;6)(q35.2;p22.3)dn, detectada em paciente com suspeita clínica de SRS, a interrupção de um gene, pela quebra no cromossomo 6, pode ser a causa do quadro clínico; alternativamente, a translocação pode ter impactado a regulação de genes de desenvolvimento localizados próximos aos pontos de quebra. A análise de expressão em sangue periférico mostrou que os níveis de cDNA do gene, interrompido pelo ponto de quebra da translocação, estavam reduzidos à metade. Além de sinais típicos da SRS, a paciente apresentava algumas características clínicas sugestivas de displasia cleidocraniana. Assim, a translocação t(5;6) pode ter alterado a interação de genes de desenvolvimento e seus elementos reguladores, levando à desregulação de sua expressão espaço-temporal, e resultando num fenótipo atípico, com características sobrepostas de mais de uma síndrome genética

Sequencing and functional analysis of the genome of a nematode egg-parasitic fungus, Pochonia chlamydosporia

Pochonia chlamydosporia is a worldwide-distributed soil fungus with a great capacity to infect and destroy the eggs and kill females of plant-parasitic nematodes. Additionally, it has the ability to colonize endophytically roots of economically-important crop plants, thereby promoting their growth and eliciting plant defenses. This multitrophic behavior makes P. chlamydosporia a potentially useful tool for sustainable agriculture approaches. We sequenced and assembled ∼41 Mb of P. chlamydosporia genomic DNA and predicted 12,122 gene models, of which many were homologous to genes of fungal pathogens of invertebrates and fungal plant pathogens. Predicted genes (65%) were functionally annotated according to Gene Ontology, and 16% of them found to share homology with genes in the Pathogen Host Interactions (PHI) database. The genome of this fungus is highly enriched in genes encoding hydrolytic enzymes, such as proteases, glycoside hydrolases and carbohydrate esterases. We used RNA-Seq technology in order to identify the genes expressed during endophytic behavior of P. chlamydosporia when colonizing barley roots. Functional annotation of these genes showed that hydrolytic enzymes and transporters are expressed during endophytism. This structural and functional analysis of the P. chlamydosporia genome provides a starting point for understanding the molecular mechanisms involved in the multitrophic lifestyle of this fungus. The genomic information provided here should also prove useful for enhancing the capabilities of this fungus as a biocontrol agent of plant-parasitic nematodes and as a plant growth-promoting organism.

Endophytic colonization of barley (Hordeum vulgare) roots by the nematophagous fungus Pochonia chlamydosporia reveals plant growth promotion and a general defense and stress transcriptomic response

Plant crop yields are negatively conditioned by a large set of biotic and abiotic factors. An alternative to mitigate these adverse effects is the use of fungal biological control agents and endophytes. The egg-parasitic fungus Pochonia chlamydosporia has been traditionally studied because of its potential as a biological control agent of plant-parasitic nematodes. This fungus can also act as an endophyte in monocot and dicot plants, and has been shown to promote plant growth in different agronomic crops. An Affymetrix 22K Barley GeneChip was used in this work to analyze the barley root transcriptomic response to P. chlamydosporia root colonization. Functional gene ontology (GO) and gene set enrichment analyses showed that genes involved in stress response were enriched in the barley transcriptome under endophytism. An 87.5 % of the probesets identified within the abiotic stress response group encoded heat shock proteins. Additionally, we found in our transcriptomic analysis an up-regulation of genes implicated in the biosynthesis of plant hormones, such as auxin, ethylene and jasmonic acid. Along with these, we detected induction of brassinosteroid insensitive 1-associated receptor kinase 1 (BR1) and other genes related to effector-triggered immunity (ETI) and pattern-triggered immunity (PTI). Our study supports at the molecular level the growth-promoting effect observed in plants endophytically colonized by P. chlamydosporia, which opens the door to further studies addressing the capacity of this fungus to mitigate the negative effects of biotic and abiotic factors on plant crops.

Systematic expression profiling of the mouse transcriptome using RIKEN cDNA mircoarrays

The number of known mRNA transcripts in the mouse has been greatly expanded by the RIKEN Mouse Gene Encyclopedia project. Validation of their reproducible expression in a tissue is an important contribution to the study of functional genomics. In this report, we determine the expression profile of 57,931 clones on 20 mouse tissues using cDNA microarrays. Of these 57,931 clones, 22,928 clones correspond to the FANTOM2 clone set. The set represents 20,234 transcriptional units (TUs) out of 33,409 TUs in the FANTOM2 set. We identified 7206 separate clones that satisfied stringent criteria for tissue-specific expression. Gene Ontology terms were assigned for these 7206 clones, and the proportion of 'molecular function' ontology for each tissue-specific clone was examined. These data will provide insights into the function of each tissue. Tissue-specific gene expression profiles obtained using our cDNA microarrays were also compared with the data extracted from the GNF Expression Atlas based on Affymetrix microarrays. One major outcome of the RIKEN transcriptome analysis is the identification of numerous nonprotein-coding mRNAs. The expression profile was also used to obtain evidence of expression for putative noncoding RNAs. In addition, 1926 clones (70%) of 2768 clones that were categorized as unknown EST, and 1969 (58%) clones of 3388 clones that were categorized as unclassifiable were also shown to be reproducibly expressed.

ASD: the Alternative Splicing Database

Alternative splicing is widespread in mammalian gene expression, and variant splice patterns are often specific to different stages of development, particular tissues or a disease state. There is a need to systematically collect data on alternatively spliced exons, introns and splice isoforms, and to annotate this data. The Alternative Splicing Database consortium has been addressing this need, and is committed to maintaining and developing a value-added database of alternative splice events, and of experimentally verified regulatory mechanisms that mediate splice variants. In this paper we present two of the products from this project: namely, a database of computationally delineated alternative splice events as seen in alignments of EST/cDNA sequences with genome sequences, and a database of alternatively spliced exons collected from literature. The reported splice events are from nine different organisms and are annotated for various biological features including expression states and cross-species conservation. The data are presented on our ASD web pages (http://www.ebi.ac.uk/asd).