Introduction
Autism spectrum disorder (ASD) affects approximately 7.6 per 1 000 children worldwide (Baxter et al., 2015; Elsabbagh et al., 2012; Fombonne et al., 2016; Siu et al., 2016). The core symptoms of ASD are deficits in social communication and restricted interests and stereotyped and repetitive moments (Christiansz, Gray, Taffe, & Tonge, 2016; Kent et al., 2013; Lee, Thomas, & Lee, 2015). Despite advances in the investigation of ASD, there is still no clear evidence concerning the etiology of these disorders (Muhle, Reed, Stratigos, & Veenstra-VanderWeele, 2018). Some risk factors for ASD have attracted considerable attention, like advanced paternal age (APA) (Croen, Najjar, Fireman, & Grether, 2007; Glasson et al., 2004; Larsson et al., 2005; Lauritsen, Pedersen, & Mortensen, 2005; Maimburg & Vaeth, 2006; Reichenberg, Gross, Sandin, & Susser, 2010). Advanced paternal age is known as a risk factor for ASD (Bilder, Pinborough-Zimmerman, Miller, & McMahon, 2009; Buizer-Voskamp et al., 2011; Burstyn, Sithole, & Zwaigenbaum, 2010; D’Onofrio et al., 2014; Durkin et al., 2008; Frans et al., 2013; Golding, Steer, & Pembrey, 2010; Grether, Anderson, Croen, Smith, & Windham, 2009; Hultman, Sandin, Levine, Lichtenstein, & Reichenberg, 2011; Idring et al., 2014; Lampi et al., 2013; Lundstrom et al., 2010; Mamidala et al., 2013; Parner et al., 2012; Pinborough-Zimmerman et al., 2011; Sasanfar et al., 2010; Shelton, Tancredi, & Hertz-Picciotto, 2010; Sun, Allison, Auyeung, Baron-Cohen, & Brayne, 2014; van Balkom et al., 2012; Williams, Helmer, Duncan, Peat, & Mellis, 2008; Wu et al., 2017; Zhang et al., 2010). Nevertheless, the biological mechanism underlying this association remains unknown. Meanwhile, some hypotheses have been established to explore the molecular mechanism underlying the association of advanced paternal age and ASD. These hypotheses may be divided in two (but they would not be mutually exclusive): one that involves sequence changes in the DNA (de novo mutations) and the other that does not involve this DNA sequence changes (epigenetic changes) (de Kluiver, Buizer-Voskamp, Dolan, & Boomsma, 2017; Janecka et al., 2017). Supporting the hypothesis of sequence DNA changes in APA, it has been proposed that whilst male age increases, the rate of genetic mutations during spermatogenesis also increases, affecting the probability of inheriting de novo mutations to the off-spring which could contribute to an increased morbidity (Buwe, Guttenbach, & Schmid, 2005; Crow, 2000; Flatscher-Bader et al., 2011; Kong et al., 2012; Sanders et al., 2011; Szatmari et al., 2007). In the case of the epigenetic alterations hypothesis (no DNA sequence alterations), establishes that aging affects the chromatin structure and DNA methylated patterns of sexual cells, which could lead to inherited gene expression alterations in the newly formed zygote derivate from altered aging-related epigenetic changes (Jenkins, Aston, Pflueger, Cairns, & Carrell, 2014; Kimura, Yoshizaki, & Osumi, 2018). Contrasting this hypothesis of aging-related changes, in some analysis it has been reported that the difference between the conceptional parental age is a risk factor for ASD, principally by younger mothers (Sandin et al., 2016). Where this difference has to be explored in other populations.
However, most studies are centered in Scandinavian or Asian populations, leaving behind other populations, like, Latin Americans and African Americans. Most research includes children living outside their countries of birth (principally in the United States) and the possible effects of this risk factor in the individuals from these populations residing in their countries of birth, is scarce (Becerra et al., 2014). A recent meta-analysis of APA and ASD, screening almost 67 000 individuals with ASD, concluded that APA is an ASD-risk factor influenced possibly by the residing country (Wu et al., 2017). In this study, Wu et al. reported that the effect of APA in North America is lower than the effect in children residing in Europe and is even lower than the effect of children residing in Asia. In the case of Mexican children, there are no studies analyzing the effect of APA as a risk factor for ASD. This study aims to perform an exploratory analysis of APA as a risk factor for ASD in children residing in Mexico City.
Method
Design of the study
The study design was a case-control, observational, cross-sectional, with a non-probabilistic, and consecutive sampling.
Participants
The total sample included (N = 1 068 children). We considered children with ASD diagnosis (n = 162) as cases and neurotypical children (NT, n = 906) as controls.
Procedure
The structure of the sample evaluation was performed following the same criteria previously reported for a co-validation of an autism instrument (Albores-Gallo et al., 2012). Children with ASD (case group) were recruited from the outpatient psychiatric services of the Hospital Psiquiátrico Infantil Juan N. Navarro from March 2014 to November 2016. Neurotypical children (control group) were from a community-based sample recruited from nurseries located in four different districts of Mexico City. The four kindergartens belong to the DIF system and located in Tlalpan, which receives children from parents who have no social security. This characteristic makes this population very similar to the one that receives medical care in the Hospital Psiquiátrico Infantil. Male and female children, between 18 to 72 months of age were recruited. We excluded children with known comorbid severe chronic diseases, such as diabetes, cancer, or sensory impairments such as deafness, blindness, or a genetic syndrome associated with autism such as tuberous sclerosis, Rett syndrome, or fragile X. All children were screened for ASD with M-CHAT scale (Albores-Gallo et al., 2012). NT children were excluded if they failed more than two critical items or more than three items in the M-CHAT scale. In the case of children in the ASD group, all the children were evaluated by an expert child psychiatry and a senior board-certified child psychiatrist. The ASD diagnosis was based on DSM-IV criteria and confirmed by the ADI-R instrument (Lord, Rutter, & Le Couteur, 1994). For this study, we analyzed the following variables: age of the parents at the moment of the birth of the child (conceptional age).
Statistical analysis
In all the manuscript, we will refer to the conceptional maternal or the conceptional paternal age as only paternal age or maternal age. To explore the normality distribution of the parental age, we used a Kolmogorov-Smirnov test. To compare the differences in parental age between ASD and NT groups, we performed a Wilcoxon test. We categorized the individuals based in the distribution of the overall parental age, based in intervals of 10 years. We used as reference for the next models, were the overall mean, median and mode of the parental age was found.
A series of multivariate logistic regression models were performed to estimate the odds ratios (OR) and 95% confidence intervals (CI 95%) to analyze the paternal age and maternal age as a risk factor for ASD. All the models were adjusted for the opposite parental age (i.e, if we analyzed the parental age, we fitted the models with maternal age as a covariate), age, and gender of the child (covariables). A p less than .05 was considered statistically significant. We fitted three different independent models.
In the first model we measured the associations between the independent variable of maternal or paternal age (i.e., in two independent models, one for maternal age and one for paternal age) and the dependent variable, the ASD case status (all models adjusted by covariables).
In the second model, we also measured the association between the maternal or paternal age and the ASD case status, but in this analysis we categorized the individuals based in the conception maternal and paternal age. Five categories: lower than 20, 20 to 29, 30 to 39, 40 to 49, and higher than 50. In all the models, the OR and CI 95% was calculated, the category where the mean, median, and mode overall sample, parental age was used as the reference category.
Additionally, in the third model, we performed an analysis to search if parental age differences as a risk factor for ASD. For this analysis, we calculated the difference of the paternal age subtracted from the maternal age (paternal age-maternal age) and coded as a continuous variable. Once we performed this parental age difference calculation we categorized the children in four categories: children with differences in parental age higher than 10 years (with higher maternal age); children with difference in parental ages > zero and < 10 (with higher maternal age); children without differences in parental age; children with difference in parental ages > zero and < 10 (with higher paternal age); and children with differences in parental age higher than 10 years (with higher paternal age). Again, multivariate logistic regressions were done to calculate OR and CI 95%. We considered the group without difference between parental age as the reference group for these parameters calculations. All models were fitted using R (R Development Core Team, 2008).
Ethical considerations
The ethics and research committee approved the study of the Hospital Psiquiátrico Infantil Juan N. Navarro (approval number = II1/01/0304/Tc). To be included in the study, written informed consent and assent was obtained from each study participant through the primary caregiver of each infant.
Results
The mean age of the children with ASD (M = 4.40 SD ± 1.11) was slightly lower, but not statistically significant compared to the mean age of the children diagnosed with ASD (M = 4.48 SD ± 1.13) (Table 1). The proportion of males/females in children diagnosed with ASD was higher than those children with NT, being this relationship of three boys to one girl in the children diagnosed with ASD.
ASD (n = 162) | NT (n = 906) | ||
---|---|---|---|
Sex, n (%)a | p < .00c | ||
Male | 123 (75.93) | 480 (52.98) | |
Female | 39 (24.07) | 426 (47.02) | |
Age, M (SD)b | 4.40 (1.11) | 4.48 (1.13) | p = .41 |
Maternal age | 27.60 (5.85) | 26.73 (5.97) | p = .13 |
Paternal age | 31.54 (8.67) | 29.17 (8.67) | p < .001 |
M-CHAT, M (SD) | 6.66 (4.21) | 3.27 (2.19) | p < .0001 |
Notes:
a For sex the total individuals (n) and the percentage inside the parenthesis (%) is reported.
b For age and M-CHAT the mean (M) and the standard deviation (SD) inside the parenthesis is reported.
c For sex a p value is reported for a chi-square test, meanwhile, for age and M-CHAT a Wilcoxon test p value is reported.
Maternal age was slightly higher in mothers of children diagnosed with ASD (M = 27.60 SD ± 5.85), but not statistically significant different from the maternal age of the mothers with NT children (M = 26.73 SD ± 5.97, p = .13). Meanwhile, paternal age was higher in children diagnosed with ASD (M = 31.54 SD ± 8.67) compared to the paternal age of NT children (M = 29.17 SD ± 8.67), being this difference statistically significant (p < .001). In the first regression models, maternal age did not reach statistical significance as a risk factor for ASD (OR = 1.02, CI 95% [.99, 1.06], p = .07). In contrast, paternal age increased risk for ASD (OR = 1.04, CI 95% [1.02, 1.07], p < .001). Also, in the stratified analysis (five categories of parental age), advanced paternal age increased the risk to ASD, exclusively when the age of the father was between 40 to 49 years (OR = 2.57, CI 95% [1.47, 4.49], p < .001) (Table 2). Also, a reduced parental age, principally between 20 to 29 years, reduced the risk to ASD. In contrast, in this stratified analysis, the advanced maternal age did not increase the risk for ASD.
Categories | Paternal age | Maternal age |
---|---|---|
Lower than 20 years | 1.28 [.61 – 2.54] p = .32c | .34 [.14, .69] p < .01 |
20 to 29 yearsa | .65 [.44 – .98] p = .04 d | Reference |
30 to 39 yearsb | Reference | 1.16 [.78, 1.71] p = .45 |
40 to 49 years | 2.57 [1.47 – 4.49] p < .001 | 1.14 [.40, 2.85] p = .79 |
Higher than 50 years | 2.79 [.53 – 12.53] p = .19 | NAe |
Notes:
a The category where the overall mean, median and mode of the paternal age (29.53, 29, 29 years) in the sample was used as reference for OR and CI 95% calculations (≥ 30, < 40).
b The category where the overall mean, median and mode of the paternal age (26.86, 26, 25 years) in the sample was used as reference for OR and CI 95% calculations (≥ 20, < 30).
c The reported values are OR [CI 95%] p.
d In bold are the statistical significant results (values of p less than .05).
e No children of the overall sample have a mother with maternal age higher or equal than 50 years.
In the third model, the difference in parental age was tested as a risk factor for ASD. In this analysis, when the father had a higher age than the mother, the risk for ASD increased for a difference between 10 years (OR = 3.37, CI 95% [1.65, 7.21], p = .001) (Table 3). Also, when the mother had a higher age than the father, this results in a higher risk, but this model reaches only nominal statistical significance.
Higher Maternal Age | Higher Paternal Age | |||
---|---|---|---|---|
Higher or equal 10 years | Higher than 0 and lower 10 years | Higher or equal 10 yearsa | Higher than 0 and lower 10 years | |
.00 [.00 – NA] p = .99 | 1.97 [.99, 4.11] p = .06 | 3.37 [1.65, 7.21] p = .001 | 1.15 [.63, 2.23] p = .66 |
Notes: Individuals within the group with no difference in parental ages were used as a reference group for OR and CI 95% calculation.
a In bold letters, the statistically significant values are observed (values of p less than .05). The reported values are OR [CI 95%] p.
Discussion and conclusion
In this study, we analyzed the advanced parental age as a risk for ASD in children residing in Mexico City. The results showed a small increased risk for ASD when the father had an advanced parental age, and higher when the age difference between parents was larger than 10 years. However, the risk of APA in our sample (OR of 1.04) is lower compared to the result reported in a recent worldwide meta-analysis (OR of 1.41) (Wu et al., 2017). Comparing our preliminary estimate with each geographical area reported in the meta-analysis, this is almost the half of the one reported for Asia and Oceania (OR of 2.56), and lower than Europe (OR of 1.44). This estimate is even lower than the results reported for the North America geographical area (which has the lower risk effect of advanced paternal age in ASD, OR: 1.44). The studies included in this meta-analysis, located in North America, were only recruited in the United States, which is the principal country with migration from Mexico. In this sense, the effect of APA for ASD could be inferred to be the same for the Mexican population residing in Mexico or in the United States. However, epidemiological studies done with Hispanic population residing in the US (included Mexican) suggest that ethnicity might influence the prevalence of ASD. For example, mothers with Hispanic ancestry were less likely to have a child with ASD. Moreover, when the mother is an immigrant, the likelihood of having a child with ASD is reduced (Durkin et al., 2008; Liptak et al., 2008; Pedersen et al., 2012; Pinborough-Zimmerman et al., 2011; Schieve et al., 2012). This could mean that the effect of ASD risk factors in Hispanic residing individuals could be lower than the effects in individuals residing outside their country of birth; nevertheless, this lower probability of ASD remains to be explored in further studies. Another point to be considered in this reduced of likelihood of having a child with ASD in Mexican immigrants could be the reduced level of migration that ASD individuals could have and in consequence reduce the prevalence of ASD in Mexican population residing in the United States. In contrast, in other studies exploring the immigrant status and ASD, an increase has been reported for other populations (Gillberg, Steffenburg, Börjesson, & Andersson, 1987; Gillberg, Schaumann, & Gillberg, 1995), but there are still no plausible hypotheses about these differences. That could be a point to be explored about what role could have, if any, the Mexican ethnicity in the estimates for ASD, and the relationship with APA.
In our analysis, we found that higher differences in parental ages increase the risks for ASD, which could be reinforcing the effect that both parents age could be playing in this relationship. In Mexico, there are only a few studies exploring the dynamics of parental conceptional ages, could prompt future studies in order to explore how this dynamic could alter the risk for ASD and correlate with other sociodemographic variables. One limitation of our study was the collection of the sample from only one city of Mexico, which could only give preliminary data about the effect of advanced paternal age as a risk factor for ASD in the Mexican population. Another limitation is that we sampled children to be in the same age range for children diagnosed with ASD and NT children. This sampling could lead to some stochastic errors in the calculation of the paternal age mean differences between cases and controls. These errors could be increased as an effect of the differences in the sample sizes of cases and controls, we have less than 200 ASD diagnosed children. The limitation in the sample size could impact the differences found in paternal age between ASD diagnosed (M = 31.54 SD ± 8.67) and TD children (M = 29.17 SD ± 8.67). Nevertheless, further studies in larger sample sizes are required.
Even when our results are preliminary, we point out some interesting evidence. Our analysis is the first in estimating the effect of advanced parental age in the risk for ASD in a Mexican children sample, and we found that advanced paternal age ASD could have a lower effect in the risk for ASD in our population; still, analysis on larger samples sizes are required. Also, analysis from other sources of data like epigenetic or genetic factors could help to clarify this difference.