LncRNA co-expression network analysis reveals novel biomarkers for pancreatic cancer.

High mortality and low survival rates for pancreatic ductal adenocarcinoma (PDAC) mainly result from the delay in diagnosis and treatment. Therefore, there is an urgent need to identify early PDAC biomarkers and new therapeutic targets. In this study, we applied a commonly used systems biology approach, the weighted gene co-expression network analysis (WGCNA), on lncRNA expression data. Eleven lncRNAs, namely A2M-AS1, DLEU2, LINC01133, LINC00675, MIR155HG, SLC25A25-AS1, LINC01857, LOC642852 (LINC00205), ITGB2-AS1, TSPOAP1-AS1 and PSMB8-AS1 have been identified and validated on an independent PDAC expression dataset. Furthermore, we characterized them by functional and pathway enrichment analysis and identified which lncRNAs showed differential expression, differential promoter methylation levels and copy number alterations between normal and PDAC samples. Finally, we also performed a survival analysis and identified A2M-AS1, LINC01133, LINC00205 and TSPOAP1-AS1 as prognostic biomarkers for PDAC. Interestingly, although only a few cancer-associated lncRNAs have been functionally characterized, LINC00675 and LINC01133 lncRNAs have already been demonstrated to be involved in PDAC development and progression. Therefore, our results provide new potential diagnostic/prognostic biomarkers and therapeutic targets for PDAC that deserve to be further investigated. Moreover, these lncRNAs may improve the understanding about molecular pathogenesis of PDAC.

ExportAid: database of RNA elements regulating nuclear RNA export in mammals.

MOTIVATION: Regulation of nuclear mRNA export or retention is carried out by RNA elements but the mechanism is not yet well understood. To understand the mRNA export process, it is important to collect all the involved RNA elements and their trans-acting factors. RESULTS: By hand-curated literature screening we collected, in ExportAid database, experimentally assessed data about RNA elements regulating nuclear export or retention of endogenous, heterologous or artificial RNAs in mammalian cells. This database could help to understand the RNA export language and to study the possible export efficiency alterations owing to mutations or polymorphisms. Currently, ExportAid stores 235 and 96 RNA elements, respectively, increasing and decreasing export efficiency, and 98 neutral assessed sequences. AVAILABILITY AND IMPLEMENTATION: Freely accessible without registration at http://www.introni.it/ExportAid/ExportAid.html. Database and web interface are implemented in Perl, MySQL, Apache and JavaScript with all major browsers supported.

A possible S-glutathionylation of specific proteins by glyoxalase II: An in vitro and in silico study.

Glyoxalase II, the second of 2 enzymes in the glyoxalase system, is a hydroxyacylglutathione hydrolase that catalyses the hydrolysis of S-d-lactoylglutathione to form d-lactic acid and glutathione, which is released from the active site. The tripeptide glutathione is the major sulfhydryl antioxidant and has been shown to control several functions, including S-glutathionylation of proteins. S-Glutathionylation is a way for the cells to store reduced glutathione during oxidative stress, or to protect protein thiol groups from irreversible oxidation, and few enzymes involved in protein S-glutathionylation have been found to date. In this work, the enzyme glyoxalase II and its substrate S-d-lactoylglutathione were incubated with malate dehydrogenase or with actin, resulting in a glutathionylation reaction. Glyoxalase II was also submitted to docking studies. Computational data presented a high propensity of the enzyme to interact with malate dehydrogenase or actin through its catalytic site and further in silico investigation showed a high folding stability of glyoxalase II toward its own reaction product glutathione both protonated and unprotonated. This study suggests that glyoxalase II, through a specific interaction of its catalytic site with target proteins, could be able to perform a rapid and specific protein S-glutathionylation using its natural substrate S-d-lactoylglutathione. SIGNIFICANCE: This article reports for the first time a possible additional role of Glo2 that, after interacting with a target protein, is able to promote S-glutathionylation using its natural substrate SLG, a glutathione derived compound. In this perspective, Glo2 can play a new important regulatory role inS-glutathionylation, acquiring further significance in cellular post-translational modifications of proteins.

KRAS mutation status is associated with specific pattern of genes expression in pancreatic adenocarcinoma.

AIMS: To evaluate potential differences at a molecular level between KRAS mutant tumors (MT) and KRAS wild-type (WT) pancreatic tumors and the biological and prognostic significance of different KRAS mutations. MATERIALS & METHODS: Expression of a panel of 29 genes was analyzed in KRAS WT and MT tumors. Effects of KRAS mutation and gene expression levels were assessed on patients' survival. RESULTS: MUC6 (p = 0.009), HGF (p = 0.011), VEGFR-2 (p = 0.020) and VEGFB (p = 0.026) were significantly more expressed and SMAD4 was less suppressed (p = 0.003) in WT KRAS. Contrariwise, SHH (p = 0.012) and IHH (p = 0.031) were more expressed in MT KRAS patients. No OS difference was found between WT and MT KRAS tumors. CONCLUSION: KRAS mutation status seems to identify two different subtypes of pancreatic ductal adenocarcinoma with similar outcome but distinct molecular features and probably different therapeutic targets.

Lgr5 expression, cancer stem cells and pancreatic cancer: results from biological and computational analyses.

AIMS: To determine the relationship between Lgr5 and other stemness markers and pathologic features in pancreatic ductal adenocarcinoma (PDAC) samples. MATERIALS & METHODS: In 69 samples, Lgr5 was analyzed by qRT-PCR together with a panel of 29 genes. Bioinformatic analysis was carried out to identify a possible pathway regulating Lgr5 expression in PDAC. RESULTS: Lgr5 expression was not associated with the expression of tested cancer stem cell markers. Moreover, it was not an independent predictor of survival neither at univariate analysis (p = 0.21) nor at multivariate analysis (p = 0.225). CONCLUSION: Based on the lack of correlation between Lgr5 and tested cancer stem cell markers, Lgr5 does not seem to be a potential stemness marker or prognostic factor in PDAC.

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.

SpliceAid 2: a database of human splicing factors expression data and RNA target motifs.

Splicing is the most frequently altered biological process by mutations within gene regions. Information for splicing is recognized by several factors that bind pre-mRNA sequence and, through coordinated interaction, yield mature transcripts. Some in silico methods have been developed to predict if a mutation leads to aberrant splicing patterns. We previously created SpliceAid tool that is able to minimize false positive predictions because it adopts strictly experimental RNA target motifs bound by splicing proteins in humans. In order to improve prediction accuracy and better understand the splicing outcome, the tissue specificity of each splicing regulatory factor has to be taken into account. Here, we have developed SpliceAid 2 by adding the expression data related to the splicing factors extracted from the main proteomic and transcriptomic databases, true 5' and 3' splice sites, polypyrimidine tracts, and branch point sequences. The new version collects 2,220 target sites of 62 human splicing proteins and their expression data in 320 tissues per cell. SpliceAid 2 can be useful to foresee the splicing pattern alteration, to guide the identification of the molecular effect due to the mutations and to understand the tissue-specific alternative splicing. SpliceAid 2 is freely accessible at www.introni.it/spliceaid.html.

ADH4 intronic variations are associated with alcohol dependence: results from an Italian case-control association study.

OBJECTIVES: This study investigated the involvement of ADH4 gene polymorphisms in the susceptibility to alcohol use disorders. METHODS: Thirty-eight single-nucleotide polymorphisms (SNPs) in and around the ADH4 gene were investigated in 136 Italian alcoholics and 276 healthy controls. A new approach based on a bioinformatic method selected 26 SNPs that may affect the splicing sites, destroying or creating binding sites of splicing regulatory proteins. RESULTS: Case-control comparisons for allele and genotype frequencies showed that ADH4 SNPs were associated with alcohol dependence but not with alcohol abuse. The association signal was strongest for rs1009145, rs13148577 (both P=0.0008) and rs7689753 (P=0.0007), whose minor alleles were predicted to alter the target protein sequences involved in mRNA splicing. A pairwise linkage disequilibrium analysis showed that all SNPs except five were located in a single haplotype block. Six haplotype tag SNPs were selected to infer haplotypes and to estimate their frequency distributions. A logistic regression analysis confirmed the association between ADH4 variants and alcohol dependence when sex, age, years of education, marital status and the allele genotype, haplotype and diplotype data of the six haplotype tag SNP were considered. Haplotype ATAAAT, which contained the minor allele of rs10009145 and the major allele of rs7689753, increased the risk of alcohol dependence, whereas haplotype GGGGAT, bearing the major allele of rs10009145 and the minor allele of rs7689753, protected against it. Again, there was no evidence of an association with alcohol abuse. CONCLUSION: These data suggest that ADH4 intronic variants play a role in alcohol dependence susceptibility in Italian populations. Functional studies are needed to establish the role of the genetic variations that seem to affect the splicing mechanism.

TRAM (Transcriptome Mapper): database-driven creation and analysis of transcriptome maps from multiple sources.

BACKGROUND: Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input. RESULTS: TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified. CONCLUSIONS: TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.

Bioinformatic analyses to select phenotype affecting polymorphisms in HTR2C gene.

OBJECTIVE: Single nucleotide polymorphisms (SNPs) in serotonin related genes influence mental disorders, responses to pharmacological and psychotherapeutic treatments. In planning association studies, researchers that want to investigate new SNPs have to select some among a large number of candidates. Our aim is to guide researchers in the selection of the most likely phenotype affecting polymorphisms. Here, we studied serotonin receptor 2C (HTR2C) SNPs because, till now, only relatively few of about 2000 are investigated. METHODS: We used the most updated and assessed bioinformatic tools to predict which variations can give rise to biological effects among 2450 HTR2C SNPs. RESULTS: We suggest 48 SNPs that are worth considering in future association studies in the field of psychiatry, psychology and pharmacogenomics. Moreover, our analyses point out the biological level probably affected, such as transcription, splicing, miRNA regulation and protein structure, thus allowing to suggest future molecular investigations. CONCLUSIONS: Although few association studies are available in literature, their results are in agreement with our predictions, showing that our selection methods can help to guide future association studies.

SpliceAid: a database of experimental RNA target motifs bound by splicing proteins in humans.

UNLABELLED: The correct post-transcriptional RNA processing is finely regulated by RNA-binding proteins. Unfortunately, there is little experimental information on target RNA sequences of RNA-binding proteins and moreover such experimentally derived target sequences are annotated in a compact form by the score matrices that overestimate the number of possible recognized sequences. We carried out an exhaustive hand curated literature search to create a database, SpliceAid, collecting all the experimentally assessed target RNA sequences that are bound by splicing proteins in humans. We built a web resource, database driven, to easy query SpliceAid and give back the results by an accurate and dynamic graphic representation. AVAILABILITY: SpliceAid database is freely accessible at http://www.introni.it/splicing.html.

An improved in silico selection of phenotype affecting polymorphisms in SLC6A4, HTR1A and HTR2A genes.

OBJECTIVE: Among the experimentally assessed DNA variations in serotonin related genes, some influence physiological expression of personality and mental disorders, others alter the responses to pharmacological and/or psychotherapeutic treatments. Because of the huge number of polymorphisms lying in genes and of the great length of time necessary to perform association studies, a selection of the variations being studied is a necessary and crucial step. METHODS: In this work we used the most updated and assessed bioinformatic tools to predict the phenotype affecting polymorphisms of the human HTR1A, HTR2A and SLC6A4 serotonin related genes. Moreover, we carried out a literature search to collect information about the recent association studies to compare it versus our prediction data. RESULTS: Gene polymorphism analysis indicated the variations that are worth considering in the association studies in the field of psychiatry, psychology and pharmacogenomics. The literature revision allowed to show both the few well and the most not enough investigated polymorphisms. CONCLUSIONS: Our data can be useful to select polymorphisms for new association studies, especially those not yet investigated that can be related to behaviour, mental disorders and individual treatment response.

Sequential or Concomitant Inhibition of Cyclin-Dependent Kinase 4/6 Before mTOR Pathway in Hormone-Positive HER2 Negative Breast Cancer: Biological Insights and Clinical Implications.

About 75% of all breast cancers are hormone receptor-positive (HR+). However, the efficacy of endocrine therapy is limited due to the high rate of either pre-existing or acquired resistance. In this work we reconstructed the pathways around estrogen receptor (ER), mTOR, and cyclin D in order to compare the effects of CDK4/6 and PI3K/AKT/mTOR inhibitors. A positive feedback loop links mTOR and ER that support each other. We subsequently considered whether a combined or sequential inhibition of CDK4/6 and PI3K/AKT/mTOR could ensure better results. Studies indicate that inhibition of CDK4/6 activates mTOR as an escape mechanism to ensure cell proliferation. In literature, the little evidence dealing with this topic suggests that pre-treatment with mTOR pathway inhibitors could prevent or delay the onset of CDK4/6 inhibitor resistance. Additional studies are needed in order to find biomarkers that can identify patients who will develop this resistance and in whom the sensitivity to CDK4/6 inhibitors can be restored.

CXCL12 and Its Isoforms: Different Roles in Pancreatic Cancer?

CXCL12 is a chemokine that acts through CXCR4 and ACKR3 receptors and plays a physiological role in embryogenesis and haematopoiesis. It has an important role also in tumor development, since it is released by stromal cells of tumor microenvironment and alters the behavior of cancer cells. Many studies investigated the roles of CXCL12 in order to understand if it has an anti- or protumor role. In particular, it seems to promote tumor invasion, proliferation, angiogenesis, epithelial to mesenchymal transition (EMT), and metastasis in pancreatic cancer. Nevertheless, some evidence shows opposite functions; therefore research on CXCL12 is still ongoing. These discrepancies could be due to the presence of at least six CXCL12 splicing isoforms, each with different roles. Interestingly, three out of six variants have the highest levels of expression in the pancreas. Here, we report the current knowledge about the functions of this chemokine and then focus on pancreatic cancer. Moreover, we discuss the methods applied in recent studies in order to understand if they took into account the existence of the CXCL12 isoforms.

Different Cardiotoxicity of Palbociclib and Ribociclib in Breast Cancer: Gene Expression and Pharmacological Data Analyses, Biological Basis, and Therapeutic Implications.

Breast cancer is the most frequent tumor in women. The recent advent of cyclin-dependent kinase (CDK) 4/6 inhibitors palbociclib and ribociclib has represented a major step forward for patients with hormone receptor-positive breast cancer. These two agents have showed similar efficacy in terms of breast cancer outcome but different cardiotoxic effects. In particular, ribociclib, but not palbociclib, has been associated with QT interval prolongation, and the mechanisms underlying this event are still unclear. In order to clarify such difference, we matched the candidate genes associated with QT interval prolongation with genes whose expression is altered following palbociclib or ribociclib treatment. We also investigated whether pharmacokinetic and pharmacodynamic characteristics, such as IC50 (hERG) [concentration of drug producing 50% inhibition (human ether-à-go-go related gene)] and maximum concentration (Cmax), could justify the different effects on QT interval prolongation. Our results show that ribociclib, but not palbociclib, could act by down-regulating the expression of KCNH2 (encoding for potassium channel hERG) and up-regulating SCN5A and SNTA1 (encoding for sodium channels Nav1.5 and syntrophin-α1, respectively), three genes associated with long QT syndrome. Consistent with the cardiotoxicity induced by ribociclib, its IC50 (hERG)/free concentration (Cmax free) ratio is closer to the safety threshold than that of palbociclib. In summary, we hypothesize that the different cardiotoxicity associated with ribociclib and palbociclib could be due to the alteration of potassium and sodium channels induced by ribociclib. A better comprehension of the mechanisms of cardiac channelopathies and drug-induced QT interval prolongation will be fundamental to avoid serious and potentially lethal adverse events and, as a consequence, optimize the management of breast cancer patients.

Glutathione compartmentalization and its role in glutathionylation and other regulatory processes of cellular pathways.

Glutathione is considered the major non-protein low molecular weight modulator of redox processes and the most important thiol reducing agent of the cell. The biosynthesis of glutathione occurs in the cytosol from its constituent amino acids, but this tripeptide is also present in the most important cellular districts, such as mitochondria, nucleus, and endoplasmic reticulum, thus playing a central role in several metabolic pathways and cytoprotection mechanisms. Indeed, glutathione is involved in the modulation of various cellular processes and, not by chance, it is a ubiquitous determinant for redox signaling, xenobiotic detoxification, and regulation of cell cycle and death programs. The balance between its concentration and redox state is due to a complex series of interactions between biosynthesis, utilization, degradation, and transport. All these factors are of great importance to understand the significance of cellular redox balance and its relationship with physiological responses and pathological conditions. The purpose of this review is to give an overview on glutathione cellular compartmentalization. Information on its subcellular distribution provides a deeper understanding of glutathione-dependent processes and reflects the importance of compartmentalization in the regulation of specific cellular pathways. © 2018 BioFactors, 45(2):152-168, 2019.

Clinical impact of different exosomes' protein expression in pancreatic ductal carcinoma patients treated with standard first line palliative chemotherapy.

INTRODUCTION: Pancreatic ductal adenocarcinoma is associated to dismal prognosis despite the use of palliative chemotherapy, partly due to the lack of knowledge of biological processes underlying disease progression. Exosomes have been identified as biomarkers sources in different cancer types. Aim of the study was to analyse the contents of circulating exosomes in patients with pancreatic cancer who received palliative chemotherapy. PATIENTS AND METHODS: Patients were submitted to blood sample collection before chemotherapy (T0) and after 3 months (T3). We quantified by an ELISA-based technique specific proteins of cancer-derived exosomes (CD44,CD44v6,EpCAM,CD9,CD81,Tspan8,Integrin α6,Integrin ß4,CD24,CXCR4). We correlated the baseline levels of these factors and changes between T3 and T0 and survival outcomes. Survival analyses were performed by Kaplan-Meier method. Correlation was assessed by log-rank test and level of statistical significance was set at 0.05. Multivariate analysis was performed by logistic regression analysis. RESULTS: Nineteen patients were enrolled. EpCAM T0 levels and increased EpCAM levels from T0 to T3 were those mostly associated with differences in survival. Patients having higher EpCAM had median progression free survival (PFS) of 3.18vs7.31 months (HR:2.82,95%CI:1.03-7.73,p = 0.01). Overall survival (OS) was shorter for patients having higher EpCAM (5.83vs16.45 months,HR:6.16,95%CI:1.93-19.58,p = 0.0001) and also response rates (RR) were worse (20%vs87%,p = 0.015). EpCAM increase during treatment was associated with better median PFS (2.88vs7.31 months,HR:0.24,95%CI:0.04-1.22,p = 0.003). OS was also better (8.75vs11.04 months, HR:0.77,95%CI:0.21-2.73,p = 0.66) and RR were 60%vs20% (p = 0.28). Among clinical factors that might determine changes on PFS and OS, only ECOG PS was associated to significantly worse PFS and OS (p = 0.0137and<0.001 respectively).Multivariate analysis confirmed EpCAM T0 levels and EpCAM T0/T3 changes as independent prognostic factors for PFS. CONCLUSIONS: Pancreatic cancer patients exosomes express EpCAM, whose levels change during treatment. This represents a useful prognostic factor and also suggests that future treatment modalities who target EpCAM should be tested in pancreatic cancer patients selected by exosome EpCAM expression.

KRIT1 Loss-Of-Function Associated with Cerebral Cavernous Malformation Disease Leads to Enhanced S-Glutathionylation of Distinct Structural and Regulatory Proteins.

Loss-of-function mutations in the KRIT1 gene are associated with the pathogenesis of cerebral cavernous malformations (CCMs), a major cerebrovascular disease still awaiting therapies. Accumulating evidence demonstrates that KRIT1 plays an important role in major redox-sensitive mechanisms, including transcriptional pathways and autophagy, which play major roles in cellular homeostasis and defense against oxidative stress, raising the possibility that KRIT1 loss has pleiotropic effects on multiple redox-sensitive systems. Using previously established cellular models, we found that KRIT1 loss-of-function affects the glutathione (GSH) redox system, causing a significant decrease in total GSH levels and increase in oxidized glutathione disulfide (GSSG), with a consequent deficit in the GSH/GSSG redox ratio and GSH-mediated antioxidant capacity. Redox proteomic analyses showed that these effects are associated with increased S-glutathionylation of distinct proteins involved in adaptive responses to oxidative stress, including redox-sensitive chaperonins, metabolic enzymes, and cytoskeletal proteins, suggesting a novel molecular signature of KRIT1 loss-of-function. Besides providing further insights into the emerging pleiotropic functions of KRIT1, these findings point definitively to KRIT1 as a major player in redox biology, shedding new light on the mechanistic relationship between KRIT1 loss-of-function and enhanced cell sensitivity to oxidative stress, which may eventually lead to cellular dysfunctions and CCM disease pathogenesis.

Emerging Biomarkers in Bladder Cancer Identified by Network Analysis of Transcriptomic Data.

Bladder cancer is a very common malignancy. Although new treatment strategies have been developed, the identification of new therapeutic targets and reliable diagnostic/prognostic biomarkers for bladder cancer remains a priority. Generally, they are found among differentially expressed genes between patients and healthy subjects or among patients with different tumor stages. However, the classical approach includes processing these data taking into consideration only the expression of each single gene regardless of the expression of other genes. These complex gene interaction networks can be revealed by a recently developed systems biology approach called Weighted Gene Co-expression Network Analysis (WGCNA). It takes into account the expression of all genes assessed in an experiment in order to reveal the clusters of co-expressed genes (modules) that, very probably, are also co-regulated. If some genes are co-expressed in controls but not in pathological samples, it can be hypothesized that a regulatory mechanism was altered and that it could be the cause or the effect of the disease. Therefore, genes within these modules could play a role in cancer and thus be considered as potential therapeutic targets or diagnostic/prognostic biomarkers. Here, we have reviewed all the studies where WGCNA has been applied to gene expression data from bladder cancer patients. We have shown the importance of this new approach in identifying candidate biomarkers and therapeutic targets. They include both genes and miRNAs and some of them have already been identified in the literature to have a role in bladder cancer initiation, progression, metastasis, and patient survival.

Possible role of nucleotide correlations between human exon junctions.

There is ample evidence that prediction of human splice sites can be refined by analyzing the nucleotides surrounding splice sites. This could mean that exon nucleotides over splice sites harbour information for the splicing process in addition to the coding information to specify aminoacids. We analyzed the correlations among the nucleotides lying at the end and at the beginning of all the consecutive human exons to seek relationships among the nucleotides. We have divided the sequences taking into account the phase of interruption. Even though exon sequences are involved in the coding function, we found phase-dependent, specific correlations in the area of exon junctions. These regularities do not give rise to specific motifs, but rather to a phase-specific nucleotide context that could contribute to define the splice site or aid the splicing machinery to join the exon ends. Results provide further evidence that accurate selection of human splice sites likely requires the contribution of exon regulatory sequences.