Gender | MMSE | Age | |
AD | 10/15 | 20.48 ± 10.22 | 75.99 ± 4.28 |
MCI | 90/115 | 28.13 ± 1.734 | 71.48 ± 7.57 |
HC | 31/34 | 28.94 ± 1.25 | 75.17 ± 5.86 |
IL-10-induced DC (DC10) secrete IL-10, thereby promoting aeroallergen tolerance by activation of CD25+Foxp3+ Treg, while retinoic acid-induced DC (DC-RA) foster food allergen tolerance via IL-27-dependent induction of Foxp3− Treg. In some respects, these outcomes reflect those seen with aero- or food allergen-presenting lung and intestinal DC, respectively. Herein we asked whether these DCreg would also be functional in their reciprocal settings. DC-RA expressed lower levels of CCR5, CCR9, and CD103 than DC10, but higher levels of CD40, CD86, MHC II, TGF-β, IL-27 and Aldh1A2. DC-RA were also more effective in suppressing OVA-specific T cell proliferation in vitro (p ≤ 0.05). Co-culture of DC10 or DC-RA with OVA-specific T cells activated a regulatory phenotype therein, with the DC-10-induced Treg being modestly more suppressive in vitro than the DC-RA-induced Treg. We previously reported on lung trafficking of DC10, but DC-RA also traffic through the lungs and mediastinal lymph nodes before accumulating in the mesenteric lymph nodes. Both DCreg populations equally reduced airway hyperresponsiveness, Th2 responses (Th2 cytokines, eosinophilia) to airway allergen challenge, and plasma IgE/IgG1 levels in an OVA-asthma model. Similarly, they were equally effective in our OVA-food allergy model, where they reduced clinical scores (p ≤ 0.001), mast cell activation (p ≤ 0.05) and Th2 cytokine responses to allergen challenge, as well as OVA-specific IgE/IgG1 levels (p ≤ 0.01). Taken together, this data indicates that while DC10 and DC-RA employ distinct operative mechanisms, they were both able to induce tolerance in lung- and gut-associated allergic disease.
Citation: Chunyan Li, Wojciech Dawicki, Xiaobei Zhang, Chris Rudulier, John R. Gordon. IL-10- and retinoic acid-induced regulatory dendritic cells are therapeutically equivalent in mouse models of asthma and food allergy[J]. AIMS Allergy and Immunology, 2021, 5(2): 73-91. doi: 10.3934/Allergy.2021007
[1] | Wei Kong, Feifan Xu, Shuaiqun Wang, Kai Wei, Gen Wen, Yaling Yu . Application of orthogonal sparse joint non-negative matrix factorization based on connectivity in Alzheimer's disease research. Mathematical Biosciences and Engineering, 2023, 20(6): 9923-9947. doi: 10.3934/mbe.2023435 |
[2] | Zhi Yang, Kang Li, Haitao Gan, Zhongwei Huang, Ming Shi, Ran Zhou . An Alzheimer's Disease classification network based on MRI utilizing diffusion maps for multi-scale feature fusion in graph convolution. Mathematical Biosciences and Engineering, 2024, 21(1): 1554-1572. doi: 10.3934/mbe.2024067 |
[3] | Fang Niu, Zongwei Liu, Peidong Liu, Hongrui Pan, Jiaxue Bi, Peng Li, Guangze Luo, Yonghui Chen, Xiaoxing Zhang, Xiangchen Dai . Identification of novel genetic biomarkers and treatment targets for arteriosclerosis-related abdominal aortic aneurysm using bioinformatic tools. Mathematical Biosciences and Engineering, 2021, 18(6): 9761-9774. doi: 10.3934/mbe.2021478 |
[4] | Yin Zhang, Xintong Wu, Jingwen Sun, Kecen Yue, Shuangshuang Lu, Bingjian Wang, Wenjia Liu, Haifeng Shi, Ling Zou . Exploring changes in brain function in IBD patients using SPCCA: a study of simultaneous EEG-fMRI. Mathematical Biosciences and Engineering, 2024, 21(2): 2646-2670. doi: 10.3934/mbe.2024117 |
[5] | Xinlin Liu, Viktor Krylov, Su Jun, Natalya Volkova, Anatoliy Sachenko, Galina Shcherbakova, Jacek Woloszyn . Segmentation and identification of spectral and statistical textures for computer medical diagnostics in dermatology. Mathematical Biosciences and Engineering, 2022, 19(7): 6923-6939. doi: 10.3934/mbe.2022326 |
[6] | Boyang Wang, Wenyu Zhang . MARnet: multi-scale adaptive residual neural network for chest X-ray images recognition of lung diseases. Mathematical Biosciences and Engineering, 2022, 19(1): 331-350. doi: 10.3934/mbe.2022017 |
[7] | Atsushi Kawaguchi . Network-based diagnostic probability estimation from resting-state functional magnetic resonance imaging. Mathematical Biosciences and Engineering, 2023, 20(10): 17702-17725. doi: 10.3934/mbe.2023787 |
[8] | Marios G. Krokidis, Themis P. Exarchos, Panagiotis Vlamos . Data-driven biomarker analysis using computational omics approaches to assess neurodegenerative disease progression. Mathematical Biosciences and Engineering, 2021, 18(2): 1813-1832. doi: 10.3934/mbe.2021094 |
[9] | Shi Liu, Kaiyang Li, Yaoying Wang, Tianyou Zhu, Jiwei Li, Zhenyu Chen . Knowledge graph embedding by fusing multimodal content via cross-modal learning. Mathematical Biosciences and Engineering, 2023, 20(8): 14180-14200. doi: 10.3934/mbe.2023634 |
[10] | Guanghua Fu, Qingjuan Wei, Yongsheng Yang . Bearing fault diagnosis with parallel CNN and LSTM. Mathematical Biosciences and Engineering, 2024, 21(2): 2385-2406. doi: 10.3934/mbe.2024105 |
IL-10-induced DC (DC10) secrete IL-10, thereby promoting aeroallergen tolerance by activation of CD25+Foxp3+ Treg, while retinoic acid-induced DC (DC-RA) foster food allergen tolerance via IL-27-dependent induction of Foxp3− Treg. In some respects, these outcomes reflect those seen with aero- or food allergen-presenting lung and intestinal DC, respectively. Herein we asked whether these DCreg would also be functional in their reciprocal settings. DC-RA expressed lower levels of CCR5, CCR9, and CD103 than DC10, but higher levels of CD40, CD86, MHC II, TGF-β, IL-27 and Aldh1A2. DC-RA were also more effective in suppressing OVA-specific T cell proliferation in vitro (p ≤ 0.05). Co-culture of DC10 or DC-RA with OVA-specific T cells activated a regulatory phenotype therein, with the DC-10-induced Treg being modestly more suppressive in vitro than the DC-RA-induced Treg. We previously reported on lung trafficking of DC10, but DC-RA also traffic through the lungs and mediastinal lymph nodes before accumulating in the mesenteric lymph nodes. Both DCreg populations equally reduced airway hyperresponsiveness, Th2 responses (Th2 cytokines, eosinophilia) to airway allergen challenge, and plasma IgE/IgG1 levels in an OVA-asthma model. Similarly, they were equally effective in our OVA-food allergy model, where they reduced clinical scores (p ≤ 0.001), mast cell activation (p ≤ 0.05) and Th2 cytokine responses to allergen challenge, as well as OVA-specific IgE/IgG1 levels (p ≤ 0.01). Taken together, this data indicates that while DC10 and DC-RA employ distinct operative mechanisms, they were both able to induce tolerance in lung- and gut-associated allergic disease.
retinaldehyde 1A2;
carboxyfluoryl succinylate ester;
CC subfamily chemokine receptor 9;
IL-10-skewed tolerogenic dendritic cells;
LPS-matured immunostimulatory dendritic cell;
retinoic acid-skewed dendritic cells;
fetal bovine serum;
the transcription factor forkhead box P3;
lymphocyte activation gene 3;
lipopolysaccaride;
mouse mast cell protease-1;
ovalbumin;
transgenic mice expressing an OVA-specific T cell receptor;
regulatory T cell
Alzheimer's disease (AD) is a neurodegenerative disease caused by many factors, and its incidence is increasing yearly [1]. The change of gene expression in genetics involves gene variation. However, the brain structure and function of AD patients are also different from those of the control group. Image genetics explores changes in brain structure and function from the perspective of genetic variation. Through a correlation analysis of imaging and genetics, we can explore the potential markers of AD.
Machine learning algorithms have been widely used in various bioinformatics fields, such as miRNA-disease relationship prediction. Ha and Park [2] proposed the metric learning for predicting miRNA-disease association (MLMD) algorithm, which can reveal not only novel miRNAs associated with disease, but also miRNA-miRNA and disease-disease similarities. Moreover, they proposed the matrix factorization with disease similarity constraint (MDMF) algorithm based on matrix factorization, which incorporates disease similarity constraints to improve the prediction performance [3]. In addition, they proposed the simple yet effective computational framework (SMAP) algorithm to predict the relationship between and accurately predict the association between miRNA-diseases, which combines miRNA functional similarity, disease semantic similarity and Gaussian interaction spectrum kernel similarity [4]. Recently, they combined a deep neural network to propose the node2vec-based neural collaborative filtering for predicting miRNA-disease association (NCMD) algorithm, which uses Node2vec to understand the low-dimensional vector representation of miRNA and disease, and combines the linear ability of generalized matrix factorization and the nonlinear ability of multi-layer perceptron [5]. They tested and confirmed the effectiveness of the algorithm on three datasets of breast cancer, lung cancer and pancreatic cancer. In the search for biomarkers of AD, many scholars have proposed algorithms related to the diagnosis and prediction of AD. Park et al. [6] proposed a deep learning model to integrate large-scale gene expression and DNA methylation data to predict AD. This method is superior to traditional machine learning algorithms in that it uses typical dimensionality reduction methods and improves the accuracy of prediction. Wang et al. created the multi-task sparse canonical correlation analysis and regression (MT-SCCAR) model, combined with the annual total score of depression level, the clinical dementia rating scale, the functional activity questionnaire, and the neuropsychiatric symptom questionnaire. These four clinical data were used as compensation information and embedded in the algorithm by a linear regression. They confirmed the superiority and robustness of the algorithm on real and simulated data [7].
Canonical correlation analysis (CCA) is an algorithm to obtain the maximum correlation between two kinds of data. However, it is not suitable for an association analysis of high-dimensional data. For this reason, some scholars put forward a sparse canonical correlation analysis (SCCA) algorithm [8]. Based on CCA, the SCCA algorithm assists the CCA algorithm in feature selection in high-dimensional features through l1 norm constraints. However, because the l1 norm constraints only considers sparsity at the individual level, it is only partially applicable to image data. Lesions in different brain regions may play a role at the same time; therefore, it is necessary to add structural constraints to the SCCA algorithm. Du et al. [9] proposed a graph-guided pairwise group lasso (GGL) and applied it to image data. GGL can be used in a data-driven mode that does not provide prior knowledge. It thinks that each group consists of only two variables, and they will be extracted with similar or equal weights. They found that the performance of this algorithm with actual data is due to other competitive algorithms. However, they did not consider the diagnostic information of the subjects. Previous studies showed that the addition of diagnostic information could effectively improve the correlation analysis performance of the algorithm [10,11]. In addition, there may be feature redundancy in imaging and genetic data, and orthogonal constraints can be added to the algorithm by linear programming.
Thus, it was suggested to integrate structural magnetic resonance imaging (sMRI) and the gene expression data of AD patients and its control group using an orthogonal structured SCCA algorithm based on diagnostic information fusion. Specifically, after preprocessing sMRI data, we extracted the gray matter volume of each region of interest (ROI) as a feature. Then, we picked the expression of differentially expressed genes (DEGs) between the sick group and the control group as characteristics from the gene expression data. Based on the SCCA algorithm, this method added GGL constraints on images and orthogonal constraints on two kinds of data, which could improve the performance of the association analysis and prevent the influence of feature redundancy on the results. The experimental results showed that this algorithm was superior to other CCA-based algorithms and had a stronger correlation analysis ability. Top ROIs and top genes with biological and diagnostic significance can be obtained. The selected top biomarkers can provide a reference for the diagnosis of AD and drug target discovery.
The SCCA algorithm adds sparse constraints to the CCA algorithm. Given n samples, p ROIs and q genes, sMRI data can be expressed as X∈Rn×p, and gene expression data can be expressed as Y∈Rn×q. The objective of the SCCA algorithm is to adjust the typical correlation weights u∈Rp×l and v∈Rq×l to maximize the correlation between Xu and Yv, and its objective function is shown in Formula (1):
minu,v −uTXTYv+λ1||u||1+λ2||v||1 |
s.t.∥Xu∥2≤1,∥Yv∥2≤1 | (1) |
where λ1 and λ2 control the sparsity of u and v, respectively.
The graph-guided fusion lasso differs from the conventional group lasso in that it does not rely on prior knowledge; however, the graph-guided fusion lasso will introduce estimation bias. GGL uses group lasso and graph-guided fusion lasso. It can be defined as follows:
ΩGGL(u)=∑(i,j)∈E √u2i+u2j | (2) |
where E is the edge set of the graph which highly related features are connected.
This paper presented an OSSCCA-DIF algorithm, which uses GGL as a structural constraint based on the SCCA algorithm. Orthogonal constraints were used to prevent the redundant characteristics of image genetics data from affecting the results. In addition, we added the diagnostic information of samples as prior knowledge to the algorithm to improve its correlation performance. The objective function of the OSSCCA-DIF algorithm is given as follows:
minu,v −uTXTYv+λ1||u||1+λ2||v||1+XTz+β1∥uuT−I∥22+β2∥vvT−I∥22+μΩGGL(u) |
s.t.∥Xu∥2≤1,∥Yv∥2≤1 | (3) |
where z∈Rn×1 is used to store the diagnostic information of the sample, β1 and β2 control the orthogonal constraint strength of u and v, respectively, and γ is applied to control the strength of the GGL constraint.
For the optimization of Formula (3), the Lagrange multiplier method can be used to solve the partial derivatives of the weight u of the ROI and the weight v of the gene, respectively. Firstly, u is regarded as a constant term. Then, the objective function can be rewritten as Eq (4):
L(u,v)=−uTXTYv+β1∥uuT−I∥22+XTz+β2∥vvT−I∥22+∥Xu∥2+∥Yv∥2+μΩGGL(u)+λ1||u||1+λ2||v||1+XTz. | (4) |
For Eq (4), take the derivative of v and make it zero, and we can get Eq (5):
VTYv-YTXu+2β2(vvT−I)v+λ2Dvv+2YTYv = 0. | (5) |
The iterative solution formula of v can be written as Formula (6):
v = (YTY+λ2Dv+XTz+2YTY+2β2(vvT−I)v)-1(YTXu). | (6) |
Similarly, for u, if v is regarded as a constant term, then the solution formula of u can be obtained, as shown in Eq (7):
u = (XTX+λ1Du+2XTX+μGu+2β1(uuT−I)u)-1(XTYv). | (7) |
In Eq (7), Du and Dv are diagonal matrices, and their diagonal elements in the i th row can be expressed as 1|u1i|(i=1,⋯,p) and 1|v1i|(i=1,⋯,q), respectively. Gu is also a diagonal matrix, and the diagonal elements in the i th row can be expressed as 1√u21(i−1)+u21i+1√u21i+u21(i+1)(i=2,⋯,p−1).
We obtained 296 samples of AD, mild cognitive impairment (MCI), and healthy controls (HC) from the Alzheimer's Disease Neuroimaging Disease (ADNI) database (https://adni.loni.usc.edu/). Table 1 provides the statistical information for each set of samples. We collected sMRI and gene expression data from these samples. For sMRI, we first calibrated the head movement using the DiffusionKit software and then segmented the images using the Statistical Parametric Mapping (SPM) software package of the Computational Anatomy Toolbox (CAT) toolkit of Matlab software and divided them into gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF). In this paper, we used a GM volume of 140 brain regions divided by the brain as the ROI. For the gene expression data, we utilized the limma algorithm to examine the differential expression and obtained 962 DEGs. We extracted the expression of DEGs as the genetic feature of the sample.
Gender | MMSE | Age | |
AD | 10/15 | 20.48 ± 10.22 | 75.99 ± 4.28 |
MCI | 90/115 | 28.13 ± 1.734 | 71.48 ± 7.57 |
HC | 31/34 | 28.94 ± 1.25 | 75.17 ± 5.86 |
In order to obtain the best result, this paper took the canonical correlation coefficient (CCC) as the performance measure of the algorithm, and the solution method of the CCC was as follows. Using the grid search method [0.01 0.1 1], the parameters of super parameters (β1, β2, λ1, λ2, and μ) of the OSSCCA-DIF algorithm were selected for the real data sets (Figure 1). Finally, the results of the sixth parameter selection were taken as the best parameters (β1=0.01, β2=0.01, λ1=0.01, λ2=0.1 and μ=1). We detail the optimal parameter information in Table S1 of the Supplementary material. It can be seen from the figure that the CCC obtained by different parameter combinations is quite different. The proposed algorithm is sensitive to the parameters and thus affects the stability of the model. Therefore, it is helpful to select the optimal results by appropriately enlarging the selection range of parameters in the practical application. Figure 2 shows the heat map of weights u and v. Specifically, the abscissa of the graph represents either an ROI or a gene. The distribution of different colors represents the weight of the ROIs or genes; the darker the color, the higher the weight. After taking the absolute values of U and V, we give the names and weight information of the top 10 ROIs and the top 10 genes with the greatest weight in Table 2. We will discuss the biological significance of these ROIs and genes in detail in the discussion section. Figure 3 is a visual display of the top 10 ROI. Figure 4 displays the enrichment analysis results of the top 10 genes. We will discuss the relationship between these channels and AD in detail in the discussion section.
ROI | Weights | Gene | Weights |
rEnt | 0.06482 | ADAM17 | 0.006038 |
lAmy | 0.000438 | TBX4 | 0.00569 |
rAmy | 1.98 × 10−06 | ZC2HC1C | 0.005581 |
lFroOpe | 2.94 × 10−07 | CCND1 | 0.005188 |
lMidOccGy | 8.01 × 10−08 | POM121L12 | 0.004765 |
rInfTemGy | 3.23 × 10−08 | HIST1H2BM | 0.004729 |
rParHipGy | 3.10 × 10−08 | ALDH3A1 | 0.004531 |
lHip | 2.67 × 10−08 | ANO5 | 0.004486 |
rTemPo | 2.38 × 10−08 | SPDYE4 | 0.004477 |
rFusGy | 2.34 × 10−08 | HERC2P9 | 0.00443 |
To prove that the proposed algorithm has a good correlation analysis ability, we compared the CCC of the proposed algorithm with the SCCA-FGL, SCCA, and CCA algorithms (Table 3), and the CCC of the proposed algorithm was higher than the other three algorithms. In addition, we present the Pearson correlation heat map of the top 10 ROIs and top 10 genes in Figure 5. Among them, lHip and ZC2HC1C reached the maximum positive correlation (0.4551), while lMidOccGy and POM121L12 reached the maximum negative correlation (−0.3227).
Algorithm | CCC |
OSSCCA-DIF | 0.3033 |
SCCA_FGL | 0.2357 |
SCCA | 0.1563 |
CCA | 0.0427 |
Note: rEnt = Right Entorhinal Area, lAmy = Left Amygdala, rAmy = Right Amygdala, lFroOpe = Left Frontal Operculum, lMidOccGy = Left Middle Occipital Gyrus, rInfTemGy = Right Inferior Temporal Gyrus, rParHipGy = Right Parahippocampus Gyrus, lHip = Left Hippocampus, rTemPo = Right Temporal Pole, rFusGy = Right Fusiform Gyrus.
In this section, we used the Receiver Operating Characteristic (ROC) curve to verify the diagnostic performance of the top markers (Figure 6). ROC curves have been widely used in the biomedical field [12,13]. Based on the logistic regression algorithm, we constructed the diagnosis model using the top 10 ROIs and top 10 genes. Among them, the AUC of the diagnosis model constructed by the ROIs reached 0.898. The AUC of the diagnosis model created by the genes reached 0.853. In addition, we present the diagnostic models constructed using the top 10 ROIs and the top 10 genes from several other algorithms in Figure S1 of the Supplementary material. The AUC of the top 10 ROIs selected by the proposed algorithm was higher than that of several other algorithms. The AUC of the top 10 genes selected was slightly lower than that of the CCA algorithm. In addition, we present the details of the ROC curves of the diagnostic models constructed by the top ROIs/genes selected by the OSSCCA-DIF, SCCA-FGL, SCCA, and CCA algorithms in Table S2 and Table S3 of the Supplementary material.
In order to further verify the effectiveness of the algorithm, we constructed two synthetic data sets (X1∈Rn1×p1, Y1∈Rn1×q1, X2∈Rn2×p2, X2∈Rn2×q2), and generated two sets of weights (u1∈Rp1×1,v1∈Rq1×1,u2∈Rp2×1, v2∈Rq2×1). Here, n1=100, n2=500, p1 = 300, q1=500, and p2=800. In addition, we generated the variable z∈Rn×1. Xk and Yk can be generated by (xij)k N(ziukj,l∗σk). Here, k=1,2. σk and l denote the variance and noise level of the noise, respectively. We set l to 10 and present the CCCS of several algorithms on two synthetic datasets in Table 4. As can be seen from the table, the CCC of the proposed algorithm is larger than that of the other algorithms on both datasets.
Algorithm | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
OSSCCA-DIF | 0.2733 | 0.2320 |
SCCA_FGL | 0.0372 | 0.1119 |
SCCA | 0.1510 | 0.0836 |
CCA | 0.1692 | 0.0424 |
Additionally, we conducted ablation experiments on the proposed OSSCCA-DIF algorithm. Specifically, by removing each part of the OSSCCA-DIF algorithm either individually or in pairs (except for sparse constraints), we compare the CCC between the real and simulation datasets under the same parameter conditions (Table 5). The objective functions to be compared (scenarios 1–6) are given below, as shown in Eqs (8)–(13):
L(u,v)=−uTXTYv+∥Xu∥2+∥Yv∥2+μΩGGL(u)+λ1||u||1+λ2||v||1, | (8) |
L(u,v)=−uTXTYv+β1∥uuT−I∥22+β2∥vvT−I∥22+∥Xu∥2+∥Yv∥2+λ1||u||1+λ2||v||1, | (9) |
L(u,v)=−uTXTYv+∥Xu∥2+∥Yv∥2+λ1||u||1+λ2||v||1+XTz, | (10) |
L(u,v)=−uTXTYv+β1∥uuT−I∥22+β2∥vvT−I∥22+∥Xu∥2+∥Yv∥2+μΩGGL(u)+λ1||u||1+λ2||v||1, | (11) |
L(u,v)=−uTXTYv+β1∥uuT−I∥22+β2∥vvT−I∥22+∥Xu∥2+∥Yv∥2+λ1||u||1+λ2||v||1+XTz, | (12) |
L(u,v)=−uTXTYv+μΩGGL(u)+λ1||u||1+λ2||v||1+XTz+∥Xu∥2+∥Yv∥2. | (13) |
Case | CCC (real dataset 1) | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
Case 1 | 0.1427 | 0.1236 | 0.1198 |
Case 2 | 0.1542 | 0.1468 | 0.1333 |
Case 3 | 0.1642 | 0.1539 | 0.1248 |
Case 4 | 0.1344 | 0.1326 | 0.1164 |
Case 5 | 0.1529 | 0.1432 | 0.1409 |
Case 6 | 0.1243 | 0.1142 | 0.1022 |
OSCCA-DIF | 0.3033 | 0.2733 | 0.2320 |
As a neurodegenerative disease, there is no effective treatment for AD. Alongside aging, the incidence of AD is increasing year by year. Image genetics can mine disease-related markers by integrating image genomics and genetic data through a series of correlation analysis algorithms. Therefore, this paper proposed an OSSCCA-DIF algorithm. Based on the SCCA algorithm, this algorithm added orthogonal constraints on weight vectors u and v and GGL structural constraints on image feature weight u. In addition, the diagnostic information of the sample was added to the algorithm. The experimental results showed that by integrating sMRI and ROI data in real data, the performance of this algorithm was better than other CCA-based algorithms.
Most of the top ROIs mined by the proposed algorithm have proven to be closely related to AD. First, our algorithm determined that rEnt was the brain region with the most significant weight, and the weight value reached 0.06482. The entorhinal cortex (EC) is unanimously considered to be the earliest pathological structure of AD [14,15]. Thaker et al. [16] analyzed the relationship between the thickness of the olfactory cortex (sMRI) and pathological changes (autopsy) in 50 AD patients and found that the thickness of the olfactory cortex may be related to the severity of AD. The experiment confirmed that, compared with the control group, the volume and average thickness of the right inner olfactory cortex in AD patients were lower, and the expression level of lncRNA BACE1-AS in plasma exosomes isolated from AD patients was significantly increased. Therefore, Wang et al. [17] proposed that the level of BACE1-AS in peripheral blood exosomes should be combined with the volume and thickness of the right entorhinal cortex as a potential biomarker of AD. Second, our algorithm identified that lAmy and rAmy were top ROIs. These two brain regions are both sides of the amygdala. As an important structure of emotional learning and memory, the amygdala is related to a series of mental diseases such as AD. The MRI volume of the amygdala may be related to the severity of dementia in AD, and it shows neuron loss and atrophy in AD patients [18,19,20]. The amygdala has an excellent diagnostic value for sMRI of AD [21]. Finally, our algorithm found that the top brain regions (rParHipGy and lHip) also proved to be closely related to AD. In the memory system of the human brain, it is essential to connect the posterior cingulate cortex with the hippocampus, either directly or indirectly through the parahippocampal gyrus. These brain regions all play a vital role in the early progression of AD [22]. In the experiment evaluating the correlation between brain metabolism and the orientation of AD, it was found that its improved orientation performance was related to the more significant brain metabolism in brain regions such as the left middle occipital gyrus, and the higher CERAD identification score was more related to the metabolic activity in the left medial temporal lobe regions (including the hippocampus, parahippocampal gyrus, and left fusiform gyrus) [23].
The top genes (ADAM17, CCND1, HIST1H2BM, ALDH3A1) determined by our algorithm are proven to be either directly or indirectly related to AD. It is common knowledge that a feature of AD is the accumulation of extracellular Amyloid-β (Aβ) plaques and neurofibrillary tangles composed of tau in neurons [24,25,26]. Among them, Aβ is a protein hydrolysate separated from its precursor, namely the amyloid precursor protein (APP), by β- and γ-secretase, and tau is a microtubule-associated protein involved in microtubule stability. The main manifestations of AD patients are decreased memory, attention, spatial orientation, language ability, and olfactory function, which are all related to the deposition of tau protein and APP. It has been proven that ADAM17 is a potential therapeutic target for AD because ADAM17 can be used as an α-secretase regulating APP, thereby affecting the production of Aβ. Additionally, protease encoded by ADAM17 plays a role in the shedding of tumor necrosis factor-α (TNF-α). As a key pro-inflammatory cytokine in inflammation, TNF-α's signal transduction aggravates Aβ and tau pathology in vivo [27]. In order to explore the role of propionic acid (PPA) in the pathogenesis of AD, Aliashrafi et al. [28] selected 284 genes related to PPA and AD and identified CCND1 as an important hub gene, bottleneck gene, and seed gene through a network analysis and an Molecular Complex Detection (MCODE) analysis. Zeng et al. [29] also determined that CCND1 can be the core goal of AD treatment. H2BC14 (HIST1H2BM) is the core component of the nucleosome. The nucleosome assembly protein 1-like 5(NAP1L5) is downregulated in the brain tissue of AD patients, and the overexpression of NAP1L5 can alleviate APP metabolism and Tau phosphorylation [30]. ALDH3A1, a protein-coding gene, belongs to A1, a member of the aldehyde dehydrogenase 3 family. In the enrichment analysis of the top 10 genes, we also found that aldehyde dehydrogenase is closely related to AD. The relationship between other genes and AD needs further study.
Through the enrichment analysis of the top 10 genes, we found that many significant pathways were related to the occurrence and development of AD. McKibben and Rhoades [31] studied the role of the proline-rich region (PRR) in regulating the interaction between Tau and soluble tubulin. Aldehyde dehydrogenase can balance the amine metabolism of neurodegenerative diseases such as AD [32]. Tao et al. [33] indicated that aldehyde dehydrogenase-2 could be a potential target for AD treatment. Rapamycin can block G1/S conversion between AD patients and normal controls. Experiments have confirmed that compared with the control group, AD patients who used rapamycin still progressed to the late cell cycle [34]. Cyclin dependent kinase 5 (Cdk5) can also be used as a potential therapeutic target for AD [35].
Finally, we also determined the diagnostic significance of the top markers for AD. We verified the AUC of each top ROI and gene in the test set and found that the AUCs of all ROIs and genes were more significant than 0.5. In addition, through the logistic regression algorithm, we found that the AUC of the diagnosis model constructed by the top 10 ROIs and top 10 genes reached 0.898 and 0.853, respectively, and were both within reasonable confidence intervals.
In this paper, we discussed the risk brain regions and risk genes closely related to the occurrence and development of AD by studying the relationship between sMRI and bivariate variables of gene expression data. The proposed OSSCCA-DIF algorithm has advantages in correlation analysis performance and biomarker selection. However, most algorithms based on the SCCA algorithm assume that the image genetic data is linear; however, this assumption is not necessarily valid in real data. Therefore, in future research, we will try to introduce a deep-learning algorithm to make up for this defect.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors declare there is no conflict of interest.
The data used in this paper came the AD Neuroimaging Initiative (ADNI) database (https://adni.loni.usc.edu).
Hyperparameter | Value |
β1 | 0.01 |
β2 | 0.01 |
λ1 | 0.01 |
λ2 | 0.1 |
μ | 1 |
Random seed | 1 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.898 | 0.815–0.912 | 0.86 | 0.812 | 0.673 |
SCCA-FGL | 0.85 | 0.741–0.960 | 0.791 | 0.875 | 0.666 |
SCCA | 0.712 | 0.562–0.863 | 0.605 | 0.812 | 0.417 |
CCA | 0.89 | 0.807–0.972 | 0.767 | 0.875 | 0.642 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.853 | 0.729–0.978 | 0.791 | 0.875 | 0.666 |
SCCA-FGL | 0.763 | 0.635–0.891 | 0.581 | 0.875 | 0.456 |
SCCA | 0.760 | 0.628–0.893 | 0.907 | 0.5 | 0.407 |
CCA | 0.865 | 0.770–0.959 | 0.791 | 0.812 | 0.603 |
[1] |
Gordon JR, Ma Y, Churchman L, et al. (2014) Regulatory dendritic cells for immunotherapy in immunologic diseases. Front Immunol 5: 7. doi: 10.3389/fimmu.2014.00007
![]() |
[2] |
Weiner HL, da Cunha AP, Quintana F, et al. (2011) Oral tolerance. Immunol Rev 241: 241-259. doi: 10.1111/j.1600-065X.2011.01017.x
![]() |
[3] |
Dawicki W, Li C, Town J, et al. (2017) Therapeutic reversal of food allergen sensitivity by mature retinoic acid-differentiated dendritic cell induction of LAG3+CD49b−Foxp3− regulatory T cells. J Allergy Clin Immunol 139: 1608-1620. doi: 10.1016/j.jaci.2016.07.042
![]() |
[4] |
Huang H, Dawicki W, Lu M, et al. (2013) Regulatory dendritic cell expression of MHCII and IL-10 are jointly requisite for induction of tolerance in a murine model of OVA-asthma. Allergy 68: 1126-1135. doi: 10.1111/all.12273
![]() |
[5] |
Koya T, Matsuda H, Takeda K, et al. (2007) IL-10-treated dendritic cells decrease airway hyperresponsiveness and airway inflammation in mice. J Allergy Clin Immunol 119: 1241-1250. doi: 10.1016/j.jaci.2007.01.039
![]() |
[6] |
Lu M, Dawicki W, Zhang X, et al. (2011) Therapeutic induction of tolerance by IL-10-differentiated dendritic cells in a mouse model of house dust mite-asthma. Allergy 66: 612-620. doi: 10.1111/j.1398-9995.2010.02526.x
![]() |
[7] |
Nayyar A, Dawicki W, Huang H, et al. (2012) Induction of prolonged asthma tolerance by IL-10-differentiated dendritic cells: differential impact on airway hyperresponsiveness and the Th2 immunoinflammatory response. J Immunol 189: 72-79. doi: 10.4049/jimmunol.1103286
![]() |
[8] |
Akbari O, DeKruyff RH, Umetsu DT (2001) Pulmonary dendritic cells producing IL-10 mediate tolerance induced by respiratory exposure to antigen. Nat Immunol 2: 725-731. doi: 10.1038/90667
![]() |
[9] |
Huang H, Dawicki W, Zhang X, et al. (2010) Tolerogenic dendritic cells induce CD4+ CD25hiFoxp3+ regulatory T cell differentiation from CD4+CD25−/loFoxp3− effector T cells. J Immunol 185: 5003-5010. doi: 10.4049/jimmunol.0903446
![]() |
[10] |
Bakdash G, Vogelpoel LTC, Van Capel TM, et al. (2015) Retinoic acid primes human dendritic cells to induce gut-homing, IL-10-producing regulatory T cells. Mucosal Immunol 8: 265-278. doi: 10.1038/mi.2014.64
![]() |
[11] |
Lutz MB, Kukutsch N, Ogilvie AL, et al. (1999) An advanced culture method for generating large quantities of highly pure dendritic cells from mouse bone marrow. J Immunol Methods 223: 77-92. doi: 10.1016/S0022-1759(98)00204-X
![]() |
[12] | Nayyar A (2009) Therapeutic immunomodulation of allergic lung disease using regulatory dendritic cells in a mouse model of asthma [PhD's thesis] Canada: University of Saskatchewan. |
[13] |
Gordon JR, Li F, Nayyar A, et al. (2005) CD8 alpha+, but not CD8 alpha−, dendritic cells tolerize Th2 responses via contact-dependent and -independent mechanisms, and reverse airway hyperresponsiveness, Th2, and eosinophil responses in a mouse model of asthma. J Immunol 175: 1516-1522. doi: 10.4049/jimmunol.175.3.1516
![]() |
[14] |
Schneider AM, Li F, Zhang X, et al. (2001) Induction of pulmonary allergen-specific IgA responses or airway hyperresponsiveness in the absence of allergic lung disease following sensitization with limiting doses of ovalbumin-alum. Cell Immunol 212: 101-109. doi: 10.1006/cimm.2001.1854
![]() |
[15] |
Glaab T, Daser A, Braun A, et al. (2001) Tidal midexpiratory flow as a measure of airway hyperresponsiveness in allergic mice. Am J Physiol Lung Cell Mol Physiol 280: 565-573. doi: 10.1152/ajplung.2001.280.3.L565
![]() |
[16] |
Vijayaraghavan R, Schaper M, Thompson R, et al. (1994) Computer assisted recognition and quantitation of the effects of airborne chemicals acting at different areas of the respiratory tract in mice. Arch Toxicol 68: 490-499. doi: 10.1007/s002040050101
![]() |
[17] |
Li X, Yang A, Huang H, et al. (2010) Induction of Type 2 T helper cell allergen tolerance by IL-10-differentiated regulatory dendritic cells. Am J Resp Cell Molec Biol 42: 190-199. doi: 10.1165/rcmb.2009-0023OC
![]() |
[18] |
Campbell I (2007) Chi-squared and Fisher–Irwin tests of two-by-two tables with small sample recommendations. Stat Med 26: 3661-3675. doi: 10.1002/sim.2832
![]() |
[19] |
Huang H, Ma Y, Dawicki W, et al. (2013) Comparison of induced versus natural regulatory T cells of the same TCR specificity for induction of tolerance to an environmental antigen. J Immunol 191: 1136-1143. doi: 10.4049/jimmunol.1201899
![]() |
[20] |
Dai Z, Zhang S, Xie Q, et al. (2014) Natural CD8+CD122+ T cells are more potent in suppression of allograft rejection than CD4+CD25+ regulatory T cells. Am J Transplant 14: 39-48. doi: 10.1111/ajt.12515
![]() |
[21] | Yoffey JM, Courtice FC (1970) Lymphatics, Lymph and the Lymphomyeloid Complex London: Academic Press. |
[22] |
Ruane D, Brane L, Reis BS, et al. (2013) Lung dendritic cells induce migration of protective T cells to the gastrointestinal tract. J Exp Med 210: 1871-1888. doi: 10.1084/jem.20122762
![]() |
[23] |
Lopez-Pacheco C, Soldevila G, Du Pont G, et al. (2016) CCR9 is a key regulator of early phases of allergic airway inflammation. Mediators Inflamm 2016: 3635809. doi: 10.1155/2016/3635809
![]() |
[24] |
Yurchenko E, Shio MT, Huang TC, et al. (2012) Inflammation-driven reprogramming of CD4+Foxp3+ regulatory T cells into pathogenic Th1/Th17 T effectors is abrogated by mTOR inhibition in vivo. PLoS One 7: e35572. doi: 10.1371/journal.pone.0035572
![]() |
[25] | Ness S, Lin S, Gordon JR (2021) Regulatory dendritic cells, T cell tolerance and dendritic cell therapy for immunologic disease. Front Immunol 12: 633436. |
[26] |
Giannoukakis N, Phillips B, Finegold D, et al. (2011) Phase I (safety) study of autologous tolerogenic dendritic cells in type 1 diabetic patients. Diabetes Care 34: 2026-2032. doi: 10.2337/dc11-0472
![]() |
[27] |
Bell GM, Anderson AE, Diboll J, et al. (2017) Autologous tolerogenic dendritic cells for rheumatoid and inflammatory arthritis. Ann Rheum Dis 76: 227-234. doi: 10.1136/annrheumdis-2015-208456
![]() |
[28] |
Jauregui-Amezaga A, Cabezón R, Ramírez-Morros A, et al. (2015) Intraperitoneal administration of autologous tolerogenic dendritic cells for refractory Crohn's disease: a phase I study. J Crohns Colitis 9: 1071-1078. doi: 10.1093/ecco-jcc/jjv144
![]() |
[29] |
Zubizarreta I, Florez-Grau G, Vila G, et al. (2019) Immune tolerance in multiple sclerosis and neuromyelitis optica with peptide-loaded tolerogenic dendritic cells in a phase 1b trial. Proc Natl Acad Sci USA 116: 8463-8470. doi: 10.1073/pnas.1820039116
![]() |
[30] | Macedo C, Tran LM, Zahorchak AF, et al. (2020) Donor-derived regulatory dendritic cell infusion results in host cell cross-dressing and T cell subset changes in prospective living donor liver transplant recipients. Am J Transplant In press. |
[31] |
Thomson AW, Ezzelarab MB (2020) Generation and functional assessment of nonhuman primate regulatory dendritic cells and their therapeutic efficacy in renal transplantation. Cell Immunol 351: 104087. doi: 10.1016/j.cellimm.2020.104087
![]() |
1. | Uddalak Mitra, Shafiq Ul Rehman, ML-Powered Handwriting Analysis for Early Detection of Alzheimer’s Disease, 2024, 12, 2169-3536, 69031, 10.1109/ACCESS.2024.3401104 | |
2. | Shafiq Ul Rehman, Uddalak Mitra, Cognitive Handwriting Insights for Alzheimer’s Diagnosis: A Hybrid Framework, 2025, 16, 2078-2489, 249, 10.3390/info16030249 | |
3. | Simone Treccarichi, Mirella Vinci, Antonino Musumeci, Rosanna Galati Rando, Carla Papa, Salvatore Saccone, Concetta Federico, Pinella Failla, Martino Ruggieri, Francesco Calì, Agata Polizzi, Andrea Praticò, Investigating the Role of the Zinc Finger Protein ZC2HC1C on Autism Spectrum Disorder Susceptibility, 2025, 61, 1648-9144, 574, 10.3390/medicina61040574 |
Gender | MMSE | Age | |
AD | 10/15 | 20.48 ± 10.22 | 75.99 ± 4.28 |
MCI | 90/115 | 28.13 ± 1.734 | 71.48 ± 7.57 |
HC | 31/34 | 28.94 ± 1.25 | 75.17 ± 5.86 |
ROI | Weights | Gene | Weights |
rEnt | 0.06482 | ADAM17 | 0.006038 |
lAmy | 0.000438 | TBX4 | 0.00569 |
rAmy | 1.98 × 10−06 | ZC2HC1C | 0.005581 |
lFroOpe | 2.94 × 10−07 | CCND1 | 0.005188 |
lMidOccGy | 8.01 × 10−08 | POM121L12 | 0.004765 |
rInfTemGy | 3.23 × 10−08 | HIST1H2BM | 0.004729 |
rParHipGy | 3.10 × 10−08 | ALDH3A1 | 0.004531 |
lHip | 2.67 × 10−08 | ANO5 | 0.004486 |
rTemPo | 2.38 × 10−08 | SPDYE4 | 0.004477 |
rFusGy | 2.34 × 10−08 | HERC2P9 | 0.00443 |
Algorithm | CCC |
OSSCCA-DIF | 0.3033 |
SCCA_FGL | 0.2357 |
SCCA | 0.1563 |
CCA | 0.0427 |
Algorithm | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
OSSCCA-DIF | 0.2733 | 0.2320 |
SCCA_FGL | 0.0372 | 0.1119 |
SCCA | 0.1510 | 0.0836 |
CCA | 0.1692 | 0.0424 |
Case | CCC (real dataset 1) | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
Case 1 | 0.1427 | 0.1236 | 0.1198 |
Case 2 | 0.1542 | 0.1468 | 0.1333 |
Case 3 | 0.1642 | 0.1539 | 0.1248 |
Case 4 | 0.1344 | 0.1326 | 0.1164 |
Case 5 | 0.1529 | 0.1432 | 0.1409 |
Case 6 | 0.1243 | 0.1142 | 0.1022 |
OSCCA-DIF | 0.3033 | 0.2733 | 0.2320 |
Hyperparameter | Value |
β1 | 0.01 |
β2 | 0.01 |
λ1 | 0.01 |
λ2 | 0.1 |
μ | 1 |
Random seed | 1 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.898 | 0.815–0.912 | 0.86 | 0.812 | 0.673 |
SCCA-FGL | 0.85 | 0.741–0.960 | 0.791 | 0.875 | 0.666 |
SCCA | 0.712 | 0.562–0.863 | 0.605 | 0.812 | 0.417 |
CCA | 0.89 | 0.807–0.972 | 0.767 | 0.875 | 0.642 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.853 | 0.729–0.978 | 0.791 | 0.875 | 0.666 |
SCCA-FGL | 0.763 | 0.635–0.891 | 0.581 | 0.875 | 0.456 |
SCCA | 0.760 | 0.628–0.893 | 0.907 | 0.5 | 0.407 |
CCA | 0.865 | 0.770–0.959 | 0.791 | 0.812 | 0.603 |
Gender | MMSE | Age | |
AD | 10/15 | 20.48 ± 10.22 | 75.99 ± 4.28 |
MCI | 90/115 | 28.13 ± 1.734 | 71.48 ± 7.57 |
HC | 31/34 | 28.94 ± 1.25 | 75.17 ± 5.86 |
ROI | Weights | Gene | Weights |
rEnt | 0.06482 | ADAM17 | 0.006038 |
lAmy | 0.000438 | TBX4 | 0.00569 |
rAmy | 1.98 × 10−06 | ZC2HC1C | 0.005581 |
lFroOpe | 2.94 × 10−07 | CCND1 | 0.005188 |
lMidOccGy | 8.01 × 10−08 | POM121L12 | 0.004765 |
rInfTemGy | 3.23 × 10−08 | HIST1H2BM | 0.004729 |
rParHipGy | 3.10 × 10−08 | ALDH3A1 | 0.004531 |
lHip | 2.67 × 10−08 | ANO5 | 0.004486 |
rTemPo | 2.38 × 10−08 | SPDYE4 | 0.004477 |
rFusGy | 2.34 × 10−08 | HERC2P9 | 0.00443 |
Algorithm | CCC |
OSSCCA-DIF | 0.3033 |
SCCA_FGL | 0.2357 |
SCCA | 0.1563 |
CCA | 0.0427 |
Algorithm | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
OSSCCA-DIF | 0.2733 | 0.2320 |
SCCA_FGL | 0.0372 | 0.1119 |
SCCA | 0.1510 | 0.0836 |
CCA | 0.1692 | 0.0424 |
Case | CCC (real dataset 1) | CCC (synthetic dataset 1) | CCC (synthetic dataset 2) |
Case 1 | 0.1427 | 0.1236 | 0.1198 |
Case 2 | 0.1542 | 0.1468 | 0.1333 |
Case 3 | 0.1642 | 0.1539 | 0.1248 |
Case 4 | 0.1344 | 0.1326 | 0.1164 |
Case 5 | 0.1529 | 0.1432 | 0.1409 |
Case 6 | 0.1243 | 0.1142 | 0.1022 |
OSCCA-DIF | 0.3033 | 0.2733 | 0.2320 |
Hyperparameter | Value |
β1 | 0.01 |
β2 | 0.01 |
λ1 | 0.01 |
λ2 | 0.1 |
μ | 1 |
Random seed | 1 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.898 | 0.815–0.912 | 0.86 | 0.812 | 0.673 |
SCCA-FGL | 0.85 | 0.741–0.960 | 0.791 | 0.875 | 0.666 |
SCCA | 0.712 | 0.562–0.863 | 0.605 | 0.812 | 0.417 |
CCA | 0.89 | 0.807–0.972 | 0.767 | 0.875 | 0.642 |
Algorithm | Area under curve | Confidence interval | Sensitivity | Specificity | Youden index |
OSSCA-DIF | 0.853 | 0.729–0.978 | 0.791 | 0.875 | 0.666 |
SCCA-FGL | 0.763 | 0.635–0.891 | 0.581 | 0.875 | 0.456 |
SCCA | 0.760 | 0.628–0.893 | 0.907 | 0.5 | 0.407 |
CCA | 0.865 | 0.770–0.959 | 0.791 | 0.812 | 0.603 |