BALB/c and C57BL/6 Mice Differ in Polyreactive IgA Abundance, which Impacts the Generation of Antigen-Specific IgA and Microbiota Diversity.

The interrelationship between IgAs and microbiota diversity is still unclear. Here we show that BALB/c mice had higher abundance and diversity of IgAs than C57BL/6 mice and that this correlated with increased microbiota diversity. We show that polyreactive IgAs mediated the entrance of non-invasive bacteria to Peyer's patches, independently of CX3CR1(+) phagocytes. This allowed the induction of bacteria-specific IgA and the establishment of a positive feedback loop of IgA production. Cohousing of mice or fecal transplantation had little or no influence on IgA production and had only partial impact on microbiota composition. Germ-free BALB/c, but not C57BL/6, mice already had polyreactive IgAs that influenced microbiota diversity and selection after colonization. Together, these data suggest that genetic predisposition to produce polyreactive IgAs has a strong impact on the generation of antigen-specific IgAs and the selection and maintenance of microbiota diversity.

Next generation sequencing of SARS-CoV-2 genomes: challenges, applications and opportunities.

Various next generation sequencing (NGS) based strategies have been successfully used in the recent past for tracing origins and understanding the evolution of infectious agents, investigating the spread and transmission chains of outbreaks, as well as facilitating the development of effective and rapid molecular diagnostic tests and contributing to the hunt for treatments and vaccines. The ongoing COVID-19 pandemic poses one of the greatest global threats in modern history and has already caused severe social and economic costs. The development of efficient and rapid sequencing methods to reconstruct the genomic sequence of SARS-CoV-2, the etiological agent of COVID-19, has been fundamental for the design of diagnostic molecular tests and to devise effective measures and strategies to mitigate the diffusion of the pandemic. Diverse approaches and sequencing methods can, as testified by the number of available sequences, be applied to SARS-CoV-2 genomes. However, each technology and sequencing approach has its own advantages and limitations. In the current review, we will provide a brief, but hopefully comprehensive, account of currently available platforms and methodological approaches for the sequencing of SARS-CoV-2 genomes. We also present an outline of current repositories and databases that provide access to SARS-CoV-2 genomic data and associated metadata. Finally, we offer general advice and guidelines for the appropriate sharing and deposition of SARS-CoV-2 data and metadata, and suggest that more efficient and standardized integration of current and future SARS-CoV-2-related data would greatly facilitate the struggle against this new pathogen. We hope that our 'vademecum' for the production and handling of SARS-CoV-2-related sequencing data, will contribute to this objective.

YAP contributes to DNA methylation remodeling upon mouse embryonic stem cell differentiation.

The Yes-associated protein (YAP), one of the major effectors of the Hippo pathway together with its related protein WW-domain-containing transcription regulator 1 (WWTR1; also known as TAZ), mediates a range of cellular processes from proliferation and death to morphogenesis. YAP and WW-domain-containing transcription regulator 1 (WWTR1; also known as TAZ) regulate a large number of target genes, acting as coactivators of DNA-binding transcription factors or as negative regulators of transcription by interacting with the nucleosome remodeling and histone deacetylase complexes. YAP is expressed in self-renewing embryonic stem cells (ESCs), although it is still debated whether it plays any crucial roles in the control of either stemness or differentiation. Here we show that the transient downregulation of YAP in mouse ESCs perturbs cellular homeostasis, leading to the inability to differentiate properly. Bisulfite genomic sequencing revealed that this transient knockdown caused a genome-wide alteration of the DNA methylation remodeling that takes place during the early steps of differentiation, suggesting that the phenotype we observed might be due to the dysregulation of some of the mechanisms involved in regulation of ESC exit from pluripotency. By gene expression analysis, we identified two molecules that could have a role in the altered genome-wide methylation profile: the long noncoding RNA ephemeron, whose rapid upregulation is crucial for the transition of ESCs into epiblast, and the methyltransferase-like protein Dnmt3l, which, during the embryo development, cooperates with Dnmt3a and Dnmt3b to contribute to the de novo DNA methylation that governs early steps of ESC differentiation. These data suggest a new role for YAP in the governance of the epigenetic dynamics of exit from pluripotency.

VINYL: Variant prIoritizatioN bY survivaL analysis.

MOTIVATION: Clinical applications of genome re-sequencing technologies typically generate large amounts of data that need to be carefully annotated and interpreted to identify genetic variants potentially associated with pathological conditions. In this context, accurate and reproducible methods for the functional annotation and prioritization of genetic variants are of fundamental importance. RESULTS: In this article, we present VINYL, a flexible and fully automated system for the functional annotation and prioritization of genetic variants. Extensive analyses of both real and simulated datasets suggest that VINYL can identify clinically relevant genetic variants in a more accurate manner compared to equivalent state of the art methods, allowing a more rapid and effective prioritization of genetic variants in different experimental settings. As such we believe that VINYL can establish itself as a valuable tool to assist healthcare operators and researchers in clinical genomics investigations. AVAILABILITY AND IMPLEMENTATION: VINYL is available at http://beaconlab.it/VINYL and https://github.com/matteo14c/VINYL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNentropy: an entropy-based tool for the detection of significant variation of gene expression across multiple RNA-Seq experiments.

RNA sequencing (RNA-Seq) has become the experimental standard in transcriptome studies. While most of the bioinformatic pipelines for the analysis of RNA-Seq data and the identification of significant changes in transcript abundance are based on the comparison of two conditions, it is common practice to perform several experiments in parallel (e.g. from different individuals, developmental stages, tissues), for the identification of genes showing a significant variation of expression across all the conditions studied. In this work we present RNentropy, a methodology based on information theory devised for this task, which given expression estimates from any number of RNA-Seq samples and conditions identifies genes or transcripts with a significant variation of expression across all the conditions studied, together with the samples in which they are over- or under-expressed. To show the capabilities offered by our methodology, we applied it to different RNA-Seq datasets: 48 biological replicates of two different yeast conditions; samples extracted from six human tissues of three individuals; seven different mouse brain cell types; human liver samples from six individuals. Results, and their comparison to different state of the art bioinformatic methods, show that RNentropy can provide a quick and in depth analysis of significant changes in gene expression profiles over any number of conditions.

REDIportal: a comprehensive database of A-to-I RNA editing events in humans.

RNA editing by A-to-I deamination is the prominent co-/post-transcriptional modification in humans. It is carried out by ADAR enzymes and contributes to both transcriptomic and proteomic expansion. RNA editing has pivotal cellular effects and its deregulation has been linked to a variety of human disorders including neurological and neurodegenerative diseases and cancer. Despite its biological relevance, many physiological and functional aspects of RNA editing are yet elusive. Here, we present REDIportal, available online at http://srv00.recas.ba.infn.it/atlas/, the largest and comprehensive collection of RNA editing in humans including more than 4.5 millions of A-to-I events detected in 55 body sites from thousands of RNAseq experiments. REDIportal embeds RADAR database and represents the first editing resource designed to answer functional questions, enabling the inspection and browsing of editing levels in a variety of human samples, tissues and body sites. In contrast with previous RNA editing databases, REDIportal comprises its own browser (JBrowse) that allows users to explore A-to-I changes in their genomic context, empathizing repetitive elements in which RNA editing is prominent.

Reduced levels of protein recoding by A-to-I RNA editing in Alzheimer's disease.

Adenosine to inosine (A-to-I) RNA editing, catalyzed by the ADAR enzyme family, acts on dsRNA structures within pre-mRNA molecules. Editing of the coding part of the mRNA may lead to recoding, amino acid substitution in the resulting protein, possibly modifying its biochemical and biophysical properties. Altered RNA editing patterns have been observed in various neurological pathologies. Here, we present a comprehensive study of recoding by RNA editing in Alzheimer's disease (AD), the most common cause of irreversible dementia. We have used a targeted resequencing approach supplemented by a microfluidic-based high-throughput PCR coupled with next-generation sequencing to accurately quantify A-to-I RNA editing levels in a preselected set of target sites, mostly located within the coding sequence of synaptic genes. Overall, editing levels decreased in AD patients' brain tissues, mainly in the hippocampus and to a lesser degree in the temporal and frontal lobes. Differential RNA editing levels were observed in 35 target sites within 22 genes. These results may shed light on a possible association between the neurodegenerative processes typical for AD and deficient RNA editing.

Profile of microbial communities on carbonate stones of the medieval church of San Leonardo di Siponto (Italy) by Illumina-based deep sequencing.

Comprehensive studies of the biodiversity of the microbial epilithic community on monuments may provide critical insights for clarifying factors involved in the colonization processes. We carried out a high-throughput investigation of the communities colonizing the medieval church of San Leonardo di Siponto (Italy) by Illumina-based deep sequencing. The metagenomic analysis of sequences revealed the presence of Archaea, Bacteria, and Eukarya. Bacteria were Actinobacteria, Proteobacteria, Bacteroidetes, Cyanobacteria, Chloroflexi, Firmicutes and Candidatus Saccharibacteria. The predominant phylum was Actinobacteria, with the orders Actynomycetales and Rubrobacteriales, represented by the genera Pseudokineococcus, Sporichthya, Blastococcus, Arthrobacter, Geodermatophilus, Friedmanniella, Modestobacter, and Rubrobacter, respectively. Cyanobacteria sequences showing strong similarity with an uncultured bacterium sequence were identified. The presence of the green algae Oocystaceae and Trebuxiaceae was revealed. The microbial diversity was explored at qualitative and quantitative levels, evaluating the richness (the number of operational taxonomic units (OTUs)) and the abundance of reads associated with each OTU. The rarefaction curves approached saturation, suggesting that the majority of OTUs were recovered. The results highlighted a structured community, showing low diversity, made up of extremophile organisms adapted to desiccation and UV radiation. Notably, the microbiome appeared to be composed not only of microorganisms possibly involved in biodeterioration but also of carbonatogenic bacteria, such as those belonging to the genus Arthrobacter, which could be useful in bioconservation. Our investigation demonstrated that molecular tools, and in particular the easy-to-run next-generation sequencing, are powerful to perform a microbiological diagnosis in order to plan restoration and protection strategies.

RAP: RNA-Seq Analysis Pipeline, a new cloud-based NGS web application.

BACKGROUND: The study of RNA has been dramatically improved by the introduction of Next Generation Sequencing platforms allowing massive and cheap sequencing of selected RNA fractions, also providing information on strand orientation (RNA-Seq). The complexity of transcriptomes and of their regulative pathways make RNA-Seq one of most complex field of NGS applications, addressing several aspects of the expression process (e.g. identification and quantification of expressed genes and transcripts, alternative splicing and polyadenylation, fusion genes and trans-splicing, post-transcriptional events, etc.). METHODS: In order to provide researchers with an effective and friendly resource for analyzing RNA-Seq data, we present here RAP (RNA-Seq Analysis Pipeline), a cloud computing web application implementing a complete but modular analysis workflow. This pipeline integrates both state-of-the-art bioinformatics tools for RNA-Seq analysis and in-house developed scripts to offer to the user a comprehensive strategy for data analysis. RAP is able to perform quality checks (adopting FastQC and NGS QC Toolkit), identify and quantify expressed genes and transcripts (with Tophat, Cufflinks and HTSeq), detect alternative splicing events (using SpliceTrap) and chimeric transcripts (with ChimeraScan). This pipeline is also able to identify splicing junctions and constitutive or alternative polyadenylation sites (implementing custom analysis modules) and call for statistically significant differences in genes and transcripts expression, splicing pattern and polyadenylation site usage (using Cuffdiff2 and DESeq). RESULTS: Through a user friendly web interface, the RAP workflow can be suitably customized by the user and it is automatically executed on our cloud computing environment. This strategy allows to access to bioinformatics tools and computational resources without specific bioinformatics and IT skills. RAP provides a set of tabular and graphical results that can be helpful to browse, filter and export analyzed data, according to the user needs.

SpliceAid-F: a database of human splicing factors and their RNA-binding sites.

A comprehensive knowledge of all the factors involved in splicing, both proteins and RNAs, and of their interaction network is crucial for reaching a better understanding of this process and its functions. A large part of relevant information is buried in the literature or collected in various different databases. By hand-curated screenings of literature and databases, we retrieved experimentally validated data on 71 human RNA-binding splicing regulatory proteins and organized them into a database called 'SpliceAid-F' (http://www.caspur.it/SpliceAidF/). For each splicing factor (SF), the database reports its functional domains, its protein and chemical interactors and its expression data. Furthermore, we collected experimentally validated RNA-SF interactions, including relevant information on the RNA-binding sites, such as the genes where these sites lie, their genomic coordinates, the splicing effects, the experimental procedures used, as well as the corresponding bibliographic references. We also collected information from experiments showing no RNA-SF binding, at least in the assayed conditions. In total, SpliceAid-F contains 4227 interactions, 2590 RNA-binding sites and 1141 'no-binding' sites, including information on cellular contexts and conditions where binding was tested. The data collected in SpliceAid-F can provide significant information to explain an observed splicing pattern as well as the effect of mutations in functional regulatory elements.

Endometrial Cancer: A Pilot Study of the Tissue Microbiota.

BACKGROUND: The endometrium remains a difficult tissue for the analysis of microbiota, mainly due to the low bacterial presence and the sampling procedures. Among its pathologies, endometrial cancer has not yet been completely investigated for its relationship with microbiota composition. In this work, we report on possible correlations between endometrial microbiota dysbiosis and endometrial cancer. METHODS: Women with endometrial cancer at various stages of tumor progression were enrolled together with women with a benign polymyomatous uterus as the control. Analyses were performed using biopsies collected at two specific endometrial sites during the surgery. This study adopted two approaches: the absolute quantification of the bacterial load, using droplet digital PCR (ddPCR), and the analysis of the bacterial composition, using a deep metabarcoding NGS procedure. RESULTS: ddPCR provided the first-ever assessment of the absolute quantification of bacterial DNA in the endometrium, confirming a generally low microbial abundance. Metabarcoding analysis revealed a different microbiota distribution in the two endometrial sites, regardless of pathology, accompanied by an overall higher prevalence of pathogenic bacterial genera in cancerous tissues. CONCLUSIONS: These results pave the way for future studies aimed at identifying potential biomarkers and gaining a deeper understanding of the role of bacteria associated with tumors.

A platform independent RNA-Seq protocol for the detection of transcriptome complexity.

BACKGROUND: Recent studies have demonstrated an unexpected complexity of transcription in eukaryotes. The majority of the genome is transcribed and only a little fraction of these transcripts is annotated as protein coding genes and their splice variants. Indeed, most transcripts are the result of antisense, overlapping and non-coding RNA expression. In this frame, one of the key aims of high throughput transcriptome sequencing is the detection of all RNA species present in the cell and the first crucial step for RNA-seq users is represented by the choice of the strategy for cDNA library construction. The protocols developed so far provide the utilization of the entire library for a single sequencing run with a specific platform. RESULTS: We set up a unique protocol to generate and amplify a strand-specific cDNA library representative of all RNA species that may be implemented with all major platforms currently available on the market (Roche 454, Illumina, ABI/SOLiD). Our method is reproducible, fast, easy-to-perform and even allows to start from low input total RNA. Furthermore, we provide a suitable bioinformatics tool for the analysis of the sequences produced following this protocol. CONCLUSION: We tested the efficiency of our strategy, showing that our method is platform-independent, thus allowing the simultaneous analysis of the same sample with different NGS technologies, and providing an accurate quantitative and qualitative portrait of complex whole transcriptomes.

Cumulus Cell Transcriptome after Cumulus-Oocyte Complex Exposure to Nanomolar Cadmium in an In Vitro Animal Model of Prepubertal and Adult Age.

Cadmium (Cd), a highly toxic pollutant, impairs oocyte fertilization, through oxidative damage on cumulus cells (CCs). This study analysed the transcriptomic profile of CCs of cumulus-oocyte complexes (COCs) from adult and prepubertal sheep, exposed to Cd nanomolar concentration during in vitro maturation. In both age-groups, CCs of matured oocytes underwent RNA-seq, data analysis and validation. Differentially expressed genes (DEGs) were identified in adult (n = 99 DEGs) and prepubertal (n = 18 DEGs) CCs upon Cd exposure. Transcriptomes of adult CCs clustered separately between Cd-exposed and control samples, whereas prepubertal ones did not as observed by Principal Component Analysis. The transcriptomic signature of Cd-induced CC toxicity was identified by gene annotation and literature search. Genes associated with previous studies on ovarian functions and/or Cd effects were confirmed and new genes were identified, thus implementing the knowledge on their involvement in such processes. Enrichment and validation analysis showed that, in adult CCs, Cd acted as endocrine disruptor on DEGs involved in hormone biosynthesis, cumulus expansion, regulation of cell signalling, growth and differentiation and oocyte maturation, whereas in prepubertal CCs, Cd affected DEGs involved in CC development and viability and CC-oocyte communications. In conclusion, these DEGs could be used as valuable non-invasive biomarkers for oocyte competence.

Cyclic AMP induces reversible EPAC1 condensates that regulate histone transcription.

The second messenger cyclic AMP regulates many nuclear processes including transcription, pre-mRNA splicing and mitosis. While most functions are attributed to protein kinase A, accumulating evidence suggests that not all nuclear cyclic AMP-dependent effects are mediated by this kinase, implying that other effectors may be involved. Here we explore the nuclear roles of Exchange Protein Activated by cyclic AMP 1. We find that it enters the nucleus where forms reversible biomolecular condensates in response to cyclic AMP. This phenomenon depends on intrinsically disordered regions present at its amino-terminus and is independent of protein kinase A. Finally, we demonstrate that nuclear Exchange Protein Activated by cyclic AMP 1 condensates assemble at genomic loci on chromosome 6 in the proximity of Histone Locus Bodies and promote the transcription of a histone gene cluster. Collectively, our data reveal an unexpected mechanism through which cyclic AMP contributes to nuclear spatial compartmentalization and promotes the transcription of specific genes.

Yeast as a Model to Unravel New BRCA2 Functions in Cell Metabolism.

Mutations in BRCA2 gene increase the risk for breast cancer and for other cancer types, including pancreatic and prostate cancer. Since its first identification as an oncosupressor in 1995, the best-characterized function of BRCA2 is in the repair of DNA double-strand breaks (DSBs) by homologous recombination. BRCA2 directly interacts with both RAD51 and single-stranded DNA, mediating loading of RAD51 recombinase to sites of single-stranded DNA. In the absence of an efficient homologous recombination pathway, DSBs accumulate resulting in genome instability, thus supporting tumorigenesis. Yet the precise mechanism by which BRCA2 exerts its tumor suppressor function remains unclear. BRCA2 has also been involved in other biological functions including protection of telomere integrity and stalled replication forks, cell cycle progression, transcriptional control and mitophagy. Recently, we and others have reported a role of BRCA2 in modulating cell death programs through a molecular mechanism conserved in yeast and mammals. Here we hypothesize that BRCA2 is a multifunctional protein which exerts specific functions depending on cell stress response pathway. Based on a differential RNA sequencing analysis carried out on yeast cells either growing or undergoing a regulated cell death process, either in the absence or in the presence of BRCA2, we suggest that BRCA2 causes central carbon metabolism reprogramming in response to death stimuli and encourage further investigation on the role of metabolic reprogramming in BRCA2 oncosuppressive function.

Natural and after colon washing fecal samples: the two sides of the coin for investigating the human gut microbiome.

To date several studies address the important role of gut microbiome and its interplay with the human host in the health and disease status. However, the selection of a universal sampling matrix representative of the microbial biodiversity associated with the gastrointestinal (GI) tract, is still challenging. Here we present a study in which, through a deep metabarcoding analysis of the 16S rRNA gene, we compared two sampling matrices, feces (F) and colon washing feces (CWF), in order to evaluate their relative effectiveness and accuracy in representing the complexity of the human gut microbiome. A cohort of 30 volunteers was recruited and paired F and CWF samples were collected from each subject. Alpha diversity analysis confirmed a slightly higher biodiversity of CWF compared to F matched samples. Likewise, beta diversity analysis proved that paired F and CWF microbiomes were quite similar in the same individual, but remarkable inter-individual variability occurred among the microbiomes of all participants. Taxonomic analysis in matched samples was carried out to investigate the intra and inter individual/s variability. Firmicutes, Bacteroidota, Proteobacteria and Actinobacteriota were the main phyla in both F and CWF samples. At genus level, Bacteirodetes was the most abundant in F and CWF samples, followed by Faecalibacterium, Blautia and Escherichia-Shigella. Our study highlights an inter-individual variability greater than intra-individual variability for paired F and CWF samples. Indeed, an overall higher similarity was observed across matched F and CWF samples, suggesting, as expected, a remarkable overlap between the microbiomes inferred using the matched F and CWF samples. Notably, absolute quantification of total 16S rDNA by droplet digital PCR (ddPCR) revealed comparable overall microbial load between paired F and CWF samples. We report here the first comparative study on fecal and colon washing fecal samples for investigating the human gut microbiome and show that both types of samples may be used equally for the study of the gut microbiome. The presented results suggest that the combined use of both types of sampling matrices could represent a suitable choice to obtain a more complete overview of the human gut microbiota for addressing different biological and clinical questions.

Identification and functional characterization of two new transcriptional variants of the human p63 gene.

p63 belongs to a family of transcription factors, which, while demonstrating striking conservation of functional domains, regulate distinct biological functions. Its principal role is in the regulation of epithelial commitment, differentiation and maintenance programs, during embryogenesis and in adult tissues. The p63 gene has a complex transcriptional pattern, producing two subclasses of N-terminal isoforms (TA and DeltaN) which are alternatively spliced at the C-terminus. Here, we report the identification of two new C-terminus p63 variants, we named p63 delta and epsilon, that increase from 6 to 10 the number of the p63 isoforms. Expression analysis of all p63 variants demonstrates a tissue/cell-type-specific nature of p63 alternative transcript expression, probably related to their different cellular functions. We demonstrate that the new p63 variants as DeltaN isoforms are active as transcription factors as they have nuclear localization and can modulate the expression of p63 target genes. Moreover, we report that, like DeltaNp63alpha, DeltaNp63delta and epsilon sustain cellular proliferation and that their expression decreases during keratinocyte differentiation, suggesting their involvement in this process. Taken together, our results demonstrate the existence of novel p63 proteins whose expression should be considered in future studies on the roles of p63 in the regulation of cellular functions.

Accurate detection and quantification of SARS-CoV-2 genomic and subgenomic mRNAs by ddPCR and meta-transcriptomics analysis.

SARS-CoV-2 replication requires the synthesis of a set of structural proteins expressed through discontinuous transcription of ten subgenomic mRNAs (sgmRNAs). Here, we have fine-tuned droplet digital PCR (ddPCR) assays to accurately detect and quantify SARS-CoV-2 genomic ORF1ab and sgmRNAs for the nucleocapsid (N) and spike (S) proteins. We analyzed 166 RNA samples from anonymized SARS-CoV-2 positive subjects and we observed a recurrent and characteristic pattern of sgmRNAs expression in relation to the total viral RNA content. Additionally, expression profiles of sgmRNAs, as determined by meta-transcriptomics sequencing of a subset of 110 RNA samples, were highly correlated with those obtained by ddPCR. By providing a comprehensive and dynamic snapshot of the levels of SARS-CoV-2 sgmRNAs in infected individuals, our results may contribute a better understanding of the dynamics of transcription and expression of the genome of SARS-CoV-2 and facilitate the development of more accurate molecular diagnostic tools for the stratification of COVID-19 patients.

Campylobacter vulpis sp. nov. isolated from wild red foxes.

During a sampling of wild red foxes (Vulpes vulpes) for the detection of Epsilonproteobacteria, 14 strains were isolated from the caecal contents of 14 epidemiologically-unrelated animals. A genus-specific PCR indicated that the isolates belonged to the genus Campylobacter. Based on the results of a species-specific PCR, the isolates were initially identified as C. upsaliensis. However, multi-locus sequence typing (MLST) revealed that the isolates were significantly different from the C. upsaliensis present in the MLST database. A polyphasic study, including conventional biochemical and tolerance characteristics, morphology by transmission electron microscopy (TEM), MALDI-TOF analysis, and genetic comparisons based on partial 16S rDNA and atpA gene sequences, was undertaken. Finally, the complete genome sequence of the type strain 251/13T and the draft genome sequences of the other isolates were determined. Average nucleotide identity, average amino acid identity and in silico DNA-DNA hybridization analyses confirmed that the isolates represent a novel taxon for which the name Campylobacter vulpis sp. nov. is proposed, with isolate 251/13T (=CCUG 70587T = LMG 30110T) as the type strain. In order to allow a rapid discrimination of C. vulpis from the closely-related C. upsaliensis, a specific PCR test was designed, based on atpA gene sequences.

Identification of tumor-associated cassette exons in human cancer through EST-based computational prediction and experimental validation.

BACKGROUND: Many evidences report that alternative splicing, the mechanism which produces mRNAs and proteins with different structures and functions from the same gene, is altered in cancer cells. Thus, the identification and characterization of cancer-specific splice variants may give large impulse to the discovery of novel diagnostic and prognostic tumour biomarkers, as well as of new targets for more selective and effective therapies. RESULTS: We present here a genome-wide analysis of the alternative splicing pattern of human genes through a computational analysis of normal and cancer-specific ESTs from seventeen anatomical groups, using data available in AspicDB, a database resource for the analysis of alternative splicing in human. By using a statistical methodology, normal and cancer-specific genes, splice sites and cassette exons were predicted in silico. The condition association of some of the novel normal/tumoral cassette exons was experimentally verified by RT-qPCR assays in the same anatomical system where they were predicted. Remarkably, the presence in vivo of the predicted alternative transcripts, specific for the nervous system, was confirmed in patients affected by glioblastoma. CONCLUSION: This study presents a novel computational methodology for the identification of tumor-associated transcript variants to be used as cancer molecular biomarkers, provides its experimental validation, and reports specific biomarkers for glioblastoma.