Abstract
Background: Cerebral palsy (CP) is a syndrome of non-progressive motor dysfunction caused by early brain development injury. Recent evidence has shown that immunological abnormalities are associated with an increased risk of CP.
Methods: We recruited 782 children with CP as the case group and 770 healthy children as the control group. The association between IL-23R single nucleotide polymorphisms (SNPs; namely, rs10889657, rs6682925, rs1884444, rs17375018, rs1004819, rs11805303, and rs10889677) and CP was studied by using a case–control method and SHEsis online software. Subgroup analysis based on complications and clinical subtypes was also carried out.
Results: There were differences in the allele and genotype frequencies between CP cases and controls at the rs11805303 and rs10889677 SNPs (Pallele = 0.014 and 0.048, respectively; Pgenotype = 0.023 and 0.008, respectively), and the difference in genotype frequency of rs10889677 remained significant after Bonferroni correction (Pgenotype = 0.048). Subgroup analysis revealed a more significant association of rs10889677 with CP accompanied by global developmental delay (Pgenotype = 0.024 after correction) and neonatal encephalopathy (Pgenotype = 0.024 after correction).
Conclusion: The present results showed a significant association between IL-23R and CP, suggesting that IL-23R may play a potential role in CP pathogenesis.
Introduction
Cerebral palsy (CP) is a central motor disorder syndrome that manifests with abnormal muscle tension and motor function (Cheng et al., 2011; Wu et al., 2011). Individuals with CP often exhibit sensory, perceptual, cognitive, communication, and behavioral disorders, as well as epilepsy and secondary musculoskeletal problems (Beliakoff and Sun, 2006). Although perinatal medicine has developed rapidly in recent years, epidemiological studies have shown that the incidence of CP has remained stable at 2–3.5 children out of every 1,000 children (Pennington et al., 2013). CP has been a prominent disease among children's disabilities for a long time due to a lack of effective treatments (Tator et al., 2007; Tatla et al., 2013). CP seriously affects the quality of life of individuals and also brings a heavy financial burden to families and society (Duval et al., 2003; Yunus and Lima, 2009). The identification of its etiology and pathogenesis is essential to the prevention and control of CP.
CP is caused by non-progressive brain damage during the development of the fetus or infant, which can be divided into congenital and acquired damage. Congenital non-progressive injuries are caused by prenatal developmental defects, such as genetic defects, developmental defects, malformations, intrauterine infections, etc. Acquired injuries are caused by postpartum acquired factors, such as premature labor, asphyxia, hypoxic–ischemic encephalopathy (HIE), low birth weight (LOW), and pathogenic jaundice (Sachdev et al., 2001; Jacobsson and Hagberg, 2004). Increasing evidence now indicates that genetic factors are likely to play an important role in CP pathogenesis. In general, CP is regarded as the result of the combined effects of multiple genes and environmental factors. Moreover, CP has been reported to have multiple susceptibility genes, including IL-6, NOS1, OLIG2, ATG5, and ATG7 (Xu et al., 2017; Yu et al., 2018; Sun et al., 2019; Xia et al., 2019).
A great deal of evidence suggests that neuroinflammation has been found to participate in, modulate, and even induce the pathological process of immature brain injury and various cytokines have been associated with CP and neurodevelopmental disability. Previous studies have found that immature brain injury induced by secondary inflammation is one of the important pathological mechanisms of hypoxic–ischemic brain injury (Elovitz et al., 2011; Albertsson et al., 2014; Marshall and Plotkin, 2019). Abnormal activation of cytokines can cause brain damage, and fetuses are more likely to produce inflammatory reactions or brain damage after being affected by inflammatory cytokines due to immature brain development, thus further affecting the normal development of the brain. At present, many inflammatory cytokines have been reported to be significantly associated with CP or neurodevelopmental disorders, such as IL-6, IL-8, IL-10, and IL-17 (Strle et al., 2001; Chiricozzi et al., 2011; Chen et al., 2013; Magalhaes et al., 2019).
Interleukin-23 (IL-23), also known as p19, is a member of the IL-12 heterodimer cytokine family (Oppmann et al., 2000), which is mainly produced by activated dendritic cells, macrophages, and monocytes. IL-23 plays an important role in the regulation of tissue homeostasis and congenital or adaptive immunity. IL-23 is involved in the pathogenesis of many chronic inflammatory diseases, such as psoriasis, arthritis, inflammatory bowel disease, and multiple sclerosis (MS). Drugs targeting IL-23 have been used in clinical research regarding immunologic diseases (Wiekowski et al., 2001; Schon and Erpenbeck, 2018). IL-23 binds to the IL-23 receptor (IL-23R) through its N-terminal immunoglobulin domain, which activates downstream signaling pathways and exerts biological functions. Human IL-23 receptors (IL-23R) are mainly expressed in activated memory T cells, natural killer (NK) cells, and intrinsic immune cells (ILCs). Its extracellular domain contains a signal sequence, one N-terminal Ig-like domain, and two cytokine receptor domains. In a genome-wide association study in 2006, IL-23R was significantly associated with Crohn's disease, an inflammatory bowel disease. The A allele of rs11209026, a low-frequency IL-23R variant, was negatively correlated with Crohn's disease (Bloch et al., 2018). At the same time, some studies have shown that IL-23R mutation is significantly correlated with inflammatory demyelinating diseases, such as MS (Li et al., 2016).
Based on the above information, we speculate that IL-23R may be associated with susceptibility to CP, but no relevant studies have been reported thus far. Then, we used a case–control study to explore the possible association of IL-23R with CP, which will provide genetic evidence for evaluating the role of IL-23R in the etiology of CP and its related potential mechanisms.
Materials and Methods
Participants
In this study, 782 children with CP and 771 healthy controls were recruited from the centers for CP rehabilitation and Child Healthcare Departments in the Third Affiliated Hospital of Zhengzhou University, Zhengzhou Children's Hospital. This study was approved by the ethics committee of Zhengzhou University. The guardians of these participants provided written informed consent. The case group comprised 542 males (69.3%) and 240 females (30.7%), and the mean age was 18.5 ± 15.4 months. The control group comprised 771 healthy children, including 515 males (66.8%) and 256 females (33.2%), and the mean age was 19.3 ± 16.8 months (Table 1).
Table 1
| Characteristic | CP cases (n = 782) | Controls (n = 771) |
|---|---|---|
| Sex (male:female) | 542:240 | 515:256 |
| Preterm (<37 weeks) | 47 | 10 |
| <2,500 g | 40 | 2 |
| Birth asphyxia | 234 | 13 |
| Type of CP | ||
| Spastic CP | 522 | NA |
| CP with quadriplegia | 284 | NA |
| CP with diplegia | 126 | NA |
| Complications | ||
| CP with PVL | 67 | NA |
| CP with HIE | 108 | NA |
| CP with GDD | 299 | NA |
| Type of CP | 310 | NA |
| Maternal factors | ||
| PROM | 71 | 26 |
| TPL | 58 | 0 |
| PIH | 26 | 7 |
Clinical characteristics of all participants.
CP, cerebral palsy; PVL, periventricular leukomalacia; HIE, hypoxic–ischemic encephalopathy; GDD, global developmental delay; PROM, premature rupture of membrane; TPL, threatened premature labor; PIH, pregnancy-induced hypertension.
CP Diagnosis, Classification, and Exclusion Criteria
In the case group, we excluded children diagnosed with congenital metabolic diseases and myopathy as well as children with a family history of nervous system diseases. Pediatric rehabilitation specialists confirmed the CP diagnosis using standard criteria related to non-progressive disorders of movement control and posture (Rosenbaum et al., 2007). All participants received a detailed clinical evaluation with comprehensive pre-test counseling.
The available clinical information included demographic variables, such as sex, gestational age, mode of delivery, singletons, and twins, as well as the known risk factors [such as pregnancy-induced hypertension (PIH), perinatal asphyxia, and threatened premature labor], CP complications [such as global developmental delay (GDD)], and neonatal complications (such as HIE).
GDD diagnosis is limited to individuals under the age of 5 years old when the clinical severity level cannot be reliably assessed during early childhood. GDD is diagnosed when an individual fails to meet the expected developmental milestones in several areas of intellectual functioning and applies to individuals who are unable to undergo systematic assessments of intellectual functioning, including children who are too young to participate in standardized testing. Neonatal encephalopathy (NE) is a clinical syndrome that includes HIE, intracranial hemorrhage, various metabolic disorders, neurodegenerative diseases, and so on; its diagnosis requires at least two senior neonatologists.
Genotyping and Statistical Analysis
Peripheral blood samples were obtained from the subjects for genomic DNA extraction. According to the single nucleotide polymorphism (SNP) location in IL-23R, a minor allele frequency (MAF) >0.1, and potential function, we selected seven SNPs (rs10889657, rs6682925, rs1884444, rs17375018, rs1004819, rs11805303, and rs10889677, Figure 1) as candidates and genotyped them by the MassARRAY system. Shanghai Perchant Biotechnology Co., Ltd. synthesized primers and probes.
Figure 1
We performed statistical analysis with SHEsis, an online program (http://analysis.bio-x.cn/) that can test Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium (LD) and calculate allele frequencies, genotype frequencies, and haplotype frequencies for each SNP locus in the case group and control group. The P values were two-tailed, and we considered P < 0.05 to be significant. We also calculated the odds ratio (OR) and its 95% confidence interval (CI). We employed the Bonferroni correction to account for multiple testing on each individual SNP and haplotype. We used the G*power 3.1 software to evaluate the statistical efficacy.
Results
Overall Analysis
By performing power calculations, the sample size utilized in the present study has >85% power to detect a significant association (α < 0.05) when using an effect size index of 0.1. The genotype distributions of rs17375018 among control subjects showed a significant deviation from HWE (P = 0.011); therefore, we only analyzed the remaining six SNPs, namely, rs10889657, rs6682925, rs1884444, rs1004819, rs11805303, and rs10889677. There were two linked LD blocks with coefficient D′ value more than 0.8 (Figure 2). Block 1 included rs10889657, rs6682925, and rs1884444, whereas block 2 comprised rs1004819, rs11805303, and rs10889677.
Figure 2
For all subjects, the allele frequencies of rs11805303 (P = 0.014) and rs10889677 (P = 0.048) and the genotype frequencies of rs10889657 (P = 0.028), rs11805303 (P = 0.023), and rs10889677 (P = 0.008) were significantly different between CP patients and controls. After Bonferroni correction, only the rs10889677 AA genotype frequency was significantly more enriched in CP children than in controls (OR = 1.178, 95% CI = 1.002–1.386, Pc = 0.048) (Table 2).
Table 2
| Group | Allele frequency | p value | OR [95% CI] | Genotype frequency | p-value | H-W | |||
|---|---|---|---|---|---|---|---|---|---|
| rs10889657 | C | T | C/C | C/T | T/T | ||||
| CP | 347 (0.267) | 953 (0.733) | 0.09 | 0.861 [0.725–1.023] | 59 (0.091) | 229 (0.352) | 362 (0.557) | 0.028a | 0.011 |
| Control | 372 (0.297) | 880 (0.703) | 53 (0.085) | 266(0.425) | 307 (0.490) | 0.665 | |||
| rs6682925 | C | T | C/C | C/T | T/T | ||||
| CP | 625 (0.422) | 855 (0.578) | 0.406 | 0.940 [0.813–1.088] | 133 (0.180) | 359 (0.485) | 248 (0.335) | 0.681 | 0.876 |
| Control | 634 (0.437) | 827 (0.563) | 139 (0.189) | 365 (0.497) | 231 (0.314) | 0.807 | |||
| rs1884444 | G | T | G/G | G/T | T/T | ||||
| CP | 517 (0.344) | 985 (0.656) | 0.169 | 0.900 [0.776–1.046] | 91 (0.121) | 335 (0.446) | 325 (0.433) | 0.275 | 0.744 |
| Control | 556 (0.368) | 954 (0.632) | 97 (0.128) | 362 (0.479) | 296 (0.392) | 0.402 | |||
| rs1004819 | A | G | A/A | A/G | G/G | ||||
| CP | 532 (0.566) | 408 (0.434) | 0.05 | 1.194 [1.000–1.425] | 168 (0.357) | 196 (0.417) | 106 (0.226) | 0.063 | 0.001 |
| Control | 544 (0.522) | 498 (0.478) | 150 (0.288) | 244 (0.468) | 127 (0.244) | 0.16 | |||
| rs11805303 | C | T | C/C | C/T | T/T | ||||
| CP | 628 (0.415) | 884 (0.585) | 0.014b | 0.835 [0.723–0.964] | 150 (0.198) | 328 (0.434) | 278 (0.368) | 0.023c | 0.003 |
| Control | 696 (0.460) | 818 (0.540) | 167 (0.221) | 362 (0.478) | 228 (0.301) | 0.304 | |||
| rs10889677 | A | C | A/A | A/C | C/C | ||||
| CP | 1121 (0.750) | 373 (0.250) | 0.048d | 1.178 [1.002–1.386] | 429 (0.574) | 263 (0.352) | 55 (0.074) | 0.008e | 0.099 |
| Control | 1079 (0.718) | 423 (0.282) | 378 (0.503) | 323 (0.430) | 50 (0.067) | 0.085 | |||
Allele and genotype frequency analysis of SNPs of IL-23R.
OR, odds ratio; CI, confidence interval.The P values corrected by Bonferroni are
0.168,
0.084,
0.138,
0.288,
0.048.
Haplotype analysis is a powerful strategy to determine whether or not the above-mentioned CP-associated SNPs have a greater effect when analyzed together. Hence, we performed a haplotype analysis of rs11805303 and rs10889677 SNPs. The haplotypes CC (P = 0.036) and TA (P = 0.013) presented a significant association with CP; the positive association of TA with CP was significant even after Bonferroni correction (P = 0.039). Furthermore, there was a statistically significant global effect of the haplotype (P = 0.039) (Table 3).
Table 3
| Haplotype | Case (frequency) | Control (frequency) | P value | OR (95% CI) |
|---|---|---|---|---|
| CA | 250.91 (0.171) | 271.19 (0.183) | 0.413244 | 0.924 (0.765–1.117) |
| CC | 356.09 (0.242) | 410.81 (0.277) | 0.036346 | 0.838 (0.711–0.989) |
| TA | 850.09 (0.578) | 793.81 (0.536) | 0.013013 | 1.203 (1.040–1.392) |
| Fisher's P value | 0.039 | |||
| Pearson's P value | 0.039 |
Haplotype analysis of rs11805303 and rs10889677 SNPs.
Subgroup Analysis
CP is a highly heterogeneous condition that likely has multiple etiologies. It is reasonable to speculate that the broad clinical spectrum of CP can be at least partially attributed to considerable genetic heterogeneity. Therefore, we conducted subgroup analysis according to clinical phenotypes. There was a significant association of CP + GDD with rs10889657, rs1004819, rs11805303, and rs10889677, and the association of rs10889677 met the Bonferroni correction cut-off for multiple testing (OR = 1.293, 95% CI = 1.034–1.61, P = 0.024, Table 4).
Table 4
| Group | Allele frequency | P value | OR (95% CI) | Genotype frequency | P value | Group | Allele frequency | P value | OR (95% CI) |
|---|---|---|---|---|---|---|---|---|---|
| rs10889657 | C | T | C/C | C/T | T/T | ||||
| CP | 132 (0.253) | 390 (0.747) | 0.06 | 0.801 (0.635–1.009) | 23 (0.088) | 86 (0.330) | 152 (0.582) | 0.026a | 0.038 |
| Control | 372 (0.297) | 880 (0.703) | 53 (0.085) | 266 (0.425) | 307 (0.490) | 0.665 | |||
| rs6682925 | C | T | C/C | C/T | T/T | ||||
| CP | 239 (0.424) | 325 (0.576) | 0.578 | 0.946 (0.777–1.151) | 54 (0.191) | 131 (0.465) | 97 (0.344) | 0.608 | 0.412 |
| Control | 634 (0.437) | 827 (0.563) | 139 (0.189) | 365 (0.497) | 231 (0.314) | 0.807 | |||
| rs184444 | G | T | G/G | G/T | T/T | ||||
| CP | 196 (0.343) | 376 (0.657) | 0.279 | 0.894 (0.731–1.094) | 38 (0.133) | 120 (0.420) | 128 (0.448) | 0.201 | 0.246 |
| Control | 556 (0.368) | 954 (0.632) | 97 (0.128) | 362 (0.479) | 296 (0.392) | 0.402 | |||
| rs1004819 | A | G | A/A | A/G | G/G | ||||
| CP | 219 (0.592) | 151 (0.408) | 0.021b | 1.328 (1.044–1.688) | 71 (0.384) | 77 (0.416) | 37 (0.200) | 0.051 | 0.06 |
| Control | 544 (0.522) | 498 (0.478) | 150 (0.288) | 244 (0.468) | 127 (0.244) | 0.16 | |||
| rs11805303 | C | T | C/C | C/T | T/T | ||||
| CP | 232 (0.403) | 344 (0.597) | 0.019c | 0.793 (0.652–0.963) | 53 (0.184) | 126 (0.438) | 109 (0.378) | 0.051 | 0.124 |
| Control | 696 (0.460) | 818 (0.540) | 167 (0.221) | 362 (0.478) | 228 (0.301) | 0.304 | |||
| rs10889677 | A | C | A/A | A/C | C/C | ||||
| CP | 439 (0.767) | 133 (0.233) | 0.024d | 1.293 (1.034–1.619) | 174 (0.608) | 91 (0.318) | 21 (0.073) | 0.004e | 0.067 |
| Control | 1,079 (0.718) | 423 (0.282) | 378 (0.503) | 323 (0.430) | 50 (0.067) | 0.084 |
Allele and genotype frequencies of IL-23R in CP with GDD and the control group.
P value after Bonferroni
0.077,
0.126,
0.114,
0.144,
0.024.
Furthermore, we analyzed the associations between IL-23R and CP subtypes with various risk factors. There were significant differences in the rs11805303 allele frequency and the rs10889677 genotype frequency between controls and CP patients with NE; the association of rs10889677 with CP + NE remained significant even after Bonferroni correction (OR = 1.176, 95% CI = 0.947–1.460, P = 0.024, Table 5). There were no significant differences in either allele or genotype frequencies of the other SNPs in the other CP subgroups defined by clinical phenotypes and available risk factors (Supplementary Table 1).
Table 5
| Group | Allele frequency | P value | OR (95% CI) | Genotype frequency | P value | Group | Allele frequency | P value | OR (95% CI) |
|---|---|---|---|---|---|---|---|---|---|
| rs10889657 | C | T | C/C | C/T | T/T | ||||
| CP | 146 (0.271) | 392 (0.729) | 0.271 | 0.881 (0.703–1.104) | 25 (0.093) | 96 (0.357) | 148 (0.550) | 0.163 | 0.109 |
| Control | 372 (0.297) | 880 (0.703) | 53 (0.085) | 266 (0.425) | 307 (0.490) | 0.665 | |||
| rs6682925 | C | T | C/C | C/T | T/T | ||||
| CP | 257 (0.425) | 347 (0.575) | 0.619 | 0.952 (0.787–1.154) | 48 (0.159) | 161 (0.533) | 93 (0.308) | 0.434 | 0.116 |
| Control | 634 (0.437) | 827 (0.563) | 139 (0.189) | 365 (0.497) | 231 (0.314) | 0.807 | |||
| rs184444 | G | T | G/G | G/T | T/T | ||||
| CP | 219 (0.369) | 375 (0.631) | 0.984 | 1.002 (0.823–1.220) | 36 (0.121) | 140 (0.495) | 114 (0.384) | 0.891 | 0.276 |
| Control | 556 (0.368) | 954 (0.632) | 97 (0.128) | 362 (0.479) | 296 (0.392) | 0.402 | |||
| rs1004819 | A | G | A/A | A/G | G/G | ||||
| CP | 201 (0.552) | 163 (0.448) | 0.322 | 1.128 (0.888–1.434) | 62 (0.341) | 77 (0.423) | 43 (0.236) | 0.391 | 0.051 |
| Control | 544 (0.522) | 498 (0.478) | 150 (0.288) | 244 (0.468) | 127 (0.244) | 0.16 | |||
| rs11805303 | C | T | C/C | C/T | T/T | ||||
| CP | 244 (0.405) | 358 (0.595) | 0.023a | 0.801 (0.661–0.970) | 52 (0.173) | 140 (0.465) | 109 (0.362) | 0.083 | 0.541 |
| Control | 696 (0.460) | 818 (0.540) | 167 (0.221) | 362 (0.478) | 228 (0.301) | 0.304 | |||
| rs10889677 | A | C | A/A | A/C | C/C | ||||
| CP | 450 (0.750) | 150 (0.250) | 0.141 | 1.176 (0.947–1.460) | 177 (0.590) | 96 (0.320) | 27 (0.090) | 0.004b | 0.011 |
| Control | 1,079 (0.718) | 423 (0.282) | 378 (0.503) | 323 (0.430) | 50 (0.067) | 0.084 |
Allele and genotype frequencies of IL-23R in CP with NE and the control group.
P value after Bonferroni
0.138,
0.024.
Discussion
CP is a heterogeneous condition with multiple causes (Badawi and Keogh, 2013). The etiology in an individual patient is often multifactorial. These known CP causes, such as periventricular leukomalacia (PVL), NE, infarct, and premature delivery, account for only a minority of the total cases (Cowan et al., 2003; Yildiz et al., 2012; Colver et al., 2014; Chang et al., 2015). A single severe adverse event can be sufficient to cause CP, but much more often it is not a single cause, rather multiple concurrent risk factors that precede CP (Hankins and Speer, 2003; Djukic et al., 2009). Secondary neuroinflammation is associated with many CP risk factors. Findings from animal and clinical studies suggested that persistent neuroinflammation might prevent regeneration or exacerbate brain damage (Elovitz et al., 2011). Altered inflammation is one of the common causes of CP (Chiricozzi et al., 2011; Chen et al., 2013; Xia et al., 2018). Although the exact cause of CP is largely unknown, it is thought to be due to a combination of an altered fetal inflammatory response and primary brain damages.
Abnormal inflammation is one of the important pathological causes of immature brain injury (Djukic et al., 2009; Du et al., 2011), which is involved in the pathogenesis of central nervous system diseases, such as epilepsy, Parkinson's disease, cerebral ischemia, and hemorrhage, and may easily lead to secondary brain insult. During cerebral ischemia–reperfusion injury, inflammatory cells, such as microglia, astrocytes, and leukocytes, are activated. The activated inflammatory cells synthesize and secrete inflammatory mediators, phenomena that, in turn, further activate inflammatory cells and aggravate brain injury (Moon et al., 2009; Albertsson et al., 2018).
Some inflammatory cytokines are involved in immune responses, which cause brain injury through an inflammatory mechanism (McAdams and Juul, 2012). Several inflammatory cytokines have been identified as being involved in CP or other neurodevelopmental disorders, such as IL-6, IL-8, and IL-17 (Chiricozzi et al., 2011; Chen et al., 2013). Our previous studies found that rs1800795 (G174C), located in the IL-6 promoter region, was significantly associated with CP, and that the risk of spastic hemiplegia and quadriplegia in carriers of the rs1800795 C allele also significantly increased (Djukic et al., 2009; Wu et al., 2009). Central nervous system inflammation is often characterized by microglia activation, and active microglia will mediate neurotoxicity by secreting inflammatory cytokines, proteins, or other bioactive substances, resulting in secondary brain injury (Tang and Le, 2016).
A recent study demonstrated that inflammatory factors are related to microglia activation (Zhao et al., 2020). IL-23 mainly acts as a pro-inflammatory cytokine and has potential anti-tumor and anti-infection effects. IL-23R can activate Janus kinase (JAK). IL-23 bound IL-23R can active downstream JAKs and phosphorylate the signal transducer and activator of transcription (STAT) binding site in the intracellular region of the receptor. STAT dimerizes and is phosphorylated by JAKs. The phosphorylated STAT dimers enter the nucleus and act on downstream target genes. IL-23R has been identified to associate with multiple diseases, including alopecia areata (rs10889677) and nephropathy (rs10805303) (Figure 3) (Safrany et al., 2011; Yu et al., 2012; Huang et al., 2016; Qin et al., 2016; Poomarimuthu et al., 2018; Sode et al., 2018; Zhong et al., 2018; Kramer et al., 2019; Loures et al., 2019; Ruyssen-Witrand et al., 2019; Tabatabaei-Panah et al., 2020).
Figure 3
Our results showed that IL-23R rs10889677 increases susceptibility to CP at the overall level and in some subgroups. These findings suggest that IL-23R is a potential susceptibility gene for CP. Furthermore, studies have shown that rs10889677 is related with different diseases in different races, such as Crohn's disease in both Jewish and non-Jewish populations, ankylosing spondylitis in Caucasian patients, Graves' disease in North Americans, and ulcerative colitis in Europeans (Duerr et al., 2006; Wellcome Trust Case Control et al., 2007; Brown, 2008; Huber et al., 2008; Silverberg et al., 2009). Therefore, we have sufficient reason to assume the positive association of IL-23R gene with the CP etiology. Whether IL-23R, as a gene associated with inflammatory bowel disease, can also affect the brain development of children by affecting intestinal flora remains to be seen.
Given that TNF-α and IL-6 are significantly increased after brain injury (Leviton and Dammann, 2004; Xie et al., 2014), we hypothesize that IL-23R, like other inflammatory factors TNF-α and IL-6, may cause brain damage and lead to CP through the following steps. (1) Increased inflammatory cytokines promote the release of nitric oxide synthase and free radicals and excitatory amino acids, which have toxic effects on neurons, especially developing brain tissue. (2) They will trigger a systemic inflammatory response that leads to brain damage when undergoing intrauterine infection. (3) Endothelial cell damage can cause thrombosis. Inflammatory factors will activate platelets, lead to their aggregation, activate coagulation factors, and damage white matter neurons. (4) The damage will increase the permeability of the blood–brain barrier, thereby allowing peripheral bacteria and inflammatory factors to enter the brain, aggravating brain damage. (5) They promote the release of prostaglandins and other substances, resulting in pregnancy's advance labor, leading to an increased risk of CP.
Our study has some limitations. First, this is a study based on a single gene for susceptibility to CP. Given the genetic heterogeneity and gene–gene interaction involved in the CP etiology, other candidate genes that are part of the IL-23R signaling pathway need to be analyzed together. Second, we were unable to measure IL-23R protein expression in the brains of the subjects in the current study; future studies are encouraged to examine the inflammatory cytokine alteration in the brain. Third, although our study demonstrated an association between the IL-23R rs10889677 SNP and CP, further functional and replicated studies are necessary to verify the association of IL-23R with CP, which is of great significance to identify the CP etiology and pathogenesis.
In summary, a significant association between IL-23R and CP was firstly detected in Han Chinese, suggesting that the IL-23R gene has a significant effect on the risk of CP, especially in subjects with GDD or NE. The inflammatory response and cytokine cascade are likely to play a role in the occurrence and development of CP. This result needs to be further validated with well-designed studies with large sample sizes and in other populations. We should also pay attention to the possibility of increased risk of CP if the fetus is found to carry the IL-23R risk genotypes before or after delivery.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Zhengzhou University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
QX and CZ conceived and designed the study. YX, JS, LX, QS, CG, and DZ recruited subjects and sorted out clinical information. YW and DB performed all of the laboratory work. YW and YF performed all data and statistical analyses. YF drafted the manuscript, and QX, CZ, XW, YW, and YX revised the manuscript critically for important intellectual content. YQ and JS provided data, developed models, reviewed results, and provided guidance on the methods. All authors contributed and critically reviewed the final version of the manuscript. All authors have read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Plan for Stem Cell and Transformation Research (2017YFA0104202), the Department of Science and Technology of Henan Province, China (171100310200), the Health Department of Henan Province (SB201901055), the Shanghai Municipal Commission of Science and Technology Research Project (19JC1411000), the Capacity-Building Planning Program for Shanghai Women and Children's Health Service, the collaborative innovation center project construction for Shanghai women and children's health, the Fourth Round of the Shanghai Three-year Action Plan on Public Health Discipline and Talent Program: Women and Children's Health (No. 15GWZK0401), the Swedish Research Council (2015-04780, 2018-02667), the Swedish Governmental grants to scientists working in health care (ALFGBG-717791), the Swedish Brain Foundation (FO2018-0034), the VINNMER-Marie Curie (VINNOVA, 2015-04780), and the National Natural Science Foundation of China (81801499).
Acknowledgments
We thank all participants involved in this study and all authors.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewer HB-P declared a past co-authorship with the authors XW and CZ to the handling editor.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2020.590098/full#supplementary-material
References
1
AlbertssonA. M.BiD.DuanL.ZhangX.LeavenworthJ. W.QiaoL.et al. (2014). The immune response after hypoxia-ischemia in a mouse model of preterm brain injury. J. Neuroinflammation11:153. 10.1186/s12974-014-0153-z
2
AlbertssonA. M.ZhangX.VontellR.BiD.BronsonR. T.SupramaniamV.et al. (2018). gammadelta T Cells Contribute to injury in the developing brain. Am. J. Pathol.188, 757–767. 10.1016/j.ajpath.2017.11.012
3
BadawiN.KeoghJ. M. (2013). Causal pathways in cerebral palsy. J. Paediatr. Child Health49, 5–8. 10.1111/jpc.12068
4
BeliakoffJ.SunZ. (2006). Zimp7 and Zimp10, two novel PIAS-like proteins, function as androgen receptor coregulators. Nucl. Recept. Signal4:e017. 10.1621/nrs.04017
5
BlochY.BouchareychasL.MerceronR.SkladanowskaK.Van den BosscheL.DetryS.et al. (2018). Structural activation of pro-inflammatory human Cytokine IL-23 by Cognate IL-23 receptor enables recruitment of the shared receptor IL-12Rbeta1. Immunity48, 45–58 e46. 10.1016/j.immuni.2017.12.008
6
BrownM. A. (2008). Breakthroughs in genetic studies of ankylosing spondylitis. Rheumatology47, 132–137. 10.1093/rheumatology/kem269
7
ChangM. J.MaH. I.LuT. H. (2015). Estimating the prevalence of cerebral palsy in Taiwan: a comparison of different case definitions. Res. Dev. Disabil.36C, 207–212. 10.1016/j.ridd.2014.10.001
8
ChenM.LiT.LinS.BiD.ZhuD.ShangQ.et al. (2013). Association of interleukin 6 gene polymorphisms with genetic susceptibilities to spastic tetraplegia in males: a case-control study. Cytokine61, 826–830. 10.1016/j.cyto.2013.01.011
9
ChengX.LiT.WangH.ZhuD.MaC.MaB.et al. (2011). Methylenetetrahydrofolate reductase gene polymorphisms and cerebral palsy in Chinese infants. J. Hum. Genet.56, 17–21. 10.1038/jhg.2010.127
10
ChiricozziA.Guttman-YasskyE.Suarez-FarinasM.NogralesK. E.TianS.CardinaleI.et al. (2011). Integrative responses to IL-17 and TNF-alpha in human keratinocytes account for key inflammatory pathogenic circuits in psoriasis. J. Invest. Dermatol.131, 677–687. 10.1038/jid.2010.340
11
ColverA.FairhurstC.PharoahP. O. (2014). Cerebral palsy. Lancet383, 1240–1249. 10.1016/S0140-6736(13)61835-8
12
CowanF.RutherfordM.GroenendaalF.EkenP.MercuriE.BydderG. M.et al. (2003). Origin and timing of brain lesions in term infants with neonatal encephalopathy. Lancet361, 736–742. 10.1016/S0140-6736(03)12658-X
13
DjukicM.GibsonC. S.MaclennanA. H.GoldwaterP. N.HaanE. A.McMichaelG.et al. (2009). Genetic susceptibility to viral exposure may increase the risk of cerebral palsy. Aust. N. Z. J. Obstet. Gynaecol.49, 247–253. 10.1111/j.1479-828X.2009.00999.x
14
DuX.FleissB.LiH.D'AngeloB.SunY.ZhuC.et al. (2011). Systemic stimulation of TLR2 impairs neonatal mouse brain development. PLoS ONE6:e19583. 10.1371/journal.pone.0019583
15
DuerrR. H.TaylorK. D.BrantS. R.RiouxJ. D.SilverbergM. S.DalyM. J.et al. (2006). A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science314, 1461–1463. 10.1126/science.1135245
16
DuvalD.DuvalG.KedingerC.PochO.BoeufH. (2003). The ‘PINIT' motif, of a newly identified conserved domain of the PIAS protein family, is essential for nuclear retention of PIAS3L. FEBS Lett.554, 111–118. 10.1016/S0014-5793(03)01116-5
17
ElovitzM. A.BrownA. G.BreenK.AntonL.MaubertM.BurdI. (2011). Intrauterine inflammation, insufficient to induce parturition, still evokes fetal and neonatal brain injury. Int. J. Dev. Neurosci.29, 663–671. 10.1016/j.ijdevneu.2011.02.011
18
HankinsG. D.SpeerM. (2003). Defining the pathogenesis and pathophysiology of neonatal encephalopathy and cerebral palsy. Obstet. Gynecol.102, 628–636. 10.1016/S0029-7844(03)00574-X
19
HuangJ.YangY.ZhouF.LiangZ.KangM.KuangY.et al. (2016). Meta-analysis of the IL23R and IL12B polymorphisms in multiple sclerosis. Int. J. Neurosci.126, 205–212. 10.3109/00207454.2015.1007508
20
HuberA. K.JacobsonE. M.JazdzewskiK.ConcepcionE. S.TomerY. (2008). Interleukin (IL)-23 receptor is a major susceptibility gene for Graves' ophthalmopathy: the IL-23/T-helper 17 axis extends to thyroid autoimmunity. J. Clin. Endocrinol. Metab.93, 1077–1081. 10.1210/jc.2007-2190
21
JacobssonB.HagbergG. (2004). Antenatal risk factors for cerebral palsy. Best Pract. Res. Clin. Obstet. Gynaecol.18, 425–436. 10.1016/j.bpobgyn.2004.02.011
22
KramerM.HasanreisogluM.WeissS.KumovaD.Schaap-FoglerM.Guntekin-ErgunS.et al. (2019). Single-nucleotide polymorphisms in IL23R-IL12RB2 (rs1495965) are highly prevalent in patients with Behcet's Uveitis and vary between populations. Ocul. Immunol. Inflamm.27, 766–773. 10.1080/09273948.2018.1467463
23
LevitonA.DammannO. (2004). Coagulation, inflammation, and the risk of neonatal white matter damage. Pediatr. Res.55, 541–545. 10.1203/01.PDR.0000121197.24154.82
24
LiF. F.ZhuX. D.YanP.JinM. H.YueH.ZhangQ.et al. (2016). Characterization of variations in IL23A and IL23R genes: possible roles in multiple sclerosis and other neuroinflammatory demyelinating diseases. Aging8, 2734–2746. 10.18632/aging.101058
25
LouresM. A. R.AlvesH. V.de MoraesA. G.SantosT. D. S.LaraF. F.NevesJ. S. F.et al. (2019). Association of TNF, IL12, and IL23 gene polymorphisms and psoriatic arthritis: meta-analysis. Expert Rev. Clin. Immunol15, 303–313. 10.1080/1744666X.2019.1564039
26
MagalhaesR. C.MoreiraJ. M.LauarA. O.da SilvaA. A. S.TeixeiraA. L.. (2019). Inflammatory biomarkers in children with cerebral palsy: a systematic review. Res. Dev. Disabil.95:103508. 10.1016/j.ridd.2019.103508
27
MarshallH. S.PlotkinS. (2019). The changing epidemiology of mumps in a high vaccination era. Lancet Infect. Dis.19, 118–119. 10.1016/S1473-3099(18)30541-3
28
McAdamsR. M.JuulS. E. (2012). The role of cytokines and inflammatory cells in perinatal brain injury. Neurol. Res. Int.2012:561494. 10.1155/2012/561494
29
MoonJ. B.LeeC. H.ParkC. W.ChoJ. H.HwangI. K.YooK. Y.et al. (2009). Neuronal degeneration and microglial activation in the ischemic dentate gyrus of the gerbil. J. Vet. Med. Sci.71, 1381–1386. 10.1292/jvms.001381
30
OppmannB.LesleyR.BlomB.TimansJ. C.XuY.HunteB.et al. (2000). Novel p19 protein engages IL-12p40 to form a cytokine, IL-23, with biological activities similar as well as distinct from IL-12. Immunity13, 715–725. 10.1016/S1074-7613(00)00070-4
31
PenningtonL.VirellaD.MjoenT.da Graca AndradaM.MurrayJ.ColverA.et al. (2013). Development of The Viking Speech Scale to classify the speech of children with cerebral palsy. Res. Dev. Disabil.34, 3202–3210. 10.1016/j.ridd.2013.06.035
32
PoomarimuthuM.ElangoS.SolomonP. R.SoundrapandianS.MariakuttikanJ. (2018). Association of IL17 and IL23R gene polymorphisms with rheumatic heart disease in South Indian population. Immunol. Invest.47, 754–764. 10.1080/08820139.2018.1493053
33
QinX.XuJ.WuZ.SunF.ChenH.ZhengW.et al. (2016). Association study of rs924080 and rs11209032 polymorphisms of IL23R-IL12RB2 in a Northern Chinese Han population with Behcet's disease. Hum. Immunol.77, 1284–1290. 10.1016/j.humimm.2016.09.006
34
RosenbaumP.PanethN.LevitonA.GoldsteinM.BaxM.DamianoD.et al. (2007). A report: the definition and classification of cerebral palsy April 2006. Dev. Med. Child Neurol. Suppl.109, 8–14. 10.1111/j.1469-8749.2007.tb12610.x
35
Ruyssen-WitrandA.LuxembourgerC.CantagrelA.NigonD.ClaudepierreP.DegboeY.et al. (2019). Association between IL23R and ERAP1 polymorphisms and sacroiliac or spinal MRI inflammation in spondyloarthritis: DESIR cohort data. Arthritis Res. Ther.21:22. 10.1186/s13075-018-1807-5
36
SachdevS.BruhnL.SieberH.PichlerA.MelchiorF.GrosschedlR. (2001). PIASy, a nuclear matrix-associated SUMO E3 ligase, represses LEF1 activity by sequestration into nuclear bodies. Genes Dev.15, 3088–3103. 10.1101/gad.944801
37
SafranyE.SzellM.CsongeiV.JaromiL.SipekyC.SzaboT.et al. (2011). Polymorphisms of the IL23R gene are associated with psoriasis but not with immunoglobulin A nephropathy in a Hungarian population. Inflammation34, 603–608. 10.1007/s10753-010-9268-2
38
SchonM. P.ErpenbeckL. (2018). The Interleukin-23/Interleukin-17 axis links adaptive and innate immunity in psoriasis. Front. Immunol.9:1323. 10.3389/fimmu.2018.01323
39
SilverbergM. S.ChoJ. H.RiouxJ. D.McGovernD. P.WuJ.AnneseV.et al. (2009). Ulcerative colitis-risk loci on chromosomes 1p36 and 12q15 found by genome-wide association study. Nat. Genet.41, 216–220. 10.1038/ng.275
40
SodeJ.BankS.VogelU.AndersenP. S.SorensenS. B.BojesenA. B.et al. (2018). Genetically determined high activities of the TNF-alpha, IL23/IL17, and NFkB pathways were associated with increased risk of ankylosing spondylitis. BMC Med. Genet.19:165. 10.1186/s12881-018-0680-z
41
StrleK.ZhouJ. H.ShenW. H.BroussardS. R.JohnsonR. W.FreundG. G.et al. (2001). Interleukin-10 in the brain. Crit. Rev. Immunol.21, 427–449. 10.1615/CritRevImmunol.v21.i5.20
42
SunL.XiaL.WangM.ZhuD.WangY.BiD.et al. (2019). Variants of the OLIG2 gene are associated with cerebral palsy in Chinese Han infants with hypoxic-ischemic encephalopathy. Neuromolecular. Med.21, 75–84. 10.1007/s12017-018-8510-1
43
Tabatabaei-PanahP. S.MoravvejH.DelpasandS.JafariM.SepehriS.AbgoonR.et al. (2020). IL12B and IL23R polymorphisms are associated with alopecia areata. Genes Immun.21, 203–210. 10.1038/s41435-020-0100-1
44
TangY.LeW. (2016). Differential roles of M1 and M2 microglia in neurodegenerative diseases. Mol. Neurobiol.53, 1181–1194. 10.1007/s12035-014-9070-5
45
TatlaS. K.SauveK.Virji-BabulN.HolstiL.ButlerC.Van Der LoosH. F. (2013). Evidence for outcomes of motivational rehabilitation interventions for children and adolescents with cerebral palsy: an American Academy for Cerebral Palsy and Developmental Medicine systematic review. Dev. Med. Child Neurol.55, 593–601. 10.1111/dmcn.12147
46
TatorC.BrayG.MorinD. (2007). The CBANCH report–the burden of neurological diseases, disorders, and injuries in Canada. Can. J. Neurol. Sci.34, 268–269. 10.1017/S0317167100006673
47
Wellcome Trust Case ControlC.Australo-Anglo-American SpondylitisC.BurtonP. R.ClaytonD. G.CardonL. R.CraddockN.et al. (2007). Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat. Genet.39, 1329–1337. 10.1038/ng.2007.17
48
WiekowskiM. T.LeachM. W.EvansE. W.SullivanL.ChenS. C.VassilevaG.et al. (2001). Ubiquitous transgenic expression of the IL-23 subunit p19 induces multiorgan inflammation, runting, infertility, and premature death. J. Immunol.166, 7563–7570. 10.4049/jimmunol.166.12.7563
49
WuY. W.CroenL. A.TorresA. R.Van De WaterJ.GretherJ. K.HsuN. N. (2009). Interleukin-6 genotype and risk for cerebral palsy in term and near-term infants. Ann. Neurol.66, 663–670. 10.1002/ana.21766
50
WuY. W.CroenL. A.VanderwerfA.GelfandA. A.TorresA. R. (2011). Candidate genes and risk for CP: a population-based study. Pediatr. Res.70, 642–646. 10.1203/PDR.0b013e31823240dd
51
XiaL.ChenM.BiD.SongJ.ZhangX.WangY.et al. (2018). Combined analysis of Interleukin-10 gene polymorphisms and protein expression in children with cerebral palsy. Front. Neurol.9:182. 10.3389/fneur.2018.00182
52
XiaL.XuJ.SongJ.XuY.ZhangB.GaoC.et al. (2019). Autophagy-related gene 7 polymorphisms and cerebral palsy in Chinese infants. Front. Cell Neurosci.13:494. 10.3389/fncel.2019.00494
53
XieC.ZhouK.WangX.BlomgrenK.ZhuC. (2014). Therapeutic benefits of delayed lithium administration in the neonatal rat after cerebral hypoxia-ischemia. PLoS ONE9:e107192. 10.1371/journal.pone.0107192
54
XuJ.XiaL.ShangQ.DuJ.ZhuD.WangY.et al. (2017). A variant of the autophagy-related 5 gene is associated with child cerebral palsy. Front. Cell Neurosci.11:407. 10.3389/fncel.2017.00407
55
YildizE. P.TatliB.EkiciB.EraslanE.AydinliN.CaliskanM.et al. (2012). Evaluation of etiologic and prognostic factors in neonatal convulsions. Pediatr. Neurol.47, 186–192. 10.1016/j.pediatrneurol.2012.05.015
56
YuP.ShenF.ZhangX.CaoR.ZhaoX.LiuP.et al. (2012). Association of single nucleotide polymorphisms of IL23R and IL17 with ulcerative colitis risk in a Chinese Han population. PLoS ONE7:e44380. 10.1371/journal.pone.0044380
57
YuT.XiaL.BiD.WangY.ShangQ.ZhuD.et al. (2018). Association of NOS1 gene polymorphisms with cerebral palsy in a Han Chinese population: a case-control study. BMC Med Genomics11:56. 10.1186/s12920-018-0374-6
58
YunusA. A.LimaC. D. (2009). Structure of the Siz/PIAS SUMO E3 ligase Siz1 and determinants required for SUMO modification of PCNA. Mol. Cell35, 669–682. 10.1016/j.molcel.2009.07.013
59
ZhaoY.GanY.XuG.YinG.LiuD. (2020). MSCs-derived exosomes attenuate acute brain injury and inhibit microglial inflammation by reversing CysLT2R-ERK1/2 mediated microglia M1 polarization. Neurochem. Res.45, 1180–1190. 10.1007/s11064-020-02998-0
60
ZhongL.WangW.SongH. (2018). Complex role of IL-23R polymorphisms on ankylosing spondylitis: a meta-analysis. Expert Rev. Clin. Immunol.14, 635–643. 10.1080/1744666X.2018.1491308
Summary
Keywords
cerebral palsy, inflammatory cytokines, interleukin, gene polymorphism, IL23R
Citation
Wang Y, Xu Y, Fan Y, Bi D, Song J, Xia L, Shang Q, Gao C, Zhang X, Zhu D, Qiao Y, Su Y, Wang X, Zhu C and Xing Q (2020) The Association Study of IL-23R Polymorphisms With Cerebral Palsy in Chinese Population. Front. Neurosci. 14:590098. doi: 10.3389/fnins.2020.590098
Received
31 July 2020
Accepted
14 October 2020
Published
25 November 2020
Volume
14 - 2020
Edited by
Patricia Pelufo Silveira, McGill University, Canada
Reviewed by
Hilla Ben-Pazi, Shaare Zedek Medical Center, Israel; Claudia Vianna Maurer-Morelli, State University of Campinas, Brazil
Updates
Copyright
© 2020 Wang, Xu, Fan, Bi, Song, Xia, Shang, Gao, Zhang, Zhu, Qiao, Su, Wang, Zhu and Xing.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Changlian Zhu zhuc@zzu.edu.cn; Changlian.zhu@neuro.gu.seQinghe Xing xingqinghe@hotmail.com
This article was submitted to Neurogenomics, a section of the journal Frontiers in Neuroscience
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.