Loss of hepatic Mboat7 leads to liver fibrosis.

OBJECTIVE: The rs641738C>T variant located near the membrane-bound O-acyltransferase domain containing 7 (MBOAT7) locus is associated with fibrosis in liver diseases, including non-alcoholic fatty liver disease (NAFLD), alcohol-related liver disease, hepatitis B and C. We aim to understand the mechanism by which the rs641738C>T variant contributes to pathogenesis of NAFLD. DESIGN: Mice with hepatocyte-specific deletion of MBOAT7 (Mboat7Δhep) were generated and livers were characterised by histology, flow cytometry, qPCR, RNA sequencing and lipidomics. We analysed the association of rs641738C>T genotype with liver inflammation and fibrosis in 846 NAFLD patients and obtained genotype-specific liver lipidomes from 280 human biopsies. RESULTS: Allelic imbalance analysis of heterozygous human liver samples pointed to lower expression of the MBOAT7 transcript on the rs641738C>T haplotype. Mboat7Δhep mice showed spontaneous steatosis characterised by increased hepatic cholesterol ester content after 10 weeks. After 6 weeks on a high fat, methionine-low, choline-deficient diet, mice developed increased hepatic fibrosis as measured by picrosirius staining (p<0.05), hydroxyproline content (p<0.05) and transcriptomics, while the inflammatory cell populations and inflammatory mediators were minimally affected. In a human biopsied NAFLD cohort, MBOAT7 rs641738C>T was associated with fibrosis (p=0.004) independent of the presence of histological inflammation. Liver lipidomes of Mboat7Δhep mice and human rs641738TT carriers with fibrosis showed increased total lysophosphatidylinositol levels. The altered lysophosphatidylinositol and phosphatidylinositol subspecies in MBOAT7Δhep livers and human rs641738TT carriers were similar. CONCLUSION: Mboat7 deficiency in mice and human points to an inflammation-independent pathway of liver fibrosis that may be mediated by lipid signalling and a potentially targetable treatment option in NAFLD.

Tissue-, sex-, and age-specific DNA methylation of rat glucocorticoid receptor gene promoter and insulin-like growth factor 2 imprinting control region.

Tissue-, sex-, and age-specific epigenetic modifications such as DNA methylation are largely unknown. Changes in DNA methylation of the glucocorticoid receptor gene (NR3C1) and imprinting control region (ICR) of IGF2 and H19 genes during the lifespan are particularly interesting since these genes are susceptible to epigenetic modifications by prenatal stress or malnutrition. They are important regulators of development and aging. Methylation changes of NR3C1 affect glucocorticoid receptor expression, which is associated with stress sensitivity and stress-related diseases predominantly occurring during aging. Methylation changes of IGF2/H19 affect growth trajectory and nutrient use with risk of metabolic syndrome. Using a locus-specific approach, we characterized DNA methylation patterns of different Nr3c1 promoters and Igf2/H19 ICR in seven tissues of rats at 3, 9, and 24 mo of age. We found a complex pattern of locus-, tissue-, sex-, and age-specific DNA methylation. Tissue-specific methylation was most prominent at the shores of the Nr3c1 CpG island (CGI). Sex-specific differences in methylation peaked at 9 mo. During aging, Nr3c1 predominantly displayed hypomethylation mainly in females and at shores, whereas hypermethylation occurred within the CGI. Igf2/H19 ICR exhibited age-related hypomethylation occurring mainly in males. Methylation patterns of Nr3c1 in the skin correlated with those in the cortex, hippocampus, and hypothalamus. Skin may serve as proxy for methylation changes in central parts of the hypothalamic-pituitary-adrenal axis and hence for vulnerability to stress- and age-associated diseases. Thus, we provide in-depth insight into the complex DNA methylation changes of rat Nr3c1 and Igf2/H19 during aging that are tissue and sex specific.

Physiological state co-regulates thousands of mammalian mRNA splicing events at tandem splice sites and alternative exons.

Thousands of tandem alternative splice sites (TASS) give rise to mRNA insertion/deletion variants with small size differences. Recent work has concentrated on the question of biological relevance in general, and the physiological regulation of TASS in particular. We have quantitatively studied 11 representative TASS cases in comparison to one mutually exclusive exon case and two cassette exons (CEs) using a panel of human and mouse tissues, as well as cultured cell lines. Tissues show small but significant differences in TASS isoform ratios, with a variance 4- to 20-fold lower than seen for CEs. Remarkably, in cultured cells, all studied alternative splicing (AS) cases showed a cell-density-dependent shift of isoform ratios with similar time series profiles. A respective genome-wide co-regulation of TASS splicing was shown by next-generation mRNA sequencing data. Moreover, data from human and mouse organs indicate that this co-regulation of TASS occurs in vivo, with brain showing the strongest difference to other organs. Together, the results indicate a physiological AS regulation mechanism that functions almost independently from the splice site context and sequence.

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1.

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.

Genome-wide search for novel human uORFs and N-terminal protein extensions using ribosomal footprinting.

So far, the annotation of translation initiation sites (TISs) has been based mostly upon bioinformatics rather than experimental evidence. We adapted ribosomal footprinting to puromycin-treated cells to generate a transcriptome-wide map of TISs in a human monocytic cell line. A neural network was trained on the ribosomal footprints observed at previously annotated AUG translation initiation codons (TICs), and used for the ab initio prediction of TISs in 5062 transcripts with sufficient sequence coverage. Functional interpretation suggested 2994 novel upstream open reading frames (uORFs) in the 5' UTR, 1406 uORFs overlapping with the coding sequence, and 546 N-terminal protein extensions. The TIS detection method was validated on the basis of previously published alternative TISs and uORFs. Among primates, TICs in newly annotated TISs were significantly more conserved than control codons, both for AUGs and near-cognate codons. The transcriptome-wide map of novel candidate TISs derived as part of the study will shed further light on the way in which human proteome diversity is influenced by alternative translation initiation and regulation.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND: The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS: We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION: MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

Genetic and functional identification of the likely causative variant for cholesterol gallstone disease at the ABCG5/8 lithogenic locus.

UNLABELLED: The sterolin locus (ABCG5/ABCG8) confers susceptibility for cholesterol gallstone disease in humans. Both the responsible variant and the molecular mechanism causing an increased incidence of gallstones in these patients have as yet not been identified. Genetic mapping utilized patient samples from Germany (2,808 cases, 2,089 controls), Chile (680 cases, 442 controls), Denmark (366 cases, 766 controls), India (247 cases, 224 controls), and China (280 cases, 244 controls). Analysis of allelic imbalance in complementary DNA (cDNA) samples from human liver (n = 22) was performed using pyrosequencing. Transiently transfected HEK293 cells were used for [(3) H]-cholesterol export assays, analysis of protein expression, and localization of allelic constructs. Through fine mapping in German and Chilean samples, an ∼250 kB disease-associated interval could be defined for this locus. Lack of allelic imbalance or allelic splicing of the ABCG5 and ABCG8 transcripts in human liver limited the search to coding single nucleotide polymorphisms. Subsequent mutation detection and genotyping yielded two disease-associated variants: ABCG5-R50C (P = 4.94 × 10(-9) ) and ABCG8-D19H (P = 1.74 × 10(-10) ) in high pairwise linkage disequilibrium (r(2) = 0.95). [(3) H]-cholesterol export assays of allelic constructs harboring these genetic candidate variants demonstrated increased transport activity (3.2-fold, P = 0.003) only for the ABCG8-19H variant, which was also superior in nested logistic regression models in German (P = 0.018), Chilean (P = 0.030), and Chinese (P = 0.040) patient samples. CONCLUSION: This variant thus provides a molecular basis for biliary cholesterol hypersecretion as the mechanism for cholesterol gallstone formation, thereby drawing a link between "postgenomic" and "pregenomic" pathophysiological knowledge about this common complex disorder. (HEPATOLOGY 2012).

Sarcoidosis is associated with a truncating splice site mutation in BTNL2.

Sarcoidosis is a polygenic immune disorder with predominant manifestation in the lung. Genome-wide linkage analysis previously indicated that the extended major histocompatibility locus on chromosome 6p was linked to susceptibility to sarcoidosis. Here, we carried out a systematic three-stage SNP scan of 16.4 Mb on chromosome 6p21 in as many as 947 independent cases of familial and sporadic sarcoidosis and found that a 15-kb segment of the gene butyrophilin-like 2 (BTNL2) was associated with the disease. The primary disease-associated variant (rs2076530; P(TDT) = 3 x 10(-6), P(case-control) = 1.1 x 10(-8); replication P(TDT) = 0.0018, P(case-control) = 1.8 x 10(-6)) represents a risk factor that is independent of variation in HLA-DRB1. BTNL2 is a member of the immunoglobulin superfamily and has been implicated as a costimulatory molecule involved in T-cell activation on the basis of its homology to B7-1. The G --> A transition constituting rs2076530 leads to the use of a cryptic splice site located 4 bp upstream of the affected wild-type donor site. Transcripts of the risk-associated allele have a premature stop in the spliced mRNA. The resulting protein lacks the C-terminal IgC domain and transmembrane helix, thereby disrupting the membrane localization of the protein, as shown in experiments using green fluorescent protein and V5 fusion proteins.

Haplotypes of IL-12Rß1 impact on the clinical phenotype of hidradenitis suppurativa.

Antigen presentation in chronic skin disorders is mediated through the interleukin (IL)-12/IL-23 pathway and, hence, through the IL-12 receptor. Recent evidence suggesting dysregulated antigen presentation in skin lesions of hidradenitis suppurativa (HS) led to investigate the role of single nucleotide polymorphisms (SNPs) of the gene IL-12RB1 coding for the IL12-Rß1 receptor subunit. Genomic DNA was isolated from 139 patients and 113 healthy controls; nine SNPs in the transcribed region of IL12RB1 were genotyped. No significant differences of genotype and allele frequencies were found between the two groups. Two common haplotypes were recognized, namely h1 and h2. Carriage of h2 related with minor frequency alleles was associated with a greater risk for the acquisition of Hurley III disease stage and with the involvement of a greater number of skin areas. Patients with the h1 haplotype presented disease at an older age. This is the genetic study enrolling the largest number of patients with HS to date. Although SNPs of IL12RB1 do not seem to convey genetic predisposition, they are associated directly with the phenotype of the disease.

Widespread occurrence of alternative splicing at NAGNAG acceptors contributes to proteome plasticity.

Splice acceptors with the genomic NAGNAG motif may cause NAG insertion-deletions in transcripts, occur in 30% of human genes and are functional in at least 5% of human genes. We found five significant biases indicating that their distribution is nonrandom and that they are evolutionarily conserved and tissue-specific. Because of their subtle effects on mRNA and protein structures, these splice acceptors are often overlooked or underestimated, but they may have a great impact on biology and disease.

Statistical inference of allelic imbalance from transcriptome data.

Next-generation sequencing and the availability of high-density genotyping arrays have facilitated an analysis of somatic and meiotic mutations at unprecedented level, but drawing sensible conclusions about the functional relevance of the detected variants still remains a formidable challenge. In this context, the study of allelic imbalance in intermediate RNA phenotypes may prove a useful means to elucidate the likely effects of DNA variants of unknown significance. We developed a statistical framework for the assessment of allelic imbalance in next-generation transcriptome sequencing (RNA-seq) data that requires neither an expression reference nor the underlying nuclear genotype(s), and that allows for allele miscalls. Using extensive simulation as well as publicly available whole-transcriptome data from European-descent individuals in HapMap, we explored the power of our approach in terms of both genotype inference and allelic imbalance assessment under a wide range of practically relevant scenarios. In so doing, we verified a superior performance of our methodology, particularly at low sequencing coverage, compared to the more simplistic approach of completely ignoring allele miscalls. Because the proposed framework can be used to assess somatic mutations and allelic imbalance in one and the same set of RNA-seq data, it will be particularly useful for the analysis of somatic genetic variation in cancer studies.

Comprehensive assessment of sequence variation within the copy number variable defensin cluster on 8p23 by target enriched in-depth 454 sequencing.

BACKGROUND: In highly copy number variable (CNV) regions such as the human defensin gene locus, comprehensive assessment of sequence variations is challenging. PCR approaches are practically restricted to tiny fractions, and next-generation sequencing (NGS) approaches of whole individual genomes e.g. by the 1000 Genomes Project is confined by an affordable sequence depth. Combining target enrichment with NGS may represent a feasible approach. RESULTS: As a proof of principle, we enriched a ~850 kb section comprising the CNV defensin gene cluster DEFB, the invariable DEFA part and 11 control regions from two genomes by sequence capture and sequenced it by 454 technology. 6,651 differences to the human reference genome were found. Comparison to HapMap genotypes revealed sensitivities and specificities in the range of 94% to 99% for the identification of variations.Using error probabilities for rigorous filtering revealed 2,886 unique single nucleotide variations (SNVs) including 358 putative novel ones. DEFB CN determinations by haplotype ratios were in agreement with alternative methods. CONCLUSION: Although currently labor extensive and having high costs, target enriched NGS provides a powerful tool for the comprehensive assessment of SNVs in highly polymorphic CNV regions of individual genomes. Furthermore, it reveals considerable amounts of putative novel variations and simultaneously allows CN estimation.

Polymorphisms in the 3'-untranslated region of the CDH1 gene are a risk factor for primary gastric diffuse large B-cell lymphoma.

BACKGROUND: Primary gastric B-cell lymphomas arise from mucosa-associated lymphatic tissue (MALT) in patients with chronic Helicobacter pylori infection. We investigated whether germline variants in the CDH1 gene, coding for E-cadherin, genetically predispose patients to primary gastric B-cell lymphoma. DESIGN AND METHODS: Single marker analyses of the CDH1 gene were conducted in patients with primary gastric B-cell lymphoma (n=144), in patients with primary gastric high-grade lymphoma (n=61), and in healthy blood donors (n=361). Twelve single nucleotide polymorphisms were genotyped by TaqMan(®) technology. Allelic imbalance was tested by pyrosequencing and clone direct sequencing of heterozygote genomic and cDNA. Mutation detection was conducted around the poly-A signal of the CDH1 3'-untranslated region. The influence of the 3'-untranslated region on protein translation was determined by a luciferase reporter assay. RESULTS: Single marker analyses identified two single nucleotide polymorphisms in strong linkage disequilibrium located in the CDH1 3'-untranslated region. One of them was significantly associated with primary gastric diffuse large B-cell lymphomas after correction for multiple testing and this association was confirmed in an independent sample set. Patients homozygous for the rare T allele (rs1801026) had a 4.9-fold increased risk (95% CI: 1.5-15.9) of developing primary gastric diffuse large B-cell lymphoma. Allelic imbalance and reporter gene assays indicated a putative influence on mRNA stability and/or translational efficacy. CONCLUSIONS: We identified variants in CDH1 as the first potential genetic risk factors for the development of primary gastric diffuse large B-cell lymphomas. One of the potentially causative variants affects allelic CDH1 expression. These findings support the hypothesis that besides somatic alterations of B-cells, germline variants in the CDH1 gene contribute to a predisposition to the development of primary gastric diffuse large B-cell lymphomas.

Accurate prediction of NAGNAG alternative splicing.

Alternative splicing (AS) involving NAGNAG tandem acceptors is an evolutionarily widespread class of AS. Recent predictions of alternative acceptor usage reported better results for acceptors separated by larger distances, than for NAGNAGs. To improve the latter, we aimed at the use of Bayesian networks (BN), and extensive experimental validation of the predictions. Using carefully constructed training and test datasets, a balanced sensitivity and specificity of >or=92% was achieved. A BN trained on the combined dataset was then used to make predictions, and 81% (38/47) of the experimentally tested predictions were verified. Using a BN learned on human data on six other genomes, we show that while the performance for the vertebrate genomes matches that achieved on human data, there is a slight drop for Drosophila and worm. Lastly, using the prediction accuracy according to experimental validation, we estimate the number of yet undiscovered alternative NAGNAGs. State of the art classifiers can produce highly accurate prediction of AS at NAGNAGs, indicating that we have identified the major features of the 'NAGNAG-splicing code' within the splice site and its immediate neighborhood. Our results suggest that the mechanism behind NAGNAG AS is simple, stochastic, and conserved among vertebrates and beyond.

Epigenetic Modifications Associated with Maternal Anxiety during Pregnancy and Children's Behavioral Measures.

Epigenetic changes are associated with altered behavior and neuropsychiatric disorders and they modify the trajectory of aging. Maternal anxiety during pregnancy is a common environmental challenge for the fetus, causing changes in DNA methylation. Here, we determined the mediating role of DNA methylation and the moderating role of offspring sex on the association between maternal anxiety and children's behavioral measures. In 83 mother-child dyads, maternal anxiety was assessed in each trimester of pregnancy when the child was four years of age. Children's behavioral measures and children's buccal DNA methylation levels (NR3C1, IGF2/H19 ICR, and LINE1) were examined. Higher maternal anxiety during the third trimester was associated with more methylation levels of the NR3C1. Moderating effects of sex on the association between maternal anxiety and methylation were found for IGF2/H19 and LINE1 CpGs. Mediation analysis showed that methylation of NR3C1 could buffer the effects of maternal anxiety on children's behavioral measures, but this effect did not remain significant after controlling for covariates. In conclusion, our data support an association between maternal anxiety during pregnancy and DNA methylation. The results also underscore the importance of sex differences and timing effects. However, DNA methylation as underlying mechanism of the effect of maternal anxiety during pregnancy on offspring's behavioral measures was not supported.

TassDB2 - A comprehensive database of subtle alternative splicing events.

BACKGROUND: Subtle alternative splicing events involving tandem splice sites separated by a short (2-12 nucleotides) distance are frequent and evolutionarily widespread in eukaryotes, and a major contributor to the complexity of transcriptomes and proteomes. However, these events have been either omitted altogether in databases on alternative splicing, or only the cases of experimentally confirmed alternative splicing have been reported. Thus, a database which covers all confirmed cases of subtle alternative splicing as well as the numerous putative tandem splice sites (which might be confirmed once more transcript data becomes available), and allows to search for tandem splice sites with specific features and download the results, is a valuable resource for targeted experimental studies and large-scale bioinformatics analyses of tandem splice sites. Towards this goal we recently set up TassDB (Tandem Splice Site DataBase, version 1), which stores data about alternative splicing events at tandem splice sites separated by 3 nt in eight species. DESCRIPTION: We have substantially revised and extended TassDB. The currently available version 2 contains extensive information about tandem splice sites separated by 2-12 nt for the human and mouse transcriptomes including data on the conservation of the tandem motifs in five vertebrates. TassDB2 offers a user-friendly interface to search for specific genes or for genes containing tandem splice sites with specific features as well as the possibility to download result datasets. For example, users can search for cases of alternative splicing where the proportion of EST/mRNA evidence supporting the minor isoform exceeds a specific threshold, or where the difference in splice site scores is specified by the user. The predicted impact of each event on the protein is also reported, along with information about being a putative target for the nonsense-mediated decay (NMD) pathway. Links are provided to the UCSC genome browser and other external resources. CONCLUSION: TassDB2, available via http://www.tassdb.info, provides comprehensive resources for researchers interested in both targeted experimental studies and large-scale bioinformatics analyses of short distance tandem splice sites.

Haplotyping and copy number estimation of the highly polymorphic human beta-defensin locus on 8p23 by 454 amplicon sequencing.

BACKGROUND: The beta-defensin gene cluster (DEFB) at chromosome 8p23.1 is one of the most copy number (CN) variable regions of the human genome. Whereas individual DEFB CNs have been suggested as independent genetic risk factors for several diseases (e.g. psoriasis and Crohn's disease), the role of multisite sequence variations (MSV) is less well understood and to date has only been reported for prostate cancer. Simultaneous assessment of MSVs and CNs can be achieved by PCR, cloning and Sanger sequencing, however, these methods are labour and cost intensive as well as prone to methodological bias introduced by bacterial cloning. Here, we demonstrate that amplicon sequencing of pooled individual PCR products by the 454 technology allows in-depth determination of MSV haplotypes and estimation of DEFB CNs in parallel. RESULTS: Six PCR products spread over approximately 87 kb of DEFB and harbouring 24 known MSVs were amplified from 11 DNA samples, pooled and sequenced on a Roche 454 GS FLX sequencer. From approximately 142,000 reads, approximately 120,000 haplotype calls (HC) were inferred that identified 22 haplotypes ranging from 2 to 7 per amplicon. In addition to the 24 known MSVs, two additional sequence variations were detected. Minimal CNs were estimated from the ratio of HCs and compared to absolute CNs determined by alternative methods. Concordance in CNs was found for 7 samples, the CNs differed by one in 2 samples and the estimated minimal CN was half of the absolute in one sample. For 7 samples and 2 amplicons, the 454 haplotyping results were compared to those by cloning/Sanger sequencing. Intrinsic problems related to chimera formation during PCR and differences between haplotyping by 454 and cloning/Sanger sequencing are discussed. CONCLUSION: Deep amplicon sequencing using the 454 technology yield thousands of HCs per amplicon for an affordable price and may represent an effective method for parallel haplotyping and CN estimation in small to medium-sized cohorts. The obtained haplotypes represent a valuable resource to facilitate further studies of the biomedical impact of highly CN variable loci such as the beta-defensin locus.

Assessing the fraction of short-distance tandem splice sites under purifying selection.

Many alternative splice events result in subtle mRNA changes, and most of them occur at short-distance tandem donor and acceptor sites. The splicing mechanism of such tandem sites likely involves the stochastic selection of either splice site. While tandem splice events are frequent, it is unknown how many are functionally important. Here, we use phylogenetic conservation to address this question, focusing on tandems with a distance of 3-9 nucleotides. We show that previous contradicting results on whether alternative or constitutive tandem motifs are more conserved between species can be explained by a statistical paradox (Simpson's paradox). Applying methods that take biases into account, we found higher conservation of alternative tandems in mouse, dog, and even chicken, zebrafish, and Fugu genomes. We estimated a lower bound for the number of alternative sites that are under purifying (negative) selection. While the absolute number of conserved tandem motifs decreases with the evolutionary distance, the fraction under selection increases. Interestingly, a number of frameshifting tandems are under selection, suggesting a role in regulating mRNA and protein levels via nonsense-mediated decay (NMD). An analysis of the intronic flanks shows that purifying selection also acts on the intronic sequence. We propose that stochastic splice site selection can be an advantageous mechanism that allows constant splice variant ratios in situations where a deviation in this ratio is deleterious.

Systematic evaluation of the effect of common SNPs on pre-mRNA splicing.

The evolutionary and biomedical importance of differential mRNA splicing is well established. Numerous studies have assessed patterns of differential splicing in different genes and correlated these patterns to the genotypes for adjacent single-nucleotide polymorphisms (SNPs). Here, we have chosen a reverse approach and screened dbSNP for common SNPs at either canonical splice sites or exonic splice enhancers (ESEs) that would be classified as putatively splicing-relevant by bioinformatic tools. The 223 candidate SNPs retrieved from dbSNP were experimentally tested using a previously established panel of 92 matching DNAs and cDNAs. For each SNP, 16 cDNAs providing a balanced representation of the genotypes at the respective SNP were investigated by nested RT-PCR and subsequent sequencing. Putative allele-dependent splicing was verified by the cloning of PCR products. The positive predictive value of the bioinformatics tools turned out to be low, ranging from 0% for ESEfinder to 9% (in the case of acceptor-site SNPs) for a recently reported neural network. The results highlight the need for a better understanding of the sequence characteristics of functional splice-sites to improve our ability to predict in silico the splicing relevance of empirically observed DNA sequence variants.

TassDB: a database of alternative tandem splice sites.

Subtle alternative splice events at tandem splice sites are frequent in eukaryotes and substantially increase the complexity of transcriptomes and proteomes. We have developed a relational database, TassDB (TAndem Splice Site DataBase), which stores extensive data about alternative splice events at GYNGYN donors and NAGNAG acceptors. These splice events are of subtle nature since they mostly result in the insertion/deletion of a single amino acid or the substitution of one amino acid by two others. Currently, TassDB contains 114 554 tandem splice sites of eight species, 5209 of which have EST/mRNA evidence for alternative splicing. In addition, human SNPs that affect NAGNAG acceptors are annotated. The database provides a user-friendly interface to search for specific genes or for genes containing tandem splice sites with specific features as well as the possibility to download large datasets. This database should facilitate further experimental studies and large-scale bioinformatics analyses of tandem splice sites. The database is available at http://helios.informatik.uni-freiburg.de/TassDB/.