Objective
The purpose of this study was to evaluate the association between fetal inflammation and coagulation gene single-nucleotide polymorphisms (SNPs) and neurodevelopmental delay at age 2 years.
Study Design
We conducted a case-controlled secondary analysis of a randomized trial of single- vs multiple-course corticosteroids. Multiplex assay assessed 46 SNPs. Cases had mental developmental and/or psychomotor delay at age 2 years. Control subjects had normal neurodevelopment.
Results
One hundred twenty-five cases and 147 control subjects were analyzed. Allele frequencies were different between cases and control subjects for interleukin (IL)1β-511 ( P = .009), IL4R-148 ( P = .03), IL6-174 ( P = .02), and IL6-176 ( P = .007). Genotype frequencies were different for IL1β-511 ( P = .03) and IL6-174 ( P = .04). Results for IL1β-511, IL4R-148, and IL6-176 remained significant after logistic regression analysis. IL1β-511 and IL6-176 minor alleles were associated with increased risk of neurodevelopmental delay (odds ratio, 3.1; 95% confidence interval [CI], 1.2–8.2 and 2.2; 95% CI, 1.2–3.9, respectively). IL4R-148 minor allele was protective (odds ratio, 0.6; 95% CI, 0.4–0.9).
Conclusion
Fetal SNPs in IL1β, IL-4R, and IL-6 may be associated with neurodevelopmental delay at age 2 years.
There is accumulating evidence to support the hypothesis that neurodevelopmental outcomes after preterm birth are influenced by both genetic and environmental factors. The fetal inflammatory response syndrome is a multisystem disorder that has been associated with preterm birth, cerebral white-matter damage (periventricular leukomalacia), and cerebral palsy (CP). Fetal proinflammatory cytokines, interleukin (IL) 6, IL1-β, and tumor necrosis factor-α have been associated with fetal/neonatal central nervous system injury. Further, polymorphisms in genes that are related to cytokine function (tumor necrosis factor-α, mannose-binding lectin, IL8, and lymphotoxin-α) have been associated with an increased risk of CP in recent studies.
Perinatal thromboembolic events are also recognized increasingly as a cause of neurologic disabilities such as CP, epilepsy, and other cognitive abnormalities. Polymorphisms in genes that are related to thrombosis or thrombolysis (factor V Leiden, prothrombin gene mutation G20210A, factor VII, plasminogen activator inhibitor I, and methylenetetrahydrofolate reductase) have been associated with an increased risk of CP. The interaction between the inflammation and coagulation pathways may be important to the mechanism of fetal central nervous system injury. Perinatal arterial ischemic stroke (a common cause of hemiplegic CP) has been linked to several independent maternal risk factors, which include chorioamnionitis. There is compelling evidence that activated coagulation factors contribute to neonatal white-matter damage through enhanced inflammation, rather than occlusion of cerebral blood vessels.
Although previous studies suggest that certain single nucleotide polymorphisms (SNPs) in inflammation and coagulation pathway genes increase vulnerability to CP, it is uncertain whether SNPs in 1 or both of these pathways are associated with other forms of neurodevelopmental delay. We hypothesized that children with neurodevelopmental delay at age 2 years, when compared with children with normal neurodevelopment, have differences in inflammation and/or coagulation pathways gene SNPs.
Materials and Methods
Subjects
The children who were evaluated were enrolled in the National Institute of Child Health and Human Development Maternal-Fetal Medicine Units (MFMU) Network, which is a randomized, double-masked, placebo-controlled, multicenter clinical trial of single- vs multiple-course antenatal corticosteroids. The primary aim of this trial was to assess the clinical efficacy and safety of repeated doses of antenatal corticosteroids in pregnancies that were at risk for preterm birth. The details of the trial, which was conducted between 2000 and 2003, have been reported previously.
A secondary aim of the MFMU Network antenatal corticosteroid trial was the correlation of steroid regimen with neurodevelopmental outcomes as assessed by the Bayley Scales of Infant Development, second edition, at age 2 years. The Bayley Scales consist of mental and psychomotor developmental indices (MDI and PDI). A Bayley score of 85 is 1 SD below the mean and consistent with mild neurodevelopmental delay. A Bayley score of 70 is 2 SD below the mean and consistent with more significant delay. Placental samples and fetal cord serum were collected in a subset of subjects who were enrolled in the MFMU Network antenatal corticosteroid trial per study protocol.
Our secondary analysis aimed to evaluate the association between SNPs in fetal inflammation and coagulation pathway genes and neurodevelopmental delay at age 2 years. Inclusion criteria consisted of (1) maternal enrollment in the MFMU Network antenatal corticosteroid trial, (2) Bayley scores at age 2 years, and (3) DNA extracted from placental samples or fetal cord serum. Cases were children with mental and/or psychomotor delay that had been defined by a Bayley MDI and/or PDI score of <85 at age 2 years. CP cases (as defined in the primary trial) were excluded from the analysis to test our hypothesis of an association between inflammation and coagulation gene SNPs and neurodevelopmental delay at age 2 years in the absence of CP. Control subjects were children with normal neurodevelopment that had been defined by Bayley MDI and PDI score of ≥85.
The study group is illustrated in the Figure . From the original steroid trial, 583 children were discharged alive. Of these, 4 children died between discharge and follow-up evaluation, and 23 children were not available because they were born at a center that discontinued participation in the MFMU Network in 2001. Of the remaining 556 children who were available for follow-up evaluation, 459 children (83%) had Bayley scores at age 2-3 years. After exclusion of CP cases (n = 7), 452 children were eligible for participation in our analysis. There were 40 twin pairs in the original steroid trial. For this study, 1 twin from each pair (n = 40) was excluded randomly to avoid the issue of including related individuals in the analysis. Of the remaining 412 infants, DNA from placental tissues and/or fetal serum was available for 272 infants (125 cases and 147 control subjects).
After being reviewed, this study was determined to be exempt from institutional review board approval procedures because of the deidentification of data and study samples before this analysis.
DNA Extraction
DNA was extracted from placental samples and/or umbilical cord serum that had been collected at the time of delivery. Placental samples were either fresh (n = 145) or formalin-fixed and paraffin-embedded (n = 92). DNA was extracted from approximately 0.7 g of fresh placental tissue with the PureGene DNA Purification System (Qiagen, Valencia, CA), as per the manufacturer’s protocols. DNA from the formalin-fixed paraffin-embedded samples was extracted from 3 to 6 5-μm sections that were obtained from the paraffin block and isolated with the PureGene DNA Purification System, as per the manufacturer’s protocols. DNA was also extracted from umbilical cord serum samples (n = 190). Briefly, 250-1000 μL of cord serum was centrifuged at 10,000 g for 10 minutes to pellet any cells in the sera. The supernatant was discarded, and DNA was extracted with the same procedure used for placental samples.
Selection of SNPs for genotyping
We selected 46 SNPs in 27 inflammation and coagulation genes based on the available literature and hypothesized causal pathways. We also included SNPs that have been reported to be associated with CP. The SNPs that were included in our custom multiplex assay are shown in Table 1 .
Gene | Symbol | Polymorphism | RS no. |
---|---|---|---|
Lymphotoxin-α (tumor necrosis factor-β) | LTA | thr26asn(C/A) | rs1041981 a |
Mannose-binding lectin | MBL2 | –550(G/C) | rs11003125 |
gly54asp(A/G) | rs1800450 | ||
gly57glu(A/G) | rs1800451 a | ||
arg52cys(T/C) | rs5030737 a | ||
–221(G/C) | rs7096206 | ||
Tumor necrosis factor-α | TNF | –308(G/A) | rs1800629 |
–376(G/A) | rs1800750 a | ||
–238(G/A) | rs361525 | ||
Tumor necrosis factor receptor superfamily member 6 | TNFRSF6 | –670(G/A) | rs1800682 |
Interleukin 1-α | IL1A | –889(C/T) | rs1800587 |
Interleukin 1-β | IL1B | +3962(C/T) | rs12621220 |
–511(A/G) | rs16944 | ||
Interleukin 1 type 1 receptor | IL1RN | +970(C/T) | rs3917365 |
Interleukin 4 | IL4 | –590(C/T) | rs2234648 |
582(A/G) | rs35648164 | ||
Interleukin 4 receptor | IL4R | gln576arg(A/G) | rs1801275 |
ile50val(A/G) | rs1805010 | ||
Interleukin 5 receptor-α | IL5RA | –80(G/A) | rs2290608 |
Interleukin 6 | IL6 | –174(G/C) | rs1800795 |
–176(G/C) | rs2234683 | ||
Interleukin 9 | IL9 | 553(A/C) | rs2066760 |
thr113met(T/C) | rs2069885 a | ||
Interleukin 10 | IL10 | –819(C/T) | rs1800871 a |
–1082(G/A) | rs1800896 | ||
Interleukin 13 | IL13 | arg130gln(T/C) | rs20541 |
Factor II | F2 | 20210(G/A) | rs1799963 |
Factor V Leiden | F5 | arg506gln(A/G) | rs6025 |
arg353gln(T/C) | rs6046 a | ||
Factor XIII | F13 | val34leu(T/G) | rs5985 |
Methylenetetrahydrofolate reductase | MTHFR | 1298(A/C) | rs1801131 |
677(C/T) | rs1801133 | ||
Thrombomodulin | THBD | 1418(C/T) | rs1042579 a |
Plasminogen activator inhibitor 1 | PAI1 | +11053(G/T) | rs7242 a |
Fibrinogen-β polypeptide chain | FGB | –455(G/A) | rs1800790 a |
Thrombin activatable fibrinolysis inhibitor | TAFI | 1040(C/T) | rs1926447 |
Platelet glycoprotein Ia (integrin, α2) | ITGA2 | 873(G/A) | rs1062535 a |
β-adrenergic receptor | ADRB2 | arg16gly(A/G) | rs1042713 |
gln27glu(C/G) | rs1042714 | ||
Nitric oxide synthase 3, endothelial | NOS3 | –690(C/T) | rs3918226 |
glu298asp(T/G) | rs1799983 | ||
–922(A/G) | rs1800779 | ||
Selectin P | SELP | val640leu(T/G) | rs6133 |
thr756pro(C/A) | rs6136 | ||
Transforming growth factor-β 1 | TGFB1 | –800(G/A) | rs1800468 |
–509(C/T) | rs1800469 |
a Polymerase chain reaction amplification was not successful.
Genotyping
The 48-Plex GenomeLab SNPStream Genotyping System and accompanying automated genotype calling software (Beckman Coulter, Fullerton, CA) were used for genotyping. Polymerase chain reaction (PCR) assays and extension primers for these SNPs were designed with the use of Beckman Coulter’s Autoprimer multiplex primer design engine ( www.Autoprimer.com ). PCR and extension reactions were performed according to the manufacturer’s instructions. Control subjects were included appropriately in each 384-well plate. Researchers and laboratory personnel were blinded to the case/control status of the biologic samples.
When DNA was available from both placenta and serum, the placental samples were used preferentially secondary to higher quality DNA. In subjects with a placental tissue genotyping success rate of <85%, the serum samples were also genotyped with a goal of increasing the number of SNP results. In 25 subjects with both placental tissue and sera genotyped, manufacturer’s automated allele calling resulted in at least 1 discordant SNP between the 2 samples. Comparison of data quality allowed a correct genotype to be easily identified in most cases. If the subject’s correct genotype could not be identified easily, no genotype was reported for that SNP. In all but 1 case, the serum data were questionable because of the low intensity of the genotype signal caused by the poor quality and low concentration of the extracted DNA. In 3 subjects, the data for both serum and placental samples was strong, but many SNPs were discordant. These 3 subjects were removed from the analysis. Because the serum samples tended to have lower DNA quality, the automated genotype results were individually reviewed. If data quality did not allow for a clear genotype determination for a given serum sample, no genotype was reported.
Statistics
Demographic and clinical characteristics of cases and control subjects were compared with the use of the χ 2 or Fisher’s exact test for categoric variables and the Wilcoxon Rank Sum test for continuous variables. Differences in allele and genotype frequencies between cases and control subjects were tested for each SNP with the use of the χ 2 or Fisher’s exact tests, as appropriate. SNPs that were found to be significantly related to an MDI and/or PDI score of <85 in the univariate analyses ( P < .05) were analyzed further by multiple logistic regression modeling. Genotypes were included as covariates in a series of regression models that assumed an additive, dominant, or recessive genetic pattern. The additive model assumes that having 2 copies of minor allele has twice the effect of having 1 copy; the dominant model assumes that having at least 1 copy of the minor allele is sufficient for disease, and the recessive model assumes that 2 copies of the minor allele are needed for disease. Other covariates included gestational age at delivery, small-for-gestational-age status (based on a birthweight <10th percentile of published standards), sex, steroid treatment group, race, smoking, and maternal level of education. Of the 3 ways of examining genotype, the regression model that had the highest likelihood score was considered to be the best-fitting model for the respective SNP. The sample size and the rarity of some genotypes and phenotypes limited further modeling of interactions.
Tests for Hardy-Weinberg equilibrium (χ 2 or Fisher’s exact) were performed for control subjects for each SNP that was significant in the univariate analysis. Because this was an exploratory study, no adjustments were made for multiple comparisons, and all comparisons are reported. A probability value of < .05 was considered to be statistically significant. All calculations were performed with SAS statistical software (SAS Institute, Inc, Cary, NC).
Results
The demographic and clinical characteristics of those children who were included vs excluded from our study based on the availability of biologic samples for DNA analysis are shown in Table 2 . Subjects who were excluded from our study were born at an earlier gestational age ( P = .01), had a lower birthweight ( P = .01), and were more likely to have received multiple courses of antenatal corticosteroids ( P = .007). The demographic and clinical characteristics of our study subjects are shown in Table 3 . Significant differences between cases and control subjects were observed for race/ethnicity ( P < .001) and maternal education ( P = .01).
Characteristic | Children who were included | Children who were excluded a | P value |
---|---|---|---|
Subjects, n | 272 | 180 | |
Sex, n (%) | .51 | ||
Male | 143 (52.6) | 89 (49.4) | |
Female | 129 (47.4) | 91 (50.6) | |
Race/ethnicity, n (%) b | .59 | ||
White | 92 (33.8) | 37 (35.2) | |
Black | 104 (38.2) | 44 (41.9) | |
Other | 76 (27.9) | 24 (22.9) | |
Gestational age at delivery, wk b , c | 35.2 ± 3.8 | 34.2 ± 3.7 | .01 |
Gestational age at delivery, n (%) b | .03 | ||
≥37 wk | 108 (39.7) | 27 (25.7) | |
32-36 wk | 111 (40.8) | 49 (46.7) | |
<32 wk | 53 (19.5) | 29 (27.6) | |
Birthweight, g c | 2396.3 ± 778.7 | 2059.4 ± 721.4 | .01 |
Neurodevelopment, n (%) | |||
Mental and/or psychomotor developmental indices <85 | 125 (46.0) | 86 (47.8) | .94 |
Mental and/or psychomotor developmental indices <70 | 50 (18.4) | 37 (20.6) | .86 |
Chorioamnionitis, n (%) b | 8 (2.9) | 3 (2.9) | 1.00 |
Maternal tobacco use, n (%) b | 55 (20.2) | 25 (23.8) | .44 |
Antenatal multiple-course steroids, n (%) b | 126 (46.3) | 65 (61.9) | .007 |
Maternal education, y b , c | 12.3 ± 2.7 | 12.3 ± 2.2 | .55 |