Search

In the last two decades, many detailed full transcriptomic studies on complex biological samples have been published and included in large gene expression repositories. These studies primarily provide a bulk expression signal for each sample, including multiple cell-types mixed within the global signal. The cellular heterogeneity in these mixtures does not allow the activity of specific genes in specific cell types to be identified. Therefore, inferring relative cellular composition is a very powerful tool to achieve a more accurate molecular profiling of complex biological samples. In recent decades, computational techniques have been developed to solve this problem by applying deconvolution methods, designed to decompose cell mixtures into their cellular components and calculate the relative proportions of these elements. Some of them only calculate the cell proportions (supervised methods), while other deconvolution algorithms can also identify the gene signatures specific for each cell type (unsupervised methods). In these work, five deconvolution methods (CIBERSORT, FARDEEP, DECONICA, LINSEED and ABIS) were implemented and used to analyze blood and immune cells, and also cancer cells, in complex mixture samples (using three bulk expression datasets). Our study provides three analytical tools (corrplots, cell-signature plots and bar-mixture plots) that allow a thorough comparative analysis of the cell mixture data. The work indicates that CIBERSORT is a robust method optimized for the identification of immune cell-types, but not as efficient in the identification of cancer cells. We also found that LINSEED is a very powerful unsupervised method that provides precise and specific gene signatures for each of the main immune cell types tested: neutrophils and monocytes (of the myeloid lineage), B-cells, NK cells and T-cells (of the lymphoid lineage), and also for cancer cells.

Human Proteome Project (HPP) presents a systematic characterization of the protein landscape under different conditions using several complementary-omic techniques (LC-MS/MS proteomics, affinity proteomics, transcriptomics, etc.). In the present study, using a B-cell lymphoma cell line as a model, comprehensive integration of RNA-Seq transcriptomics, MS/MS, and antibody-based affinity proteomics (combined with size-exclusion chromatography) (SEC-MAP) were performed to uncover correlations that could provide insights into protein dynamics at the intracellular level. Here, 5672 unique proteins were systematically identified by MS/MS analysis and subcellular protein extraction strategies (neXtProt release 2020-21, MS/MS data are available via ProteomeXchange with identifier PXD003939). Moreover, RNA deep sequencing analysis of this lymphoma B-cell line identified 19,518 expressed genes and 5707 protein coding genes (mapped to neXtProt). Among these data sets, 162 relevant proteins (targeted by 206 antibodies) were systematically analyzed by the SEC-MAP approach, providing information about PTMs, isoforms, protein complexes, and subcellular localization. Finally, a bioinformatic pipeline has been designed and developed for orthogonal integration of these high-content proteomics and transcriptomics datasets, which might be useful for comprehensive and global characterization of intracellular protein profiles.

Background: Identification of biomarkers associated with the prognosis of different cancer subtypes is critical to achieve better therapeutic assistance. In colorectal cancer (CRC) the discovery of stable and consistent survival markers remains a challenge due to the high heterogeneity of this class of tumors. In this work, we identified a new set of gene markers for CRC associated to prognosis and risk using a large unified cohort of patients with transcriptomic profiles and survival information., Results: We built an integrated dataset with 1273 human colorectal samples, which provides a homogeneous robust framework to analyse genome-wide expression and survival data. Using this dataset we identified two sets of genes that are candidate prognostic markers for CRC in stages III and IV, showing either up-regulation correlated with poor prognosis or up-regulation correlated with good prognosis. The top 10 up-regulated genes found as survival markers of poor prognosis (i.e. low survival) were: DCBLD2, PTPN14, LAMP5, TM4SF1, NPR3, LEMD1, LCA5, CSGALNACT2, SLC2A3 and GADD45B. The stability and robustness of the gene survival markers was assessed by cross-validation, and the best-ranked genes were also validated with two external independent cohorts: one of microarrays with 482 samples; another of RNA-seq with 269 samples. Up-regulation of the top genes was also proved in a comparison with normal colorectal tissue samples. Finally, the set of top 100 genes that showed overexpression correlated with low survival was used to build a CRC risk predictor applying a multivariate Cox proportional hazards regression analysis. This risk predictor yielded an optimal separation of the individual patients of the cohort according to their survival, with a p-value of 8.25e-14 and Hazard Ratio 2.14 (95% CI: 1.75-2.61)., Conclusions: The results presented in this work provide a solid rationale for the prognostic utility of a new set of genes in CRC, demonstrating their potential to predict colorectal tumor progression and evolution towards poor survival stages. Our study does not provide a fixed gene signature for prognosis and risk prediction, but instead proposes a robust set of genes ranked according to their predictive power that can be selected for additional tests with other CRC clinical cohorts.

Background: In the study of complex diseases using genome-wide expression data from clinical samples, a difficult case is the identification and mapping of the gene signatures associated to the stages that occur in the progression of a disease. The stages usually correspond to different subtypes or classes of the disease, and the difficulty to identify them often comes from patient heterogeneity and sample variability that can hide the biomedical relevant changes that characterize each stage, making standard differential analysis inadequate or inefficient., Results: We propose a methodology to study diseases or disease stages ordered in a sequential manner (e.g. from early stages with good prognosis to more acute or serious stages associated to poor prognosis). The methodology is applied to diseases that have been studied obtaining genome-wide expression profiling of cohorts of patients at different stages. The approach allows searching for consistent expression patterns along the progression of the disease through two major steps: (i) identifying genes with increasing or decreasing trends in the progression of the disease; (ii) clustering the increasing/decreasing gene expression patterns using an unsupervised approach to reveal whether there are consistent patterns and find genes altered at specific disease stages. The first step is carried out using Gamma rank correlation to identify genes whose expression correlates with a categorical variable that represents the stages of the disease. The second step is done using a Self Organizing Map (SOM) to cluster the genes according to their progressive profiles and identify specific patterns. Both steps are done after normalization of the genomic data to allow the integration of multiple independent datasets. In order to validate the results and evaluate their consistency and biological relevance, the methodology is applied to datasets of three different diseases: myelodysplastic syndrome, colorectal cancer and Alzheimer's disease. A software script written in R, named genediseasePatterns, is provided to allow the use and application of the methodology., Conclusion: The method presented allows the analysis of the progression of complex and heterogeneous diseases that can be divided in pathological stages. It identifies gene groups whose expression patterns change along the advance of the disease, and it can be applied to different types of genomic data studying cohorts of patients in different states.

Lipid metabolism plays an essential role in carcinogenesis due to the requirements of tumoral cells to sustain increased structural, energetic and biosynthetic precursor demands for cell proliferation. We investigated the association between expression of lipid metabolism-related genes and clinical outcome in intermediate-stage colon cancer patients with the aim of identifying a metabolic profile associated with greater malignancy and increased risk of relapse. Expression profile of 70 lipid metabolism-related genes was determined in 77 patients with stage II colon cancer. Cox regression analyses using c-index methodology was applied to identify a metabolic-related signature associated to prognosis. The metabolic signature was further confirmed in two independent validation sets of 120 patients and additionally, in a group of 264 patients from a public database. The combined analysis of these 4 genes, ABCA1, ACSL1, AGPAT1 and SCD, constitutes a metabolic-signature (ColoLipidGene) able to accurately stratify stage II colon cancer patients with 5-fold higher risk of relapse with strong statistical power in the four independent groups of patients. The identification of a group of 4 genes that predict survival in intermediate-stage colon cancer patients allows delineation of a high-risk group that may benefit from adjuvant therapy, and avoids the toxic and unnecessary chemotherapy in patients classified as low-risk group.

Background: Despite the large increase of transcriptomic studies that look for gene signatures on diseases, there is still a need for integrative approaches that obtain separation of multiple pathological states providing robust selection of gene markers for each disease subtype and information about the possible links or relations between those genes., Results: We present a network-oriented and data-driven bioinformatic approach that searches for association of genes and diseases based on the analysis of genome-wide expression data derived from microarrays or RNA-Seq studies. The approach aims to (i) identify gene sets associated to different pathological states analysed together; (ii) identify a minimum subset within these genes that unequivocally differentiates and classifies the compared disease subtypes; (iii) provide a measurement of the discriminant power of these genes and (iv) identify links between the genes that characterise each of the disease subtypes. This bioinformatic approach is implemented in an R package, named geNetClassifier, available as an open access tool in Bioconductor. To illustrate the performance of the tool, we applied it to two independent datasets: 250 samples from patients with four major leukemia subtypes analysed using expression arrays; another leukemia dataset analysed with RNA-Seq that includes a subtype also present in the previous set. The results show the selection of key deregulated genes recently reported in the literature and assigned to the leukemia subtypes studied. We also show, using these independent datasets, the selection of similar genes in a network built for the same disease subtype., Conclusions: The construction of gene networks related to specific disease subtypes that include parameters such as gene-to-gene association, gene disease specificity and gene discriminant power can be very useful to draw gene-disease maps and to unravel the molecular features that characterize specific pathological states. The application of the bioinformatic tool here presented shows a neat way to achieve such molecular characterization of the diseases using genome-wide expression data.

Bone metastasis of lung adenocarcinoma (AC) is a frequent complication of advanced disease. The purpose of this study was to identify key mediators conferring robust prometastatic activity with clinical significance. We isolated highly metastatic subpopulations (HMS) using a previously described in vivo model of lung AC bone metastasis. We performed transcriptomic profiling of HMS and stringent bioinformatics filtering. Functional validation was assessed by overexpression and lentiviral silencing of single, double and triple combination in vivo and in vitro. We identified HDAC4, PITX1 and ROBO1 that decreased bone metastatic ability after their simultaneous abrogation. These effects were solely linked to defects in osseous colonization. The molecular mechanisms related to bone colonization were mediated by non-cell autonomous effects that include the following: (1) a marked decrease in osteoclastogenic activity in vitro and in vivo, an effect associated with reduced pro-osteoclastogenic cytokines IL-11 and PTHrP expression levels, as well as decreased in vitro expression of stromal rankl in conditions mimicking tumor-stromal interactions; (2) an abrogated response to TGF-β signaling by decreased phosphorylation and levels of Smad2/3 in tumor cells and (3) an impaired metalloproteolytic activity in vitro. Interestingly, coexpression of HDAC4 and PITX1 conferred high prometastatic activity in vivo. Further, levels of both genes correlated with patients at higher risk of metastasis in a clinical lung AC data set and with a poorer clinical outcome. These findings provide functional and clinical evidence that this metastatic subset is an important determinant of osseous colonization. These data suggest novel therapeutic targets to effectively block lung AC bone metastasis.

Background: Accurate analysis of whole-gene expression and individual-exon expression is essential to characterize different transcript isoforms and identify alternative splicing events in human genes. One of the omic technologies widely used in many studies on human samples are the exon-specific expression microarray platforms., Results: Since there are not many validated comparative analyses to identify specific splicing events using data derived from these types of platforms, we have developed an algorithm (called ESLiM) to detect significant changes in exon use, and applied it to a reference dataset of 270 human genes that show alternative expression in different tissues. We compared the results with three other methodological approaches and provided the R source code to be applied elsewhere. The genes positively detected by these analyses also provide a verified subset of human genes that present tissue-regulated isoforms. Furthermore, we performed a validation analysis on human patient samples comparing two different subtypes of acute myeloid leukemia (AML) and we experimentally validated the splicing in several selected genes that showed exons with highly significant signal change., Conclusions: The comparative analyses with other methods using a fair set of human genes that show alternative splicing and the validation on clinical samples demonstrate that the proposed novel algorithm is a reliable tool for detecting differential splicing in exon-level expression data.

Despite evidence about the implication of the bone marrow (BM) stromal microenvironment in multiple myeloma (MM) cell growth and survival, little is known about the effects of myelomatous cells on BM stromal cells. Mesenchymal stromal cells (MSCs) from healthy donors (dMSCs) or myeloma patients (pMSCs) were co-cultured with the myeloma cell line MM.1S, and the transcriptomic profile of MSCs induced by this interaction was analyzed. Deregulated genes after co-culture common to both d/pMSCs revealed functional involvement in tumor microenvironment cross-talk, myeloma growth induction and drug resistance, angiogenesis and signals for osteoclast activation and osteoblast inhibition. Additional genes induced by co-culture were exclusively deregulated in pMSCs and predominantly associated to RNA processing, the ubiquitine-proteasome pathway, cell cycle regulation, cellular stress and non-canonical Wnt signaling. The upregulated expression of five genes after co-culture (CXCL1, CXCL5 and CXCL6 in d/pMSCs, and Neuregulin 3 and Norrie disease protein exclusively in pMSCs) was confirmed, and functional in vitro assays revealed putative roles in MM pathophysiology. The transcriptomic profile of pMSCs co-cultured with myeloma cells may better reflect that of MSCs in the BM of myeloma patients, and provides new molecular insights to the contribution of these cells to MM pathophysiology and to myeloma bone disease.

Gene expression profiling signatures may be used to classify the subtypes of Myelodysplastic syndrome (MDS) patients. However, there are few reports on the global methylation status in MDS. The integration of genome-wide epigenetic regulatory marks with gene expression levels would provide additional information regarding the biological differences between MDS and healthy controls. Gene expression and methylation status were measured using high-density microarrays. A total of 552 differentially methylated CpG loci were identified as being present in low-risk MDS; hypermethylated genes were more frequent than hypomethylated genes. In addition, mRNA expression profiling identified 1005 genes that significantly differed between low-risk MDS and the control group. Integrative analysis of the epigenetic and expression profiles revealed that 66.7% of the hypermethylated genes were underexpressed in low-risk MDS cases. Gene network analysis revealed molecular mechanisms associated with the low-risk MDS group, including altered apoptosis pathways. The two key apoptotic genes BCL2 and ETS1 were identified as silenced genes. In addition, the immune response and micro RNA biogenesis were affected by the hypermethylation and underexpression of IL27RA and DICER1. Our integrative analysis revealed that aberrant epigenetic regulation is a hallmark of low-risk MDS patients and could have a central role in these diseases.

Most genetic studies in colorectal carcinomas have focused on those abnormalities that are acquired by primary tumors, particularly in the transition from adenoma to carcinoma, whereas few studies have compared the genetic abnormalities of primary versus paired metastatic samples. In this study, we used high-density 500K single-nucleotide polymorphism arrays to map the overall genetic changes present in liver metastases (n=20) from untreated colorectal carcinoma patients studied at diagnosis versus their paired primary tumors (n=20). MLH1, MSH2 and MSH6 gene expression was measured in parallel by immunohistochemistry. Overall, metastatic tumors systematically contained those genetic abnormalities observed in the primary tumor sample from the same subject. However, liver metastases from many cases (up to 8 out of 20) showed acquisition of genetic aberrations that were not found in their paired primary tumors. These new metastatic aberrations mainly consisted of (1) an increased frequency of genetic lesions of chromosomes that have been associated with metastatic colorectal carcinoma (1p, 7p, 8q, 13q, 17p, 18q, 20q) and, more interestingly, (2) acquisition of new chromosomal abnormalities (eg, losses of chromosomes 4 and 10q and gains of chromosomes 5p and 6p). These genetic changes acquired by metastatic tumors may be associated with either the metastatic process and/or adaption of metastatic cells to the liver microenvironment. Further studies in larger series of patients are necessary to dissect the specific role of each of the altered genes and chromosomal regions in the metastatic spread of colorectal tumors.

Background: Genome-wide expression studies have developed exponentially in recent years as a result of extensive use of microarray technology. However, expression signals are typically calculated using the assignment of "probesets" to genes, without addressing the problem of "gene" definition or proper consideration of the location of the measuring probes in the context of the currently known genomes and transcriptomes. Moreover, as our knowledge of metazoan genomes improves, the number of both protein-coding and noncoding genes, as well as their associated isoforms, continues to increase. Consequently, there is a need for new databases that combine genomic and transcriptomic information and provide updated mapping of expression probes to current genomic annotations., Results: GATExplorer (Genomic and Transcriptomic Explorer) is a database and web platform that integrates a gene loci browser with nucleotide level mappings of oligo probes from expression microarrays. It allows interactive exploration of gene loci, transcripts and exons of human, mouse and rat genomes, and shows the specific location of all mappable Affymetrix microarray probes and their respective expression levels in a broad set of biological samples. The web site allows visualization of probes in their genomic context together with any associated protein-coding or noncoding transcripts. In the case of all-exon arrays, this provides a means by which the expression of the individual exons within a gene can be compared, thereby facilitating the identification and analysis of alternatively spliced exons. The application integrates data from four major source databases: Ensembl, RNAdb, Affymetrix and GeneAtlas; and it provides the users with a series of files and packages (R CDFs) to analyze particular query expression datasets. The maps cover both the widely used Affymetrix GeneChip microarrays based on 3' expression (e.g. human HG U133 series) and the all-exon expression microarrays (Gene 1.0 and Exon 1.0)., Conclusions: GATExplorer is an integrated database that combines genomic/transcriptomic visualization with nucleotide-level probe mapping. By considering expression at the nucleotide level rather than the gene level, it shows that the arrays detect expression signals from entities that most researchers do not contemplate or discriminate. This approach provides the means to undertake a higher resolution analysis of microarray data and potentially extract considerably more detailed and biologically accurate information from existing and future microarray experiments.

Specific microRNA (miRNA) signatures have been associated with different cytogenetic subtypes in acute leukemias. This finding prompted us to investigate potential associations between genetic abnormalities in multiple myeloma (MM) and singular miRNA expression profiles. Moreover, global gene expression profiling was also analyzed to find correlated miRNA gene expression and select miRNA target genes that show such correlation. For this purpose, we analyzed the expression level of 365 miRNAs and the gene expression profiling in 60 newly diagnosed MM patients, selected to represent the most relevant recurrent genetic abnormalities. Supervised analysis showed significantly deregulated miRNAs in the different cytogenetic subtypes as compared with normal PC. It is interesting to note that miR-1 and miR-133a clustered on the same chromosomal loci, were specifically overexpressed in the cases with t(14;16). The analysis of the relationship between miRNA expression and their respective target genes showed a conserved inverse correlation between several miRNAs deregulated in MM cells and CCND2 expression level. These results illustrate, for the first time, that miRNA expression pattern in MM is associated with genetic abnormalities, and that the correlation of the expression profile of miRNA and their putative mRNA targets is useful to find statistically significant protein-coding genes in MM pathogenesis associated with changes in specific miRNAs.

Control of transgene expression from long-term expression vectors can be achieved with inducible and regulated promoters. The two most commonly used inducible systems employ doxycycline or mifepristone as the drug activating a silent trans-activator, which is expressed from a constitutive promoter. We evaluated the alterations provoked by constitutive expression in the liver of rtTA2(S)-M2 (rtTA2; second-generation reverse tetracycline-controlled trans-activator) and GLp65, which are the trans-activators of the doxycyline- and mifepristone-inducible systems, respectively. To this end we performed transcriptomic analysis of mice expressing these trans-activators in the liver over 1 month. rtTA2 expression induced alterations in a few genes (69 gene probesets; false discovery rate [FDR], approximately 0.05), whereas GLp65 caused more numerous changes (1059 gene probe-sets, an FDR of approximately 0.05). However, only 20 and 53 of the genes from the rtTA2 and GLp65 groups, respectively, showed changes (R-fold >or= 3). Functional assignments indicate that alterations were mild and of little general significance. Few additional transcriptomic changes were observed when expressing trans-activators in the presence of inducer drugs; most were due to the drugs themselves. These results and the absence of toxicity observed in treated animals indicate that the two inducible systems are well tolerated and have little impact on the liver transcriptome profile. The milder alterations found with the use of rtTA2 suggest that this system is possibly safer for gene therapy applications.

Background: Analysis of gene expression data using genome-wide microarrays is a technique often used in genomic studies to find coexpression patterns and locate groups of co-transcribed genes. However, most studies done at global "omic" scale are not focused on human samples and when they correspond to human very often include heterogeneous datasets, mixing normal with disease-altered samples. Moreover, the technical noise present in genome-wide expression microarrays is another well reported problem that many times is not addressed with robust statistical methods, and the estimation of errors in the data is not provided., Methodology/principal Findings: Human genome-wide expression data from a controlled set of normal-healthy tissues is used to build a confident human gene coexpression network avoiding both pathological and technical noise. To achieve this we describe a new method that combines several statistical and computational strategies: robust normalization and expression signal calculation; correlation coefficients obtained by parametric and non-parametric methods; random cross-validations; and estimation of the statistical accuracy and coverage of the data. All these methods provide a series of coexpression datasets where the level of error is measured and can be tuned. To define the errors, the rates of true positives are calculated by assignment to biological pathways. The results provide a confident human gene coexpression network that includes 3327 gene-nodes and 15841 coexpression-links and a comparative analysis shows good improvement over previously published datasets. Further functional analysis of a subset core network, validated by two independent methods, shows coherent biological modules that share common transcription factors. The network reveals a map of coexpression clusters organized in well defined functional constellations. Two major regions in this network correspond to genes involved in nuclear and mitochondrial metabolism and investigations on their functional assignment indicate that more than 60% are house-keeping and essential genes. The network displays new non-described gene associations and it allows the placement in a functional context of some unknown non-assigned genes based on their interactions with known gene families., Conclusions/significance: The identification of stable and reliable human gene to gene coexpression networks is essential to unravel the interactions and functional correlations between human genes at an omic scale. This work contributes to this aim, and we are making available for the scientific community the validated human gene coexpression networks obtained, to allow further analyses on the network or on some specific gene associations. The data are available free online at http://bioinfow.dep.usal.es/coexpression/.

Motivation: Alteration of gene expression often results in up- or down-regulated genes and the most common analysis strategies look for such differentially expressed genes. However, molecular disease mechanisms typically constitute abnormalities in the regulation of genes producing strong alterations in the expression levels. The search for such deregulation states in the genomic expression profiles will help to identify disease-altered genes better., Results: We have developed an algorithm that searches for the genes which present a significant alteration in the variability of their expression profiles, by comparing an altered state with a control state. The algorithm provides groups of genes and assigns a statistical measure of significance to each group of genes selected. The method also includes a prefilter tool to select genes with a threshold of differential expression that can be set by the user ad casum. The method is evaluated using an experimental set of microarrays of human control and cancer samples from patients with acute promyelocytic leukemia.