A sketch of known and novel MYCN-associated miRNA networks in neuroblastoma

Neuroblastoma (NB) originates from neural crest-derived precursors and represents the most common childhood extracranial solid tumour. MicroRNAs (miRNAs), a class of small non-coding RNAs that participate in a wide variety of biological processes by regulating gene expression, appear to play an essential role within the NB context. High-throughput next generation sequencing (NGS) was applied to study the miRNA transcriptome in a cohort of NB tumours with and without MYCN-amplification (MNA and MNnA, respectively) and in dorsal root ganglia (DRG), as a control. Out of the 128 miRNAs differentially expressed in the NB vs. DRG comparison, 47 were expressed at higher levels, while 81 were expressed at lower levels in the NB tumours. We also found that 23 miRNAs were differentially expressed in NB with or without MYCN-amplification, with 17 miRNAs being upregulated and 6 being downregulated in the MNA subtypes. Functional annotation analysis of the target genes of these differentially expressed miRNAs demonstrated that many mRNAs were involved in cancer-related pathways, such as DNA-repair and apoptosis as well as FGFR and EGFR signalling. In particular, we found that miR-628-3p negatively affects MYCN gene expression. Furthermore, we identified a novel miRNA candidate with variable expression in MNA vs. MNnA tumours, whose putative target genes are implicated in the mTOR pathway. The present study provides further insight into the molecular mechanisms that correlate miRNA dysregulation to NB development and progression.

In-depth characterization of breast cancer tumor-promoting cell transcriptome by RNA sequencing and microarrays

Numerous studies have reported the existence of tumor-promoting cells (TPC) with self-renewal potential and a relevant role in drug resistance. However, pathways and modifications involved in the maintenance of such tumor subpopulations are still only partially understood. Sequencing-based approaches offer the opportunity for a detailed study of TPC including their transcriptome modulation. Using microarrays and RNA sequencing approaches, we compared the transcriptional profiles of parental MCF7 breast cancer cells with MCF7-derived TPC (i.e. MCFS). Data were explored using different bioinformatic approaches, and major findings were experimentally validated. The different analytical pipelines (Lifescope and Cufflinks based) yielded similar although not identical results. RNA sequencing data partially overlapped microarray results and displayed a higher dynamic range, although overall the two approaches concordantly predicted pathway modifications. Several biological functions were altered in TPC, ranging from production of inflammatory cytokines (i.e., IL-8 and MCP-1) to proliferation and response to steroid hormones. More than 300 non-coding RNAs were defined as differentially expressed, and 2,471 potential splicing events were identified. A consensus signature of genes up-regulated in TPC was derived and was found to be significantly associated with insensitivity to fulvestrant in a public breast cancer patient dataset. Overall, we obtained a detailed portrait of the transcriptome of a breast cancer TPC line, highlighted the role of non-coding RNAs and differential splicing, and identified a gene signature with a potential as a context-specific biomarker in patients receiving endocrine treatment.

BITS 2015: the annual meeting of the Italian Society of Bioinformatics

This preface introduces the content of the BioMed Central journal Supplements related to the BITS 2015 meeting, held in Milan, Italy, from the 3^th to the 5^th of June, 2015.

Histone Modifications in a Mouse Model of Early Adversities and Panic Disorder: Role for Asic1 and Neurodevelopmental Genes

Hyperventilation following transient, CO₂-induced acidosis is ubiquitous in mammals and heritable. In humans, respiratory and emotional hypersensitivity to CO₂ marks separation anxiety and panic disorders, and is enhanced by early-life adversities. Mice exposed to the repeated cross-fostering paradigm (RCF) of interference with maternal environment show heightened separation anxiety and hyperventilation to 6% CO₂-enriched air. Gene-environment interactions affect CO₂ hypersensitivity in both humans and mice. We therefore hypothesised that epigenetic modifications and increased expression of genes involved in pH-detection could explain these relationships. Medullae oblongata of RCF- and normally-reared female outbred mice were assessed by ChIP-seq for H3Ac, H3K4me3, H3K27me3 histone modifications, and by SAGE for differential gene expression. Integration of multiple experiments by network analysis revealed an active component of 148 genes pointing to the mTOR signalling pathway and nociception. Among these genes, Asic1 showed heightened mRNA expression, coherent with RCF-mice’s respiratory hypersensitivity to CO₂ and altered nociception. Functional enrichment and mRNA transcript analyses yielded a consistent picture of enhancement for several genes affecting chemoception, neurodevelopment, and emotionality. Particularly, results with Asic1 support recent human findings with panic and CO₂ responses, and provide new perspectives on how early adversities and genes interplay to affect key components of panic and related disorders.

Preface: BITS2014, the annual meeting of the Italian Society of Bioinformatics

This Preface introduces the content of the BioMed Central journal Supplements related to BITS2014 meeting, held in Rome, Italy, from the 26th to the 28th of February, 2014.

DOT1L-mediated H3K79me2 modification critically regulates gene expression during cardiomyocyte differentiation

Epigenetic changes on DNA and chromatin are implicated in cell differentiation and organogenesis. For the heart, distinct histone methylation profiles were recently linked to stage-specific gene expression programs during cardiac differentiation in vitro. However, the enzymes catalyzing these modifications and the genes regulated by them remain poorly defined. We therefore decided to identify the epigenetic enzymes that are potentially involved in cardiomyogenesis by analyzing the expression profile of the 85 genes encoding the epigenetic-related proteins in mouse cardiomyocytes (CMs), and then study how they affect gene expression during differentiation and maturation of this cell type. We show here with gene expression screening of epigenetic enzymes that the highly expressed H3 methyltransferase disruptor of telomeric silencing 1-like (DOT1L) drives a transitional pattern of di-methylation on H3 lysine 79 (H3K79) in CMs at different stages of differentiation in vitro and in vivo. Through a genome-wide chromatin-immunoprecipitation DNA-sequencing approach, we found H3K79me2 enriched at genes expressed during cardiac differentiation. Moreover, knockdown of Dot1L affected the expression of H3K79me2-enriched genes. Our results demonstrate that histone methylation, and in particular DOT1L-mediated H3K79me2 modification, drives cardiomyogenesis through the definition of a specific transcriptional landscape.

Deep Sequencing the microRNA profile in rhabdomyosarcoma reveals down-regulation of miR-378 family members

BACKGROUND

Rhabdomyosarcoma (RMS) is a highly malignant tumour accounting for nearly half of soft tissue sarcomas in children. MicroRNAs (miRNAs) represent a class of short, non-coding, regulatory RNAs which play a critical role in different cellular processes. Altered miRNA levels have been reported in human cancers, including RMS.

METHODS

Using deep sequencing technology, a total of 685 miRNAs were investigated in a group of alveolar RMSs (ARMSs), embryonal RMSs (ERMSs) as well as in normal skeletal muscle (NSM). Q-PCR, MTT, cytofluorimetry, migration assay, western blot and immunofluorescence experiments were carried out to determine the role of miR-378a-3p in cancer cell growth, apoptosis, migration and differentiation. Bioinformatics pipelines were used for miRNA target prediction and clustering analysis.

RESULTS

Ninety-seven miRNAs were significantly deregulated in ARMS and ERMS when compared to NSM. MiR-378 family members were dramatically decreased in RMS tumour tissue and cell lines. Interestingly, members of the miR-378 family presented as a possible target the insulin-like growth factor receptor 1 (IGF1R), a key signalling molecule in RMS. MiR-378a-3p over-expression in an RMS-derived cell line suppressed IGF1R expression and affected phosphorylated-Akt protein levels. Ectopic expression of miR-378a-3p caused significant changes in apoptosis, cell migration, cytoskeleton organization as well as a modulation of the muscular markers MyoD1, MyoR, desmin and MyHC. In addition, DNA demethylation by 5-aza-2'-deoxycytidine (5-aza-dC) was able to up-regulate miR-378a-3p levels with a concomitant induction of apoptosis, decrease in cell viability and cell cycle arrest in G2-phase. Cells treated with 5-aza-dC clearly changed their morphology and expressed moderate levels of MyHC.

CONCLUSIONS

MiR-378a-3p may function as a tumour suppressor in RMS and the restoration of its expression would be of therapeutic benefit in RMS. Furthermore, the role of epigenetic modifications in RMS deserves further investigations.

Whole transcriptome RNA sequencing data from blood leukocytes derived from Parkinson's disease patients prior to and following deep brain stimulation treatment

Recent evidence demonstrates the power of RNA sequencing (RNA-Seq) for identifying valuable and urgently needed blood biomarkers and advancing both early and accurate detection of neurological diseases, and in particular Parkinson's disease (PD). RNA sequencing technology enables non-biased, high throughput, probe-independent inspection of expression data and high coverage and both quantification of global transcript levels as well as the detection of expressed exons and junctions given a sufficient sequencing depth (coverage). However, the analysis of sequencing data frequently presents a bottleneck. Tools for quantification of alternative splicing from sequenced libraries hardly exist at the present time, and methods that support multiple sequencing platforms are especially lacking. Here, we describe in details a whole RNA-Seq transcriptome dataset produced from PD patient's blood leukocytes. The samples were taken prior to, and following deep brain stimulation (DBS) treatment while being on stimulation and following 1 h of complete electrical stimulation cessation and from healthy control volunteers. We describe in detail the methodology applied for analyzing the RNA-Seq data including differential expression of long noncoding RNAs (lncRNAs). We also provide details of the corresponding analysis of in-depth splice isoform data from junction and exon reads, with the use of the software AltAnalyze. Both the RNA-Seq raw (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42608) and analyzed data (https://www.synapse.org/#!Synapse:syn2805267) may be found valuable towards detection of novel blood biomarkers for PD.

Emerging bioinformatics approaches for analysis of NGS-derived coding and non-coding RNAs in neurodegenerative diseases

Neurodegenerative diseases in general and specifically late-onset Alzheimer’s disease (LOAD) involve a genetically complex and largely obscure ensemble of causative and risk factors accompanied by complex feedback responses. The advent of “high-throughput” transcriptome investigation technologies such as microarray and deep sequencing is increasingly being combined with sophisticated statistical and bioinformatics analysis methods complemented by knowledge-based approaches such as Bayesian Networks or network and graph analyses. Together, such “integrative” studies are beginning to identify co-regulated gene networks linked with biological pathways and potentially modulating disease predisposition, outcome, and progression. Specifically, bioinformatics analyses of integrated microarray and genotyping data in cases and controls reveal changes in gene expression of both protein-coding and small and long regulatory RNAs; highlight relevant quantitative transcriptional differences between LOAD and non-demented control brains and demonstrate reconfiguration of functionally meaningful molecular interaction structures in LOAD. These may be measured as changes in connectivity in “hub nodes” of relevant gene networks (Zhang etal., 2013). We illustrate here the open analytical questions in the transcriptome investigation of neurodegenerative disease studies, proposing “ad hoc” strategies for the evaluation of differential gene expression and hints for a simple analysis of the non-coding RNA (ncRNA) part of such datasets. We then survey the emerging role of long ncRNAs (lncRNAs) in the healthy and diseased brain transcriptome and describe the main current methods for computational modeling of gene networks. We propose accessible modular and pathway-oriented methods and guidelines for bioinformatics investigations of whole transcriptome next generation sequencing datasets. We finally present methods and databases for functional interpretations of lncRNAs and propose a simple heuristic approach to visualize and represent physical and functional interactions of the coding and non-coding components of the transcriptome. Integrating in a functional and integrated vision coding and ncRNA analyses is of utmost importance for current and future analyses of neurodegenerative transcriptomes.

Long non-coding RNA and alternative splicing modulations in Parkinson's leukocytes identified by RNA sequencing

The continuously prolonged human lifespan is accompanied by increase in neurodegenerative diseases incidence, calling for the development of inexpensive blood-based diagnostics. Analyzing blood cell transcripts by RNA-Seq is a robust means to identify novel biomarkers that rapidly becomes a commonplace. However, there is lack of tools to discover novel exons, junctions and splicing events and to precisely and sensitively assess differential splicing through RNA-Seq data analysis and across RNA-Seq platforms. Here, we present a new and comprehensive computational workflow for whole-transcriptome RNA-Seq analysis, using an updated version of the software AltAnalyze, to identify both known and novel high-confidence alternative splicing events, and to integrate them with both protein-domains and microRNA binding annotations. We applied the novel workflow on RNA-Seq data from Parkinson's disease (PD) patients' leukocytes pre- and post- Deep Brain Stimulation (DBS) treatment and compared to healthy controls. Disease-mediated changes included decreased usage of alternative promoters and N-termini, 5′-end variations and mutually-exclusive exons. The PD regulated FUS and HNRNP A/B included prion-like domains regulated regions. We also present here a workflow to identify and analyze long non-coding RNAs (lncRNAs) via RNA-Seq data. We identified reduced lncRNA expression and selective PD-induced changes in 13 of over 6,000 detected leukocyte lncRNAs, four of which were inversely altered post-DBS. These included the U1 spliceosomal lncRNA and RP11-462G22.1, each entailing sequence complementarity to numerous microRNAs. Analysis of RNA-Seq from PD and unaffected controls brains revealed over 7,000 brain-expressed lncRNAs, of which 3,495 were co-expressed in the leukocytes including U1, which showed both leukocyte and brain increases. Furthermore, qRT-PCR validations confirmed these co-increases in PD leukocytes and two brain regions, the amygdala and substantia-nigra, compared to controls. This novel workflow allows deep multi-level inspection of RNA-Seq datasets and provides a comprehensive new resource for understanding disease transcriptome modifications in PD and other neurodegenerative diseases.

Long non-coding RNAs (lncRNAs) comprise a novel, fascinating class of RNAs with largely unknown biological functions. Parkinson's-disease (PD) is the most frequent motor disorder, and Deep-brain-stimulation (DBS) treatment alleviates the symptoms, but early disease biomarkers are still unknown and new future genetic interference targets are urgently needed. Using RNA-sequencing technology and a novel computational workflow for in-depth exploration of whole-transcriptome RNA-seq datasets, we detected and analyzed lncRNAs in sequenced libraries from PD patients' leukocytes pre and post-treatment and the brain, adding this full profile resource of over 7,000 lncRNAs to the few human tissues-derived lncRNA datasets that are currently available. Our study includes sample-specific database construction, detecting disease-derived changes in known and novel lncRNAs, exons and junctions and predicting corresponding changes in Polyadenylation choices, protein domains and miRNA binding sites. We report widespread transcript structure variations at the splice junction and exons levels, including novel exons and junctions and alteration of lncRNAs followed by experimental validation in PD leukocytes and two PD brain regions compared with controls. Our results suggest lncRNAs involvement in neurodegenerative diseases, and specifically PD. This comprehensive workflow will be of use to the increasing number of laboratories producing RNA-Seq data in a wide range of biomedical studies.

Transient exposure to low levels of insecticide affects metabolic networks of honeybee larvae

The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L(-1)) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators.

Estrogen-dependent dynamic profile of eNOS-DNA associations in prostate cancer

In previous work we have documented the nuclear translocation of endothelial NOS (eNOS) and its participation in combinatorial complexes with Estrogen Receptor Beta (ERβ) and Hypoxia Inducible Factors (HIFs) that determine localized chromatin remodeling in response to estrogen (E2) and hypoxia stimuli, resulting in transcriptional regulation of genes associated with adverse prognosis in prostate cancer (PCa). To explore the role of nuclear eNOS in the acquisition of aggressive phenotype in PCa, we performed ChIP-Sequencing on chromatin-associated eNOS from cells from a primary tumor with poor outcome and from metastatic LNCaP cells. We found that: 1. the eNOS-bound regions (peaks) are widely distributed across the genome encompassing multiple transcription factors binding sites, including Estrogen Response Elements. 2. E2 increased the number of peaks, indicating hormone-dependent eNOS re-localization. 3. Peak distribution was similar with/without E2 with ≈ 55% of them in extragenic DNA regions and an intriguing involvement of the 5′ domain of several miRs deregulated in PCa. Numerous potentially novel eNOS-targeted genes have been identified suggesting that eNOS participates in the regulation of large gene sets. The parallel finding of downregulation of a cluster of miRs, including miR-34a, in PCa cells associated with poor outcome led us to unveil a molecular link between eNOS and SIRT1, an epigenetic regulator of aging and tumorigenicity, negatively regulated by miR-34a and in turn activating eNOS. E2 potentiates miR-34a downregulation thus enhancing SIRT1 expression, depicting a novel eNOS/SIRT1 interplay fine-tuned by E2-activated ER signaling, and suggesting that eNOS may play an important role in aggressive PCa.

Deep-sequencing of endothelial cells exposed to hypoxia reveals the complexity of known and novel microRNAs

In order to understand the role of microRNAs (miRNAs) in vascular physiopathology, we took advantage of deep-sequencing techniques to accurately and comprehensively profile the entire miRNA population expressed by endothelial cells exposed to hypoxia. SOLiD sequencing of small RNAs derived from human umbilical vein endothelial cells (HUVECs) exposed to 1% O₂ or normoxia for 24 h yielded more than 22 million reads per library. A customized bioinformatic pipeline identified more than 400 annotated microRNA/microRNA* species with a broad abundance range: miR-21 and miR-126 totaled almost 40% of all miRNAs. A complex repertoire of isomiRs was found, displaying also 5' variations, potentially affecting target recognition. High-stringency bioinformatic analysis identified microRNA candidates, whose predicted pre-miRNAs folded into a stable hairpin. Validation of a subset by qPCR identified 18 high-confidence novel miRNAs as detectable in independent HUVEC cultures and associated to the RISC complex. The expression of two novel miRNAs was significantly down-modulated by hypoxia, while miR-210 was significantly induced. Gene ontology analysis of their predicted targets revealed a significant association to hypoxia-inducible factor signaling, cardiovascular diseases, and cancer. Overexpression of the novel miRNAs in hypoxic endothelial cells affected cell growth and confirmed the biological relevance of their down-modulation. In conclusion, deep-sequencing accurately profiled known, variant, and novel microRNAs expressed by endothelial cells in normoxia and hypoxia.

RSSsite: a reference database and prediction tool for the identification of cryptic Recombination Signal Sequences in human and murine genomes

Recombination signal sequences (RSSs) flanking V, D and J gene segments are recognized and cut by the VDJ recombinase during development of B and T lymphocytes. All RSSs are composed of seven conserved nucleotides, followed by a spacer (containing either 12 ± 1 or 23 ± 1 poorly conserved nucleotides) and a conserved nonamer. Errors in V(D)J recombination, including cleavage of cryptic RSS outside the immunoglobulin and T cell receptor loci, are associated with oncogenic translocations observed in some lymphoid malignancies. We present in this paper the RSSsite web server, which is available from the address http://www.itb.cnr.it/rss. RSSsite consists of a web-accessible database, RSSdb, for the identification of pre-computed potential RSSs, and of the related search tool, DnaGrab, which allows the scoring of potential RSSs in user-supplied sequences. This latter algorithm makes use of probability models, which can be recasted to Bayesian network, taking into account correlations between groups of positions of a sequence, developed starting from specific reference sets of RSSs. In validation laboratory experiments, we selected 33 predicted cryptic RSSs (cRSSs) from 11 chromosomal regions outside the immunoglobulin and TCR loci for functional testing.

In silico prediction and experimental validation of natural antisense transcripts in two cancer-associated regions of human chromosome 6

Antisense transcription has long been recognized as a mechanism involved in the regulation of gene expression. Therefore, several human diseases associated with abnormal patterns of gene expression might display antisense RNA-mediated pathogenetic mechanisms. Such issue could be particularly relevant for cancer pathogenesis, since deregulated gene expression has long been established as a hallmark of cancer cells. Herein, we report on a bioinformatic search for antisense transcription in two cancer-associated regions of human chromosome 6 (6q21 and 6q27). Natural antisense transcripts (NATs) for several genes in both genomic regions were predicted in silico and subsequently validated by strand-specific RT-PCR. Detailed experimental validation by quantitative real-time RT-PCR of five putative cancer related sense-antisense transcript pairs revealed a single candidate tumor suppressor gene (RPS6KA2) whose expression levels display marked cancer-related changes that are likely mediated by its antisense RNA in a breast cancer cell line model.

An Ariadne's thread to the identification and annotation of noncoding RNAs in eukaryotes

Non-protein coding RNAs (ncRNAs) have emerged as a vast and heterogeneous portion of eukaryotic transcriptomes. Several ncRNA families, either short (<200 nucleotides, nt) or long (>200 nt), have been described and implicated in a variety of biological processes, from translation to gene expression regulation and nuclear trafficking. Most probably, other families are still to be discovered. Computational methods for ncRNA research require different approaches from the ones normally used in the prediction of protein-coding genes. Indeed, primary sequence alone is often insufficient to infer ncRNA functionality, whereas secondary structure and local conservation of portions of the transcript could provide useful information for both the prediction and the functional annotation of ncRNAs. Here we present an overview of computational methods and bioinformatics resources currently available for studying ncRNA genes, introducing the common themes as well as the different approaches required for long and short ncRNA identification and annotation.

A transcriptional sketch of a primary human breast cancer by 454 deep sequencing

Background

The cancer transcriptome is difficult to explore due to the heterogeneity of quantitative and qualitative changes in gene expression linked to the disease status. An increasing number of "unconventional" transcripts, such as novel isoforms, non-coding RNAs, somatic gene fusions and deletions have been associated with the tumoral state. Massively parallel sequencing techniques provide a framework for exploring the transcriptional complexity inherent to cancer with a limited laboratory and financial effort. We developed a deep sequencing and bioinformatics analysis protocol to investigate the molecular composition of a breast cancer poly(A)⁺ transcriptome. This method utilizes a cDNA library normalization step to diminish the representation of highly expressed transcripts and biology-oriented bioinformatic analyses to facilitate detection of rare and novel transcripts.

Results

We analyzed over 132,000 Roche 454 high-confidence deep sequencing reads from a primary human lobular breast cancer tissue specimen, and detected a range of unusual transcriptional events that were subsequently validated by RT-PCR in additional eight primary human breast cancer samples. We identified and validated one deletion, two novel ncRNAs (one intergenic and one intragenic), ten previously unknown or rare transcript isoforms and a novel gene fusion specific to a single primary tissue sample. We also explored the non-protein-coding portion of the breast cancer transcriptome, identifying thousands of novel non-coding transcripts and more than three hundred reads corresponding to the non-coding RNA MALAT1, which is highly expressed in many human carcinomas.

Conclusion

Our results demonstrate that combining 454 deep sequencing with a normalization step and careful bioinformatic analysis facilitates the discovery and quantification of rare transcripts or ncRNAs, and can be used as a qualitative tool to characterize transcriptome complexity, revealing many hitherto unknown transcripts, splice isoforms, gene fusion events and ncRNAs, even at a relatively low sequence sampling.

Non-random retention of protein-coding overlapping genes in Metazoa

Background

Although the overlap of transcriptional units occurs frequently in eukaryotic genomes, its evolutionary and biological significance remains largely unclear. Here we report a comparative analysis of overlaps between genes coding for well-annotated proteins in five metazoan genomes (human, mouse, zebrafish, fruit fly and worm).

Results

For all analyzed species the observed number of overlapping genes is always lower than expected assuming functional neutrality, suggesting that gene overlap is negatively selected. The comparison to the random distribution also shows that retained overlaps do not exhibit random features: antiparallel overlaps are significantly enriched, while overlaps lying on the same strand and those involving coding sequences are highly underrepresented. We confirm that overlap is mostly species-specific and provide evidence that it frequently originates through the acquisition of terminal, non-coding exons. Finally, we show that overlapping genes tend to be significantly co-expressed in a breast cancer cDNA library obtained by 454 deep sequencing, and that different overlap types display different patterns of reciprocal expression.

Conclusion

Our data suggest that overlap between protein-coding genes is selected against in Metazoa. However, when retained it may be used as a species-specific mechanism for the reciprocal regulation of neighboring genes. The tendency of overlaps to involve non-coding regions of the genes leads to the speculation that the advantages achieved by an overlapping arrangement may be optimized by evolving regulatory non-coding transcripts.

Splicy: a web-based tool for the prediction of possible alternative splicing events from Affymetrix probeset data

Background

The Affymetrix™ technology is nowadays a well-established method for the analysis of gene expression profiles in cancer research studies. However, changes in gene expression levels are not the only way to link genes and disease. The existence of gene isoforms specifically linked with cancer or apoptosis is increasingly found in literature. Hence it is of great interest to associate the results of a gene expression study with updated evidences on the transcript structure and its possible variants.

Results

We present here a web-based software tool, Splicy, whose primary task is to retrieve data on the mapping of Affymetrix™ probes to single exons of gene transcripts and displaying graphically this information projected on the gene physical structure.

Starting from a list of Affymetrix™ probesets the program produces a series of graphical displays, each relative to a transcript associated with the gene targeted by a given probe. The information on the transcript-by-transcript and exon-by-exon mapping of probe pairs can be retrieved both graphically and in the form of tab-separated files. The mapping of single probes to NCBI RefSeq or EMBL cDNAs is handled by the ISREC mapping tables used in the CleanEx Expression Reference Database Project. We currently maintain these mappings for most popular human and mouse Affymetrix™ chips, and Splicy can be queried for matches with human and mouse NCBI RefSeq or EMBL cDNAs.

Conclusion

Splicy generates probeset annotations and images describing the relation between the single probes and intron/exon structure of the target transcript in all its known variants. We think that Splicy will be useful for giving to the researcher a clearer picture of the possible transcript variants linked with a given gene and an additional view on the interpretation of microarray experiment data. Splicy is publicly available and has been realized in the framework of a bioinformatics grant from the Italian Cancer Research Association.

Identification of alternatively spliced dab1 isoforms in zebrafish

We have investigated the genomic organization, the occurrence of alternative splicing and the differential expression of the zebrafish disabled1 (dab1) gene. Dab1 is a key effector of the Reelin pathway, which regulates neuronal migration during brain development in vertebrates. The coding region of the zebrafish dab1 gene spans over 600 kb of genomic DNA and is composed of 15 exons. Alternative splicing in a region enriched for tyrosine residues generates at least three different isoforms. These isoforms are developmentally regulated and show differential tissue expression. Comparison with mouse and human data shows an overall conservation of the genomic organization with different alternative splicing events generating species-specific isoforms. Because these alternative splicing events give rise to isoforms with different numbers of phosphorylateable tyrosines, we speculate that alternative splicing of the dab1 gene in zebrafish and in other vertebrates regulates the nature of the cellular response to the Reelin signal.

NemaFootPrinter: a web based software for the identification of conserved non-coding genome sequence regions between C. elegans and C. briggsae

Background

NemaFootPrinter (Nematode Transcription Factor Scan Through Philogenetic Footprinting) is a web-based software for interactive identification of conserved, non-exonic DNA segments in the genomes of C. elegans and C. briggsae. It has been implemented according to the following project specifications:

a) Automated identification of orthologous gene pairs.

b) Interactive selection of the boundaries of the genes to be compared.

c) Pairwise sequence comparison with a range of different methods.

d) Identification of putative transcription factor binding sites on conserved, non-exonic DNA segments.

Results

Starting from a C. elegans or C. briggsae gene name or identifier, the software identifies the putative ortholog (if any), based on information derived from public nematode genome annotation databases. The investigator can then retrieve the genome DNA sequences of the two orthologous genes; visualize graphically the genes' intron/exon structure and the surrounding DNA regions; select, through an interactive graphical user interface, subsequences of the two gene regions. Using a bioinformatics toolbox (Blast2seq, Dotmatcher, Ssearch and connection to the rVista database) the investigator is able at the end of the procedure to identify and analyze significant sequences similarities, detecting the presence of transcription factor binding sites corresponding to the conserved segments. The software automatically masks exons.

Discussion

This software is intended as a practical and intuitive tool for the researchers interested in the identification of non-exonic conserved sequence segments between C. elegans and C. briggsae. These sequences may contain regulatory transcriptional elements since they are conserved between two related, but rapidly evolving genomes. This software also highlights the power of genome annotation databases when they are conceived as an open resource and the possibilities offered by seamless integration of different web services via the http protocol.

Availability: the program is freely available at

Overview of BITS2005, the Second Annual Meeting of the Italian Bioinformatics Society

The BITS2005 Conference brought together about 200 Italian scientists working in the field of Bioinformatics, students in Biology, Computer Science and Bioinformatics on March 17–19 2005, in Milan. This Editorial provides a brief overview of the Conference topics and introduces the peer-reviewed manuscripts accepted for publication in this Supplement.

DG-CST (Disease Gene Conserved Sequence Tags), a database of human-mouse conserved elements associated to disease genes

The identification and study of evolutionarily conserved genomic sequences that surround disease-related genes is a valuable tool to gain insight into the functional role of these genes and to better elucidate the pathogenetic mechanisms of disease. We created the DG-CST (Disease Gene Conserved Sequence Tags) database for the identification and detailed annotation of human–mouse conserved genomic sequences that are localized within or in the vicinity of human disease-related genes. CSTs are defined as sequences that show at least 70% identity between human and mouse over a length of at least 100 bp. The database contains CST data relative to over 1088 genes responsible for monogenetic human genetic diseases or involved in the susceptibility to multifactorial/polygenic diseases. DG-CST is accessible via the internet at http://dgcst.ceinge.unina.it/ and may be searched using both simple and complex queries. A graphic browser allows direct visualization of the CSTs and related annotations within the context of the relative gene and its transcripts.

The meaning of it all: web-based resources for large-scale functional annotation and visualization of DNA microarray data

The vast amount of unstructured data emerging from the various genome projects has led to the development of a number of web-based tools designed to annotate genes with biological information. Here we discuss a selection of these tools with regards to their scope, limitations and ease of use.

Microarray data analysis and mining

DNA microarray is an innovative technology for obtaining information on gene function. Because it is a high-throughput method, computational tools are essential in data analysis and mining to extract the knowledge from experimental results. Filtering procedures and statistical approaches are frequently combined to identify differentially expressed genes. However, obtaining a list of differentially expressed genes is only the starting point because an important step is the integration of differential expression profiles in a biological context, which is a hot topic in data mining. In this chapter an integrated approach of filtering and statistical validation to select trustable differentially expressed genes is described together with a brief introduction on data mining focusing on the classification of co-regulated genes on the basis of their biological function.

AntiHunter: searching BLAST output for EST antisense transcripts

AntiHunter is a new web-based tool for the identification of expressed sequence tag (EST) antisense transcripts from BLAST output. In order to perform an analysis, user is required to input a genomic sequence plus an associated list of transcript names and coordinates of the genomic region (i.e. genome annotation). After masking the repeated regions (if any), program will perform a BLASTN search of the input sequence versus the selected EST database, reporting by Email the EST entries that reveal a putative antisense transcript with respect to the user supplied list.

The DNA sequence quality machine at IFOM: a simple Web-based tool for quantitative assessment of sequencing reactions

DNA sequence quality is a factor of paramount importance in the world of modern genetic and genomics. Both the sequencing of Human Genome in the "post-draft" era [NHGRI Standard for quality of Human Genomic Sequences, Rev. 7 July (2002) where http://www.nhgri.nih.gov/Grant_info/Funding/ Statements/RFA/quality_standard.html is the HTTP address] and recent "high-throughput" approaches to genetic investigation such as SAGE [Velculescu, V.E., Zhang, L., Vogelstein, B. et al. (1995) "Serial analysis of gene expression", Science 270, 484-487] need a reliable, standardized measure of the quality of a sequencing reaction. The increasing importance of SNP studies also requires a stronger quality control on sequencing reactions by the final user. We propose here a simple, web-based tool for integrated sequence quality evaluation, high quality region quantitative value calculation and chromatogram display. This software is aimed at the small to medium DNA sequence laboratory or to the single researcher, interested in getting a quantitative measure of the sequence quality, browsing the chromatogram and checking the quality values base by base. The program is freely available from the IFOM bioinformatics web Server at http://bio.ifom-firc.it/Phred20/index.html.

Diagnostic role and prognostic significance of a simplified immunophenotypic classification of mature B cell chronic lymphoid leukemias

We verified the diagnostic and prognostic role of a simplified immunophenotypic classification (IC) in a series of 258 patients (M/F: 1.4; median age: 64 years; median follow-up: 64 months; 75 deaths) with mature B cell lymphoid leukemias (MBC-LL) for whom no histopathological diagnosis was available because of minimal or no lymph node involvement. The IC was based on the reactivity of three pivotal immunophenotypic markers: CD5, CD23 and SIg intensity. On the basis of different expression patterns, we identified four diagnostic clusters (C) characterized by distinct clinico-biological features and different prognoses: C1 (149 patients) identified most classical B cell chronic lymphocytic leukemias (CLL-type cluster; SIg(dim)/CD5+/CD23+); C2, 38 patients whose clinico-hematological characteristics were intermediate between C1 and C3 (CLL-variant cluster; SIg(bright)/CD5+/CD23+/-or SIg(dim)/CD5-/-/CD23 indifferent); C3 (16 patients) most situations consistent with mantle cell lymphoma in leukemic phase (MCL-type cluster; SIg(bright)/CD5+/CD23-); and C4, 55 cases, most of whom were consistent with leukemic phase lymphoplasmacytic/splenic marginal zone lymphomas (LP/S-type cluster; SIg(bright)/CD5-/+/CD23 indifferent). At univariate survival analysis, prognosis worsened from C1 to C4, C2 and C3 (P = 0.0001), and this was maintained at multivariate analysis (P = 0.006), together with CD11c expression (P = 0.0043), age at diagnosis (cut-off 70 years; P = 0.0008) and platelet count (cut-off 140 x 10(9)/l; P = 0.0034). Besides recognising the two well-known situations of classic B-CLL and MCL, our IC identified situations with distinct prognostic and/or clinical behaviors.

Automated DNA chip annotation tables at IFOM: the importance of synchronisation and cross-referencing of sequence databases

The increasing popularity of DNA chip technology for the study of gene expression is producing, for each experiment, a sizable quantity of numerical data to analyse and an accompanying large number of gene identifiers that should be associated with the relevant biological annotation. We describe here a website at IFOM (FIRC Institute of Molecular Oncology) where we release regularly updated annotation tables for the most used Affymetrix oligonucleotide DNA chips and for the whole Research Genetics 46K clone collection for cDNA arrays. These tables are synchronised with every new release of the mouse and human UniGene databases (NCBI; National Center for Biotechnology Information), allowing fast and easy preliminary annotation of DNA array experiments. We also report some comparative evidence about the importance of biological database synchronisation and cross-references in the process of generating annotation tables for DNA chips.

The meaning of it all: web-based resources for large-scale functional annotation and visualization of DNA microarray data

The vast amount of unstructured data emerging from the various genome projects has led to the development of a number of web-based tools designed to annotate genes with biological information. Here we discuss a selection of these tools with regards to their scope, limitations and ease of use.

A global analysis of Caenorhabditis elegans operons

The nematode worm Caenorhabditis elegans and its relatives are unique among animals in having operons. Operons are regulated multigene transcription units, in which polycistronic pre-messenger RNA (pre-mRNA coding for multiple peptides) is processed to monocistronic mRNAs. This occurs by 3' end formation and trans-splicing using the specialized SL2 small nuclear ribonucleoprotein particle for downstream mRNAs. Previously, the correlation between downstream location in an operon and SL2 trans-splicing has been strong, but anecdotal. Although only 28 operons have been reported, the complete sequence of the C. elegans genome reveals numerous gene clusters. To determine how many of these clusters represent operons, we probed full-genome microarrays for SL2-containing mRNAs. We found significant enrichment for about 1,200 genes, including most of a group of several hundred genes represented by complementary DNAs that contain SL2 sequence. Analysis of their genomic arrangements indicates that >90% are downstream genes, falling in 790 distinct operons. Our evidence indicates that the genome contains at least 1,000 operons, 2 8 genes long, that contain about 15% of all C. elegans genes. Numerous examples of co-transcription of genes encoding functionally related proteins are evident. Inspection of the operon list should reveal previously unknown functional relationships.

Modeling molecular networks: a systems biology approach to gene function

A report on the European Science Foundation Workshop on Modeling of Molecular Networks, Granada, Spain, 11-14 June 2002.

The tripartite motif family identifies cell compartments

A functional genomic approach, based on systematic data gathering, was used to characterize a family of proteins containing a tripartite motif (TRIM). A total of 37 TRIM genes/proteins were studied, 21 of which were novel. The results demonstrate that TRIM proteins share a common function: by means of homo-multimerization they identify specific cell compartments.

OPP-EBV-CAD regimen as salvage treatment in advanced refractory or resistant multiple myeloma

With the aim of developing an effective therapy for heavily pretreated refractory MM outpatients, we evaluated the OPPEBVCAD regimen, a Hodgkin's disease-derived protocol that includes many drugs effective in MM administered in a sequential schedule. Twenty-two pts aged 42-72 years, with symptomatic highly-pretreated refractory (18 cases), or primary resistant MM (four cases. including two pts with plasma cell leukemia-PCL) received this therapy every 28 days (2-4 cycles, followed by a maintenance program). Therapeutic response (Chronic Leukemia-Myeloma Task Force criteria) and performance status (PS) and pain (W.H.O.) were evaluated. All of the pts were evaluable for response. There were 9 (40%) objective responses (OR: stabilization of blood counts and bone lesions, serum calcium normalization, 50% or more reduction in the concentration of serum monoclonal component (MC), 90% reduction in Bence-Jones proteinuria), 8 (36%) partial responses (PR: 25-50% reduction in serum MC), 1 no response or stable disease (NR), and 4 (18%) cases of progressive disease (PD). OR plus PR were 77%. Of the 4 primary resistant tumors (2 PCL and 2 MM), 2 achieved PR, 1 OR (a PCL case) and 1 progressed. Median survival was 15 months for responding pts (OR plus PR) and 4.5 months for non-responders (NR plus PD). PS and pain improved in 15 pts and did not change in 9. The most frequent side effects were cytopenias, with one drug related infective death. The OPPEBVCAD regimen proved to be an effective therapy for refractory relapsing or primary resistant MM: in responders (two-thirds of the pts), survival was prolonged by about 10 months. Its efficacy in anthracycline-treated pts, as well as the feasibility of using it on an outpatient basis without any continuous drug infusions, make this regimen a promising third line salvage therapy.

Soluble cytokine levels correlate with the activity and clinical stage of Hodgkin's disease at diagnosis

The serum levels of some cytokines seem to correlate with outcome in Hodgkin's disease (HD) and may be helpful in formulating new and better prognostic systems. The aim of this study was to analyse the correlations between the serum levels of different cytokines and the clinico-hematological features suggestive of a worse prognosis. The study involved 31 pts with a "de novo" diagnosis of HD (median age: 30 yrs; M/F: 13/18; stage I/II vs III/IV: 19/12; B symptoms: 12; bulky disease and extranodal disease: 9). The serum levels of sCD30, TNFalpha, TNF receptor I and II, IL6, IL6 receptor, IL10, sICAM-1 were evaluated at diagnosis, and correlated with gender, age (< or =/> 30), stage (I-II vs III-IV), systemic symptoms, bulky disease, ESR (</> or = 40), serum copper (< or =/> 170 microg/dL), WBC counts (< or =/> 15x10(9)/L), prognostic scores (PS) according to Hasenclever (</> or = 2), and the presence of extranodal disease. Stages III/IV were associated with significantly higher levels of sCD30 and TNF-RI (p=0.03), systemic symptoms with significantly higher levels of sCD30, TNFalpha, IL6, TNF-RI (p=0.027, 0.03, 0.0005, 0.002), bulky disease with TNF-RI (p=0.03), high ESR with IL6 and TNF-RI (p=0.0011, 0.0001), high WBC counts with sCD30, IL6, TNF-RI (p=0.03, 0.002, 0.01), high serum copper with sCD30 and IL6 (p=0.05, 0.0004), higher PS with sCD30, IL6, TNF-RI (p=0.002, 0.0003, 0.005), extranodal disease with TNFalpha and IL6 (p=0.05, 0.01). It was possible to define cut-off levels for some cytokines (sCD30 > 33.15 U/mL, TNFalpha > 29.71 pg/mL, IL6 > 12.43 pg/mL, TNF-RI > 3.23 ng/mL, IL6-R > 57 ng/mL) that significantly correlate with systemic symptoms, higher disease stages, ESR, serum copper, WBC counts and PS. Our study shows that high sCD30, TNFalpha, IL6 and TNF-RI levels are associated with advanced disease or a worse prognostic score. In the context of multiparametric HD staging, cytokine evaluation may be useful for identifying candidates for more intensive therapies.

Role of retinoic acid receptor overexpression in sensitivity to fenretinide and tumorigenicity of human ovarian carcinoma cells

The role of retinoic acid receptor (RAR) expression in sensitivity to N-(4-hydroxyphenyl)retinamide (4HPR or fenretinide) as well as on the tumorigenicity of human ovarian carcinoma cells was examined. Two human ovarian cancer cell lines, A2780 and IGROV-1, with a 10-fold difference in sensitivity to 4HPR were chosen to study RAR involvement in the response to 4HPR. To determine which RAR was effective, RARalpha, beta and gamma were individually overexpressed in A2780 cells, which are the most sensitive to 4HPR. Sensitivity to 4HPR was increased in RARbeta-overexpressing clones, whereas it was slightly decreased in RARalpha transfectants (which had diminished RARbeta expression) and was unchanged in clones transfected with RARgamma. IGROV-1 cells, which are RARbeta negative, were transfected with RARbeta. Surprisingly, none of the obtained IGROV-1 RARbeta transfectants expressed RARbeta protein, in spite of RARbeta mRNA transcription. All clones were similar to the parental IGROV-1 cells in their sensitivity to 4HPR. Treatment with a pharmacologically achievable concentration of 4HPR (1 microM) led to a rapid 2-fold increase in RARbeta mRNA levels in A2780 cells, but it did not induce RARbeta expression in IGROV-1 cells. Analysis of the tumorigenicity of A2780-transfected clones revealed that overexpression of RARalpha was associated with a significant reduction in tumor takes (50% and 67%, respectively, vs. 96% for the parent line) and with a reduced growth rate. Oncogenicity was clearly decreased in only 1 of the 2 RARbeta-overexpressing clones (33% takes) and was unchanged in the 2 clones with increased RARgamma expression. Our results demonstrate that basal expression and 4HPR inducibility of RARbeta play a role in mediating 4HPR response in ovarian cancer cells. The findings of reduced oncogenicity of clones overexpressing RARalpha and of one clone overexpressing RARbeta indicate that RARalpha and RARbeta might have a tumor-suppressive effect in ovarian tumors.

TargetFinder: searching annotated sequence databases for target genes of transcription factors

TargetFinder is a new software tool to search a database of annotated sequences for transcription factor binding sites located in context with other important transcription regulatory signals and regions, like the TATA element, the promoter, and so on, thereby greatly reducing the background usually associated with this kind of search.

AVAILABILITY

The TargetFinder Web service is available at http://hercules.tigem.it/TargetFinder.html

CONTACT

giovanni.lavorgna@hsr.it

Use of IgM monoclonal reagents licensed for tube tests in column agglutination technology

BACKGROUND

Red cell (RBC) phenotyping using column agglutination technology (CAT) is currently limited by the reagents formulated in the system. To overcome this limitation, it was investigated whether monoclonal IgM reagents licensed for use with tube tests produced valid results with CAT.

STUDY DESIGN AND METHODS

Commercial CAT, does not contain antisera, was used to evaluate Procedures A (40 microL of reagent and 10 microL of 4% RBCs) and B (50 microL of reagent and 50 microL of 0.8% RBCs) with or without incubation at room temperature. In Study 1, reagents were tested to determine whether potentiators inhibit the passage of antigen-negative RBCs through the column. In Study 2, CAT sensitivity was measured by the use of potency titrations to define a procedure for each reagent that matched or exceeded that of the tube method. In Study 3, the specificity of each reagent was determined in parallel with the CAT and tube tests. Typing of 1644 samples was performed.

RESULTS

Study 1: Free passage was obtained with all reagents. Study 2: Immediate-spin methods using CAT produced the same results as the tube method. Study 3: With 8048 comparisons made, discrepant results were found in 32 transfused patients and in 6 cord blood samples, mainly with Lewis reagents. With comparison of CAT and the standard tube method, complete agreement was obtained with Kell reagents, 99.9-percent agreement with Kidd reagents, and 98.9-percent and 99.4-percent agreement with Lewis reagents.

CONCLUSION

Most examined reagents seem suitable for use with CAT.

Treatment and prognosis in a series of primary extranodal lymphomas of the ocular adnexa

BACKGROUND

The aim of this study was to assess clinicopathological characteristics and outcome in a series of primary ocular adnexal lymphomas (POALs).

PATIENTS AND METHODS

Nineteen patients with localised (stage IE) POAL were followed for a median of 96 months (24-156). The diagnosis was based on surgical biopsies followed by immunohistochemistry in 16 cases or fine-needle aspiration followed by immunocytophenotypic analysis in three cases. Twelve patients were treated with local radiotherapy (RT), five with chemotherapy (CT), and two refused further therapy after apparently radical tumour removal achieved by the diagnostic excisional biopsy.

RESULTS

The histological and immunological pattern was consistent with a diagnosis of MALT-type lymphoma (11 cases), follicular center non-Hodgkin's lymphoma (three cases). a large-cell variant of Burkitt's lymphoma (one case), and large-cell transformed MALT lymphoma (one case). Low-grade lymphoma was diagnosed in the three cases which underwent fine-needle aspiration biopsy. All of the patients achieved and maintained complete remission except for those treated with surgical excision alone (two MALT conjunctival lymphoma cases): one of these relapsed locally, the other experienced the systemic spread of a transformed diffuse large-cell lymphoma and died 72 months after diagnosis. The side effects consisted of two cases of RT-related cataract after 52 and 72 months.

CONCLUSIONS

Regardless of histology, prognosis was excellent when surgery plus RT was adopted, and CT seems to be a valid alternative to RT. Surgery alone may be sub-optimal.

Poor prognosis in non-villous splenic marginal zone cell lymphoma is associated with p53 mutations

We have recently reported a series of 15 non-villous splenic marginal zone lymphoma patients, six of whom showed p53 mutations (40%). This molecular alteration did not correlate with any particular clinico-pathologic feature at diagnosis. After a median follow-up of 56 months, four cases evolved into aggressive fatal non-Hodgkin's lymphoma (NHL) and two had refractory progressive disease; interestingly, p53 mutations were demonstrated in five of these patients at diagnosis. As the patients with wild-type p53 presented responsive or indolent disease, this genetic alteration may be an early marker of aggressive transformation or refractoriness. p53 evaluation at diagnosis could be advisable in this particular subset of NHL.