Citation: Daniil Shevyrev, Valeriy Tereshchenko, Olesya Manova, Vladimir kozlov. Homeostatic proliferation as a physiological process and a risk factor for autoimmune pathology[J]. AIMS Allergy and Immunology, 2021, 5(1): 18-32. doi: 10.3934/Allergy.2021002
| [1] | Abdulwahab Ali Almazroi . Survival prediction among heart patients using machine learning techniques. Mathematical Biosciences and Engineering, 2022, 19(1): 134-145. doi: 10.3934/mbe.2022007 |
| [2] | Weidong Gao, Yibin Xu, Shengshu Li, Yujun Fu, Dongyang Zheng, Yingjia She . Obstructive sleep apnea syndrome detection based on ballistocardiogram via machine learning approach. Mathematical Biosciences and Engineering, 2019, 16(5): 5672-5686. doi: 10.3934/mbe.2019282 |
| [3] | Lili Jiang, Sirong Chen, Yuanhui Wu, Da Zhou, Lihua Duan . Prediction of coronary heart disease in gout patients using machine learning models. Mathematical Biosciences and Engineering, 2023, 20(3): 4574-4591. doi: 10.3934/mbe.2023212 |
| [4] | Jingchi Jiang, Xuehui Yu, Yi Lin, Yi Guan . PercolationDF: A percolation-based medical diagnosis framework. Mathematical Biosciences and Engineering, 2022, 19(6): 5832-5849. doi: 10.3934/mbe.2022273 |
| [5] | Abigail Ferreira, Rui Lapa, Nuno Vale . A retrospective study comparing creatinine clearance estimation using different equations on a population-based cohort. Mathematical Biosciences and Engineering, 2021, 18(5): 5680-5691. doi: 10.3934/mbe.2021287 |
| [6] | 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 |
| [7] | Michele L. Joyner, Cammey Cole Manning, Whitney Forbes, Valerie Bobola, William Frazier . Modeling Ertapenem: the impact of body mass index on distribution of the antibiotic in the body. Mathematical Biosciences and Engineering, 2019, 16(2): 713-726. doi: 10.3934/mbe.2019034 |
| [8] | Wenzhu Song, Lixia Qiu, Jianbo Qing, Wenqiang Zhi, Zhijian Zha, Xueli Hu, Zhiqi Qin, Hao Gong, Yafeng Li . Using Bayesian network model with MMHC algorithm to detect risk factors for stroke. Mathematical Biosciences and Engineering, 2022, 19(12): 13660-13674. doi: 10.3934/mbe.2022637 |
| [9] | Anastasia-Maria Leventi-Peetz, Kai Weber . Probabilistic machine learning for breast cancer classification. Mathematical Biosciences and Engineering, 2023, 20(1): 624-655. doi: 10.3934/mbe.2023029 |
| [10] | Tingxi Wen, Hanxiao Wu, Yu Du, Chuanbo Huang . Faster R-CNN with improved anchor box for cell recognition. Mathematical Biosciences and Engineering, 2020, 17(6): 7772-7786. doi: 10.3934/mbe.2020395 |
Obstructive sleep apnea (OSA) is featured by repeated upper airway narrowing or obstruction, causing a substantial reduction or cessation of airflow during sleep [1,2]. These disruptions in breathing lead to intermittent hypercapnia and hypoxemia, causing metabolic dysregulation, endothelial dysfunction, and systemic inflammatory responses [2,3]. These pathophysiologic changes have contributed to various major comorbidities [2,3,4,5]. OSA is a global health concern and economic burden, with estimates of the incidence of moderate or severe OSA reaching 6%-17% in the general adult population and up to 49% in the elderly population [4,6], with most of the cases being undiagnosed [7].
No single physical examination finding can be used to diagnose or exclusively account for OSA. The gold standard diagnostic test for OSA relies on in-lab all-night polysomnography (PSG) [7,8]. However, this method is costly, and patients require continuous sleep monitoring by medical staff [2]; thus, limited availability of PSG may cause delays in diagnosing OSA and an increased disease burden [7,9].
Several body measurements predispose patients to OSA. Body weight, body mass index (BMI), central body fat distribution, and large neck and waist circumferences are possible risk factors [2,6,10]. In addition, increased age and males are also at risk for OSA [2,6]. Despite the epidemiological reports demonstrating that age, sex, and BMI are associated with OSA severity, an integrated predicting model that uses only these three factors to assess OSA risk has not yet been developed.
Machine learning approaches have expanded the ability to design efficient predictive and diagnostic tools and have gained popularity in many areas of medical research and clinical applications [11,12]. Logistic regression (LR) is one of the computational methods for probing multivariate data to detect the independent effects of a variable after adjusting for the contributions of the others [13]. An artificial neural network (ANN) is a supervised learning algorithm with advantages in analyzing nonlinear data, enhancing data interpretation, and designing effective predictive and diagnostic models to assist different clinical tasks [11,14].
The present study investigated the clinical predictive validity of age, sex, and BMI for sleep apnea and applied LR and ANN algorithms to develop a risk-predicting model using these three parameters as predictors.
This retrospective observational study used clinical data from the medical records of participants who had undergone PSG at the Sleep Center of Shuang Ho Hospital, Taipei Medical University, Taiwan, between January 2015 and January 2020.
The inclusion criteria for this study were: (1) participants aged≥18 years and (2) who had undergone a full-night PSG. Participants were excluded if (1) they had received any OSA treatments or been admitted to the hospital for any invasive surgery for OSA, (2) they had undergone PSG for continuous positive airway pressure (CPAP) titration, (3) the total recording time was < 4 h [8], and (4) the total sleep time was < 3 h [15].
Before the PSG examination, a baseline screening questionnaire was used to assess each participant's basic information, medication, and surgical history. The participants were measured for height, weight, and BMI (kg/m2) at the time of registration.
All the participants underwent a full-night PSG recorded by qualified technicians. All the sleep signals were obtained using an Embla N7000 Recording System (Medcare Flaga, Reykjavik, Iceland) and Sandman Elite sleep diagnostic software (Natus Neurology Inc., Middleton, WI, USA).
Experienced sleep technologists scored all sleep apnea and hypopnea events following the American Academy of Sleep Medicine 2012 updated guideline [8,16]. Each participant's apnea-hypopnea index (AHI) was calculated according to the total number of apnea and hypopnea events divided by the total sleeping time. Hypopnea was defined as a 30-89% reduction in nasal airflow for at least 10s, with arousal or at least a decrease of 3% oxygen saturation. Apnea was defined as a more than 90% reduction in oral airflow for at least 10s with or without arousal or decreased saturation.
The severity of OSA was classified as none (AHI < 5), mild (AHI = 5-15), moderate (AHI = 15-30), and severe (AHI ≥ 30) [2,8]. To develop the predictive models for the risk of moderate-to-severe OSA, we divided the participants into none-and-mild (AHI < 15) and moderate-to-severe (AHI ≥ 15) OSA groups for further analysis.
Taipei Medical University-Joint Institutional Review Board (TMU-JIRB: N201911007) approved the protocol of this study. Waivers of informed consent were approved by the TMU-JIRB (TMU-JIRB: N201911007) for this retrospective study involving the secondary analysis of existing data.
All analyses were performed using STATISTICA software (TIBCO Software, Tulsa, OK, USA). Continuous data are presented as mean and standard deviation, whereas categorical variables are expressed as frequencies and percentages. We used Fisher's exact test to determine associations between two categorical variables and Pearson's chi-squared (X2) test to determine the statistical significance between the means of two continuous variables in the participants with none-to-mild and moderate-to-severe OSA. The one-way analysis of variance (ANOVA) was used to determine the statistical significance between the means of ≥ 3 independent groups. A two-tailed p-value less than 0.05 was considered statistically significant.
All the LR and ANN models were established and modeled using STATISTICA software. In the models, age and BMI were set as continuous variables, and sex was designated as a categorical variable. For our ANNs, the employed architecture was feed-forward, combined with a back-propagation algorithm with multilayer perceptron (MLP). A bias term was used for each neural layer in the MLP. We set the number of neurons in the hidden layer from 1 to 50. We initially partitioned the eligible PSG data into a training-testing set and an independent validation set, containing ~80% and ~20% of the data, respectively. The data from the training-testing set were used to determine the generalizability of the LR and ANN analyses using five-fold cross-validation (Figure 1A). All participants were stratified by age, sex, BMI, and OSA severity to maintain consistent proportions of demographic data across the training-testing and validation sets. We reported the mean AUCs for the training and testing sets of the LR and ANN models through five-fold cross-validation. The model with the best AUC in the five testing sets was identified to have the best predictive performance and used as the final model. The final models were then applied to the independent validation set, and we calculated the accuracy, sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and AUCs of the validation set derived from the confusion matrix [17] (Figure 1A). The classification threshold was set to 0.5 for all models.
During the study period, 11, 100 participants underwent PSG in the Sleep Center. Of these, 1, 665 participants were excluded from the baseline analyses for the following reasons: they had been admitted to the hospital for invasive surgery or CPAP titration (n = 247), the total recording time < 4h (n = 1, 386), the total sleep time < 3h (n = 32), missing BMI information (n = 7) or missing AHI information (n = 6) (Figure 1A). Overall, 9, 422 participants (2918 women and 6504 men, mean age: 49.2 ± 14.0 years) who had undergone PSG were enrolled in this study. Of these, 5, 920 (62.8%) had an AHI of ≥15, indicating moderate-to-severe OSA, and 3, 502 (37.2%) had an AHI < 15, indicating none-to-mild OSA.
A higher proportion of men have moderate-to-severe OSA than women (Table 1). Compared with participants with none-to-mild OSA, participants with moderate-to-severe OSA were older and had a higher BMI. In addition, participants with moderate-to-severe OSA had a shorter total sleep time, lower sleep efficiency, shorter rapid eye movement duration, and lower oxygen saturation during the PSG recording than those with none-to-mild OSA (Table 1).
| Variables | Moderate-to-severe OSA (AHI ≥ 15) |
None-to-mild OSA (AHI < 15) |
p-value |
| n | 5, 920 | 3, 502 | |
| Age | 50.8 ± 13.5 | 46.6 ± 14.4 | < 0.0001 |
| Female, n (%) | 1266 (21.4) | 1652 (47.2) | < 0.0001 |
| Height (cm) | 167.2 ± 8.5 | 164.4 ± 8.9 | < 0.0001 |
| Weight (kg) | 79.7±16.1 | 66.1 ± 13.1 | < 0.0001 |
| BMI (kg/m2) | 28.4 ± 4.8 | 24.3 ± 3.8 | < 0.0001 |
| Polysomnography | |||
| Total recording time (min) | 366.6±11.9 | 367.8±13.4 | < 0.0001 |
| Total sleep time (min) | 275.8±60.2 | 281.4±61.1 | < 0.0001 |
| Sleep efficiency (%) | 75.2±16.4 | 76.5±16.5 | 0.0002 |
| Non-REM sleep (min) | 240.6±51.7 | 242.2±51.3 | 0.133 |
| REM sleep (min) | 35.2±22.5 | 39.1±23.7 | < 0.0001 |
| AHI | 44.7±23.8 | 6.4±4.5 | < 0.0001 |
| SO2 mean (%) | 94.5±2.6 | 96.6±1.3 | < 0.0001 |
| Data are expressed as mean ± standard deviation. The p-value represents the comparison between the AHI ≥ 15 and AHI < 15 groups. | |||
DownLoad:
CSV
A multivariate analysis was applied to adjust for age, sex, and BMI to determine the independent predictors of moderate-to-severe OSA. Table 2 presents the variables' crude and adjusted odds ratio (OR) values in the LR model. After adjustment for the associated variables, the regression analysis revealed that the female sex was associated with a 0.28-fold reduced risk of moderate-to-severe OSA (95% CI: 0.25-0.32; p < 0.0001), and every unit increase in age or BMI was associated with a 1.04-fold (95%CI: 1.04-1.05; p < 0.0001) or 1.27-fold (95%CI: 1.26-1.29; p < 0.0001) increase in the likelihood of moderate-to-severe OSA, respectively.
| Unadjusted Model | Adjusted Model | |||
| Variables | OR | 95% CI | OR | 95% CI |
| Age | 1.02*† | 1.02-1.03 | 1.04*† | 1.04-1.05 |
| Female | 0.30* | 0.28-0.33 | 0.28* | 0.25-0.32 |
| BMI | 1.27*† | 1.25-1.28 | 1.27*† | 1.26-1.29 |
| The odds ratio (OR) and 95% confidence intervals (CI) represent the association between the variables and the presence of moderate-to-severe OSA (AHI ≥ 15). *p < 0.0001; †Odds ratio of per unit changes. | ||||
DownLoad:
CSV
As mentioned, age, sex, and BMI were used as input parameters for developing the LR and ANN models for predicting moderate-to-severe OSA. As shown in Figure 1A, the 7328 participants in the training-testing group were randomly partitioned into training and testing sets for five-fold cross-validation. Table 3 showed that age, sex, BMI, and OSA severity distribution remained identical across the five testing sets and the independent validation sets.
| Variable | Testing set 1 | Testing set 2 | Testing set 3 | Testing set 4 | Testing set 5 | Independent validation set | p-value |
| n | 1466 | 1465 | 1466 | 1465 | 1466 | 2094 | |
| AHI≥15, n (%) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 1322 (63.1) | 1.000 |
| Age | 49.0 ± 0.4 | 49.0 ± 0.4 | 49.9 ± 0.4 | 49.1 ± 0.4 | 49.0 ± 0.4 | 49.3 ± 0.3 | 0.503 |
| Female, n (%) | 483 (33.0) | 438 (29.9) | 472 (32.2) | 445 (30.4) | 431 (29.4) | 649 (31.0) | 0.266 |
| BMI (kg/m2) | 26.9 ± 4.9 | 26.9 ± 4.9 | 26.9 ± 5.0 | 26.9 ± 5.0 | 26.8 ± 4.8 | 26.9 ± 4.9 | 0.973 |
| The p-value represents the multiple comparisons across the five testing and the independent validation sets. | |||||||
DownLoad:
CSV
Following adequate training, the five testing sets that achieved the best prediction performance after five-fold cross-validation were ANN models with 3, 14, 25, 27, and 29 hidden neurons, respectively. The mean AUCs of the training sets were 0.808 ± 0.004 (Figure 2A) and 0.808 ± 0.015 for the testing sets (Figure 2B). Among the generated ANN models, the model containing 25 hidden neurons obtained the highest testing AUC and was selected as the final model. After being applied to the independent validation set, the ANN-based model achieved the accuracy, sensitivity, specificity, PPV, and NPV of 76.4%, 87.7%, 56.9%, 77.7%, 73.0%, with the AUC of 0.807, respectively. Results of the LR model revealed a mean AUC of 0.802 ± 0.004 for the training sets (Figure 2C) and 0.801 ± 0.015 for the testing sets (Figure 2D). After being validated to the independent validation set, the accuracy, sensitivity, specificity, PPV, and NPV were 76.3%, 87.5%, 57.0%, 77.7%, and 72. 7% with an AUC of 0.806 for the LR-based model. Figure 3 shows the ROC curve and AUC of the final ANN and LR models applied to the independent validation set.
The present study validated age, sex, and BMI as parameters for use in LR- and ANN models to estimate the risk of OSA in participants who had undergone PSG. The established LR- and ANN-based model integrating the three parameters achieved a notable prediction performance with an AUC of 0.806 and 0.807 for predicting moderate-to-severe OSA in the validation set. The relatively high AUC values indicated that the models could distinguish participants between groups of interest, namely, participants with none-to-mild and those with moderate-to-severe OSA. Patients with OSA are prone to various complications and require interventions [2,7]; our models may contribute to further understanding these three variables' fundamental role in assessing OSA risk.
OSA has been recognized as a global health and economic burden, and early identification and diagnosis of OSA is a crucial issue in preventive medicine [4,6,7]. Since standard diagnostic tests for OSA rely on expensive in-laboratory PSG, which can be challenging to schedule in high-demand situations, identifying crucial factors and developing a simple and reliable tool to estimate the risk of OSA is clinically relevant [7,9].
Anatomical and physiological differences in the oropharynx, larynx and upper respiratory tract, changes in respiratory physiology associated with different sleep structures, and differences in the distribution of adipose tissue have been suggested as reasons for sex differences in developing OSA [18]. Deterioration of the collapsibility of the pharyngeal airway and the sensitivity of the respiratory control system may contribute to the pathogenesis of OSA with aging [2,18]. Hormonal changes in menopausal women have been linked to an increased risk of OSA and a decreased sex difference in OSA with age [19,20]. The prevalence of OSA has been reported to increase 1.4 to 3.2-fold with each 10-year increment in age, and men had a 1.7 to 3.0-fold higher odds of developing OSA than women [2,6]. These findings support the rationale for using age and sex as parameters in a generic prediction model to predict OSA risk.
Obesity is a well-recognized risk factor for developing OSA, and a higher BMI is associated with greater severity of OSA for both sexes [19,21]. As already mentioned, anthropometric profiles, such as body weight, BMI, central body fat distribution, and neck and waist circumference, are predisposed to the risk of OSA [2,6,10]. In patients with OSA, occlusion of the nasal and oro-pharynx may occur during sleep when the tongue and epiglottis move backward and lean together. Individuals with obesity have excess soft tissue around the collapsible pharynx, making the narrow airway more prone to collapse [2,3,21]. In line with the findings of other studies, we observed that the participants with moderate-to-severe OSA had a higher BMI. Although other anthropometric profiles have been associated with OSA, neck and waist circumference were considered surrogates for BMI [4,21,22]. The study by Huxley et al. has also suggested that "measures of central obesity were superior to BMI as discriminators of cardiovascular risk factors, although the differences were small and unlikely to be of clinical relevance" [23]. Furthermore, the same study demonstrated that "combining BMI with any measure of central obesity did not improve the discriminatory capability of the individual measures" [23]. In this study, we proposed BMI as a simple, well-known, and readily available predictor of OSA for the prediction models.
Despite the epidemiological reports indicating that age, sex, and BMI are associated with OSA, an integrated predicting model using only these three attributes to assess OSA risk has not yet been developed. Advances in machine learning techniques have improved the ability of researchers and clinical applications to analyze massive amounts of data and develop appropriate diagnostic and predictive methods [24,25,26,27]. Several studies have applied machine learning algorithms to predict OSA risk. However, most of these studies have incorporated more complex parameters, such as sleep questionnaires [28,29,30], at-home PSG [15], cephalometric images [31], neck and waist circumference measurements [28,29,32,33], and facial-cervical measurements [28,32], into the predictive models.
In this context, our study used a large hospital-based sample of adults, allowing accurate comparison between none-to-mild and moderate-to-severe cases. The generated LR- and ANN-boosted models exhibited relatively high AUCs with adequate accuracy, sensitivity, PPV, and NPV, revealing the ability of these models to distinguish between the groups of interest, namely individuals with moderate-to-severe and none-to-mild OSA. Furthermore, our models used simple and easy-to-obtain parameters (age, sex, and BMI), which may provide a direction for developing predictive models under conditions of limited and constrained information.
Both LR and ANN are widely used models in biomedicine. This study's LR and ANN models achieved similar predictive validity in predicting OSA in patients undergoing PSG, which may be related to our data structure. ANN has the advantage of analyzing complex, multidimensional and nonlinear variables [34,35,36,37]. When the predictor-outcome relationships in the dataset are relatively simple and linear, LR using linear combinations of variables may achieve performance that approximates that of ANN [36,38]. This consideration was contextualized in the conclusion by Dreiseitl and Ohno-Machado [38] that "there is no single algorithm that performs better than all other algorithms on any given data set and application area." and also summarized by Sargent [34] that "both methods should continue to be used and explored in a complementary manner."
The present study has several limitations. First, in our cohort, the incidence of moderate-to-severe OSA was 62.8%, which is higher than the reported prevalence of OSA in the general population [4,19,20,39]. In clinical practice, patients referred for a PSG exam are usually expected to have OSA or another sleep disorder; therefore, the participant population in this study cannot be representative of the general population with a very low prevalence of OSA. Second, the data set used in this study was limited to a single center and included only Asian ethnic groups. Therefore, the trained model may not be generalizable beyond the institution or other ethnic groups [28,29]. Third, this retrospective observational study did not exclude undiagnosed medical conditions associated with OSA, such as neurological, cardiovascular, and pulmonary disorders, which may affect the model's predictive performance. Our model did not contain detailed anthropometric imaging or measurements, which may limit the model's ability to distinguish the disease-specific causes of OSA. Finally, the relatively low specificity for the models (57.0% and 56.9%) suggests a higher rate of false positives, meaning that a healthy person is more likely to receive an incorrect diagnosis of OSA using the proposed tool. As Jonas et al. [30] concluded, "There is uncertainty about the accuracy or clinical utility of all potential screening tools." a PSG to verify OSA severity and confirm the diagnosis is warranted for patients with clinically suspected sleep apnea.
Our study validated the diagnostic relevance of age, sex, and BMI for OSA. The generated LR- and ANN-based model effectively predicted OSA risk using these three simple parameters. Our results, therefore, may improve insights into OSA-related essential risk factors and provide knowledge to facilitate individualized risk stratification and formulate diagnostic decisions for sleep apnea.
This work was funded by the National Science and Technology Council (NSTC 111-2314-B-038-132-MY3) to CC Chung.
The authors declare there is no conflict of interest.
| [1] |
Sprent J, Tough DF (1994) Lymphocyte life-span and memory. Science 265: 1395-1400. doi: 10.1126/science.8073282
|
| [2] |
Jamieson BD, Douek DC, Killian S, et al. (1999) Generation of functional thymocytes in the human adult. Immunity 10: 569-575. doi: 10.1016/S1074-7613(00)80056-4
|
| [3] |
Moxham VF, Karegli J, Phillips RE, et al. (2008) Homeostatic proliferation of lymphocytes results in augmented memory-like function and accelerated allograft rejection. J Immunol 180: 3910-3918. doi: 10.4049/jimmunol.180.6.3910
|
| [4] |
den Braber I, Mugwagwa T, Vrisekoop N, et al. (2012) Maintenance of peripheral naive T cells is sustained by thymus output in mice but not humans. Immunity 36: 288-297. doi: 10.1016/j.immuni.2012.02.006
|
| [5] |
Kozlov VA (2014) Homeostatic proliferation as a basis for the inevitable formation of total immunodeficiency. Med Immunol (Russia) 16: 403-408. doi: 10.15789/1563-0625-2014-5-403-408
|
| [6] |
Dummer W, Ernst B, LeRoy E, et al. (2001) Autologous regulation of naive T cell homeostasis within the T cell compartment. J Immunol 166: 2460-2468. doi: 10.4049/jimmunol.166.4.2460
|
| [7] |
Chalasani G, Dai Z, Konieczny BT, et al. (2002) Recall and propagation of allospecific memory T cells independent of secondary lymphoid organs. P Natl Acad Sci USA 99: 6175-6180. doi: 10.1073/pnas.092596999
|
| [8] |
Miller CN, Hartigan-O'Connor DJ, Lee MS, et al. (2013) IL-7 production in murine lymphatic endothelial cells and induction in the setting of peripheral lymphopenia. Int Immunol 25: 471-483. doi: 10.1093/intimm/dxt012
|
| [9] |
Anthony SM, Rivas SC, Colpitts SL, et al. (2016) Inflammatory signals regulate IL-15 in response to lymphodepletion. J Immunol 196: 4544-4552. doi: 10.4049/jimmunol.1600219
|
| [10] |
Ge Q, Rao VP, Cho BK, et al. (2001) Dependence of lymphopenia-induced T cell proliferation on the abundance of peptide/MHC epitopes and strength of their interaction with T cell receptors. Proc Natl Acad Sci USA 98: 1728-1733. doi: 10.1073/pnas.98.4.1728
|
| [11] |
Yamaki S, Ine S, Kawabe T, et al. (2014) OX40 and IL-7 play synergistic roles in the homeostatic proliferation of effector memory CD4+ T cells. Eur J Immunol 44: 3015-3025. doi: 10.1002/eji.201444701
|
| [12] |
Li O, Zheng P, Liu Y (2004) CD24 expression on T cells is required for optimal T cell proliferation in lymphopenic host. J Exp Med 200: 1083-1089. doi: 10.1084/jem.20040779
|
| [13] |
Min B, Yamane H, Hu-Li J, et al. (2005) Spontaneous and homeostatic proliferation of CD4 T cells are regulated by different mechanisms. J Immunol 174: 6039-6044. doi: 10.4049/jimmunol.174.10.6039
|
| [14] |
Ge Q, Rao VP, Cho BK, et al. (2001) Dependence of lymphopenia-induced T cell proliferation on the abundance of peptide/MHC epitopes and strength of their interaction with T cell receptors. Proc Natl Acad Sci USA 98: 1728-1733. doi: 10.1073/pnas.98.4.1728
|
| [15] |
Theofilopoulos AN, Dummer W, Kono DH (2001) T cell homeostasis and systemic autoimmunity. J Clin Invest 108: 335-340. doi: 10.1172/JCI200112173
|
| [16] |
Kassiotis G, Zamoyska R, Stockinger B (2003) Involvement of avidity for major histocompatibility complex in homeostasis of naive and memory T cells. J Exp Med 197: 1007-1016. doi: 10.1084/jem.20021812
|
| [17] |
Kieper WC, Burghardt JT, Surh CD (2004) A role for TCR affinity in regulating naive T cell homeostasis. J Immunol 172: 40-44. doi: 10.4049/jimmunol.172.1.40
|
| [18] |
Goldrath AW, Bogatzki LY, Bevan MJ (2000) Naive T cells transiently acquire a memory-like phenotype during homeostasis-driven proliferation. J Exp Med 192: 557-564. doi: 10.1084/jem.192.4.557
|
| [19] |
Klein L, Kyewski B, Allen PM, et al. (2014) Positive and negative selection of the T cell repertoire: what thymocytes see (and don't see). Nat Rev Immunol 14: 377-391. doi: 10.1038/nri3667
|
| [20] |
Dolton G, Zervoudi E, Rius C, et al. (2018) Optimized peptide-MHC multimer protocols for detection and isolation of autoimmune T-Cells. Front Immunol 9: 1378. doi: 10.3389/fimmu.2018.01378
|
| [21] |
Enouz S, Carrié L, Merkler D, et al. (2012) Autoreactive T cells bypass negative selection and respond to self-antigen stimulation during infection. J Exp Med 209: 1769-1779. doi: 10.1084/jem.20120905
|
| [22] |
Yan J, Mamula MJ (2002) Autoreactive T cells revealed in the normal repertoire: escape from negative selection and peripheral tolerance. J Immunol 168: 3188-3194. doi: 10.4049/jimmunol.168.7.3188
|
| [23] |
Boehncke WH, Brembilla NC (2019) Autoreactive T-lymphocytes in inflammatory skin diseases. Front Immunol 10: 1198. doi: 10.3389/fimmu.2019.01198
|
| [24] |
Koehli S, Naeher D, Galati-Fournier V, et al. (2014) Optimal T-cell receptor affinity for inducing autoimmunity. Proc Natl Acad Sci USA 111: 17248-17253. doi: 10.1073/pnas.1402724111
|
| [25] |
Mueller DL (2010) Mechanisms maintaining peripheral tolerance. Nat Immunol 11: 21-27. doi: 10.1038/ni.1817
|
| [26] |
Vrisekoop N, Monteiro JP, Mandl JN, et al. (2014) Revisiting thymic positive selection and the mature T cell repertoire for antigen. Immunity 41: 181-190. doi: 10.1016/j.immuni.2014.07.007
|
| [27] |
Qi Q, Liu Y, Cheng Y, et al. (2014) Diversity and clonal selection in the human T-cell repertoire. Proc Natl Acad Sci USA 111: 13139-13144. doi: 10.1073/pnas.1409155111
|
| [28] |
Zheng SG, Wang J, Wang P, et al. (2007) IL-2 is essential for TGF-beta to convert naive CD4+CD25− cells to CD25+Foxp3+ regulatory T cells and for expansion of these cells. J Immunol 178: 2018-2027. doi: 10.4049/jimmunol.178.4.2018
|
| [29] |
Pacholczyk R, Ignatowicz H, Kraj P, et al. (2006) Origin and T cell receptor diversity of Foxp3+CD4+CD25+ T cells. Immunity 25: 249-259. doi: 10.1016/j.immuni.2006.05.016
|
| [30] |
Wyss L, Stadinski BD, King CG, et al. (2016) Affinity for self-antigen selects Treg cells with distinct functional properties. Nat Immunol 17: 1093-1101. doi: 10.1038/ni.3522
|
| [31] |
Bolotin DA, Poslavsky S, Davydov AN, et al. (2017) Antigen receptor repertoire profiling from RNA-seq data. Nat Biotechnol 35: 908-911. doi: 10.1038/nbt.3979
|
| [32] |
Jörg S, Grohme DA, Erzler M, et al. (2016) Environmental factors in autoimmune diseases and their role in multiple sclerosis. Cell Mol Life Sci 73: 4611-4622. doi: 10.1007/s00018-016-2311-1
|
| [33] |
Ting YT, Petersen J, Ramarathinam SH, et al. (2018) The interplay between citrullination and HLA-DRB1 polymorphism in shaping peptide binding hierarchies in rheumatoid arthritis. J Biol Chem 293: 3236-3251. doi: 10.1074/jbc.RA117.001013
|
| [34] |
Schulze-Koops H (2004) Lymphopenia and autoimmune diseases. Arthritis Res Ther 6: 178-180. doi: 10.1186/ar1208
|
| [35] |
Symmons DP, Farr M, Salmon M, et al. (1989) Lymphopenia in rheumatoid arthritis. J R Soc Med 82: 462-463. doi: 10.1177/014107688908200806
|
| [36] |
Jones JL, Thompson SA, Loh P, et al. (2013) Human autoimmunity after lymphocyte depletion is caused by homeostatic T-cell proliferation. Proc Natl Acad Sci USA 110: 20200-20205. doi: 10.1073/pnas.1313654110
|
| [37] |
Koetz K, Bryl E, Spickschen K, et al. (2000) T cell homeostasis in patients with rheumatoid arthritis. Proc Natl Acad Sci USA 97: 9203-9208. doi: 10.1073/pnas.97.16.9203
|
| [38] |
Virot E, Duclos A, Adelaide L, et al. (2017) Autoimmune diseases and HIV infection: A cross-sectional study. Medicine (Baltimore) 96: e5769. doi: 10.1097/MD.0000000000005769
|
| [39] | Jerne NK (1974) Towards a network theory of the immune system. Ann Immunol (Paris) 125: 373-389. |
| [40] |
Ernst B, Lee DS, Chang JM, et al. (1999) The peptide ligands mediating positive selection in the thymus control T cell survival and homeostatic proliferation in the periphery. Immunity 11: 173-181. doi: 10.1016/S1074-7613(00)80092-8
|
| [41] |
Milam AAV, Bartleson JM, Buck MD, et al. (2020) Tonic TCR signaling inversely regulates the basal metabolism of CD4+ T cells. ImmunoHorizons 4: 485-497. doi: 10.4049/immunohorizons.2000055
|
| [42] |
Hochweller K, Wabnitz GH, Samstag Y, et al. (2010) Dendritic cells control T cell tonic signaling required for responsiveness to foreign antigen. Proc Natl Acad Sci USA 107: 5931-5936. doi: 10.1073/pnas.0911877107
|
| [43] |
Stefanová I, Dorfman JR, Germain RN (2002) Self-recognition promotes the foreign antigen sensitivity of naive T lymphocytes. Nature 420: 429-434. doi: 10.1038/nature01146
|
| [44] |
Miller CN, Hartigan-O'Connor DJ, Lee MS, et al. (2013) IL-7 production in murine lymphatic endothelial cells and induction in the setting of peripheral lymphopenia. Int Immunol 25: 471-483. doi: 10.1093/intimm/dxt012
|
| [45] |
Saeed S, Revell PA (2001) Production and distribution of interleukin 15 and its receptors (IL-15Ralpha and IL-R2beta) in the implant interface tissues obtained during revision of failed total joint replacement. Int J Exp Pathol 82: 201-209. doi: 10.1111/j.1365-2613.2001.iep185.x
|
| [46] |
Regamey N, Obregon C, Ferrari-Lacraz S, et al. (2007) Airway epithelial IL-15 transforms monocytes into dendritic cells. Am J Respir Cell Mol Biol 37: 75-84. doi: 10.1165/rcmb.2006-0235OC
|
| [47] |
Park CS, Yoon SO, Armitage RJ, et al. (2004) Follicular dendritic cells produce IL-15 that enhances germinal center B cell proliferation in membrane-bound form. J Immunol 173: 6676-6683. doi: 10.4049/jimmunol.173.11.6676
|
| [48] |
Niu N, Qin X (2013) New insights into IL-7 signaling pathways during early and late T cell development. Cell Mol Immunol 10: 187-189. doi: 10.1038/cmi.2013.11
|
| [49] |
Mishra A, Sullivan L, Caligiuri MA (2014) Molecular pathways: interleukin-15 signaling in health and in cancer. Clin Cancer Res 20: 2044-2050. doi: 10.1158/1078-0432.CCR-12-3603
|
| [50] |
Tan JT, Ernst B, Kieper WC, et al. (2002) Interleukin (IL)-15 and IL-7 jointly regulate homeostatic proliferation of memory phenotype CD8+ cells but are not required for memory phenotype CD4+ cells. J Exp Med 195: 1523-1532. doi: 10.1084/jem.20020066
|
| [51] |
Do JS, Visperas A, Oh K, et al. (2012) Memory CD4 T cells induce selective expression of IL-27 in CD8+ dendritic cells and regulate homeostatic naive T cell proliferation. J Immunol 188: 230-237. doi: 10.4049/jimmunol.1101908
|
| [52] |
Lenz DC, Kurz SK, Lemmens E, et al. (2004) IL-7 regulates basal homeostatic proliferation of antiviral CD4+ T cell memory. Proc Natl Acad Sci USA 101: 9357-9362. doi: 10.1073/pnas.0400640101
|
| [53] |
Berard M, Brandt K, Bulfone-Paus S, et al. (2003) IL-15 promotes the survival of naive and memory phenotype CD8+ T cells. J Immunol 170: 5018-5026. doi: 10.4049/jimmunol.170.10.5018
|
| [54] |
Prlic M, Blazar BR, Khoruts A, et al. (2001) Homeostatic expansion occurs independently of costimulatory signals. J Immunol 167: 5664-5668. doi: 10.4049/jimmunol.167.10.5664
|
| [55] |
Do JS, Min B (2009) Differential requirements of MHC and of DCs for endogenous proliferation of different T-cell subsets in vivo. Proc Natl Acad Sci USA 106: 20394-20398. doi: 10.1073/pnas.0909954106
|
| [56] |
Hawse WF, Cattley RT (2019) T cells transduce T-cell receptor signal strength by generating different phosphatidylinositols. J Biol Chem 294: 4793-4805. doi: 10.1074/jbc.RA118.006524
|
| [57] |
Cosway EJ, Lucas B, James KD, et al. (2017) Redefining thymus medulla specialization for central tolerance. J Exp Med 214: 3183-3195. doi: 10.1084/jem.20171000
|
| [58] |
Logunova NN, Kriukova VV, Shelyakin PV, et al. (2020) MHC-II alleles shape the CDR3 repertoires of conventional and regulatory naïve CD4+ T cells. Proc Natl Acad Sci USA 117: 13659-13669. doi: 10.1073/pnas.2003170117
|
| [59] |
Kieper WC, Troy A, Burghardt JT, et al. (2005) Recent immune status determines the source of antigens that drive homeostatic T cell expansion. J Immunol 174: 3158-3163. doi: 10.4049/jimmunol.174.6.3158
|
| [60] |
Do JS, Foucras G, Kamada N, et al. (2012) Both exogenous commensal and endogenous self antigens stimulate T cell proliferation under lymphopenic conditions. Cell Immunol 272: 117-123. doi: 10.1016/j.cellimm.2011.11.002
|
| [61] |
Min B (2018) Spontaneous T cell proliferation: A physiologic process to create and maintain homeostatic balance and diversity of the immune system. Front Immunol 9: 547. doi: 10.3389/fimmu.2018.00547
|
| [62] |
Schluns KS, Kieper WC, Jameson SC, et al. (2000) Interleukin-7 mediates the homeostasis of naïve and memory CD8 T cells in vivo. Nat Immunol 1: 426-432. doi: 10.1038/80868
|
| [63] |
Ku CC, Murakami M, Sakamoto A, et al. (2000) Control of homeostasis of CD8+ memory T cells by opposing cytokines. Science 288: 675-678. doi: 10.1126/science.288.5466.675
|
| [64] |
Tan JT, Dudl E, LeRoy E, et al. (2001) IL-7 is critical for homeostatic proliferation and survival of naive T cells. Proc Natl Acad Sci USA 98: 8732-8737. doi: 10.1073/pnas.161126098
|
| [65] |
Goldrath AW, Sivakumar PV, Glaccum M, et al. (2002) Cytokine requirements for acute and Basal homeostatic proliferation of naive and memory CD8+ T cells. J Exp Med 195: 1515-1522. doi: 10.1084/jem.20020033
|
| [66] |
Troy AE, Shen H (2003) Cutting edge: homeostatic proliferation of peripheral T lymphocytes is regulated by clonal competition. J Immunol 170: 672-676. doi: 10.4049/jimmunol.170.2.672
|
| [67] |
Min B, Foucras G, Meier-Schellersheim M, et al. (2004) Spontaneous proliferation, a response of naive CD4 T cells determined by the diversity of the memory cell repertoire. Proc Natl Acad Sci USA 101: 3874-3879. doi: 10.1073/pnas.0400606101
|
| [68] |
Kawabe T, Sun SL, Fujita T, et al. (2013) Homeostatic proliferation of naive CD4+ T cells in mesenteric lymph nodes generates gut-tropic Th17 cells. J Immunol 190: 5788-5798. doi: 10.4049/jimmunol.1203111
|
| [69] |
Dummer W, Niethammer AG, Baccala R, et al. (2002) T cell homeostatic proliferation elicits effective antitumor autoimmunity. J Clin Invest 110: 185-192. doi: 10.1172/JCI0215175
|
| [70] |
Jones JL, Thompson SA, Loh P, et al. (2013) Human autoimmunity after lymphocyte depletion is caused by homeostatic T-cell proliferation. Proc Natl Acad Sci USA 110: 20200-20205. doi: 10.1073/pnas.1313654110
|
| [71] |
Koetz K, Bryl E, Spickschen K, et al. (2000) T cell homeostasis in patients with rheumatoid arthritis. Proc Natl Acad Sci USA 97: 9203-9208. doi: 10.1073/pnas.97.16.9203
|
| [72] |
Arruda LCM, Lima-Júnior JR, Clave E, et al. (2018) Homeostatic proliferation leads to telomere attrition and increased PD-1 expression after autologous hematopoietic SCT for systemic sclerosis. Bone Marrow Transpl 53: 1319-1327. doi: 10.1038/s41409-018-0162-0
|
| [73] |
O'Bryan JM, Woda M, Co M, et al. (2013) Telomere length dynamics in human memory T cells specific for viruses causing acute or latent infections. Immun Ageing 10: 37. doi: 10.1186/1742-4933-10-37
|
| [74] |
Fortner KA, Bond JP, Austin JW, et al. (2017) The molecular signature of murine T cell homeostatic proliferation reveals both inflammatory and immune inhibition patterns. J Autoimmun 82: 47-61. doi: 10.1016/j.jaut.2017.05.003
|
| [75] |
Scharer CD, Fortner KA, Dragon JA, et al. (2020) Selective DNA demethylation accompanies T cell homeostatic proliferation and gene regulation in lupus-prone lpr mice. ImmunoHorizons 4: 679-687. doi: 10.4049/immunohorizons.2000078
|
| [76] |
Sakaguchi S, Yamaguchi T, Nomura T, et al. (2008) Regulatory T cells and immune tolerance. Cell 133: 775-787. doi: 10.1016/j.cell.2008.05.009
|
| [77] |
Shevyrev D, Tereshchenko V (2020) Treg heterogeneity, function, and homeostasis. Front Immunol 10: 3100. doi: 10.3389/fimmu.2019.03100
|
| [78] |
Tang AL, Teijaro JR, Njau MN, et al. (2008) CTLA4 expression is an indicator and regulator of steady-state CD4+ FoxP3+ T cell homeostasis. J Immunol 181: 1806-1813. doi: 10.4049/jimmunol.181.3.1806
|
| [79] |
Chen Y, Shen S, Gorentla BK, et al. (2012) Murine regulatory T cells contain hyperproliferative and death-prone subsets with differential ICOS expression. J Immunol 188: 1698-1707. doi: 10.4049/jimmunol.1102448
|
| [80] |
Sojka DK, Hughson A, Fowell DJ (2009) CTLA-4 is required by CD4+CD25+ Treg to control CD4+ T-cell lymphopenia-induced proliferation. Eur J Immunol 39: 1544-1551. doi: 10.1002/eji.200838603
|
| [81] |
Bolton HA, Zhu E, Terry AM, et al. (2015) Selective Treg reconstitution during lymphopenia normalizes DC costimulation and prevents graft-versus-host disease. J Clin Invest 125: 3627-3641. doi: 10.1172/JCI76031
|
| [82] |
Winstead CJ, Fraser JM, Khoruts A (2008) Regulatory CD4+CD25+Foxp3+ T cells selectively inhibit the spontaneous form of lymphopenia-induced proliferation of naive T cells. J Immunol 180: 7305-7317. doi: 10.4049/jimmunol.180.11.7305
|
| [83] |
Ben Ahmed M, Belhadj Hmida N, Moes N, et al. (2009) IL-15 renders conventional lymphocytes resistant to suppressive functions of regulatory T cells through activation of the phosphatidylinositol 3-kinase pathway. J Immunol 182: 6763-6770. doi: 10.4049/jimmunol.0801792
|
| [84] |
Heninger AK, Theil A, Wilhelm C, et al. (2012) IL-7 abrogates suppressive activity of human CD4+CD25+FOXP3+ regulatory T cells and allows expansion of alloreactive and autoreactive T cells. J Immunol 189: 5649-5658. doi: 10.4049/jimmunol.1201286
|
Figures(2)
Daniil Shevyrev, Valeriy Tereshchenko, Olesya Manova, Vladimir kozlov. Homeostatic proliferation as a physiological process and a risk factor for autoimmune pathology[J]. AIMS Allergy and Immunology, 2021, 5(1): 18-32. doi: 10.3934/Allergy.2021002
| Variables | Moderate-to-severe OSA (AHI ≥ 15) |
None-to-mild OSA (AHI < 15) |
p-value |
| n | 5, 920 | 3, 502 | |
| Age | 50.8 ± 13.5 | 46.6 ± 14.4 | < 0.0001 |
| Female, n (%) | 1266 (21.4) | 1652 (47.2) | < 0.0001 |
| Height (cm) | 167.2 ± 8.5 | 164.4 ± 8.9 | < 0.0001 |
| Weight (kg) | 79.7±16.1 | 66.1 ± 13.1 | < 0.0001 |
| BMI (kg/m2) | 28.4 ± 4.8 | 24.3 ± 3.8 | < 0.0001 |
| Polysomnography | |||
| Total recording time (min) | 366.6±11.9 | 367.8±13.4 | < 0.0001 |
| Total sleep time (min) | 275.8±60.2 | 281.4±61.1 | < 0.0001 |
| Sleep efficiency (%) | 75.2±16.4 | 76.5±16.5 | 0.0002 |
| Non-REM sleep (min) | 240.6±51.7 | 242.2±51.3 | 0.133 |
| REM sleep (min) | 35.2±22.5 | 39.1±23.7 | < 0.0001 |
| AHI | 44.7±23.8 | 6.4±4.5 | < 0.0001 |
| SO2 mean (%) | 94.5±2.6 | 96.6±1.3 | < 0.0001 |
| Data are expressed as mean ± standard deviation. The p-value represents the comparison between the AHI ≥ 15 and AHI < 15 groups. | |||
DownLoad:
CSV
| Unadjusted Model | Adjusted Model | |||
| Variables | OR | 95% CI | OR | 95% CI |
| Age | 1.02*† | 1.02-1.03 | 1.04*† | 1.04-1.05 |
| Female | 0.30* | 0.28-0.33 | 0.28* | 0.25-0.32 |
| BMI | 1.27*† | 1.25-1.28 | 1.27*† | 1.26-1.29 |
| The odds ratio (OR) and 95% confidence intervals (CI) represent the association between the variables and the presence of moderate-to-severe OSA (AHI ≥ 15). *p < 0.0001; †Odds ratio of per unit changes. | ||||
DownLoad:
CSV
| Variable | Testing set 1 | Testing set 2 | Testing set 3 | Testing set 4 | Testing set 5 | Independent validation set | p-value |
| n | 1466 | 1465 | 1466 | 1465 | 1466 | 2094 | |
| AHI≥15, n (%) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 1322 (63.1) | 1.000 |
| Age | 49.0 ± 0.4 | 49.0 ± 0.4 | 49.9 ± 0.4 | 49.1 ± 0.4 | 49.0 ± 0.4 | 49.3 ± 0.3 | 0.503 |
| Female, n (%) | 483 (33.0) | 438 (29.9) | 472 (32.2) | 445 (30.4) | 431 (29.4) | 649 (31.0) | 0.266 |
| BMI (kg/m2) | 26.9 ± 4.9 | 26.9 ± 4.9 | 26.9 ± 5.0 | 26.9 ± 5.0 | 26.8 ± 4.8 | 26.9 ± 4.9 | 0.973 |
| The p-value represents the multiple comparisons across the five testing and the independent validation sets. | |||||||
DownLoad:
CSV
| Variables | Moderate-to-severe OSA (AHI ≥ 15) |
None-to-mild OSA (AHI < 15) |
p-value |
| n | 5, 920 | 3, 502 | |
| Age | 50.8 ± 13.5 | 46.6 ± 14.4 | < 0.0001 |
| Female, n (%) | 1266 (21.4) | 1652 (47.2) | < 0.0001 |
| Height (cm) | 167.2 ± 8.5 | 164.4 ± 8.9 | < 0.0001 |
| Weight (kg) | 79.7±16.1 | 66.1 ± 13.1 | < 0.0001 |
| BMI (kg/m2) | 28.4 ± 4.8 | 24.3 ± 3.8 | < 0.0001 |
| Polysomnography | |||
| Total recording time (min) | 366.6±11.9 | 367.8±13.4 | < 0.0001 |
| Total sleep time (min) | 275.8±60.2 | 281.4±61.1 | < 0.0001 |
| Sleep efficiency (%) | 75.2±16.4 | 76.5±16.5 | 0.0002 |
| Non-REM sleep (min) | 240.6±51.7 | 242.2±51.3 | 0.133 |
| REM sleep (min) | 35.2±22.5 | 39.1±23.7 | < 0.0001 |
| AHI | 44.7±23.8 | 6.4±4.5 | < 0.0001 |
| SO2 mean (%) | 94.5±2.6 | 96.6±1.3 | < 0.0001 |
| Data are expressed as mean ± standard deviation. The p-value represents the comparison between the AHI ≥ 15 and AHI < 15 groups. | |||
| Unadjusted Model | Adjusted Model | |||
| Variables | OR | 95% CI | OR | 95% CI |
| Age | 1.02*† | 1.02-1.03 | 1.04*† | 1.04-1.05 |
| Female | 0.30* | 0.28-0.33 | 0.28* | 0.25-0.32 |
| BMI | 1.27*† | 1.25-1.28 | 1.27*† | 1.26-1.29 |
| The odds ratio (OR) and 95% confidence intervals (CI) represent the association between the variables and the presence of moderate-to-severe OSA (AHI ≥ 15). *p < 0.0001; †Odds ratio of per unit changes. | ||||
| Variable | Testing set 1 | Testing set 2 | Testing set 3 | Testing set 4 | Testing set 5 | Independent validation set | p-value |
| n | 1466 | 1465 | 1466 | 1465 | 1466 | 2094 | |
| AHI≥15, n (%) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 919 (62.7) | 920 (62.8) | 1322 (63.1) | 1.000 |
| Age | 49.0 ± 0.4 | 49.0 ± 0.4 | 49.9 ± 0.4 | 49.1 ± 0.4 | 49.0 ± 0.4 | 49.3 ± 0.3 | 0.503 |
| Female, n (%) | 483 (33.0) | 438 (29.9) | 472 (32.2) | 445 (30.4) | 431 (29.4) | 649 (31.0) | 0.266 |
| BMI (kg/m2) | 26.9 ± 4.9 | 26.9 ± 4.9 | 26.9 ± 5.0 | 26.9 ± 5.0 | 26.8 ± 4.8 | 26.9 ± 4.9 | 0.973 |
| The p-value represents the multiple comparisons across the five testing and the independent validation sets. | |||||||
