 Research article
 Open access
 Published:
Disentangling the roles of maternal and paternal age on birth prevalence of Down syndrome and other chromosomal disorders using a Bayesian modeling approach
BMC Medical Research Methodology volume 19, Article number: 82 (2019)
Abstract
Background
Multiple neonatal and pediatric disorders have been linked to older paternal ages. Combining these findings with the evidence that many men are having children at much later ages generates considerable public health concern. The risk of paternal age has been difficult to estimate and interpret because children often have parents whose ages are similar and likely to be confounded. Epidemiologic studies often model the conditional effects of paternal age using regression models that typically treat maternal age as linear, curvilinear or as ageband categories. Each of these approaches has limitations. As an alternative, the current study measures age to the nearest year, and fits a Bayesian model in which each parent’s age is given a conditional autoregressive prior (CAR).
Methods
Data containing approximately 12,000,000 birth records were obtained from the United States Natality database for the years 2014 to 2016. Date were crosstabulated for maternal ages 15–49 years and for paternal ages 15–65 years. A Bayesian logistic model was implemented using conditional autoregressive priors for both maternal and paternal ages modeled separately and jointly for both Down syndrome and chromosomal disorders other than Down syndrome.
Results
Models with maternal and paternal ages given CAR priors were judged to be better fitting than traditional models. For Down syndrome, the approach attributed a very large risk to advancing maternal age with the effect of advancing paternal age having a very small sparing effect on birth prevalence. Maternal age was also related to the birth prevalence of chromosomal disorders other than Down syndrome while paternal age was not.
Conclusions
Advancing paternal age was not associated with an increase in risk for either Down syndrome or chromosomal disorders other than Down syndrome.
Background
A rapidly increasing list of neonatal and pediatric disorders has been linked to older paternal ages [1]. Combining these findings with the evidence that many men are having children at much later ages generates considerable public health concern [2, 3]. The best known of the many conditions that have been linked to older paternal age are stillbirths, birth defects, childhood cancers and neurodevelopmental disorders, specifically autism spectral disorders and schizophrenia [1]. However, considerable controversy exists in identifying the conditions caused by paternal age because the analysis needs to adjust for mother’s age [4,5,6,7]. There is considerable theory to explain potential associations between neonatal disorders and paternal age. It has long been known that males, with advancing age, have a nonlinear increase in germline mutations with age related to cumulative changes with the spermatagonial stem cells [8]. These effects result from agerelated changes that compromise DNA replication, DNA repair, cell cycle control, and epigenetic modifications in spermatagonial stem cells and these errors accumulate with successive mitotic divisions [9, 10] and contribute to de novo mutations, affecting genetic traits in a variety of ways [11]. Now, the evidence that paternal germline mutations are responsible for a variety of conditions is considered overwhelming and the list continues to grow [1]. In spite of the perceived magnitude of the problem, the epidemiologic search for the causative mutagens has stalled. Three reviews over a period of 18 years trace the history of the, so far, futile search for mutagens responsible for paternal germline mutations [12,13,14]. The most important difficulty appears to be confounding, especially the confounding by maternal age and, presumably by maternal exposures.
When modeling the joint effects of maternal and paternal ages, two approaches have predominated [15]. Often, age is measured to the closest year and maternal and paternal ages are modeled as linear or as curvilinear (linear and quadratic). This approach is usually inadequate because the linear and quadratic functions will often fit well over specific age ranges and fit poorly over other age ranges. Furthermore, the best fitting linear and quadratic forms will be dependent upon the scale of the model. For example, the scale for the logistic model is usually loglinear as opposed to linear. The second predominant approach has been to stratify ages into categories which can leave residual confounding within age categories [15]. As an alternative, we propose a Bayesian modeling approach that measures age to the nearest year and models each parent’s age as a conditional autoregressive (CAR) [16]. The CAR prior facilitates smoothing of agespecific parental risk estimates to the risk estimates of ages one year younger and one year older as an autoregressive function. This approach facilitates relatively precise estimation of agerelated risk especially under the condition of nonlinearity. While joint conditioning of both parental ages is arbitrarily complex, it is straightforward under a Markov Chain Monte Carlo (MCMC) implementation. Under the Bayesian MCMC implementation, the risk estimate for one parent’s age will be adjusted for the full distribution of possible effects of the other parent’s age (not just the mean of the expected risk, for example). The objective of this study was to parse the maternal and paternal age effects on Down syndrome (DS) and chromosomal disorders other than Down syndrome (CD). The novel approach should help resolve the current uncertainty on the direct effects of paternal age on these syndromes. Furthermore, an approach to parsing parenteral age effects for a wide variety of disorders will be illustrated and validated. The approach has potential to provide an advantage to the estimation of the risks of paternal age and, thus, could enable the identification of multiple disorders mediated by mutations during spermatogenesis. Such an advantage may promote the identification of specific cumulative exposures contributing to the causes of agerelated paternal risk.
Methods
Database
Data containing approximately 12,000,000 birth records were obtained from the United States Natality database for the years 2014 to 2016. In the United States, state laws require birth certificates to be completed for all births, and federal law mandates national collection and publication of births and other vital statistics data. The National Vital Statistics System, the federal compilation of these data, is the result of the cooperation between the National Center for Health Statistics (NCHS) and the states to provide access to statistical information from birth certificates. This study was evaluated by the Texas A&M Institutional Review Board (IRB) and determined to be exempt from IRB review.
Model 1 – maternal age randomwalk (CAR)
For each of DS and CD, case counts were crosstabulated by j = 35 maternal ages (15 to 49 years). For each row in the table Y_{j} was the count of cases, at birth, and n_{j}, the count of births. The counts, Y_{j} were modeled as independent Binomial distributions conditional on an unknown rate parameter (μ_{j}).
The logit of the rate parameter was then modeled as a linear function of the overall intercept and a random effect for each maternal age.
The intercept was given a flat, improper prior. The maternal prior was a minimally informative CAR or random walk prior of length 35 (ages (j) = 15 to 49). The precision of the CAR prior was specified as uniform (0,10) on the standard deviation scale.
Model 2 – paternal age randomwalk (CAR)
For each of DS and CD, case counts were crosstabulated by k = 51 paternal ages (15 to 65 years). For each row in the table Y_{k} was the count of cases, at birth, and n_{k}, the count of births. The counts, Y_{k} were modeled as independent Binomial distributions conditional on an unknown rate parameter (μ_{k}).
The logit of the rate parameter was then modeled as a linear function of the overall intercept and a random effect for paternal age.
The intercept was given a flat, improper prior. The paternal prior was a minimally informative CAR or random walk prior of length 51 (ages (k) = 15 to 65). The precision of the CAR prior was specified as uniform (0,10) on the standard deviation scale.
Model 3 – fully conditional randomwalk (CAR)
For each of DS and CD, case counts were crosstabulated by j = 35 maternal ages (15 to 49 years) and k = 51 paternal ages (15 to 65 years). For each row in the table Y_{jk} was the count of cases, at birth, and n_{jk}, the count of births. The counts, Y_{jk} were modeled as independent Binomial distributions conditional on an unknown rate parameter (μ_{jk}).
The logit of the rate parameter was then modeled as a linear function of the overall intercept and a random effect for each maternal and paternal age.
The intercept was given a flat, improper prior. The maternal prior was a minimally informative CAR or random walk prior of length 35 (ages (j) = 15 to 49). The paternal prior was a minimally informative CAR or random walk prior of length 51 (ages (k) = 15 to 65). The precision of both CAR priors was specified as uniform (0,10) on the standard deviation scale.
The implementation allowed a burnin of 5000 iterations then the next 10,000 iterations were sampled for the posterior distribution. Convergence was evaluated by observing convergence of separate chains with diverse starting values. The median, the lower 2.5% limit and the upper 97.5% limit were all drawn from the complete posterior distributions. The authors refer to the interval from the 2.5 percentile to 97.5 percentile values as the 95% Bayesian credible interval. When the lower bound of this credible interval is greater than 1, the value for the Bayesian exceedance probability would be greater than 95% which would be relatively analogous to a frequentist pvalue of less than 5% for a 2tailed test [17]. The CAR prior produces estimates of random effects that sum to zero at the scale of the log odds. For presentation purposes, the CAR estimates were transformed to odds ratios standardized to parental ages of 15 years.
The CAR models were compared to the Bayesian version of more common models including a linear model, a linear and quadratic model that we refer to as curvilinear and a model that divided age into 5year age categories. The age categories were 15–19, 20–24, 25–29, 30–34, 35–39, 40–44 and 45–49 for each parent’s age and three additional categories for fathers’ ages namely, 50–54, 55–59 and 60–65. Minimally informative Normal priors with zero mean and wide variance, specifically N(0,1000), were used for intercepts, linear, quadratic and agecategory effects. To compare final CAR models to these three models, all odds ratios were adjusted to use the overall mean risk as the baseline risk and the medians from the posterior distributions were plotted. Model fit was evaluated using the Deviance Information Criterion (DIC). [18] All models and the data are available in the Additional files 1–6.
Results
The study identified 11,943,020 births over the threeyear period. Of these births, 10,293,589 could be determined have the mother’s age within the 15 to 49year range and the father’s age belonged in the 15 to 65year range. Excluded observations included 1,642,373 births for which the father’s ages were not recorded. These observations included all births for which a father was not identified. Exclusions included 324 for which both mother and father were younger than 15 y and 49 births for which both the mother was older than 49 and father was older than 65. In observations for which the mothers’ ages were eligible, the father’s age was greater than 65 for 2539 births and less than 15 y for 450 births. When fathers’ ages were eligible, the mother’s age was greater than 49 y for 1894 births and less than 15 y for 1802. There were no exclusions for either combination of the father older than 65 y and mother younger than 15 y or father younger than 15 y and mother older than 49 y.
The crosstabulated data (n = 10,293,589) included 5390 children with Down syndrome with the diagnosis confirmed for 2273 and listed as probable for 3117 children. A diagnosis is considered “confirmed” after chromosomal evaluation and is considered “probable” when based on clinical signs at birth. Chromosomal disorders other than DS were identified in 4147 children including confirmed for 1349 children and probable for 2798 children.
The odds ratio for maternal age, unadjusted for paternal age, started to increase at approximately age 30 and then increased relatively constantly, on the log scale, up to age 45 y, where the odds ratio appeared to stop increasing. The maximum odds ratio relative to 15yearold women, was approximately 16fold (Fig. 1a). The odds ratio for paternal age, unadjusted for maternal age, started to increase at approximately age 30 and then increased relatively constantly, on the log scale, up to age 45 y, where the odds ratio appeared to stop increasing. The maximum odds ratio, relative to 15yearold men, was approximately fourfold (Fig. 1b). The plot of adjusted odds ratios for maternal age was very similar to the plot of the unadjusted odds ratios. When comparing a mother’s age 45 y to age 15 y, the median odds ratio and 95% credibility intervals was 18.9 (11.1, 32.7) for Down syndrome, when adjusted for paternal age (Fig. 1c). For paternal ages, the plot of odds ratios showed the odds ratios to be very near unity but tending to show risk sparing (i.e., odds ratio less than 1). After age 49 y, the odds ratio had a 95% credibility interval that excluded 1. When comparing age 45 y to age 15 y, the median odds ratio and 95% credible interval was 0.81 (0.60, 1.01) for Down syndrome (Fig. 1d).
The odds ratio for the effect of maternal age on CD, unadjusted for paternal age, started to increase at approximately age 30 and then increased relatively constantly, on the log scale, up to age 45 y, where the odds ratio appeared to stop increasing. The maximum odds ratio, relative to 15yearold women, was approximately sixfold (Fig. 2a). The odds ratio for paternal age, unadjusted for maternal age, started to increase at approximately age 30 and then increases relatively constantly, on the log scale, up to age 45 y where the odds ratio appeared to stop increasing. The maximum odds ratio, relative to 15yearold men, was approximately twofold (Fig. 2b). For the effect of maternal age on CD, the plot of adjusted odds ratios was very similar to the plot of the unadjusted odds ratios. When comparing age 45 y to age 15 y, the median odds ratio and 95% credibility interval was 5.8 (3.9, 9.0) for CD (Fig. 2c). For paternal age, the plot of odds ratios showed the odds ratios for CD were very precisely near unity throughout all paternal ages. When comparing age 45 y to age 15 y, the median odds ratio and 95% credible interval was 0.98 (0.74, 1.21) for chromosomal disorders other than Down syndrome (Fig. 2d).
When comparing the final CAR models with more traditional models, all models showed strong effects of maternal age on both Down syndrome and other Chromosomal disorders (Fig. 3a and c). The linear model did not fit the other 3 models well. The curvilinear model was very similar to the CAR model and 5year age category model up to age 45, where the CAR model and 5year age category model both produced an inflection point. In modeling paternal effects, all models showed agreement that the effect of increasing age is related to a decrease in the birth prevalence of Down syndrome (Fig. 3b and d). None of the four models showed an association with paternal age and birth prevalence of other chromosomal disorders. For modeling the effects of paternal age, the 5year age category model was much less smoothed than the other models, although this is shown on a very fine scale for the odds ratios. The Deviance Information Criterion (DIC) showed that the CAR model provided a far superior fit than the alternative models for the joint maternal and paternal effects (Table 1).
Discussion
The name “random walk prior” is used more often in time series analysis than “conditional autoregressive” but they are the same [16]. In timeseries, the model is favored when longterm trends vary from linearity. The current study provides ample evidence that the maternal age function for the logit of the birth prevalence is nonlinear in that the log odds does not increase uniformly each year of age. There were at least two inflection points. In comparison of the CAR model with more traditional models, the most obvious advantage of the CAR model was its ability to model multiple inflection points and rates of change for the risk. This is a result of the CAR or random walk prior being nonparametric in that there is no assumed structure among ages other than correlation among ages one year younger and one year older [16]. Nonparametric regression has been described as a preferable approach over ageband categories, fractional polynomials and spline regression but, at the time, the availability of userfriendly software was a limitation [19]. The current study used readily available software and incorporated commonly implemented and well justified prior values [20]. The OpenBUGS code and data used in the current study are provided in the supporting information. Further applications should be able to identify or resolve paternal age effects for a wide range of disorders including both childhood cancer and birth defects using existing databases.
There exists ample prior support to model maternal and paternal age effects as independent random walks for a wide variety of conditions. In females, age effects are expected to be attributed to meiosis which begins in the fetus, goes into a long period of arrested development and then is reinitiated at ovulation. Clearly, cumulative exposure would impact the arrested cells. In males, spermatagonial cells are formed by mitosis starting at puberty and, at puberty, males start a continuous process of meiosis. The sperm cells participating in fertilization began meiosis a few months before conception. While it has been reported that agerelated epigenetic changes to sperm are often caused by current agerelated exposures, [21] the evidence is overwhelming that, for mutations, age is acting as a surrogate for cumulative exposures of which some might be preventable. In the male, agerelated, cumulative exposures will have much more impact on mitosis than meiosis. Even though faulty mitosis is considered to be more relevant in aging males, faulty paternal meiosis has been reported to cause approximately 10% of Down syndrome cases [22]. There exist important needs to identify the exposures that cause nondisjunction in oocyte development which is certainly agerelated and nondisjunction in sperm development for which the role of paternal age remains unclear. From the epidemiologic perspective, the identification of the relevant risks attributable to parental ages is imperative.
For DS, the risk of maternal age did not change when controlling for paternal age. On the other hand, paternal age effects changed from very large risk to a small sparing risk when controlling for maternal age. According to a recent systematic review, a very small but statistically significant sparing effect, for paternal aging, is a novel finding [23]. In the systematic review, it was concluded that higher paternal age is probably associated with a small increase in the incidence of trisomy 21 [23]. The current study provides relatively precise risk estimates by maternal age but the risk is for birth prevalence. For Down syndrome, in the United States, both the elective termination rate and the natural loss rate are approximately 30% following diagnosis which is possible as early as 10 weeks of pregnancy and, thus, incidence (at conception) and prevalence at birth or at a time of fetal karyotyping will be very different [24]. In addition, the loss prior to 10 weeks is largely unknown but more than half normalappearing IVFproduced embryos are aneuploidy, including often Trisomy 21 [25]. It has been shown that older women are less likely to choose elective termination with a prenatal DS diagnosis [26] but the influence of paternal age on elective pregnancy termination appears to be unknown. The very small sparing risk of advancing paternal age on birth prevalence of DS could be explained by an influence of increasing paternal age to increase the likelihood of an elective pregnancy termination. This potential bias would not be present in conditions that are not diagnosed prenatally.
Conclusions
Advancing paternal age was not associated with an increase in risk for either Down syndrome or chromosomal disorders other than Down syndrome. For those who are familiar with Bayesian models, the proposed approach is simple to implement and interpret. Further applications are encouraged and supported.
Abbreviations
 CAR:

Conditional autoregressive
 CD:

Chromosomal disorders other than Down syndrome
 DS:

Down syndrome
 IRB:

Institutional Review Board
 MCMC:

Markov Chain Monte Carlo
 NCHS:

National Center for Health Statistics
References
Andersen AMN, Urhoj SK. Is advanced paternal age a health risk for the offspring? Fertil Steril. 2017;107(2):312–8.
Sigman M. Introduction: what to do with older prospective fathers: the risks of advanced paternal age. Fertil Steril. 2017;107(2):299–300.
Jennings MO, Owen RC, Keefe D, Kim ED. Management and counseling of the male with advanced paternal age. Fertil Steril. 2017;107(2):324–8.
King M, Bearman P: Advancing paternal and maternal age are both important for autism risk King and Bearman respond. Am J Public Health 2010, 100(5):773–773.
Durkin MS, Maenner MJ, Newschaffer CJ. Estimated autism risk, older reproductive age, and parameterization. Am J Public Health. 2010;100(3):389–90.
King MD, Fountain C, Dakhlallah D, Bearman PS: Estimated autism risk, older reproductive age, and parameterization response. Am J Public Health 2010, 100(3):390–390.
Reichenberg A, Cross R, Sandin S, Susser ES. Advancing paternal and maternal age are both important for autism risk. Am J Public Health. 2010;100(5):772–3.
Crow JF. The origins patterns and implications of human spontaneous mutation. Nat Rev Genet. 2000;1(1):40–7.
Kong A, Frigge ML, Masson G, Besenbacher S, Sulem P, Magnusson G, Gudjonsson SA, Sigurdsson A, Jonasdottir A, Jonasdottir A, et al. Rate of de novo mutations and the importance of father's age to disease risk. Nature. 2012;488(7412):471–5.
Pasqualotto FF, Pasqualotto EB. Aging and sperm DNA damage. Sperm Chromatin: Biological and Clinical Applications in Male Infertility and Assisted Reproduction; 2011. p. 337–49.
Goriely A, Maher G, Lim J, Taylor IB, McGowan SJ, Pfeifer S, RajpertDeMeyts E, McVean GAT, Wilkie AOM: Paternal age effect and selfish mutations. Schizophr Res 2012, 136:S4S4.
Shelby MD. Human GermCell Mutagens. Environ Mol Mutagen. 1994;23:30–4.
Wyrobek AJ, Mulvihill JJ, Wassom JS, Malling HV, Shelby MD, Lewis SE, Witt KL, Preston RJ, Perreault SD, Allen JW, et al. Assessing human germcell mutagenesis in the postgenome era: a celebration of the legacy of William Lawson (bill) Russell. Environ Mol Mutagen. 2007;48(2):71–95.
DeMarini DM. Declaring the existence of human germcell mutagens. Environ Mol Mutagen. 2012;53(3):166–72.
Reijneveld SA: Age in epidemiological analysis. J Epidemiol Commun H 2003, 57(6):397–397.
Besag J, York J, Mollie A. Bayesian imagerestoration, with 2 applications in spatial statistics. Ann I Stat Math. 1991;43(1):1–20.
Zhang JL. Comparative investigation of three Bayesian p values. Comput Stat Data An. 2014;79:277–91.
Spiegelhalter DJ, Best NG, Carlin BR, van der Linde A. Bayesian measures of model complexity and fit. J Roy Stat Soc B. 2002;64:583–616.
Greenland S. Doseresponse and trend analysis in epidemiology  alternatives to categorical analysis. Epidemiology. 1995;6(4):356–65.
Lunn D, Spiegelhalter D, Thomas A, Best N. The BUGS project: Evolution, critique and future directions. Stat Med. 2009;28(25):3049–67.
Sharma R, Agarwal A, Rohra VK, Assidi M, AbuElmagd M, Turki RF. Effects of increased paternal age on sperm quality, reproductive outcome and associated epigenetic risks to offspring. Reprod Biol Endocrin. 2015;13.
Zaragoza RV, Jacobs PA, James RS, Rogan P, Sherman S, Hassold T. Nondisjunction of human acrocentric chromosomes  studies of 432 Trisomic fetuses and Liveborns. Hum Genet. 1994;94(4):411–7.
Oldereid NB, Wennerholm UB, Pinborg A, Loft A, Laivuori H, Petzold M, Romundstad LB, SoderstromAnttila V, Bergh C. The effect of paternal factors on perinatal and paediatric outcomes: a systematic review and metaanalysis. Hum Reprod Update. 2018;24(3):320–89.
de Graaf G, Buckley F, Skotko BG. Estimates of the live births, natural losses, and elective terminations with Down syndrome in the United States. Am J Med Genet A. 2015;167(4):756–67.
Papadopoulos G, Templeton AA, Fisk N, Randall J. The frequency of chromosomeanomalies in human preimplantation embryos after Invitro fertilization. Hum Reprod. 1989;4(1):91–8.
Jacobs M, Cooper SA, McGowan R, Nelson SM, Pell JP. Pregnancy outcome following prenatal diagnosis of chromosomal anomaly: a record linkage study of 26,261 pregnancies. PLoS One. 2016;11(12).
Acknowledgements
Not applicable.
Funding
This research was funded in part by a grant from the National Institutes of Health: R03CA208021. The National Institutes of Health had no role in the design of the study, collection, analysis, interpretation of data, and in writing the manuscript.
Availability of data and materials
The data that support the findings of this study are available from the National Center for Health Statistics, [https://www.cdc.gov/nchs/data_access/vitalstatsonline.htm].
Author information
Authors and Affiliations
Contributions
JAT performed all of the contributions to the manuscript including reading and approving the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was evaluated by the Texas A&M Institutional Review Board (IRB) and determined to be exempt from IRB review.
Consent for publication
Not applicable.
Competing interests
The author declares that he has no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Model 1. Maternal Age Randomwalk (CAR). OpenBUGS code and data that can be used to repeat the analyses for Model 1. (TXT 1 kb)
Additional file 2:
Model 2. Paternal Age Randomwalk (CAR). OpenBUGS code and data that can be used to repeat the analyses for Model 2. (TXT 2 kb)
Additional file 3:
Model 3. Fully Conditional Randomwalk (CAR). OpenBUGS code and data that can be used to repeat the analyses for Model 3. (TXT 67 kb)
Additional file 4:
5 year age categories. Fully conditional modeling of maternal and paternal ages as 5 year age categories. OpenBUGS code and data that can be used to repeat the analyses for the 5 year age categories. The data are arranged so that the DIC is comparable to the other fully conditional models (TXT 37 kb)
Additional file 5:
Linear. Fully conditional modeling of maternal and paternal ages as linear. OpenBUGS code and data that can be used to repeat the analyses for maternal and paternal ages as linear. (TXT 65 kb)
Additional file 6:
Curvilinear. Fully conditional modeling of maternal and paternal ages as curvilinear. OpenBUGS code and data that can be used to repeat the analyses for maternal and paternal ages as curvilinear. (TXT 66 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Thompson, J.A. Disentangling the roles of maternal and paternal age on birth prevalence of Down syndrome and other chromosomal disorders using a Bayesian modeling approach. BMC Med Res Methodol 19, 82 (2019). https://doi.org/10.1186/s1287401907201
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287401907201