Introduction

The associations between older paternal age and autism spectrum disorders (ASD) is now well established,1, 2, 3, 4, 5 including two recent population-based studies.4, 5 A risk increase with advancing maternal age has also been shown.6 There are, however, important issues to resolve. Meta-analyses have shown considerable study heterogeneity2, 6 and it is still not clear whether paternal and maternal ages represent independent risk factors.7 Most importantly, there is limited information about the combined effect of paternal and maternal age as most studies have considered these as independent. There have been attempts to address this question4, 8 but with no clear conclusion (Supplementary Appendix A). For schizophrenia, parental age risk has been proposed through different biological mechanisms4, 9, 10 and an association with parental age differences has also been suggested.11 A large US study reported increased risk for adverse perinatal outcomes for mothers older than their partner.12 Even the largest population-based studies have been limited by size and thereby unable to reliably examine the risk from parental ages across the entire age distribution, or to properly separate the independent and combined effects of paternal and maternal ages.

Using a large population-based cohort across several geographic regions, the aim of this study was to test if advancing paternal and maternal age are independently associated with risk for ASD and to estimate the functional form of the association between parental age and ASD.

Materials and methods

The study builds on the 'International Collaboration for Autism Registry Epidemiology (ICARE)', combining population-based cohorts across several geographic regions and health systems with the purpose of studying risk factors for ASD.13 Ethics committee approval, with waiver for informed consent, was obtained by each site. The ASD prevalence varies across sites13 but time trends and rates are consistent with contemporary reports from other populations.14 For details see Schendel et al.13

Study population

The study population comprised all live-born singletons in Denmark, Norway and Sweden between 1985 and 2004. Data from Western Australia (WA) included all 1985–1999 singleton live-born births among the non-aboriginal population. Data from Israel comprised all 1993–2004 singleton live-born Jewish births with a subsequent diagnosis of ASD and a random sample of controls born in the same period. We excluded twins because they are highly correlated for ASD risk,15 thereby violating the assumption of independent observations, and they are also associated with maternal age and birth complications16 which could introduce bias.

ASD outcome, parental age and covariate information

Children were followed from birth to reported diagnosis of ASD through 2004 in WA; 2006 in Norway and 2009 in Denmark, Israel and Sweden. Sweden and Denmark provided diagnoses from medical registers whereas in Israel, Norway and WA, it was derived from government-maintained service/benefits registers. Case identification and registry reporting procedures have been described previously13 and Supplementary Appendix B provides details on the reliability and validity of reported diagnoses.

Parental ages, sex and birth year were obtained from birth or civil registers.

Statistical methods

Logistic regression was used to calculate odds-ratios as estimates of relative risk (RR). To achieve a qualitative (descriptive) interpretation of the functional form of ASD risk across parental ages we fitted splines.17, 18 To quantify the association, the RR of ASD was presented by categories of parental age. All statistical tests were performed by using the two-sided 5% level of significance and associated two-sided 95% Wald-type confidence intervals (CIs) were calculated. We did not adjust for multiplicity of statistical tests. Besides the exclusion of twins, it was not possible to identify family members in the data and correct for possible correlations between family members in the analyses.

Because data from Israel were obtained by using a case–control design with known sampling probabilities for cases and controls, rather than a birth cohort design, we included the sampling weights in the analyses.19 To adjust for the variation introduced by the sampling weights, robust standard errors were used.20

Paternal and maternal age independently

To visualize the shape and functional form of the RRs, we fitted splines to continuous paternal and maternal ages separately, adjusting for site, sex and potential confounding factors: the categorical age of the other parent (<20, 20–24, 25–29, 30–34, 35–39, 40 for mothers and <20, 20–24, 25–29, 30–34, 35–39, 40–44 and 45 for fathers), and child's birth year (4-year intervals). To quantify the risk, RRs were calculated by fitting logistic regression considering parental age categorically (mothers: <20, 20–29, 30–39, 40; fathers: <20, 20–29, 30–39, 40–49 and 50). All models had separate intercept parameters for each site, birth year (4-year intervals) and for male and female offspring.

Subgroup analyses

We further adjusted for the effect of the partner's age by analyzing parental age in sub-groups of the other parent's age. To avoid bias from uneven distribution of ASD cases in the different sub-groups we created sub-groups of maternal age with approximately equal numbers of ASD cases resulting in maternal age cut-offs of <27, 27–31 and >31 years. For each sub-group, we fitted the logistic spline models, both crude and adjusted. We used the same approach to define sub-groups of paternal age (<29, 29–34 and >34) while analyzing the maternal age.

Paternal and maternal age jointly

In a new approach enabled by the large sample size we analyzed paternal and maternal ages jointly on a bivariate continuous scale using thin-plate splines21 adjusting for sex and birth year categorically as above. The thin-plate spline is a two-dimensional smoother wherein the curvature is estimated in local regions moving over the bivariate (two-dimensional) surface of parental age. This is a descriptive statistical approach displaying the functional form of the relationship between bivariate paternal–maternal age and ASD with a qualitative interpretation rather than focusing on numerical estimates.

Finally, to determine the importance of independent and joint parental age effects on ASD risk, we fitted logistic regression models including separate covariates for the mean of the parent's ages (capturing the direct effect of ageing) and the difference in age between the parents (capturing effects not necessarily associated with ageing). We compared the goodness-of-fit of the models when including the two components together and each one separately by calculating the Akaike information criterion.

Supplementary analyses

We estimated the RR for male and female offspring separately. We repeated all analyses for autistic disorder (AD) only (excluding Israel where this information was not available).

The logistic regression model, assumes that the time at risk is the same for all subjects or that the length of follow-up does not affect the risk. If this assumption is violated, bias may be introduced. For this purpose, birth year was included in all models. We also performed a sensitivity analysis including Denmark and Sweden only, as both of these countries had information about date of diagnosis. Here we used Cox regression with age as the underlying time-scale and compared results with the logistic regression results. This approach adjusts for calendar effects as well.22

Site heterogeneity was addressed by visual inspection of the graphs of the spline predictions and by verifying pooled results through the use of models that excluded any potentially heterogeneous sites.

Data management

Statistical analyses were performed using SAS 9.3. Thin-plate splines were fitted using the gamm4, v0.1-6 package23 running R software 2.15.2.24

Results

The cohort comprised of 5 766 794 births (Denmark:1 270 229; Israel: 1 015 130; Norway: 1 121 392; Sweden: 2 023 017; and WA: 337 026). There were 30 902 (0.54%) children with ASD, and 10 128 (0.18%) with AD. Distributions of births, ASD cases, birth year and parental age are presented in Table 1 (site-specific data: Supplementary Table 1).

Table 1 Study population characteristics: covariate distributions by parental age

Independent effects of paternal and maternal age

Qualitatively, there was a monotonic increase in ASD risk with increasing paternal age (Figure 1). In contrast, except for Israel and WA, there was a U-shaped association with maternal age with the lowest risk observed shortly after the age of 30. Relative to the age of 30 there was a statistically significantly increased risk of ASD associated with advancing, as well as younger maternal age (Figure 2, Table 2).

Figure 1
figure 1

Relative risk for autism spectrum disorder with increasing paternal age, for each site separately and combined. Point estimates and pointwise two-sided 95% confidence intervals. Relative risk calculated using the parental age of 30 years as reference. Relative risk estimated by ordinary logistic regression adjusting for age of the mother (<20, 20–24, 25–29, 30–34, 35–39 and 40 years), sex (male/female), birth year (4-year intervals) and, for all sites combined, site (Denmark, Israel, Norway, Sweden and Western Australia).

PowerPoint slide

Figure 2
figure 2

Relative risk for autism spectrum disorder with increasing maternal age, for each site separately and combined. Point estimates and pointwise two-sided 95% confidence intervals. Relative risk calculated using the parental age 30 as reference. Relative risk estimated by ordinary logistic regression adjusting for age of the father (<20, 20–24, 25–29, 30–34, 35–39, 40–44 and 45 years), sex (male/female), birth year (4-year intervals) and, for all sites combined, site (Denmark, Israel, Norway, Sweden and Western Australia).

PowerPoint slide

Table 2 Relative risk for ASD and AD, with associated two-sided 95% confidence intervals by maternal and parental age

RRs in 10-year age categories are presented in Table 2. For example, relative to fathers aged 20–29 years, fathers 50 years or older had a statistically significantly increased risk for offspring with ASD (RR=1.66 95% CI:1.49–1.85). Relative to mothers aged 20–29 years, mothers younger than 20 years had a statistically significantly increased risk for offspring with ASD (RR=1.18 95% CI:1.08–1.29) (Table 2).

Similar patterns of association, but with slightly higher RRs for the highest parental ages, were evident for AD (Supplementary Figures 1 and 2, Table 2).

Sub-group analyses

An increase in ASD risk with advancing paternal age was evident in all three sub-groups of maternal age (<27, 27–31 and >31 years) (Supplementary Figure 3, right panel). For the oldest mothers (>31 years) the effect of advancing paternal age was delayed compared with mothers younger than 32 years, with no increase in risk before the father turned 40. For the two youngest sub-groups (mothers 31 years or younger), the increases RRs by advancing paternal age were similar (overlapping RRs), including a statistically significant decrease in RR when fathers were younger than 30 years.

The ASD risk by maternal age association was similarly U-shaped in each sub-group of paternal age (<29, 29–34 and >34 years) (Supplementary Figure 3, left panel), with increasing RR for both younger and older mothers, and, in each sub-group, the lowest RR for parents of similar age.

Joint effects of paternal and maternal ages

Couples of mothers 20–39 years of age and fathers 20–49 years of age generated 94.6% of all births. Among all children, 2.2% had a mother younger than 20 years and 2.4% had a mother 40 years or older; 0.9% had a father older than 50 years (Supplementary Table 2). The bivariate distribution of number of births by parental age across all sites is presented in Supplementary Figure 4.

Using couples wherein the parents were both aged 25 years as the reference, the RR of ASD as a function of the joint influence of paternal and maternal age is displayed in Figure 3. Colors ranging from light green toward the yellow and red indicate increasing RR, while dark green indicates decreasing RR compared with the 25-year-old parents. The corresponding pointwise upper and lower 95% CIs are presented in Supplementary Figure 5.

Figure 3
figure 3

Relative risk (RR) for autism spectrum disorder by paternal and maternal age jointly. RR calculated by using the maternal and paternal ages of 25 years as reference. RR estimated by ordinary logistic regression adjusting for sex (male/female), birth year (4-year intervals), and, site (Denmark, Israel, Norway, Sweden and Western Australia). Blue color indicate RR≈1, light green, to yellow and red indicate increasing RR>1.

PowerPoint slide

The joint effect of paternal and maternal ages on the RR of ASD (Figure 3) shows an inverse-shaped form compared with the age of parenting distribution (Supplementary Figure 4). That is, lowest risk corresponded to couples that generated the majority of births, specifically, 29–39-year-old fathers and 25–35-year-old mothers. RR increased in all directions from this region as the parental age difference increased. Furthermore, the RR pattern revealed that for mothers older than 40 years the RR increased in a U-shaped pattern with younger and older partner. For fathers 45 years or older the RR increased monotonically with increasing maternal age.

Illustrated in Figure 3, the highest ASD risk was evident among couples characterized by fathers older than ~45 years, independent of maternal age; fathers 35–44 years with mothers at least 10 years younger; mothers 30–39-years old with fathers at least 10 years younger. These couples corresponded to 7% of all births. In addition, among couples with fathers younger than 45 years and mothers younger than 40 years, the RR increased with increasing difference in age between the parents.

Examining whether mean parental age or parental age-difference explained most of the variation in ASD risk, there was strong support that they both were important, jointly (Supplementary Table 3).

Supplementary analyses

The RRs for ASD were similar for male and female offspring (Supplementary Figure 6). Cox regression using data from Denmark and Sweden produced RRs very similar to the RRs from the logistic regression models (Supplementary Table 4).

By visual inspection, the functional form of the RR for paternal age was similar across sites. For maternal age, results in Israel appeared to diverge from the other sites. This did not, however, bias the pooled results (Supplementary Figure 7). After excluding Israel, (a) the functional form of the RR for maternal age was similar (Supplementary Figure 7), (b) the RR associated with advancing maternal age in 10-year age categories showed similar results (data not shown) and (c) the three-dimensional figure for the bivariate effect of paternal and maternal ages was very similar (Supplementary Figure 8).

For all sites combined, the joint effect of paternal and maternal age on AD risk was qualitatively similar to that of ASD (Supplementary Figure 9).

Discussion

This study provides the strongest evidence to date supporting the hypothesis that advanced parental ages at the time of birth are independently associated with risk for ASD in the offspring; with no support for any modification by the sex of the child. There was also evidence for a combined parental age effect. The risk was highest when both parents were older, but the risk was also increased among disparately aged parents. Furthermore for all sites combined, and in each of the Scandinavian populations, there was an increased ASD risk for mothers younger than ~30 years; this was not observed in the Israeli population, which warrants further investigation.

The observed associations between parental age and ASD risk were similar in magnitude to those reported in earlier meta-analyses,2, 6 and with RRs higher for AD compared with ASD. The slight differences in RR between ASD and AD may indicate an effect of co-occurring intellectual disability (ID), more frequently reported in AD.25, 26 An earlier study showed higher RR for maternal and paternal age in children with ID.27 While we did not have data on ID in our database; data from a previously published study from Sweden28 estimated the proportion of children with ID among AD cases to 43% and only 15% among children with ASD, excluding AD.

Comparing fathers and mothers over the same age range, the RRs with advancing age were of similar magnitude. This suggests that advancing paternal age may contribute more to the risk than advancing maternal age overall, due to the longer male child-bearing potential. Human29, 30 and animal31 studies provide support for the hypothesis that de novo mutations contribute to the association between paternal age and ASD. In contrast, mechanisms mediating the effect of advancing maternal age on ASD risk have not been frequently investigated. Advancing maternal age has been associated with chromosomal changes32, 33 and genomic modifications.34

Our study, however, may provide insights into other potential etiological mechanisms underlying paternal and maternal age-associated factors in ASD risk. Our analyses showed that there was an increase in risk also when the difference between the parental ages was moderate-to-large (10 years or more) (Supplementary Figure 4). This suggests that the increase in risk is not attributable to advancing parental age per se, and that the risk increase cannot be explained solely by an accumulation of point mutations or other genomic alterations in the parents. Our data currently do not allow for further exploration of this phenomenon, but the mechanism underlying such reproductive strategies should be considered further. It is possible that the partners in these disparately aged couples represent other socio-economic (SES), genetic, and/or psychological characteristics that increase their risk of having children with ASD. Genetic shared influences (pleiotrophy) between ASD and uncommon reproductive ages may harbour potential biological mechanisms.35

The role of low SES in the association between parental age and ASD has not yet been resolved, and can potentially effect the associations in different ways. Firstly, SES can be a source for ascertainment bias (for example, in settings where access to health care is differential with respect to SES, as in the USA). In this situation, adjusting for maternal SES is warranted to control for confounding between maternal age and ASD. A study from Denmark36 examining multiple indices of SES and risk for ASD concluded that SES have little or no role in the risk of ASD, at least where access to the health-care system is equally available for all and is free of charge.

Secondly, younger parental age can have consequences effecting education, employment and other markers of SES. Therefore, while a low SES may be a risk factor for ASD, it may lie on the same causal path as younger parental age and act as a factor modifying the risk of ASD. Such a complex interplay between SES and maternal age can be accompanied by other risks for poor outcome, including maternal infection and poorer health monitoring. In this situation, adjusting for SES, as a mediator between maternal age and ASD, would be erroneous and may even introduce bias in the maternal age-ASD risk estimate.

Alternatively, low SES may be a proxy for other shared factors that may contribute to ASD risk, for example, genetic confounding. In this situation, adjustment for SES may be warranted as well. Therefore, the decision to adjust, or not, for SES deserves careful consideration in future studies of ASD. A recent publication has discussed this issue, also in the context of parental age, did not adjust for SES and reported a U-shaped relationship.4 Other studies that adjusted for SES mostly report a monotone association with maternal age.2 A Swedish study showed a U-shaped association with maternal age for ASD with co-morbid ID but not for ASD without.27

Social phobia and traits such as shyness and aloofness that may limit cross-sex interactions have been described for parents of children with ASD.37, 38 These factors could influence the age when a person establishes a relationship and has children, confounding the association between parental age and risk of ASD. A US study with nine million pregnancies showed altered risks of fetal death, preterm birth and small size for gestational age with increasing parental age differences.12 More generally, large trait differences between partners have been shown to be associated with adverse pregnancy outcomes. For example, in a large study from the UK, an increasing father–mother height difference was associated with abnormal pregnancy outcomes.39 Also, higher fertility has been seen in men with partners about 6 years younger and in women with partners ~4 years older, indicating fitness-benefits with certain combinations of parental age.40 Environmental and perinatal factors may also contribute cumulatively or at specific ages; pre- and perinatal complications are associated with ASD risk41 and a U-shaped risk profile has been reported for preterm birth and birth weight in relation to parental age.42

Study strengths include the large sample size, with prospective follow-up, in national cohorts, and along with rigorous data harmonization and quality control processes.13 The associations between parental age and ASD were similar regardless of whether clinical or service/benefits registers were used for case ascertainment. The sites' publicly financed health systems minimizes any potential bias due to differential access to health care.

The ASD prevalence observed was highest in Denmark and Sweden, somewhat lower in WA, and lowest in Israel and Norway. The range of prevalence rates in the contributing sites (12–84 per 10 000) is consistent with rates reported in epidemiological studies of ASD during the same decades.43 Rates in Denmark and Sweden are similar to those reported in newer epidemiological surveys.43 Prevalence's and time trends across sites are discussed in detail in an earlier methods paper.13

Study limitations include sparse data for extreme parental ages, despite the large sample size. Also we lack information about potentially confounding variables such as SES and parental psychiatric history.35, 44 We cannot rule out the possibility that other factors associated with parental age (for example, length of marriage or partnership, obstetric complications, gestational age and birth weight) have an important role in explaining our results. Time-trends are also possible sources of bias. Differences in year of birth and length of follow-up may explain part of the site-to-site variations. Using date of diagnosis would have allowed for a more detailed statistical adjustment for length of follow-up, but the assumptions about the data distribution used here follow the main approach used by most earlier studies.1, 2, 30, 45 Sweden and Denmark provided diagnostic information from medical registers whereas for Israel, Norway and WA, diagnostic information was derived from service/benefits registers. Service/benefit registers could potentially be biased towards cases associated with co-morbid ID. We did not have individual level information on co-morbid ID in ASD cases. Nevertheless, validation studies, from Sweden46 and Denmark47 reported rates of ID in ASD cases that were comparable to, or even somewhat lower than, reports in recent epidemiological studies identifying ASD through screening and diagnostic evaluations.48, 49 Analyses restricted to Sweden and Denmark produced RRs very similar to the RRs observed for all sites (Supplementary Table 4). Thus, while we could not completely rule out bias due to service registers potentially including a higher proportion of ASD cases with co-morbid ID, we found no evidence that this could substantially effect the overall results of the present study.

Conclusion

Shown in multiple geographic regions, increases in ASD risk was not only limited to advancing paternal or maternal age alone but also to differences in parental age including younger or older similarly aged parents as well as disparately aged parents.

These results suggest that multiple mechanisms are contributing to the association between parental age and ASD risk.