Review Special Issues

Towards defining the Mechanisms of Alzheimer’s disease based on a contextual analysis of molecular pathways

  • Received: 15 January 2016 Accepted: 17 March 2016 Published: 21 March 2016
  • Alzheimer’s disease (AD) is posing an increasingly profound problem to society. Our genuine understanding of the pathogenesis of AD is inadequate and as a consequence, diagnostic and therapeutic strategies are currently insufficient. The understandable focus of many studies is the identification of molecules with high diagnostic utility however the opportunity to obtain a further understanding of the mechanistic origins of the disease from such putative biomarkers is often overlooked. This study examines the involvement of biomarkers in AD to shed light on potential mechanisms and pathways through which they are implicated in the pathology of this devastating neurodegenerative disorder. The computational tools required to analyse ever-growing datasets in the context of AD are also discussed.

    Citation: Joanna L. Richens, Jonathan P. Bramble, Hannah L. Spencer, Fiona Cantlay, Molly Butler, Paul O’Shea. Towards defining the Mechanisms of Alzheimer’s disease based on a contextual analysis of molecular pathways[J]. AIMS Genetics, 2016, 3(1): 25-48. doi: 10.3934/genet.2016.1.25

    Related Papers:

    [1] Keqiang Li, Shangjiu Wang . Symmetry of positive solutions of a $ p $-Laplace equation with convex nonlinearites. AIMS Mathematics, 2023, 8(6): 13425-13431. doi: 10.3934/math.2023680
    [2] Keqiang Li, Shangjiu Wang, Shaoyong Li . Symmetry of large solutions for semilinear elliptic equations in a symmetric convex domain. AIMS Mathematics, 2022, 7(6): 10860-10866. doi: 10.3934/math.2022607
    [3] Zihan Yue, Wei Jiang, Boying Wu, Biao Zhang . A meshless method based on the Laplace transform for multi-term time-space fractional diffusion equation. AIMS Mathematics, 2024, 9(3): 7040-7062. doi: 10.3934/math.2024343
    [4] Yuhua Long, Sha Li . Results on homoclinic solutions of a partial difference equation involving the mean curvature operator. AIMS Mathematics, 2025, 10(3): 6429-6447. doi: 10.3934/math.2025293
    [5] Song Wang, Xiao-Bao Shu, Linxin Shu . Existence of solutions to a class of damped random impulsive differential equations under Dirichlet boundary value conditions. AIMS Mathematics, 2022, 7(5): 7685-7705. doi: 10.3934/math.2022431
    [6] Amjad Hussain, Muhammad Khubaib Zia, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar, Ilyas Khan . Lie analysis, conserved vectors, nonlinear self-adjoint classification and exact solutions of generalized $ \left(N+1\right) $-dimensional nonlinear Boussinesq equation. AIMS Mathematics, 2022, 7(7): 13139-13168. doi: 10.3934/math.2022725
    [7] Xiangmei Li, Kamran, Absar Ul Haq, Xiujun Zhang . Numerical solution of the linear time fractional Klein-Gordon equation using transform based localized RBF method and quadrature. AIMS Mathematics, 2020, 5(5): 5287-5308. doi: 10.3934/math.2020339
    [8] Dilara Siraeva . The invariant solution with blow-up for the gas dynamics equations from one-parameter three-dimensional subalgebra consisting of space, Galilean and pressure translations. AIMS Mathematics, 2024, 9(1): 89-98. doi: 10.3934/math.2024006
    [9] Khalid K. Ali, Mohamed S. Mohamed, Weam G. Alharbi, M. Maneea . Solving the time fractional q-deformed tanh-Gordon equation: A theoretical analysis using controlled Picard's transform method. AIMS Mathematics, 2024, 9(9): 24654-24676. doi: 10.3934/math.20241201
    [10] Zhilin Li, Guoping Chen, Weiwei Long, Xinyuan Pan . Variational approach to p-Laplacian fractional differential equations with instantaneous and non-instantaneous impulses. AIMS Mathematics, 2022, 7(9): 16986-17000. doi: 10.3934/math.2022933
  • Alzheimer’s disease (AD) is posing an increasingly profound problem to society. Our genuine understanding of the pathogenesis of AD is inadequate and as a consequence, diagnostic and therapeutic strategies are currently insufficient. The understandable focus of many studies is the identification of molecules with high diagnostic utility however the opportunity to obtain a further understanding of the mechanistic origins of the disease from such putative biomarkers is often overlooked. This study examines the involvement of biomarkers in AD to shed light on potential mechanisms and pathways through which they are implicated in the pathology of this devastating neurodegenerative disorder. The computational tools required to analyse ever-growing datasets in the context of AD are also discussed.


    Lung cancer is one of the most common cancers that endanger human health seriously, and approximately one million people die from it every year [1,2]. Lung adenocarcinoma (LUAD) is widely regarded as the most common histological subtype of lung cancer [3]. It is reported that the 5-year survival rate for LUAD patients in the United States is below 20% [4]. Over the past few decades, despite the progress has been made in the treatment of LUAD, the overall survival rate of this cancer is not evidently upgraded [5]. Recently, it has been proved that tremendous advance in cancer treatment has been obtained through immunotherapy, and immune checkpoint inhibitors are increasingly used to treat LUAD [6].

    Evidence has indicated that targeting programmed cell-death protein 1/programmed cell death 1 ligand 1 (PD-1/PD-L1)-related pathways is an emerging concept in cancer immunotherapy [7]. PD-1 (also called PDCD-1) is expressed on both exhausted and activated T cells, which is an immune checkpoint molecular [8]. PD-1 blockade had remarkable efficacy against extensive advanced malignancies [9], and it also serves critical role in lung cancer. For example, Kagamu et al. [10] observed that many advanced lung cancer patients achieved long-term survival without disease progression after treatment with PD-1 blockade, suggesting that PD-1 blockade therapy could evidently improve the outcome of lung cancer patients. In addition, with the development of gene expression profiles, researchers have discovered several gene expression characteristics related to survival outcomes of various cancers. The constructed prognostic model may be helpful for early diagnosis of LUAD and personalized treatment [11]. For instance, Zhang et al. [12] constructed a risk model with four genes that could predict the overall survival of LUAD patients. Moreover, Sun et al. [13] established an immune prognostic model for LUAD, which could independently identify high-risk patients with poor survival. However, the PD-1 associated risk model that was used for predicting the prognosis of LUAD has not been constructed.

    In the present study, we aimed to identify PD-1-related gene signature to predict the survival of LUAD patients. The gene expression data and clinical information related to LUAD were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. After a series of bioinformatics analysis, the genes related to both PD-1 expression level and immune cells were identified, and then Lasso Cox regression analysis was performed to obtain key genes to construct prognostic model. Furthermore, the predictive performance of this model was validated by using the data from an independent GEO and TCGA cohorts. Figure S1 shows the analysis process of the present study. This study will provide new insights into the mechanisms of PD-1 immune checkpoints and provide novel therapeutic targets for LUAD patients with poor prognosis.

    The National Center for Biotechnology information (NCBI) GEO database [14] was used to search the gene expression profiles of LUAD with "lung adenocarcinoma" as a keyword. The search criteria were as follows: 1) the species was set to "Homo sapiens" 2) the clinical parameters of the samples in dataset must include prognostic information, 3) the sample size of cancer tissues in the dataset was at least 60. Finally, four datasets (GSE26939, GSE68465, GSE72094, and GSE68571) that met the requirements were obtained and included in this analysis.

    Concretely, GSE26939 dataset included 116 LUAD samples, among which, one sample without prognostic information was removed. The platform was GPL9053 Agilent-UNC-custom-4X44K. GSE68465 dataset included 443 LUAD samples and 19 normal lung tissue samples. Of which, 442 samples with complete prognostic information were selected for further analysis. The platform was GPL96 [HG-U133A] Affymetrix Human Genome U133A Array. GSE72094 dataset included 442 LUAD samples. After filtering, 398 samples were enrolled for analysis. The platform was GPL15048 Rosetta/Merck Human RSTA Custom Affymetrix 2.0 microarray [HuRSTA_2a520709.CDF]. In addition, GSE68571 dataset included 86 LUAD samples, which were analyzed based on GPL80 [Hu6800] Affymetrix Human Full Length HuGeneFL Array platform.

    For these four datasets, standardized probe expression matrix and the corresponding annotation information were downloaded from GEO database. Then, probes that not matched to gene symbol were removed. If different probes mapped to the same gene symbol, the average value of these probes was applied as the final expression value.

    Combined with the expression level of PD-1 in each sample and the corresponding survival information, K-M survival analysis was performed by using R package survminer (Version 0.4.3) to determine the optimal cutpoint. According to the optimal cutpoint, samples were divided into PD-1 low expression group (PD-1 expression level in sample < cutpoint) and high expression group (PD-1 expression level in sample > cutpoint). Next, survival analysis and log-rank test were used to compare the prognostic value of low and high expression groups by using the R package survival (Version 2.42-6).

    The distribution of clinical factors in each sample of the PD-1 low and high expression groups was counted and plotted using R package ggstatsplot (Version: 0.6.5). In addition, the difference in continuous variables between two groups was compared using Welch's t-test, and the difference in discrete variables between two groups was calculated using Chi-square test. Clinical factors with p-value < 0.05 were regarded as statistically significant.

    R package gene set variation analysis (GSVA, Version: 1.36.2) [15] was applied to calculate the enrichment score of each Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway in each sample to obtain a scoring matrix. Then, differential pathway between PD-1 low and high expression groups was identified by using R package limma [16] (Version 3.10.3). Then, Benjamini- Hochberg method was used to perform multiple test correction and obtained adjusted p-value (adj. p-value). The adj. p-value < 0.05 was considered as the cut-off criterion of significantly enriched pathway.

    DEGs between PD-1 high expression and low expression groups were screened by using the Bayesian method provided by limma in R package. Further, Benjamini- Hochberg method was used for the adjustment of p-value. |log fold change (FC)| > 0.585 and adj. p-value < 0.05 were considered as statistically significant.

    According to the expression matrix of all genes, variance analysis was adopted to find the TOP2000 genes with higher degree of variation among samples. Then, these genes were analyzed by using R package weighted gene co-expression network analysis (WGCNA) [17] (Versions 1.61 and modules closely related to PD-1 were identified. After that, R package clusterprofiler [18] (Version: 3.8.1) was used to perform KEGG pathway enrichment analysis of genes in these modules with the following parameters: pAdjustMethod = "BH" and pvalueCutoff = "0.05". P-value < 0.05 was the threshold value of significantly enriched results.

    Subsequently, we narrowed down the range of candidate genes to obtain closely related genes. In brief, the above obtained DEGs and genes in modules were integrated, and then the overlapping genes were selected for further analysis.

    The gene sets used to label the type of each immune cell in tumor microenvironment (TME) were obtained from the study by Charoentong et al. [19]. Based on the ssGSEA algorithm in R package GSVA [20], the enrichment score of each immune cell in samples was calculated, which represented the relative abundance of each TME infiltrating cell in each sample. Then, based on the relative abundance, pheatmap package (Version: 1.0.12) was used to perform unsupervised clustering for samples, and the samples were clarified into high and low infiltrating immune groups. After that, DEGs between these two groups were screened using limma in R package. The threshold values of DEGs were |log FC| > 1 and adj. p-value < 0.05, and these genes were considered to be closely related to the immune microenvironment.

    Genes related to PD-1 expression as well as immune cells were integrated, and the shared genes were obtained for further analyses. Next, Gene ontology (GO)-biological process (BP) [21] and KEGG [22] pathway analyses of these genes were performed using DAVID [23] (Version 6.8). Count ≥ 2 and p-value < 0.05 were regarded as statistically significant.

    Combined with the expression values of genes and the overall survival time of each sample, univariate Cox regression analysis was conducted. Genes with p-value < 0.05 were regarded as prognosis-related genes and included for the next analysis.

    Then, based on these prognosis-related genes, a LASSO Cox regression was applied to identify key genes for the construction of prognostic models by using glmnet package [24] (Version 2.0-18). The 50-fold cross-validation was performed to determine the optimal value of Lasso penalty parameter. After that, multivariate Cox regression analysis of these key genes was conducted to construct prognostic risk model using the following formula: Risk score = βgene1*exprgene1 + βgene2*exprgene2 +... +βgeneN*exprgeneN (β represents correlation coefficient of prognostic gene and expr represents the expression value of the corresponding genes). According to the median value of risk score, the samples were divided into low risk group and high risk group. Next, K-M survival analysis of these two groups was performed using log-rank test.

    To determine whether the risks core model could be applied as an independent prognostic factor, the prognostic value of the risk score and other clinical parameters, including age, sex, stage, grade, marked lymphocytes, smoking status, TP53 mutation, EGFR mutation, KRAS mutation, and STK11 mutation, was evaluated with univariate and multivariate Cox regression analyses. Variables with p-value < 0.05 were regarded as independent prognostic factors in both analyses.

    We used a microarray datasets GSE72094 to validate the stability and repeatability of prognostic model. At first, the expression value of model genes in each sample of GSE72094 was extracted. According to the prognostic coefficient obtained from the previous analysis, the risk score was calculated using the same formula. Then, the samples were divided into low risk group and high risk group based on the median, followed by K-M survival analysis. In addition to the GEO dataset, we also included data from TCGA to validate the prediction performance of the model. In brief, the RNA-seq data of LUAD and survival information of corresponding samples were downloaded from the UCSC Xena platform (https://toil.xenahubs.net). In the previous study, the samples with number "-01A", that is, the primary tumor samples were selected for subsequent analysis. Finally, a total of 497 samples were obtained. Based on the expression value of seven genes in model, the risk score was calculated using the same formula. Similarly, all samples were assigned to low-risk and high-risk groups according to the median value of risk score, followed by a survival analysis.

    Moreover, we compared the predictive ability of the seven-gene prognostic model with other molecular biomarkers such as PD-1 and PD-L1 in the GSE26939 dataset. Meanwhile, a 1-, 3-, and 5-year receiver operating characteristic (ROC) curves were plotted using the R package timeROC [25], and then the area under the curve (AUC) value of the survival ROC curves was calculated.

    The clinical information of cases in four datasets is presented in Table S1. As shown in Figure 1, results showed that the survival probability of patients in the PD-1 high expression group was higher than that of low expression group. Despite there was no significant difference in survival time between PD-1 high expression group and low expression group in the GSE68465 and GSE72094 datasets, the overall trend was consistent.

    Figure 1.  K-M survival curves of patients in PD-1 high expression group and low expression group of GSE26939 (A), GSE68465 (B), GSE68571 (C), and GSE72094 (D).

    According to the analysis results in the previous step, we observed that the number of samples in the PD-1 low expression group of GSE68571 and GSE72094 was significantly lower than that in the high expression group. In order to reduce the impact of the imbalance of samples in the two groups on the results, we only selected samples in the GSE26939 and GSE68465 datasets for further analyses.

    The distribution of different clinical factors in the PD-1 high and low expression groups in these two datasets was calculated, and factors with p-value < 0.05 were obtained. In brief, in GSE26939 dataset, the age of the PD-1 high expression group was evidently lower than that of the low expression group (Figure 2A, p value = 0.004); the high expression group had more TP53 wild type (Figure 2B, p-value = 0.037); and the PD-1 high expression group had more Squamoid subtype (Figure 2C, p-value ≤ 0.001). Meanwhile, in the dataset of GSE68465, the age of patients in the PD-1 high expression group was markedly younger than that of the low expression group (Figure 2D, p-value = 0.007); the high expression group had more T1-T2 samples and less T3-T4 samples than the low expression group (Figure 2E, p-value = 0.031); and PD-1 high expression group had fewer white people than low expression group (Figure 2F, p-value ≤ 0.001).

    Figure 2.  The distribution of different clinical factors between PD-1 high expression group and low expression group in the datasets (p-value < 0.05). A: Age in the GSE26939 dataset. B: TP53 mutation status in the GSE26939 dataset. C: Subtype in the GSE26939 dataset. D: Age in the GSE68465 dataset. E: T stage in the GSE68465 dataset F: Race in the GSE68465 dataset.

    Afterward, GSVA was performed to explore difference in KEGG pathways between the PD-1 high expression and low expression groups from the GSE26939 as well as GSE68465 datasets. In total, there were 72 KEGG pathways enriched by samples in GSE26939 dataset and 137 KEGG pathways in GSE68465 dataset. The top20 pathways are presented in Figure 3. In the GSE26939 dataset, several pathways related to immunity were observed (Figure 3A), such as immunodeficiency, intestinal immune network for IgA production, and primary immunodeficiency. These pathways were increased in high expression group in comparison with that in low expression group. Meanwhile, in the GSE68465 dataset, 12 pathways were significantly higher in the PD-1 high expression group than in the low expression group, such as type II diabetes mellitus pathway and JAK-STAT-signaling pathway (Figure 3B). Notably, JAK-STAT-signaling pathway was evidently increased in the PD-1 high expression group of two GEO datasets.

    Figure 3.  The top20 results of KEGG pathways enriched in GSE26939 dataset (A) and GSE68465 dataset (B). The closer the color is to yellow, the more significant the p-value.

    Combined with the results of the previous analyses, we found that the expression level of PD-1 in GSE26939 was significantly related to prognosis of LUAD, and samples in this dataset were mainly enriched in immune-related pathways. Therefore, GSE26939 dataset was retained for further analysis. According to the threshold (|log FC| > 0.585 and adj. p-value < 0.05), 822 DEGs (709 up-regulated and 113 down-regulated genes) were screened between PD-1 high expression group and low expression group. As shown in Figure S2, DEGs could distinguish the PD-1 low expression group from the high expression group, which proved initially that these DEGs were associated with the expression level of PD-1.

    WGCNA was used to identify PD-1 expression related modules. This analysis was performed by a step-by-step method. First, soft threshold power of 8 was selected to calculate the adjacencies because it was the lowest power at which the scale-free topology fit index reached 0.85 (Figure 4A, B). Second, based on clustering and dynamic tree cut methods, the similar modules were clustered. Modules with correlation coefficient > 0.7 were selected and finally seven modules were obtained (Figure 4C). Next, the relationship between module and clinical characteristics (age, TP53 mutation status, PD-1 high- or low-expressed groups, and PD-1 expression level) was analyzed (Figure 4D). For the expression level of PD-1, magenta module, pink module, red module, and purple module represented evidently positive correlations (r > 0 and p-value < 0.05), while yellow module showed significantly negative correlation (r < 0 and p-value < 0.05). In addition, KEGG pathway enrichment analysis for genes in five modules showed that magenta module, pink module, red module, purple module, and yellow module were respectively enriched in 4, 39, 12, 52, and 8 pathways (Figure S3). For example, genes in magenta module, pink module, red module, purple module, and yellow module were mainly enriched in primary immunodeficiency, cytokine-cytokine receptor interaction, protein digestion and absorption, TNF signaling pathway, and cholesterol metabolism, respectively.

    Figure 4.  Weighted gene co-expression network construction. Analyses of network topologies for various soft-thresholding powers via scale-free fit index (A) and mean connectivity (B). C: Clustering dendrogram of genes based on topological overlapping. Different colors represent different modules, and the gray part represents genes that cannot be merged into any other modules. D: Correlation matrix for each module and clinical characteristics. The upper number represents the correlation coefficient, and the lower number in parenthesis represents the p value.

    Furthermore, genes in the five models and DEGs in PD-1 low expression vs. high expression groups were intersected, and 388 overlapping genes were obtained (Figure S4), which were regarded as candidate genes for further analysis.

    As presented in Figure S5A, the samples were divided into low immune infiltration group and high immune infiltration group according to the enrichment fraction of immune cells in each sample. We also observed that the high immune infiltration group was mainly composed of samples with high expression level of PD-1, and the low immune infiltration group was mainly contained samples with low expression level of PD-1. Furthermore, differential analysis was performed between high immune infiltration and low immune infiltration groups. A total of 467 DEGs were identified, including 445 up-regulated and 22 down-regulated genes.

    A total of common 237 genes were obtained via the intersection of 467 immune related genes and 388 PD-1 related genes. Thus, these 237 genes were considered as PD-1 related immune genes. To explore the biological function of these genes, GO_BP terms and KEGG pathways analyses were conducted. Results showed that these 237 genes were significantly enriched in 87 GO-BP terms (such as immune response and inflammatory response) and 38 KEGG pathways (such as cytokine-cytokine receptor interaction and rheumatoid arthritis). The top20 results are presented in Figure S5B.

    Based on the 237 PD-1 related immune genes, univariate Cox regression analysis was performed. A total of 47 genes that were significantly associated with overall survival of LUAD were obtained. Next, LASSO Cox regression analysis for these 47 genes was conducted to reveal the optimal modeling parameter (λ), and results showed λ = 7 was considered as optimum parameter (Figure 5A). Finally, a seven-gene (DPT, ITGAD, CLECL1, SYT13, DUSP26, AMPD1, and NELL2) signature was used to construct the prognostic model. Risk score was calculated by using the following formula: (-9.322*exprAMPD1) + (-0.169*exprCLECL1) + (-0.176*exprDPT) + (-0.029*exprDUSP26) + (-0.056*exprITGAD) + (-0.063*exprNELL2) + (0.077*exprSYT13). Based on the median of risk score, samples were divided into low risk group and high risk group. Patients in high risk group had significantly worse survival rate than those in low risk group (Figure 5B). Distributions of risk score, follow-up time, and seven-gene expression profiles are shown in Figure 5C. The expression level of SYT13 was up-regulated in the high risk group, and the expression level of the remaining six genes was down-regulated in the high risk group compared with the low risk group. In addition, K-M survival curve showed that the prognosis of patients with high expression of SYT13 was poor, and the prognosis of patients with low expression of the remaining 6 genes was worse (Figure 5D).

    Figure 5.  Prediction performances of prognostic model in test cohort (GSE26939). A: The selection of optimal modeling parameters λ in LASSO model. The two dashed lines respectively indicate two special λ values, lambda.min on the left and lambda.1se on the right. The lambda between these two values is considered appropriate. The model constructed by lambda.1se is the simplest, that is, the number of genes is small, while lambda.min has a higher accuracy rate and uses more genes. B: K-M survival curve for high risk group and low risk group. C: Risk score distribution chart and heatmap of 7 genes in each sample. D: K-M survival curve for these seven genes.

    Further, univariate and multivariate Cox regression analyses were performed to evaluate whether risk score could be used as an independent prognostic factor for LUAD. Several clinical factors were included in this analysis, including sex, age, stage, grade, marked_lymphocytes, smoking_status, and genetic mutation status (TP53, EGFR, KRAS, STK11). Univariate Cox analysis indicated that sex, age, and risk groups were significantly associated with prognosis of LUAD (Figure 6A). Next, multivariate analysis revealed that the risk score was an independent prognostic predictor (Figure 6B), suggesting that the constructed model could independently predict the prognosis of LUAD.

    Figure 6.  The forest plot of univariate (A) and multivariate (B) Cox regression analysis.

    According to the seven-gene risk model, the risk score of the 398 LUAD patients from GSE72094 was calculated using the same formula. Then, these patients were classified into low and high risk groups according to the median of risk score. K-M survival curve presented that the prognosis of patients in high risk group was significantly worse than those in low risk group (Figure 7A). The risk scores distribution and follow up of the patients in these two groups were analyzed, indicating that the high risk score group had more mortality cases (Figure 7B). Moreover, the expression level of these seven genes in the two groups was consistent with the results of above analysis (Figure 7C). These findings suggested that the constructed prognostic model was reliable and efficient.

    Figure 7.  Predictive performances of prognostic model in validation cohort. A: K-M survival curve for high and low risk groups of GSE72094 dataset. B: The distribution of risk score as well as patients' follow-up time in GSE72094. C: The heatmap of expression profile of seven genes in GSE72094. D: K-M survival of prognostic model in the TCGA cohort. E: The distribution of risk score and follow-up time of cases in the TCGA cohort. F: The heatmap of seven PD-1 related prognostic genes.

    To further determine whether the prognostic model was reliable in different populations, we executed the same formula and determined the cut-off value in the TCGA cohort. Next, patients were divided into low- or high-risk groups. Consistent with the results in the GSE72094, patients with high risk score had a lower overall survival than those with low risk score (Figure 7D). The distributions of risk score of LUAD patients as well as the relationship between risk score and follow up were visualized in Figure 7E. In addition, the expression level of SYT13 was up-regulated in the high risk group, and the level of the remaining six genes was down-regulated in the high risk group (Figure 7F), which was consistent with that in the GEO datasets.

    Next, we plotted the ROC curves for prognostic model, PD-1 and PD-L1, and then the AUC values were calculated. Results showed that AUC values for prognostic model, PD-1, and PD-L1 were respectively 0.637, 0.541, and 0.504 (Figure 8AC), indicating that the predictive ability of prognostic model was better than that of PD-1 and PD-L1. In addition, the accuracy of these three biomarkers for survival of LUAD patients 1, 3, and 5 years were evaluated. The AUC values of the ROC curves for seven-gene, PD-1, and PD-L1 predicted the 1-year overall survival were 0.764, 0.363, and 0.344, respectively (Figure 8D); AUC values of ROC predicted the 3-year survival were 0.671, 0.446, and 0.442, respectively (Figure 8E); and AUC values of ROC predicted the 5-year survival were 0.702, 0.423, and 0.418, respectively (Figure 8F). These results further confirmed that the constructed model outperformed other biomarkers with superior predictive performance.

    Figure 8.  Comparisons of prognostic model with other PD-1 related biomarkers in GSE26939. ROC curve of prognostic model (A), PD-1 (B), and PD-L1 (C). ROC curves predicted the overall survival of LUAD patients 1 (D), 3 (E), and 5 years (F).

    In the present study, a PD-1 expression-related immune prognostic model was constructed, which included seven genes (DPT, ITGAD, CLECL1, SYT13, DUSP26, AMPD1, and NELL2). All samples were divided into high-risk and low-risk groups based on the risk score for this model. Notably, the survival time of the patients in the high risk group was significantly lower than those in the low risk group. Meanwhile, this prognostic model was verified to be effective and reliable by other GEO dataset and TCGA cohort. In addition, survival analysis suggested that the prognosis of patients with higher expression level of SYT13 was poor, and the survival time of patients with lower expression level of the remaining 6 genes was worse. Our constructed PD-1-related prognostic model may be helpful for early diagnosis and treatment of LUAD in clinical practice.

    DPT, a noncollagenous extracellular matrix (ECM) protein, is involved in ECM assembly activities and cellular adhesion. Yamatoji et al. [26] indicated that DPT played a significant role in tumor metastasis and invasiveness of oral squamous cell carcinomas (OSCCs), and might be a therapeutic target for the treatment of OSCCs. Guo et al. [27] suggested that DPT could suppress the proliferation of papillary thyroid carcinoma cells via repression of MYC. Some previous studies also suggested that DPT was differentially expressed in some other cancers, such as hepatocellular carcinoma [28]. ITGAD belongs to the beta-2 integrin family of membrane glycoproteins, and is expressed in the tissue and circulating myeloid leukocytes. Evidences demonstrated that ITGAD was associated with metastatic gastric carcinoma [29] and histiocytic sarcomas [30]. However, there was no previous research describing the association between ITGAD and LUAD or lung cancer. Another gene, SYT13, encodes a member of the large synaptotagmin protein family, and it is considered as a target for the treatment of peritoneal metastasis of gastric cancer [31]. Jahn et al. [32] identified SYT13 as a liver tumor suppressor gene. Meanwhile, silencing of SYT13 promoted apoptosis and suppressed proliferation of colorectal cancer cells and thus inhibited tumor growth [33]. The study of Zhang et al. [34] suggested that SYT13 was involved in the growth, invasion, migration, and apoptosis of LUAD. Together, these studies further highlighted that SYT13 was closely associated with the development of cancer. Recently, DUSP26 has been proved to be a potential therapeutic target for human cancer [35]. Bourgonje et al. [36] indicated that lower expression level of DUSP26 was correlated with poor prognosis of gliomas and overexpression of it in E98 glioblastoma cells could reduce tumorigenicity. In addition, the findings of Yu et al. [37] suggested that DUSP26 might be regarded as an oncogene in anaplastic thyroid cancer. For the NELL2, it encodes a glycoprotein containing several von Willebrand factor C domains and epidermal growth factor (EGF)-like domains. Study in mice has shown that the protein plays a role in the growth and differentiation of nerve cells and tumorigenesis [38]. Nakamura et al. [39] indicated that NELL2 played roles in the inhibition of cell migration of renal cell carcinoma. Meanwhile, Kim et al. [40] confirmed that NELL2 was involved in cancer development by the regulation of E2F transcription factor 1. Taken together, these studies suggested that genes such as DPT, ITGAD, SYT13, DUSP26, and NELL2 were involved in pathogenesis of cancers, and were associated with the prognosis of various cancers. Nevertheless, few studies reported the roles of CLECL1 and AMPD1 in cancer especially in the LUAD. Thus, we speculated that these genes also could predict the prognosis of LUAD.

    In this analysis, we also explore the relationship between immune cells and genes. It is reported that altered cell adhesion played promoting roles in immune suppression of cancers [41], and DPT was involved in cellular adhesion [26], suggesting that DPT was associated with immunity in cancers. Integrin adhesion is required for aspects of immune function such as antigen presentation and migration in germinal centers, inflammation sites, and lymph nodes [42]. Synaptotagmin (SYT) is a large class of membrane transporters with 13 members, and is considered to be the main calcium sensor in the process of immune cell exocytosis [43]. Among these, SYT16 is reported to be correlated with immune infiltrates and a prognostic biomarker in glioma [44]. The member of dual-specificity phosphatases (DUSPs) family, such as DUSP1, DUSP2 DUSP4 DUSP10, and DUSP16 are important for immune response and metabolic homeostasis [45]. CLECL1 is involved in the regulating immunity [46], and AMPD1 plays a role in the immune function in sepsis patients [47]. These evidences indicated the indirect or direct relationship between genes in the model and immunity. Therefore, these genes may affect the prognosis of LUAD via regulating the immune microenvironment. However, the specific molecular mechanism still needs further experimental verification.

    Importantly, a model composed of seven PD-1 related immune genes (DPT, ITGAD, CLECL1, SYT13, DUSP26, AMPD1, and NELL2) was constructed, which could distinguish LUAD patients and predict prognosis. After calculation of the risk score, results of survival analysis showed that the overall survival of patients in the high-risk group was significantly worse than those in the low-risk group. This model could be used as an independent prognostic factor for LUAD patients. According to our research, the models constructed from these seven PD-1 related genes can well predict the prognosis of patients with LUAD. We think that these seven genes are biomarkers for predicting the prognosis of LUAD patients, and may become new research targets. In short, our study provided a new method for evaluating the prognosis of LUAD patients.

    Nevertheless, our results have to be interpreted in light of some limitations. First, the molecular functions of these identified genes in LUAD remain largely unknown, which require further studies in vivo and in vitro. Secondly, all genes in the constructed model were derived from the public databases (GEO and TCGA); therefore, more external samples from medical centers are needed to validate the predictive performance of the risk model.

    In summary, a PD-1 related seven-gene (DPT, ITGAD, CLECL1, SYT13, DUSP26, AMPD1, and NELL2) prognostic model was constructed, which could predict the overall survival of LUAD patients well. These seven genes help us to further understand the molecular mechanism of PD-1, and they may be potential biomarkers for the treatment of LUAD. However, the lack of experimental verification is a limitation of the present study, thus further researches are needed to verify the present results.

    This study was supported by the Shanghai Minhang District Health Commission Research Project (No. 2020MW05).

    The authors declare that they have no conflict of interest.

    [1] Strimbu K, Tavel JA (2010) What are Biomarkers? Curr Opin HIV AIDS 5: 463-466. doi: 10.1097/COH.0b013e32833ed177
    [2] Robinson DH, Toledo AH (2012) Historical development of modern anesthesia. J Invest Surg 25: 141-149. doi: 10.3109/08941939.2012.690328
    [3] Pitt D, Aubin JM (2012) Joseph Lister: father of modern surgery. Can J Surg 55: E8-9.
    [4] Hsu JL (2013) A brief history of vaccines: smallpox to the present. S D Med Spec no: 33-37.
    [5] Blevins SM, Bronze MS (2010) Robert Koch and the 'golden age' of bacteriology. Int J Infect Dis 14: e744-751. doi: 10.1016/j.ijid.2009.12.003
    [6] Biomarkers Definitions Working Group (2001) Biomarkers and surrogate endpoints: preferred definitions and conceptual framework. Clin Pharmacol Ther 69: 89-95. doi: 10.1067/mcp.2001.113989
    [7] O'Donnell MJ, Xavier D, Liu L, et al. (2010) Risk factors for ischaemic and intracerebral haemorrhagic stroke in 22 countries (the INTERSTROKE study): a case-control study. Lancet 376: 112-123. doi: 10.1016/S0140-6736(10)60834-3
    [8] Richens JL, Urbanowicz RA, Lunt EA, et al. (2009) Systems biology coupled with label-free high-throughput detection as a novel approach for diagnosis of chronic obstructive pulmonary disease. Respir Res 10: 29. doi: 10.1186/1465-9921-10-29
    [9] Song IU, Chung YA, Chung SW, et al. (2014) Early diagnosis of Alzheimer's disease and Parkinson's disease associated with dementia using cerebral perfusion SPECT. Dement Geriatr Cogn Disord 37: 276-285. doi: 10.1159/000357128
    [10] Perazella MA (2015) The Urine Sediment as a Biomarker of Kidney Disease. Am J Kidney Dis 66: 748-755. doi: 10.1053/j.ajkd.2015.02.342
    [11] Felton HT, Derrick JB, Swartz DP (1964) A Simple Immunologic Test For Pregnancy. Can Med Assoc J 91: 996-1000.
    [12] Slamon DJ, Clark GM, Wong SG, et al. (1987) Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science 235: 177-182. doi: 10.1126/science.3798106
    [13] Puglisi F, Fontanella C, Amoroso V, et al. (2015) Current challenges in HER2-positive breast cancer. Crit Rev Oncol Hematol 98: 211-221.
    [14] Anderson NL, Ptolemy AS, Rifai N (2013) The riddle of protein diagnostics: future bleak or bright? Clin Chem 59: 194-197. doi: 10.1373/clinchem.2012.184705
    [15] Mayeux R (2004) Biomarkers: potential uses and limitations. NeuroRx 1: 182-188. doi: 10.1602/neurorx.1.2.182
    [16] Anderson NL, Anderson NG (2002) The human plasma proteome: history, character, and diagnostic prospects. Mol Cell Proteomics 1: 845-867. doi: 10.1074/mcp.R200007-MCP200
    [17] Krebs HA (1950) Chemical composition of blood plasma and serum. Annu Rev Biochem 19: 409-430. doi: 10.1146/annurev.bi.19.070150.002205
    [18] Richens JL, Lunt EA, Sanger D, et al. (2009) Avoiding nonspecific interactions in studies of the plasma proteome: practical solutions to prevention of nonspecific interactions for label-free detection of low-abundance plasma proteins. J Proteome Res 8: 5103-5110. doi: 10.1021/pr900487y
    [19] Drucker E, Krapfenbauer K (2013) Pitfalls and limitations in translation from biomarker discovery to clinical utility in predictive and personalised medicine. EPMA J 4: 7. doi: 10.1186/1878-5085-4-7
    [20] Lim LS, Sherin K (2008) Screening for prostate cancer in U.S. men ACPM position statement on preventive practice. Am J Prev Med 34: 164-170. doi: 10.1016/j.amepre.2007.10.003
    [21] Prince M, Bryce R, Albanese E, et al. (2013) The global prevalence of dementia: a systematic review and metaanalysis. Alzheimers Dement 9: 63-75 e62. doi: 10.1016/j.jalz.2012.11.007
    [22] Galvin JE, Sadowsky CH (2012) Practical guidelines for the recognition and diagnosis of dementia. J Am Board Fam Med 25: 367-382. doi: 10.3122/jabfm.2012.03.100181
    [23] Bekris LM, Yu CE, Bird TD, et al. (2010) Genetics of Alzheimer disease. J Geriatr Psychiatry Neurol 23: 213-227. doi: 10.1177/0891988710383571
    [24] Wu L, Rosa-Neto P, Hsiung GY, et al. (2012) Early-onset familial Alzheimer's disease (EOFAD). Can J Neurol Sci 39: 436-445. doi: 10.1017/S0317167100013949
    [25] Ballard C, Gauthier S, Corbett A, et al. (2011) Alzheimer's disease. Lancet 377: 1019-1031. doi: 10.1016/S0140-6736(10)61349-9
    [26] Castello MA, Jeppson JD, Soriano S (2014) Moving beyond anti-amyloid therapy for the prevention and treatment of Alzheimer’s disease. BMC Neurol 14: 1-5. doi: 10.1186/1471-2377-14-1
    [27] Cedazo-Minguez A, Winblad B (2010) Biomarkers for Alzheimer's disease and other forms of dementia: clinical needs, limitations and future aspects. Exp Gerontol 45: 5-14. doi: 10.1016/j.exger.2009.09.008
    [28] Hardy J, Allsop D (1991) Amyloid deposition as the central event in the aetiology of Alzheimer's disease. Trends Pharmacol Sci 12: 383-388. doi: 10.1016/0165-6147(91)90609-V
    [29] Puzzo D, Gulisano W, Arancio O, et al. (2015) The keystone of Alzheimer pathogenesis might be sought in Abeta physiology. Neuroscience 307: 26-36. doi: 10.1016/j.neuroscience.2015.08.039
    [30] Fiandaca MS, Mapstone ME, Cheema AK, et al. (2014) The critical need for defining preclinical biomarkers in Alzheimer'sdisease. Alzheimers Dement 10: S196-S212. doi: 10.1016/j.jalz.2014.04.015
    [31] Racchi M, Govoni S (2003) The pharmacology of amyloid precursor protein processing. Exp Gerontol 38: 145-157. doi: 10.1016/S0531-5565(02)00158-4
    [32] Nhan HS, Chiang K, Koo EH (2015) The multifaceted nature of amyloid precursor protein and its proteolytic fragments: friends and foes. Acta Neuropathol 129: 1-19. doi: 10.1007/s00401-014-1347-2
    [33] Chasseigneaux S, Allinquant B (2012) Functions of Abeta, sAPPalpha and sAPPbeta : similarities and differences. J Neurochem 120 Suppl 1: 99-108. doi: 10.1111/j.1471-4159.2011.07584.x
    [34] Simic G, Babic Leko M, Wray S, et al. (2016) Tau Protein Hyperphosphorylation and Aggregation in Alzheimer's Disease and Other Tauopathies, and Possible Neuroprotective Strategies. Biomolecules 6.
    [35] Wang Y, Mandelkow E (2016) Tau in physiology and pathology. Nat Rev Neurosci 17: 22-35.
    [36] Wang JZ, Xia YY, Grundke-Iqbal I, et al. (2013) Abnormal hyperphosphorylation of tau: sites, regulation, and molecular mechanism of neurofibrillary degeneration. J Alzheimers Dis 33 Suppl 1: S123-139.
    [37] Robakis NK (2010) Are Abeta and its derivatives causative agents or innocent bystanders in AD? Neurodegener Dis 7: 32-37. doi: 10.1159/000266476
    [38] Castello MA, Soriano S (2014) On the origin of Alzheimer's disease. Trials and tribulations of the amyloid hypothesis. Ageing Res Rev 13: 10-12. doi: 10.1016/j.arr.2013.10.001
    [39] Loy CT, Schofield PR, Turner AM, et al. (2014) Genetics of dementia. Lancet 383: 828-840. doi: 10.1016/S0140-6736(13)60630-3
    [40] Arriagada PV, Growdon JH, Hedley-Whyte ET, et al. (1992) Neurofibrillary tangles but not senile plaques parallel duration and severity of Alzheimer's disease. Neurology 42: 631-639. doi: 10.1212/WNL.42.3.631
    [41] Crystal H, Dickson D, Fuld P, et al. (1988) Clinico-pathologic studies in dementia: nondemented subjects with pathologically confirmed Alzheimer's disease. Neurology 38: 1682-1687. doi: 10.1212/WNL.38.11.1682
    [42] Benilova I, Karran E, De Strooper B (2012) The toxic A[beta] oligomer and Alzheimer's disease: an emperor in need of clothes. Nat Neurosci 15: 349-357. doi: 10.1038/nn.3028
    [43] Hardy J (2009) The amyloid hypothesis for Alzheimer's disease: a critical reappraisal. J Neurochem 110: 1129-1134. doi: 10.1111/j.1471-4159.2009.06181.x
    [44] Hsia AY, Masliah E, McConlogue L, et al. (1999) Plaque-independent disruption of neural circuits in Alzheimer's disease mouse models. Proc Natl Acad Sci U S A 96: 3228-3233. doi: 10.1073/pnas.96.6.3228
    [45] Mucke L, Masliah E, Yu GQ, et al. (2000) High-level neuronal expression of abeta 1-42 in wild-type human amyloid protein precursor transgenic mice: synaptotoxicity without plaque formation. J Neurosci 20: 4050-4058.
    [46] Holmes C, Boche D, Wilkinson D, et al. (2008) Long-term effects of Abeta42 immunisation in Alzheimer's disease: follow-up of a randomised, placebo-controlled phase I trial. Lancet 372: 216-223. doi: 10.1016/S0140-6736(08)61075-2
    [47] Karran E, Mercken M, De Strooper B (2011) The amyloid cascade hypothesis for Alzheimer's disease: an appraisal for the development of therapeutics. Nat Rev Drug Discov 10: 698-712. doi: 10.1038/nrd3505
    [48] Robakis NK (2011) Mechanisms of AD neurodegeneration may be independent of Abeta and its derivatives. Neurobiol Aging 32: 372-379. doi: 10.1016/j.neurobiolaging.2010.05.022
    [49] Richens JL, Morgan K, O'Shea P (2014) Reverse engineering of Alzheimer's disease based on biomarker pathways analysis. Neurobiol Aging 35: 2029-2038. doi: 10.1016/j.neurobiolaging.2014.02.024
    [50] Morgan K (2011) The three new pathways leading to Alzheimer's disease. Neuropathol Appl Neurobiol 37: 353-357. doi: 10.1111/j.1365-2990.2011.01181.x
    [51] Baumgart M, Snyder HM, Carrillo MC, et al. (2015) Summary of the evidence on modifiable risk factors for cognitive decline and dementia: A population-based perspective. Alzheimers Dement 11: 718-726. doi: 10.1016/j.jalz.2015.05.016
    [52] Hebert LE, Scherr PA, Bienias JL, et al. (2003) Alzheimer disease in the US population: prevalence estimates using the 2000 census. Arch Neurol 60: 1119-1122. doi: 10.1001/archneur.60.8.1119
    [53] Moser VA, Pike CJ (2015) Obesity and sex interact in the regulation of Alzheimer's disease. Neurosci Biobehav Rev in press.
    [54] Li R, Singh M (2014) Sex differences in cognitive impairment and Alzheimer's disease. Front Neuroendocrinol 35: 385-403. doi: 10.1016/j.yfrne.2014.01.002
    [55] Barron AM, Pike CJ (2012) Sex hormones, aging, and Alzheimer's disease. Front Biosci (Elite Ed) 4: 976-997.
    [56] Luchsinger JA, Cheng D, Tang MX, et al. (2012) Central obesity in the elderly is related to late-onset Alzheimer disease. Alzheimer Dis Assoc Disord 26: 101-105. doi: 10.1097/WAD.0b013e318222f0d4
    [57] Emmerzaal TL, Kiliaan AJ, Gustafson DR (2015) 2003-2013: a decade of body mass index, Alzheimer's disease, and dementia. J Alzheimers Dis 43: 739-755.
    [58] Tsoi KK, Chan JY, Hirai HW, et al. (2015) Cognitive Tests to Detect Dementia: A Systematic Review and Meta-analysis. JAMA Intern Med 175: 1450-1458. doi: 10.1001/jamainternmed.2015.2152
    [59] Folstein MF, Folstein SE, McHugh PR (1975) "Mini-mental state". A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res 12: 189-198. doi: 10.1016/0022-3956(75)90026-6
    [60] Mioshi E, Dawson K, Mitchell J, et al. (2006) The Addenbrooke's Cognitive Examination Revised (ACE-R): a brief cognitive test battery for dementia screening. Int J Geriatr Psychiatry 21: 1078-1085. doi: 10.1002/gps.1610
    [61] Larner AJ, Mitchell AJ (2014) A meta-analysis of the accuracy of the Addenbrooke's Cognitive Examination (ACE) and the Addenbrooke's Cognitive Examination-Revised (ACE-R) in the detection of dementia. Int Psychogeriatr 26: 555-563. doi: 10.1017/S1041610213002329
    [62] Hampel H, Frank R, Broich K, et al. (2010) Biomarkers for Alzheimer's disease: academic, industry and regulatory perspectives. Nat Rev Drug Discov 9: 560-574. doi: 10.1038/nrd3115
    [63] Corder EH, Saunders AM, Strittmatter WJ, et al. (1993) Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science 261: 921-923. doi: 10.1126/science.8346443
    [64] Saunders AM, Strittmatter WJ, Schmechel D, et al. (1993) Association of apolipoprotein E allele epsilon 4 with late-onset familial and sporadic Alzheimer's disease. Neurology 43: 1467-1472. doi: 10.1212/WNL.43.8.1467
    [65] Xiang Y, Lam SM, Shui G (2015) What can lipidomics tell us about the pathogenesis of Alzheimer disease? Biol Chem 396: 1281-1291.
    [66] Bertram L, McQueen MB, Mullin K, et al. (2007) Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet 39: 17-23. doi: 10.1038/ng1934
    [67] Corder EH, Saunders AM, Risch NJ, et al. (1994) Protective effect of apolipoprotein E type 2 allele for late onset Alzheimer disease. Nat Genet 7: 180-184. doi: 10.1038/ng0694-180
    [68] Myers RH, Schaefer EJ, Wilson PW, et al. (1996) Apolipoprotein E epsilon4 association with dementia in a population-based study: The Framingham study. Neurology 46: 673-677. doi: 10.1212/WNL.46.3.673
    [69] Barber RC, Phillips NR, Tilson JL, et al. (2015) Can Genetic Analysis of Putative Blood Alzheimer's Disease Biomarkers Lead to Identification of Susceptibility Loci? PLoS One 10: e0142360. doi: 10.1371/journal.pone.0142360
    [70] Van Cauwenberghe C, Van Broeckhoven C, Sleegers K (2015) The genetic landscape of Alzheimer disease: clinical implications and perspectives. Genet Med in press.
    [71] Karch CM, Goate AM (2015) Alzheimer's disease risk genes and mechanisms of disease pathogenesis. Biol Psychiatry 77: 43-51. doi: 10.1016/j.biopsych.2014.05.006
    [72] Bertram L, Lill CM, Tanzi RE (2010) The genetics of Alzheimer disease: back to the future. Neuron 68: 270-281. doi: 10.1016/j.neuron.2010.10.013
    [73] Guerreiro R, Wojtas A, Bras J, et al. (2013) TREM2 variants in Alzheimer's disease. N Engl J Med 368: 117-127. doi: 10.1056/NEJMoa1211851
    [74] Harold D, Abraham R, Hollingworth P, et al. (2009) Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat Genet 41: 1088-1093. doi: 10.1038/ng.440
    [75] Hollingworth P, Harold D, Sims R, et al. (2011) Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease. Nat Genet 43: 429-435. doi: 10.1038/ng.803
    [76] Lambert JC, Heath S, Even G, et al. (2009) Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet 41: 1094-1099. doi: 10.1038/ng.439
    [77] Naj AC, Jun G, Beecham GW, et al. (2011) Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer's disease. Nat Genet 43: 436-441. doi: 10.1038/ng.801
    [78] Seshadri S, Fitzpatrick AL, Ikram MA, et al. (2010) Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA 303: 1832-1840. doi: 10.1001/jama.2010.574
    [79] Huang Y (2010) Abeta-independent roles of apolipoprotein E4 in the pathogenesis of Alzheimer's disease. Trends Mol Med 16: 287-294. doi: 10.1016/j.molmed.2010.04.004
    [80] Kim J, Yoon H, Basak J, et al. (2014) Apolipoprotein E in synaptic plasticity and Alzheimer's disease: potential cellular and molecular mechanisms. Mol Cells 37: 767-776. doi: 10.14348/molcells.2014.0248
    [81] DeMattos RB, Cirrito JR, Parsadanian M, et al. (2004) ApoE and clusterin cooperatively suppress Abeta levels and deposition: evidence that ApoE regulates extracellular Abeta metabolism in vivo. Neuron 41: 193-202. doi: 10.1016/S0896-6273(03)00850-X
    [82] Nathan BP, Bellosta S, Sanan DA, et al. (1994) Differential effects of apolipoproteins E3 and E4 on neuronal growth in vitro. Science 264: 850-852. doi: 10.1126/science.8171342
    [83] Ji Y, Gong Y, Gan W, et al. (2003) Apolipoprotein E isoform-specific regulation of dendritic spine morphology in apolipoprotein E transgenic mice and Alzheimer's disease patients. Neuroscience 122: 305-315. doi: 10.1016/j.neuroscience.2003.08.007
    [84] Lidstrom AM, Bogdanovic N, Hesse C, et al. (1998) Clusterin (apolipoprotein J) protein levels are increased in hippocampus and in frontal cortex in Alzheimer's disease. Exp Neurol 154: 511-521. doi: 10.1006/exnr.1998.6892
    [85] Sihlbom C, Davidsson P, Sjogren M, et al. (2008) Structural and quantitative comparison of cerebrospinal fluid glycoproteins in Alzheimer's disease patients and healthy individuals. Neurochem Res 33: 1332-1340. doi: 10.1007/s11064-008-9588-x
    [86] Nilselid AM, Davidsson P, Nagga K, et al. (2006) Clusterin in cerebrospinal fluid: analysis of carbohydrates and quantification of native and glycosylated forms. Neurochem Int 48: 718-728. doi: 10.1016/j.neuint.2005.12.005
    [87] Jongbloed W, van Dijk KD, Mulder SD, et al. (2015) Clusterin Levels in Plasma Predict Cognitive Decline and Progression to Alzheimer's Disease. J Alzheimers Dis 46: 1103-1110. doi: 10.3233/JAD-150036
    [88] Thambisetty M, Simmons A, Velayudhan L, et al. (2010) Association of plasma clusterin concentration with severity, pathology, and progression in Alzheimer disease. Arch Gen Psychiatry 67: 739-748. doi: 10.1001/archgenpsychiatry.2010.78
    [89] Schrijvers EM, Koudstaal PJ, Hofman A, et al. (2011) Plasma clusterin and the risk of Alzheimer disease. JAMA 305: 1322-1326. doi: 10.1001/jama.2011.381
    [90] McGeer PL, Kawamata T, Walker DG (1992) Distribution of clusterin in Alzheimer brain tissue. Brain Res 579: 337-341. doi: 10.1016/0006-8993(92)90071-G
    [91] Ghiso J, Matsubara E, Koudinov A, et al. (1993) The cerebrospinal-fluid soluble form of Alzheimer's amyloid beta is complexed to SP-40,40 (apolipoprotein J), an inhibitor of the complement membrane-attack complex. Biochem J 293 (Pt 1): 27-30.
    [92] Zlokovic BV (1996) Cerebrovascular transport of Alzheimer's amyloid beta and apolipoproteins J and E: possible anti-amyloidogenic role of the blood-brain barrier. Life Sci 59: 1483-1497. doi: 10.1016/0024-3205(96)00310-4
    [93] Giannakopoulos P, Kovari E, French LE, et al. (1998) Possible neuroprotective role of clusterin in Alzheimer's disease: a quantitative immunocytochemical study. Acta Neuropathol 95: 387-394. doi: 10.1007/s004010050815
    [94] Oda T, Wals P, Osterburg HH, et al. (1995) Clusterin (apoJ) alters the aggregation of amyloid beta-peptide (A beta 1-42) and forms slowly sedimenting A beta complexes that cause oxidative stress. Exp Neurol 136: 22-31. doi: 10.1006/exnr.1995.1080
    [95] Li X, Ma Y, Wei X, et al. (2014) Clusterin in Alzheimer's disease: a player in the biological behavior of amyloid-beta. Neurosci Bull 30: 162-168. doi: 10.1007/s12264-013-1391-2
    [96] Jiang T, Tan L, Zhu XC, et al. (2014) Upregulation of TREM2 ameliorates neuropathology and rescues spatial cognitive impairment in a transgenic mouse model of Alzheimer's disease. Neuropsychopharmacology 39: 2949-2962. doi: 10.1038/npp.2014.164
    [97] Zhao Y, Bhattacharjee S, Jones BM, et al. (2013) Regulation of TREM2 expression by an NF-small ka, CyrillicB-sensitive miRNA-34a. Neuroreport 24: 318-323. doi: 10.1097/WNR.0b013e32835fb6b0
    [98] Lue LF, Schmitz C, Walker DG (2015) What happens to microglial TREM2 in Alzheimer's disease: Immunoregulatory turned into immunopathogenic? Neuroscience 302: 138-150. doi: 10.1016/j.neuroscience.2014.09.050
    [99] Jiang T, Yu JT, Zhu XC, et al. (2013) TREM2 in Alzheimer's disease. Mol Neurobiol 48: 180-185. doi: 10.1007/s12035-013-8424-8
    [100] Frank S, Burbach GJ, Bonin M, et al. (2008) TREM2 is upregulated in amyloid plaque-associated microglia in aged APP23 transgenic mice. Glia 56: 1438-1447. doi: 10.1002/glia.20710
    [101] Paradowska-Gorycka A, Jurkowska M (2013) Structure, expression pattern and biological activity of molecular complex TREM-2/DAP12. Hum Immunol 74: 730-737. doi: 10.1016/j.humimm.2013.02.003
    [102] Paris D, Ait-Ghezala G, Bachmeier C, et al. (2014) The spleen tyrosine kinase (Syk) regulates Alzheimer amyloid-beta production and Tau hyperphosphorylation. J Biol Chem 289: 33927-33944. doi: 10.1074/jbc.M114.608091
    [103] Lebouvier T, Scales TM, Hanger DP, et al. (2008) The microtubule-associated protein tau is phosphorylated by Syk. Biochim Biophys Acta 1783: 188-192. doi: 10.1016/j.bbamcr.2007.11.005
    [104] Vafadar-Isfahani B, Ball G, Coveney C, et al. (2012) Identification of SPARC-like 1 Protein as Part of a Biomarker Panel for Alzheimer's Disease in Cerebrospinal Fluid. J Alzheimers Dis 28: 625-636.
    [105] Richens JL, Vere KA, Light RA, et al. (2014) Practical detection of a definitive biomarker panel for Alzheimer's disease; comparisons between matched plasma and cerebrospinal fluid. Int J Mol Epidemiol Genet 5: 53-70.
    [106] Mendis DB, Ivy GO, Brown IR (2000) Induction of SC1 mRNA encoding a brain extracellular matrix glycoprotein related to SPARC following lesioning of the adult rat forebrain. Neurochem Res 25: 1637-1644. doi: 10.1023/A:1026626805612
    [107] Yin GN, Lee HW, Cho JY, et al. (2009) Neuronal pentraxin receptor in cerebrospinal fluid as a potential biomarker for neurodegenerative diseases. Brain Res 1265: 158-170. doi: 10.1016/j.brainres.2009.01.058
    [108] Xu PT, Li YJ, Qin XJ, et al. (2006) Differences in apolipoprotein E3/3 and E4/4 allele-specific gene expression in hippocampus in Alzheimer disease. Neurobiol Dis 21: 256-275. doi: 10.1016/j.nbd.2005.07.004
    [109] Abdi F, Quinn JF, Jankovic J, et al. (2006) Detection of biomarkers with a multiplex quantitative proteomic platform in cerebrospinal fluid of patients with neurodegenerative disorders. J Alzheimers Dis 9: 293-348.
    [110] Weisel JW (2005) Fibrinogen and fibrin. Adv Protein Chem 70: 247-299. doi: 10.1016/S0065-3233(05)70008-5
    [111] Cortes-Canteli M, Strickland S (2009) Fibrinogen, a possible key player in Alzheimer's disease. J Thromb Haemost 7 Suppl 1: 146-150.
    [112] Cortes-Canteli M, Paul J, Norris EH, et al. (2010) Fibrinogen and beta-amyloid association alters thrombosis and fibrinolysis: a possible contributing factor to Alzheimer's disease. Neuron 66: 695-709. doi: 10.1016/j.neuron.2010.05.014
    [113] Ahn HJ, Glickman JF, Poon KL, et al. (2014) A novel Abeta-fibrinogen interaction inhibitor rescues altered thrombosis and cognitive decline in Alzheimer's disease mice. J Exp Med 211: 1049-1062. doi: 10.1084/jem.20131751
    [114] Ahn HJ, Zamolodchikov D, Cortes-Canteli M, et al. (2010) Alzheimer's disease peptide beta-amyloid interacts with fibrinogen and induces its oligomerization. Proc Natl Acad Sci U S A 107: 21812-21817. doi: 10.1073/pnas.1010373107
    [115] Mueller C, Zhou W, Vanmeter A, et al. (2010) The heme degradation pathway is a promising serum biomarker source for the early detection of Alzheimer's disease. J Alzheimers Dis 19: 1081-1091.
    [116] Langbein L, Heid HW, Moll I, et al. (1993) Molecular characterization of the body site-specific human epidermal cytokeratin 9: cDNA cloning, amino acid sequence, and tissue specificity of gene expression. Differentiation 55: 57-71. doi: 10.1111/j.1432-0436.1993.tb00033.x
    [117] Ghosh D, Lippert D, Krokhin O, et al. (2010) Defining the membrane proteome of NK cells. J Mass Spectrom 45: 1-25.
    [118] van Niel G, Raposo G, Candalh C, et al. (2001) Intestinal epithelial cells secrete exosome-like vesicles. Gastroenterology 121: 337-349. doi: 10.1053/gast.2001.26263
    [119] de Mateo S, Castillo J, Estanyol JM, et al. (2011) Proteomic characterization of the human sperm nucleus. Proteomics 11: 2714-2726. doi: 10.1002/pmic.201000799
    [120] Bohm D, Keller K, Pieter J, et al. (2012) Comparison of tear protein levels in breast cancer patients and healthy controls using a de novo proteomic approach. Oncol Rep 28: 429-438.
    [121] Kim YS, Gu BH, Choi BC, et al. (2013) Apolipoprotein A-IV as a novel gene associated with polycystic ovary syndrome. Int J Mol Med 31: 707-716.
    [122] Fu BS, Liu W, Zhang JW, et al. (2009) [Serum proteomic analysis on metastasis-associated proteins of hepatocellular carcinoma]. Nan Fang Yi Ke Da Xue Xue Bao 29: 1775-1778.
    [123] Li X, Long J, He T, et al. (2015) Integrated genomic approaches identify major pathways and upstream regulators in late onset Alzheimer's disease. Sci Rep 5: 12393. doi: 10.1038/srep12393
    [124] Femminella GD, Ferrara N, Rengo G (2015) The emerging role of microRNAs in Alzheimer's disease. Front Physiol 6: 40.
    [125] Lukiw WJ, Zhao Y, Cui JG (2008) An NF-kappaB-sensitive micro RNA-146a-mediated inflammatory circuit in Alzheimer disease and in stressed human brain cells. J Biol Chem 283: 31315-31322. doi: 10.1074/jbc.M805371200
    [126] Zhu HC, Wang LM, Wang M, et al. (2012) MicroRNA-195 downregulates Alzheimer's disease amyloid-beta production by targeting BACE1. Brain Res Bull 88: 596-601. doi: 10.1016/j.brainresbull.2012.05.018
    [127] Santa-Maria I, Alaniz ME, Renwick N, et al. (2015) Dysregulation of microRNA-219 promotes neurodegeneration through post-transcriptional regulation of tau. J Clin Invest 125: 681-686. doi: 10.1172/JCI78421
    [128] Kosicek M, Zetterberg H, Andreasen N, et al. (2012) Elevated cerebrospinal fluid sphingomyelin levels in prodromal Alzheimer's disease. Neurosci Lett 516: 302-305. doi: 10.1016/j.neulet.2012.04.019
    [129] Mulder C, Wahlund LO, Teerlink T, et al. (2003) Decreased lysophosphatidylcholine/phosphatidylcholine ratio in cerebrospinal fluid in Alzheimer's disease. J Neural Transm (Vienna) 110: 949-955. doi: 10.1007/s00702-003-0007-9
    [130] Satoi H, Tomimoto H, Ohtani R, et al. (2005) Astroglial expression of ceramide in Alzheimer's disease brains: a role during neuronal apoptosis. Neuroscience 130: 657-666. doi: 10.1016/j.neuroscience.2004.08.056
    [131] Mapstone M, Cheema AK, Fiandaca MS, et al. (2014) Plasma phospholipids identify antecedent memory impairment in older adults. Nat Med 20: 415-418. doi: 10.1038/nm.3466
    [132] Wood PL, Medicherla S, Sheikh N, et al. (2015) Targeted Lipidomics of Fontal Cortex and Plasma Diacylglycerols (DAG) in Mild Cognitive Impairment and Alzheimer's Disease: Validation of DAG Accumulation Early in the Pathophysiology of Alzheimer's Disease. J Alzheimers Dis 48: 537-546. doi: 10.3233/JAD-150336
    [133] Hu Z, Chang YC, Wang Y, et al. (2013) VisANT 4.0: Integrative network platform to connect genes, drugs, diseases and therapies. Nucleic Acids Res 41: W225-231. doi: 10.1093/nar/gkt401
    [134] Chatr-Aryamontri A, Breitkreutz BJ, Oughtred R, et al. (2015) The BioGRID interaction database: 2015 update. Nucleic Acids Res 43: D470-478. doi: 10.1093/nar/gku1204
    [135] Johnson C, Tinti M, Wood NT, et al. (2011) Visualization and biochemical analyses of the emerging mammalian 14-3-3-phosphoproteome. Mol Cell Proteomics 10: M110 005751.
    [136] Reddy A, Huang CC, Liu H, et al. (2010) Robust gene network analysis reveals alteration of the STAT5a network as a hallmark of prostate cancer. Genome Inform 24: 139-153.
    [137] Araujo DJ, Anderson AG, Berto S, et al. (2015) FoxP1 orchestration of ASD-relevant signaling pathways in the striatum. Genes Dev 29: 2081-2096. doi: 10.1101/gad.267989.115
    [138] Haas LT, Salazar SV, Kostylev MA, et al. (2015) Metabotropic glutamate receptor 5 couples cellular prion protein to intracellular signalling in Alzheimer's disease. Brain 139: 526-546.
    [139] Shahani N, Seshadri S, Jaaro-Peled H, et al. (2015) DISC1 regulates trafficking and processing of APP and Abeta generation. Mol Psychiatry 20: 874-879. doi: 10.1038/mp.2014.100
    [140] Croft D, O'Kelly G, Wu G, et al. (2011) Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res 39: D691-697. doi: 10.1093/nar/gkq1018
    [141] D'Eustachio P (2013) Pathway Databases: Making Chemical and Biological Sense of the Genomic Data Flood. Chem Biol 20: 629-635. doi: 10.1016/j.chembiol.2013.03.018
    [142] Wang X, Liotta L (2011) Clinical bioinformatics: a new emerging science. J Clin Bioinforma 1: 1. doi: 10.1186/2043-9113-1-1
    [143] Bellazzi R, Masseroli M, Murphy S, et al. (2012) Clinical Bioinformatics: challenges and opportunities. BMC Bioinform 13: 1-8. doi: 10.1186/1471-2105-13-1
    [144] Butte AJ (2008) Translational Bioinformatics: Coming of Age. J Am Med Inf Assoc 15: 709-714. doi: 10.1197/jamia.M2824
    [145] Larrañaga P, Calvo B, Santana R, et al. (2006) Machine learning in bioinformatics. Brief Bioinform 7: 86-112. doi: 10.1093/bib/bbk007
    [146] Pavlopoulos GA, Secrier M, Moschopoulos CN, et al. (2011) Using graph theory to analyze biological networks. BioData Min 4: 1-27. doi: 10.1186/1756-0381-4-1
    [147] Bullmore E, Sporns O (2009) Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci 10: 186-198. doi: 10.1038/nrn2575
    [148] Stam CJ (2014) Modern network science of neurological disorders. Nat Rev Neurosci 15: 683-695. doi: 10.1038/nrn3801
    [149] Baumgartner C, Osl M, Netzer M, et al. (2011) Bioinformatic-driven search for metabolic biomarkers in disease. J Clin Bioinforma 1: 1-10. doi: 10.1186/2043-9113-1-1
    [150] Inza I, Calvo B, Armañanzas R, et al. (2010) Machine Learning: An Indispensable Tool in Bioinformatics. In: Matthiesen R, editor. Bioinformatics Methods in Clinical Research: Humana Press. pp. 25-48.
    [151] Guyon I, Weston J, Barnhill S, et al. Gene Selection for Cancer Classification using Support Vector Machines. Mach Learn 46: 389-422.
    [152] Hsu W-C, Denq C, Chen S-S (2013) A diagnostic methodology for Alzheimer’s disease. J Clin Bioinforma 3: 1-13.
    [153] Horwitz B, Rowe JB (2011) Functional biomarkers for neurodegenerative disorders based on the network paradigm. Prog Neurobiol 95: 505-509. doi: 10.1016/j.pneurobio.2011.07.005
    [154] de Haan W, Pijnenburg YAL, Strijers RLM, et al. (2009) Functional neural network analysis in frontotemporal dementia and Alzheimer's disease using EEG and graph theory. BMC Neurosci 10: 101-101. doi: 10.1186/1471-2202-10-101
    [155] Brier MR, Thomas JB, Fagan AM, et al. (2014) Functional connectivity and graph theory in preclinical Alzheimer's disease. Neurobiol Aging 35: 757-768. doi: 10.1016/j.neurobiolaging.2013.10.081
    [156] Zhang D, Wang Y, Zhou L, et al. (2011) Multimodal classification of Alzheimer's disease and mild cognitive impairment. NeuroImage 55: 856-867. doi: 10.1016/j.neuroimage.2011.01.008
  • Reader Comments
  • © 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7428) PDF downloads(1495) Cited by(2)

Article outline

Figures and Tables

Figures(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog