Rank Aggregation Method Research Articles

Hybrid multi-criteria decision-making approach to select appropriate biomass resources for biofuel production

Biofuel generation from local biomass resources can significantly contribute to greenhouse gas mitigation and cleaner energy production. In this regard, a hybrid Multi-Criteria Decision-Making (MCDM) approach was employed to prioritize appropriate biomass resources for biofuel production. Technique for Order Preference by Similarity to Ideal Solution (TOPSIS), Additive Ratio Assessment (ARAS), and Weighted Aggregates Sum Product Assessment (WASPAS) were the employed MCDM approaches. Subsequently, ranking aggregation methods, i.e., Borda, Copeland, and Rank Mean, were applied to integrate the rankings obtained from the MCDM approaches. Guilan province of Iran was selected as a case study based on its promising potential for biofuel production from first-, second-, and third-generation biofuel resources. Initially, through an in-depth review of the literature and the use of academic professors' expert opinions, ten criteria were selected as the evaluation indices of the study: 1) creating technical side jobs, 2) preserving non-renewable energy resources, 3) relative advantage of biofuel production, 4) complexity of biofuel production process, 5) cost of the biomass conversion process, 6) biomass reusability, 7) cost of biomass supply, 8) environmental impacts of biomass accumulation, 9) adaptability of the biofuel production process to the size of biomass production units and the attitude and knowledge of the producers, and 10) energy self-sufficiency of the biomass producer. Moreover, the 11 investigated potential sources of biofuel production were rice, peanut, livestock and poultry wastes, rice waste, peanut waste, tea residues and its processing wastes, olive residues and its processing wastes, livestock and poultry slaughter and farm-raised fish wastes, municipal solid waste and sewage, forest and wood farming wastes, algae and Azolla. The results indicated that “municipal solid wastes and sewage”, “forest and wood farming wastes” and “livestock and poultry wastes” from the second-generation biofuels were identified as the most important biomass resources in the studied area.

Robust rank aggregation and cibersort algorithm applied to the identification of key genes in head and neck squamous cell cancer

<abstract> <sec><title>Objective</title> Although multiple hub genes have been identified in head and neck squamous cell cancer (HNSCC) in recent years, because of the limited sample size and inconsistent bioinformatics analysis methods, the results are not reliable. Therefore, it is urgent to use reliable algorithms to find new prognostic markers of HNSCC. </sec> <sec><title>Method</title> The Robust Rank Aggregation (RRA) method was used to integrate 8 microarray datasets of HNSCC downloaded from the Gene Expression Omnibus (GEO) database to screen differentially expressed genes (DEGs). Later, Gene Ontology (GO) functional annotation together with Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was carried out to discover functions of those discovered DEGs. According to the KEGG results, those discovered DEGs showed tight association with the occurrence and development of HNSCC. Then cibersort algorithm was used to analyze the infiltration of immune cells of HNSCC and we found that the main infiltrated immune cells were B cells, dendritic cells and macrophages. A protein-protein interaction (PPI) network was established; moreover, key modules were also constructed to select 5 hub genes from the whole network using cytoHubba. 3 hub genes showed significant relationship with prognosis for TCGA-derived HNSCC patients. </sec> <sec><title>Result</title> The potent DEGs along with hub genes were selected by the combined bioinformatic approach. <italic>AURKA</italic>, <italic>BIRC5</italic> and <italic>UBE2C</italic> genes may be the potential prognostic biomarker and therapeutic targets of HNSCC. </sec> <sec><title>Conclusions</title> The Robust Rank Aggregation method and cibersort algorithm method can accurately predict the potential prognostic biomarker and therapeutic targets of HNSCC through multiple GEO datasets. </sec> </abstract>

Identification of TYR, TYRP1, DCT and LARP7 as related biomarkers and immune infiltration characteristics of vitiligo via comprehensive strategies

ABSTRACT This study aims to explore biomarkers associated with vitiligo and analyze the pathological role of immune cell infiltration in the disease. We used the robust rank aggregation (RRA) method to integrate three vitiligo data sets downloaded from gene expression omnibus database, identify the differentially expressed genes (DEGs) and analyze the functional correlation. Then, the comprehensive strategy of combined weighted gene coexpression network analysis (WGCNA) and logical regression of the selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) machine learning algorithm are employed to screen and biomarkers associated with vitiligo. Finally, the immune cell infiltration of vitiligo was evaluated by CIBERSORT, and the correlation between biomarkers and infiltrating immune cells was analyzed. Herein, we identified 131 robust DEGs, and enrichment analysis results showed that robust DEGs and melanogenesis were closely associated with vitiligo development and progression. TYR, TYRP1, DCT and LARP7 were identified as vitiligo-related biomarkers. Immune infiltration analysis demonstrated that CD4 T Cell, CD8 T Cell, Tregs, NK cells, dendritic cells, and macrophages were involved in vitiligo’s pathogenesis. In summary, we adopted a comprehensive strategy to screen biomarkers related to vitiligo and explore the critical role of immune cell infiltration in vitiligo. Abbreviations: TYR, Tyrosinase; TYRP1, Tyrosinase-related protein-1; DCT, dopachrome tautomerase; LARP7, La ribonucleoprotein domain family, member-7; RRA, robust rank aggregation; DEGs, differentially expressed genes; WGCNA, weighted gene coexpression network analysis; LASSO, logical regression of the selection operator; SVM-RFE, support vector machine recursive feature elimination; RF, random forest; GWAS, Genome-wide association study; FasL, Fas-Fas ligand; Tregs, T-regulatory cells; NK, natural killer; GEPCs, gene expression profiling chips; GO, gene ontology; GSEA, gene set enrichment analysis; FDR, false discovery rate; AUC, area under the curve; ROC, receiver-operating characteristic; BP, biological process; CC, cellular component; MF, molecular function.

A Novel Rank Aggregation-Based Hybrid Multifilter Wrapper Feature Selection Method in Software Defect Prediction.

The high dimensionality of software metric features has long been noted as a data quality problem that affects the performance of software defect prediction (SDP) models. This drawback makes it necessary to apply feature selection (FS) algorithm(s) in SDP processes. FS approaches can be categorized into three types, namely, filter FS (FFS), wrapper FS (WFS), and hybrid FS (HFS). HFS has been established as superior because it combines the strength of both FFS and WFS methods. However, selecting the most appropriate FFS (filter rank selection problem) for HFS is a challenge because the performance of FFS methods depends on the choice of datasets and classifiers. In addition, the local optima stagnation and high computational costs of WFS due to large search spaces are inherited by the HFS method. Therefore, as a solution, this study proposes a novel rank aggregation-based hybrid multifilter wrapper feature selection (RAHMFWFS) method for the selection of relevant and irredundant features from software defect datasets. The proposed RAHMFWFS is divided into two stepwise stages. The first stage involves a rank aggregation-based multifilter feature selection (RMFFS) method that addresses the filter rank selection problem by aggregating individual rank lists from multiple filter methods, using a novel rank aggregation method to generate a single, robust, and non-disjoint rank list. In the second stage, the aggregated ranked features are further preprocessed by an enhanced wrapper feature selection (EWFS) method based on a dynamic reranking strategy that is used to guide the feature subset selection process of the HFS method. This, in turn, reduces the number of evaluation cycles while amplifying or maintaining its prediction performance. The feasibility of the proposed RAHMFWFS was demonstrated on benchmarked software defect datasets with Naïve Bayes and Decision Tree classifiers, based on accuracy, the area under the curve (AUC), and F-measure values. The experimental results showed the effectiveness of RAHMFWFS in addressing filter rank selection and local optima stagnation problems in HFS, as well as the ability to select optimal features from SDP datasets while maintaining or enhancing the performance of SDP models. To conclude, the proposed RAHMFWFS achieved good performance by improving the prediction performances of SDP models across the selected datasets, compared to existing state-of-the-arts HFS methods.

Analysis of Hub Genes Involved in Distinction Between Aged and Fetal Bone Marrow Mesenchymal Stem Cells by Robust Rank Aggregation and Multiple Functional Annotation Methods.

Stem cells from fetal tissue protect against aging and possess greater proliferative capacity than their adult counterparts. These cells can more readily expand in vitro and senesce later in culture. However, the underlying molecular mechanisms for these differences are still not fully understood. In this study, we used a robust rank aggregation (RRA) method to discover robust differentially expressed genes (DEGs) between fetal bone marrow mesenchymal stem cells (fMSCs) and aged adult bone marrow mesenchymal stem cells (aMSCs). Multiple methods, including gene set enrichment analysis (GSEA), Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for functional annotation of the robust DEGs, and the results were visualized using the R software. The hub genes and other genes with which they interacted directly were detected by protein–protein interaction (PPI) network analysis. Correlation of gene expression was measured by Pearson correlation coefficient. A total of 388 up-regulated and 289 down-regulated DEGs were identified between aMSCs and fMSCs. We found that the down-regulated genes were mainly involved in the cell cycle, telomerase activity, and stem cell proliferation. The up-regulated DEGs were associated with cell adhesion molecules, extracellular matrix (ECM)–receptor interactions, and the immune response. We screened out four hub genes, MYC, KIF20A, HLA-DRA, and HLA-DPA1, through PPI-network analysis. The MYC gene was negatively correlated with TXNIP, an age-related gene, and KIF20A was extensively involved in the cell cycle. The results suggested that MSCs derived from the bone marrow of an elderly donor present a pro-inflammatory phenotype compared with that of fMSCs, and the HLA-DRA and HLA-DPA1 genes are related to the immune response. These findings provide new insights into the differences between aMSCs and fMSCs and may suggest novel strategies for ex vivo expansion and application of adult MSCs.

Co-Occurring Mutation Clusters Predict Drug Sensitivity in Acute Myeloid Leukemia

Background The molecular origin of cancer drug resistance is not apparent for most cases of acute myeloid leukemia (AML). Clonal evolution appears to be associated with increasing drug resistance. We sought to determine whether the mutation patterns are associated with drug susceptibility in AML. Methods Seventy-two patient blood or marrow samples were enriched for CD34+ blasts by immunomagnetic bead selection. High throughput drug sensitivity screens were performed for 223 drugs after a 72-hour exposure to 8-12 customized drug concentrations (within the range of 5pM to 100µM) of each drug spanning 4-5 logs. Post exposure viability was determined using CellTiter Glo luminescent reagent. XLFit was used to analyze the data and generate dose response curves based on a standard 4-parameter logistic fit. Mutation analysis was performed by MyAML™ utilizing next generation sequencing (NGS) to analyze the 3' and 5' UTRs and exonic regions of 194 AML-associated genes and genomic breakpoints. Co-mutation and mutual exclusivity scores of gene aberrations were computed using three different statistical methods on three large publicly available datasets: 1) TCGA-AML, 2) BeatAML, and 3) data from 1540 patients (Moritz Gerstung, et al., Nat Genet. 2017). The scores for each co-mutation or mutually exclusive gene pair were then logarithmically transformed from the aggregated p-value using a robust rank aggregation method and used to construct a graph from which communities of co-mutated genes could be detected, resulting in higher mutual exclusivity between modules and higher co-occurrence within modules. Results We identified five main groups of co-mutations, including the following: 1) RUNX1 group, 2) CEBPA group, 3) NPM1 group, 4) TP53 group, and 5) RAS group (see Fig 1A). The co-mutation community in the RUNX1 group is featured with transcriptional dysregulation (RUNX1, ASXL1 and EZH2), and dysregulation in splicing (U2AF1, SRSF2 or SF3B1). The TP53 and CEBPA groups exhibit transcriptional dysregulation, and transcription factor alterations. The NPM1 and RAS groups exhibit signaling alteration, particularly the Ras/MAPK signaling pathway. Among each co-mutation community, the variant allele frequencies (VAFs) exhibited differences for specific mutations. For example, in the NPM1 group, FLT3 exhibits the lowest VAF, followed by NPM1 while DNMT3A shows the highest VAF. The five co-mutation groups are associated with different overall survival, including better survival for the CEBPA group and poorer survival for the TP53 group. Within each co-mutation group, specific associations with overall survival were detected, including better survival for NPM1-RAD21 co-mutated vs. wild type (WT), better survival for WT vs. ASXL1-RUNX1 co-mutated, worse survival for DNMT3A-FLT3 co-mutated vs. WT, and others. We also found that the co-mutation clusters were different between the de novo and relapse groups. Most importantly, we found significant correlations with drug sensitivity for the different co-occurring mutation groups. Cells from patients with mutations in the NPM1 co-mutation group show higher sensitivity to most of the drugs, including the tyrosine kinases inhibitors and PI3K-AKT-MTOR inhibitors. Cells from patients with mutations in the RAS cluster exhibited sensitivity to MEK inhibitors (Fig 1B), while cells from patients with mutations in the TP53 cluster show resistance to many drugs, including the MDM2 inhibitor AMG232 (Fig 1C). By mapping the mutated genes and drug targets into signaling pathways, we found that mutations in downstream signaling of the drug target exhibited resistance, while mutations in upstream signaling conferred higher sensitivity. We also constructed drug sensitivity prediction models based on the co-mutation groups. Conclusion We detected co-mutation groups through integrative analysis of large publicly available AML data sets. Patients with mutations in different groups of genes show different overall survival, while cells from patients with mutations in different co-mutation groups show different drug sensitivity. The co-mutation groups may not only reflect the clonal evolution history of AML, but also serve as important features for drug sensitivity prediction. Disclosures Carson: Invivoscribe, Inc: Current Employment. Patay:Invivoscribe, Inc: Current Employment. Becker:Cardiff Oncology: Research Funding; Novartis: Research Funding; SecuraBio: Research Funding; JW Pharmaceutical: Research Funding; Glycomimetics: Research Funding; Abbvie: Research Funding; Bristol Myers Squibb: Research Funding; Pfizer: Research Funding; Invivoscribe: Research Funding; Accordant Health Services/Caremark: Membership on an entity's Board of Directors or advisory committees. OffLabel Disclosure: There are ~200 drugs on the screening panel, and many would be off label for acute myeloid leukemia.

Transcriptomic profiles in Parkinson's disease.

Transcriptomics in Parkinson's disease offers insights into the pathogenesis of Parkinson's disease but obtaining brain tissue has limitations. In order to bypass this issue, we profile and compare differentially expressed genes and enriched pathways (KEGG) in two peripheral tissues (blood and skin) of 12 Parkinson's disease patients and 12 healthy controls using RNA-sequencing technique and validation with RT-qPCR. Furthermore, we compare our results to previous Parkinson's disease post mortem brain tissue and blood results using the robust rank aggregation method. The results show no overlapping differentially expressed genes or enriched pathways in blood vs. skin in our sample sets (25 vs. 1068 differentially expressed genes with an FDR ≤ 0.05; 1 vs. 9 pathways in blood and skin, respectively). A meta-analysis from previous transcriptomic sample sets using either microarrays or RNA-Seq yields a robust rank aggregation list of cortical gene expression changes with 43 differentially expressed genes; a list of substantia nigra changes with 2 differentially expressed genes and a list of blood changes with 1 differentially expressed gene being statistically significant at FDR ≤ 0.05. In cortex 1, KEGG pathway was enriched, four in substantia nigra and two in blood. None of the differentially expressed genes or pathways overlap between these tissues. When comparing our previously published skin transcription analysis, two differentially expressed genes between the cortex robust rank aggregation and skin overlap. In this study, for the first time a meta-analysis is applied on transcriptomic sample sets in Parkinson's disease. Simultaneously, it explores the notion that Parkinson's disease is not just a neuronal tissue disease by exploring peripheral tissues. The comparison of different Parkinson's disease tissues yields surprisingly few significant differentially expressed genes and pathways, suggesting that divergent gene expression profiles in distinct cell lineages, metabolic and possibly iatrogenic effects create too much transcriptomic noise for detecting significant signal. On the other hand, there are signs that point towards Parkinson's disease-specific changes in non-neuronal peripheral tissues in Parkinson's disease, indicating that Parkinson's disease might be a multisystem disorder.

Identification of Hub Genes Associated With Hepatocellular Carcinoma Using Robust Rank Aggregation Combined With Weighted Gene Co-expression Network Analysis.

BackgroundBioinformatics provides a valuable tool to explore the molecular mechanisms underlying pathogenesis of hepatocellular carcinoma (HCC). To improve prognosis of patients, identification of robust biomarkers associated with the pathogenic pathways of HCC remains an urgent research priority.MethodsWe employed the Robust Rank Aggregation method to integrate nine qualified HCC datasets from the Gene Expression Omnibus. A robust set of differentially expressed genes (DEGs) between tumor and normal tissue samples were screened. Weighted gene co-expression network analysis was applied to cluster DEGs and the key modules related to clinical traits identified. Based on network topology analysis, novel risk genes derived from key modules were mined and biological verification performed. The potential functions of these risk genes were further explored with the aid of miRNA–mRNA regulatory networks. Finally, the prognostic ability of these genes was assessed by constructing a clinical prediction model.ResultsTwo key modules showed significant association with clinical traits. In combination with protein–protein interaction analysis, 29 hub genes were identified. Among these genes, 19 from one module showed a pattern of upregulation in HCC and were associated with the tumor node metastasis stage, and 10 from the other module displayed the opposite trend. Survival analyses indicated that all these genes were significantly related to patient prognosis. Based on the miRNA-mRNA regulatory network, 29 genes strongly linked to tumor activity were identified. Notably, five of the novel risk genes, ABAT, DAO, PCK2, SLC27A2, and HAO1, have rarely been reported in previous studies. Gene set enrichment analysis for each gene revealed regulatory roles in proliferation and prognosis of HCC. Least absolute shrinkage and selection operator regression analysis further validated DAO, PCK2, and HAO1 as prognostic factors in an external HCC dataset.ConclusionAnalysis of multiple datasets combined with global network information presents a successful approach to uncover the complex biological mechanisms of HCC. More importantly, this novel integrated strategy facilitates identification of risk hub genes as candidate biomarkers for HCC, which could effectively guide clinical treatments.

Identification and Validation of Novel Biomarkers for Diagnosis and Prognosis of Hepatocellular Carcinoma

Introduction: Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related deaths worldwide due to poor survival outcome. Thus, there is an urgent need to identify effective biomarkers for early diagnosis and prognosis prediction.Methods: A total of 389 differentially expressed genes (DEGs) between HCC samples and normal were selected based on the Robust Rank Aggregation (RRA) method. We combined DEGs expression and clinical traits to construct a gene co-expression network through WGCNA. Forty hub genes were selected from the key module. Among them, YWHAB, PPAT, NOL10 were eventually identified as prognostic biomarkers using multivariate Cox regression model. Biomarkers expression pattern was investigated by informatic analysis and verified by RNA-seq of 32 patients with HCC. DiseaseMeth 2.0, MEXPRESS, and Tumor Immune Estimation Resource (TIMER) were used to assess the methylation and immune status of biomarkers. GSVA, CCK8, colony formation assay, Edu imaging kit, wound-healing assay, and xenograft tumor model were utilized to investigate the effects of biomarkers on proliferation, metastasis of HCC cells in vitro, and in vivo. The Kaplan–Meier (KM) plotter and ROC curves were used to validate the prognostic and diagnostic value of biomarker expression.Results: All the selected biomarkers were upregulated in HCC samples and higher expression levels were associated with advanced tumor stages and T grades. The regulation of YWHAB, PPAT, NOL10 promoter methylation varied in tumors, and precancerous normal tissues. Immune infiltration analysis suggested that the abnormal regulations of these biomarkers were likely attributed to B cells and dendritic cells. GSVA for these biomarkers showed their great contributions to proliferation of HCC. Specific inhibition of their expression had strong effects on tumorigenesis in vitro and in vivo. ROC and KM curves confirmed their usefulness of diagnosis and prognosis of HCC.Conclusions: These findings identified YWHAB, PPAT, and NOL10 as novel biomarkers and validated their diagnostic and prognostic value for HCC.

Comprehensive world university ranking based on ranking aggregation

Many university rankings have been proposed in recent decades. The remarkable divergence among various rankings leads to confusion for decision-makers. In this paper, we propose to generate a comprehensive world university ranking by aggregating existing individual university rankings. We present a new graph-based rank aggregation method by defining a competition graph of universities, in which each node represents a university and each directed edge represents an outranking relation between two universities. We propose to measure the quality of a university by the out-in degree ratio based on which we rank all universities. Moreover, We evaluate the effectiveness of our comprehensive world university ranking from the perspectives of normality and impartiality, respectively. It is shown that the aggregated ranking will be applied as a blend integrating all the information from individual university rankings and can efficiently eliminate the outliers and regional partiality as a “smoother”.

Identification of a Novel Signature and Construction of a Nomogram Predicting Overall Survival in Clear Cell Renal Cell Carcinoma.

BackgroundClear cell renal cell carcinoma (ccRCC) is the most common subtype of renal cell carcinoma (RCC), which accounts for majority of RCC-related deaths. It is clearly essential to further identify more novel prognostic signatures and therapeutic targets.Material and MethodsWe identified differentially expressed genes (DEGs) between ccRCC and adjacent normal tissues in GEO database using a Robust Rank Aggregation (RRA) method. An mRNA signature (mRNASig) based on DEGs was developed using Cox and LASSO analysis in the TCGA database and validated in the ICGC database. Afterward, the influence of mRNASig mRNAs on the immune microenvironment in ccRCC was explored using comprehensive bioinformatics analysis.ResultsA total of 957 robust DEGs were identified using the RRA method. mRNASig comprised CEP55, IFI44, NCF4, and TCIRG1 and was developed and validated to identify high-risk patients who had poorer prognosis than low-risk patients. A nomogram was also constructed based on mRNASig, AJCC stage, and tumor grade. The mRNASig were closely related to a variety of tumor-infiltrating lymphocytes, especially including CD8+ T cells, activated CD4+ memory T cells, regulatory T cells, activated NK cells, and resting NK cells. The mRNASig were also correlated positively with the expression of CTLA4, LAG3, PDCD1, TIGIT, and HAVCR2.ConclusionWe developed and validated mRNASig to assist clinicians in making personalized treatment decisions. Furthermore, CEP55, IFI44, NCF4, and TCIRG1 may be novel potential targets for future treatment of ccRCC.

Identification of Potential Therapeutic Targets and Immune Cell Infiltration Characteristics in Osteosarcoma Using Bioinformatics Strategy.

Osteosarcoma is one of the most aggressive malignant bone tumors worldwide. Although great advancements have been made in its treatment owing to the advent of neoadjuvant chemotherapy, the problem of lung metastasis is a major obstacle in the improvement of survival outcomes. Thus, the aim of the present study is to screen novel and key biomarkers, which may act as potential prognostic markers and therapeutic targets in osteosarcoma. We utilized the robust rank aggregation (RRA) method to integrate three osteosarcoma microarray datasets downloaded from the Gene Expression Omnibus (GEO) database, and we identified the robust differentially expressed genes (DEGs) between primary and metastatic osteosarcoma tissues. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to explore the functions of robust DEGs. The results of enrichment analysis showed that the robust DEGs were closely associated with osteosarcoma development and progression. Immune cell infiltration analysis was also conducted by CIBERSORT algorithm, and we found that macrophages are the most principal infiltrating immune cells in osteosarcoma, especially macrophages M0 and M2. Then, the protein–protein interaction network and key modules were constructed by Cytoscape, and 10 hub genes were selected by plugin cytoHubba from the whole network. The survival analysis of hub genes was also carried out based on the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. The integrated bioinformatics analysis was utilized to provide new insight into osteosarcoma development and metastasis and identified EGR1, CXCL10, MYC, and CXCR4 as potential biomarkers for prognosis of osteosarcoma.

Dysregulated microRNAs in laryngeal cancer: a comprehensive meta-analysis using a robust rank aggregation approach.

We need a reasonable method of compiling data from different studies regarding the expression of microRNA (miRNA) in laryngeal squamous cell carcinoma (LSCC). The robust rank aggregation method was used to integrate the rank lists of miRNAs from 11 studies. The enrichment analysis was performed on target genes of meta-signature miRNAs. The Cancer Genome Atlas database was used to confirm the results of meta-analysis. Three meta-signature miRNAs (miR-21-5p, miR-196a-5p and miR-145-5p) were obtained. All three miRNAs could be prognostic for LSCC. The enrichment analysis showed that these miRNAs were associated significantly with multiple cancer-related signaling pathways. The robust rank aggregation approach is an effective way to identify important miRNAs from different studies. All identified miRNAs could be candidates for LSCC diagnostic and prognostic biomarkers.

Comprehensive analysis of location-specific hub genes related to the pathogenesis of colon cancer.

The molecular mechanisms underlying colon cancer lesions at different sites are not entirely clear. Herein, we aimed to explore location-specific gene profiles related to the pathogenesis of colon cancer and to identify their function. The robust rank aggregation (RRA) method was used to integrate colon cancer microarray datasets and screen differentially expressed gene (DEG) profiles between left- and right-sided colon cancers. Then, weighted gene co-expression network analysis (WGCNA) was performed to cluster the DEGs into modules and identify hub genes. The selected hub genes were validated using The Cancer Genome Atlas dataset and clinical tissues. We assessed the association of selected hub genes with the methylation status in immune cells. In total, 905 DEGs were identified by RRA; five gene modules and 18 hub genes were related to the clinical traits of colon cancer by WGCNA. Four hub genes were selected and shown to be associated with colon cancers on different sides and distant metastasis in the validation analysis. The four hub genes showed a low methylation status, and their expression was significantly associated with methylation status. Positive correlations were observed between the four hub genes and tumor purity and among the four types of immune cells. Gene set enrichment analysis revealed that the four hub genes were mainly involved in two cancer-related pathways. In conclusion, this study identified a set of location-specific genes related to the pathogenesis of colon cancer. These four hub genes may act as novel candidate targets for the treatment of colon cancer.

Integrated bioinformatics analysis for the identification of potential key genes affecting the pathogenesis of clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (CCRCC) is a typical type of RCC with the worst prognosis among the common epithelial neoplasms of the kidney. However, its molecular pathogenesis remains unknown. Therefore, the aim of the present study was to screen for effective and potential pathogenic biomarkers of CCRCC. The gene expression profile of the GSE16441, GSE36895, GSE40435, GSE46699, GSE66270 and GSE71963 datasets were downloaded from the Gene Expression Omnibus database. First, the limma package in R language was used to identify differentially expressed genes (DEGs) in each dataset. The robust and strong DEGs were explored using the robust rank aggregation method. A total of 980 markedly robust DEGs were identified (429 upregulated and 551 downregulated). According to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, these DEGs exhibited an obvious enrichment in various cancer-related biological pathways and functions. The Search Tool for the Retrieval of Interacting Genes/Proteins database was used for the construction of a protein-protein interaction (PPI) network, the Cytoscape MCODE plug-in for module analysis and the cytoHubba plug-in to identify hub genes from the aforementioned DEGs. A total of four key modules were identified in the PPI network. A total of six hub genes, including C-X-C motif chemokine ligand 12, bradykinin receptor B2, adenylate cyclase 7, calcium sensing receptor (CASR), kininogen 1 and lysophosphatidic acid receptor 5, were identified. The DEG results of the hub genes were verified using The Cancer Genome Atlas database, and CASR was found to be significantly associated with the prognosis of patients with CCRCC. In conclusion, the present study provided new insight and potential biomarkers for the diagnosis and prognosis of CCRCC.

A robust twelve-gene signature for prognosis prediction of hepatocellular carcinoma

BackgroundThe prognosis of hepatocellular carcinoma (HCC) patients remains poor. Identifying prognostic markers to stratify HCC patients might help to improve their outcomes.MethodsSix gene expression profiles (GSE121248, GSE84402, GSE65372, GSE51401, GSE45267 and GSE14520) were obtained for differentially expressed genes (DEGs) analysis between HCC tissues and non-tumor tissues. To identify the prognostic genes and establish risk score model, univariable Cox regression survival analysis and Lasso-penalized Cox regression analysis were performed based on the integrated DEGs by robust rank aggregation method. Then Kaplan–Meier and time-dependent receiver operating characteristic (ROC) curves were generated to validate the prognostic performance of risk score in training datasets and validation datasets. Multivariable Cox regression analysis was used to identify independent prognostic factors in liver cancer. A prognostic nomogram was constructed based on The Cancer Genome Atlas (TCGA) dataset. Finally, the correlation between DNA methylation and prognosis-related genes was analyzed.ResultsA twelve-gene signature including SPP1, KIF20A, HMMR, TPX2, LAPTM4B, TTK, MAGEA6, ANX10, LECT2, CYP2C9, RDH16 and LCAT was identified, and risk score was calculated by corresponding coefficients. The risk score model showed a strong diagnosis performance to distinguish HCC from normal samples. The HCC patients were stratified into high-risk and low-risk group based on the cutoff value of risk score. The Kaplan–Meier survival curves revealed significantly favorable overall survival in groups with lower risk score (P < 0.0001). Time-dependent ROC analysis showed well prognostic performance of the twelve-gene signature, which was comparable or superior to AJCC stage at predicting 1-, 3-, and 5-year overall survival. In addition, the twelve-gene signature was independent with other clinical factors and performed better in predicting overall survival after combining with age and AJCC stage by nomogram. Moreover, most of the prognostic twelve genes were negatively correlated with DNA methylation in HCC tissues, which SPP1 and LCAT were identified as the DNA methylation-driven genes.ConclusionsWe identified a twelve-gene signature as a robust marker with great potential for clinical application in risk stratification and overall survival prediction in HCC patients.

MicroRNA expression integrated analysis and identification of novel biomarkers in small cell lung cancer: a meta-analysis.

BackgroundSmall cell lung cancer (SCLC) is an aggressive and recalcitrant cancer. In recent years, studies focused on the abnormal expression of microRNA which has proven valuable in terms of prognosis, diagnosis and treatment in SCLC. To address the limitations of independent studies data, a meta-analysis seems necessary for further exploration of microRNA as biological target and regulatory factor in SCLC.MethodsWe performed comprehensive literature retrieval in GEO database and EBI ArrayExpress database. The microRNA expression data was extracted from 4 related researches (GSE15008, GSE74190, GSE19945, GSE77380), which was obtained from GEO database. In each included study, the R. Affymetrix Expression Console’s Limma package and RMA algorithms were used to screen for raw data for gene chip quality control, standardization, log2 conversion and differential expression of the gene chip, respectively. Significant microRNA meta-signatures were identified by Robust Rank Aggregation method. Subsequently, gene ontology (GO) enrichment analysis and pathway analysis were performed using bioinformatics tools.ResultsWe found a significant microRNA meta-signature of six up-regulated (hsa-miR-182-5p, hsa-miR-96-5p, hsa-miR-7-5p, hsa-miR-301b-3p, hsa-miR-130b-3p, hsa-miR-210-3p) and four down-regulated (hsa-miR-126-3p, hsa-miR-451a, hsa-miR-145-5p, hsa-miR-486-5p) microRNA s in meta-analysis approaches. GO analysis showed that target gene of meta-signatures microRNA was mainly enriched in endosome, chordate embryonic development and transforming growth factor beta receptor. The related functional gene of microRNA meta signature synergistically targeting SCLC signaling pathway was confirmed by enrichment analysis. In particular, neurotrophin and TGF-beta signaling pathway play the most important roles in the pathway network.ConclusionsOur study identified 10 highly significant and consistently dysregulated microRNA s from 4 datasets, which offering convincing molecular targets and regulatory factors in future research of SCLC.

Is Rank Aggregation Effective in Recommender Systems? An Experimental Analysis

Recommender Systems are tools designed to help users find relevant information from the myriad of content available online. They work by actively suggesting items that are relevant to users according to their historical preferences or observed actions. Among recommender systems, top- N recommenders work by suggesting a ranking of N items that can be of interest to a user. Although a significant number of top- N recommenders have been proposed in the literature, they often disagree in their returned rankings, offering an opportunity for improving the final recommendation ranking by aggregating the outputs of different algorithms. Rank aggregation was successfully used in a significant number of areas, but only a few rank aggregation methods have been proposed in the recommender systems literature. Furthermore, there is a lack of studies regarding rankings’ characteristics and their possible impacts on the improvements achieved through rank aggregation. This work presents an extensive two-phase experimental analysis of rank aggregation in recommender systems. In the first phase, we investigate the characteristics of rankings recommended by 15 different top- N recommender algorithms regarding agreement and diversity. In the second phase, we look at the results of 19 rank aggregation methods and identify different scenarios where they perform best or worst according to the input rankings’ characteristics. Our results show that supervised rank aggregation methods provide improvements in the results of the recommended rankings in six out of seven datasets. These methods provide robustness even in the presence of a big set of weak recommendation rankings. However, in cases where there was a set of non-diverse high-quality input rankings, supervised and unsupervised algorithms produced similar results. In these cases, we can avoid the cost of the former in favor of the latter.

Identification of candidate RNA signatures in triple-negative breast cancer by the construction of a competing endogenous RNA network with integrative analyses of Gene Expression Omnibus and The Cancer Genome Atlas data.

Triple-negative breast cancer (TNBC) is a subtype of breast cancer that is characterized by aggressive and metastatic clinical characteristics and generally leads to earlier distant recurrence and poorer prognosis than other molecular subtypes. Accumulating evidence has demonstrated that long non-coding RNAs (lncRNAs) serve a crucial role in a wide variety of biological processes by interacting with microRNAs (miRNAs) as competing endogenous RNAs (ceRNAs) and, thus, affect the expression of target genes in multiple types of cancer. Seven datasets from the Gene Expression Omnibus (GEO) database, including 444 tumor and 88 healthy tissue samples, were utilized to investigate the underlying mechanisms of TNBC and identify prognostic biomarkers. Differentially expressed genes (DEGs) were further validated in The Cancer Genome Atlas database and the associations between their expression levels and clinical information were analyzed to identify prognostic values. A potential lncRNA-miRNA-mRNA ceRNA network was also constructed. Finally, 69 mRNAs from the integrated Gene Expression Omnibus datasets were identified as DEGs using the robust rank aggregation method with |log2FC|>1 and adjusted P<0.01 set as the significance cut-off levels. In addition, 29 lncRNAs, 21 miRNAs and 27 mRNAs were included in the construction of the ceRNA network. The present study elucidated the mechanisms underlying the progression of TNBC and identified novel prognostic biomarkers for TNBC.

Integrated bioinformatics analysis reveals the response of Arabidopsis roots to salt stress through multiple gene expression omnibus datasets

ABSTRACT Plant growth and crop productivity has been adversely affected by salt stress around the world. Arabidopsis root is an important model to demonstrate how the glycophyte adapts to salt stress. To better understand the salt-responsive genes in Arabidopsis root, a rank aggregation method has been applied to disclose 836 dysregulated genes including 449 upregulated and 387 down-regulated genes out of more than 20,000 genes in NaCl-treated Arabidopsis root from the Gene Expression Omnibus (GEO) database. Then, the integrated approaches of PANTHER, DAVID, KEGG, Cytoscape and STRING have been used for bioinformatic analysis. Our results showed that these salt-responsive genes were involved in various gene ontology enrichments (transcription, oxidative stress, cell wall organization), differential pathways (MAPK signaling pathway, hormone signal transduction), and complicated interactions (cross-talking with wounding stress, abscisic acid stress, heat stress, water deprivation stress, cold stress). Furthermore, the activities of the oxidative stress related gene, peroxidase, were determined under salt stress, cold stress and water deprivation stress, which validated the bioinformatic outcome of cross-tolerance to different stresses. Our work will throw new lights on the response of Arabidopsis roots to salt stress, and provide potential salt-targets to improve plant growth and increase crop output.