Skip to main content

Sperm DNA methylation epimutation biomarker for paternal offspring autism susceptibility

Abstract

Background

Autism spectrum disorder (ASD) has increased over tenfold over the past several decades and appears predominantly associated with paternal transmission. Although genetics is anticipated to be a component of ASD etiology, environmental epigenetics is now also thought to be an important factor. Epigenetic alterations, such as DNA methylation, have been correlated with ASD. The current study was designed to identify a DNA methylation signature in sperm as a potential biomarker to identify paternal offspring autism susceptibility.

Methods and results

Sperm samples were obtained from fathers that have children with or without autism, and the sperm then assessed for alterations in DNA methylation. A genome-wide analysis (> 90%) for differential DNA methylation regions (DMRs) was used to identify DMRs in the sperm of fathers (n = 13) with autistic children in comparison with those (n = 13) without ASD children. The 805 DMR genomic features such as chromosomal location, CpG density and length of the DMRs were characterized. Genes associated with the DMRs were identified and found to be linked to previously known ASD genes, as well as other neurobiology-related genes. The potential sperm DMR biomarkers/diagnostic was validated with blinded test sets (n = 8–10) of individuals with an approximately 90% accuracy.

Conclusions

Observations demonstrate a highly significant set of 805 DMRs in sperm that can potentially act as a biomarker for paternal offspring autism susceptibility. Ancestral or early-life paternal exposures that alter germline epigenetics are anticipated to be a molecular component of ASD etiology.

Introduction

Autism spectrum disorder (ASD) is a complex neurological disorder involving deficits in communication, social behaviors and stereotypic movements [1, 2]. The prevalence of ASD in 1975 was reported as 1 in 5000 and then in 2009 as 1 in 110 [3]. The American Centers for Disease Control and Prevention reported a 1 in 88 prevalence in 2012 and then a 1 in 68 in 2014. Although improved diagnosis and current awareness have played a role in this increase, particularly in the first couple decades (1975–2000), the increase in the last two decades is thought to be due to environmental and molecular factors [2,3,4]. This is supported by twin studies and numerous environmental studies. Genetic studies using genome-wide association studies (GWAS) have identified multiple genetic mutations, but they are correlated with only a small percentage of the autism patients [5]. A recent study identified sperm genetic alterations associated with offspring autism [6]. Combining genetic mutations and altered epigenetics appear to improve associations [7]. Many specific toxicants and factors have been suggested to be involved, but generally more extensive analysis is required [8]. Environmental factors are now believed to be involved in the etiology of autism. A number of molecular alterations in the genome have been correlated to the neurobiology of ASD [2]; however, the specific environmental factors, molecular processes and etiology of autism remain to be fully elucidated.

Although there are both paternal and maternal transmission of ASD, the prevalence of paternal transmission is higher in most populations. One of the main factors proposed to be involved is paternal age [9], with an increased percentage risk of 28% between 40–49 years and nearly 70% when greater than 50 years of age [4]. Increased paternal age has been associated with epigenetic DNA methylation alterations in sperm [10], including specific genes associated with autism [11, 12]. Paternal age-associated DNA methylation alterations have been shown to impact offspring health and disease susceptibility [13, 14]. Therefore, the current study controlled for age at conception and sample collection for the comparison. In addition to paternal age effects, ancestral and early-life exposures to toxicants, abnormal nutrition and stress can also impact sperm DNA methylation to potentially affect disease susceptibility of offspring [15]. The current study was designed to examine the father’s sperm epigenetics (DNA methylation) in families with or without autistic children. The hypothesis examined is that a father’s specific sperm DNA methylation alterations will correlate with offspring autism susceptibility.

Epigenetics is defined as “molecular factors and processes around DNA that regulate genome activity independent of DNA sequence and are mitotically stable.” The molecular factors and processes currently known are DNA methylation, histone modifications, chromatin structural changes, noncoding RNA and RNA methylation [15]. When the epigenetic alterations become programmed in the germ cells (sperm or egg), they have the potential to promote in subsequent generations the epigenetic transgenerational inheritance of disease and phenotypic alterations [15]. Environmental factors that promote these early-life epigenetic alterations have the ability to promote epigenetic inheritance to subsequent generations and dramatically increase disease susceptibility and prevalence [15,16,17]. The current study was designed to use an epigenome-wide association study approach and develop a potential paternal sperm biomarker for offspring autism susceptibility.

The use of specific sperm epigenetic (DNA methylation) alterations (i.e., biomarkers) could be used for a fathers (i.e., paternal) offspring autism susceptibility, and applications in an assisted reproduction setting could be considered. Although genetic tests are common in assisted reproduction and preimplantation diagnostics, epigenetic analysis is less common. Sperm DNA methylation diagnostics have been proposed for the use in assisted reproduction [18]. The availability of a sperm DNA methylation biomarker for offspring autism susceptibility would allow improved clinical management and early treatment options to be considered. An epigenome-wide association study for DNA methylation alterations in sperm from fathers with or without autistic children was used to identify potential sperm epigenetic alterations as a biomarker for paternal offspring autism susceptibility.

Results

Paternal males with children affected by autism (case) or without (control) were recruited, and paternal sperm samples were collected at the Andrology Laboratory of IVIRMA Clinic in Valencia, Spain. The sperm sample was collected upon enrollment. Thirty-six patients were enrolled, which included thirteen in the control group, thirteen in the autism case group, and eight or ten for the blinded test groups. The differences (mean ± SD) between the semen analysis for both control and case group are shown in Table 1. Observations from the groups showed no significant difference in age, fathers age at pregnancy, fathers age upon sperm collection, sperm volume, concentration, or sperm concentration between the groups. Progressive sperm motility was greater in the autism case group, with no difference in non-progressive sperm motility, as shown in Table 1a. The motile percentage was higher in the control group, and no difference was observed in the total motile sperm count. One of the control samples, IVI 14, had a very high sperm count of 396.62 million that was outside two standard deviations of the mean (2 ± SD), so the analysis was redone without this sample. When the IVI 14 sample was not used in the analysis, the total sperm number was increased in the autism case group (p < 0.02), and the total motile sperm count (time) was increased in the autism case group (p < 0.017), as well as the progressive spermatozoa (%) (p < 0.019) and immotile % (p < 0.019) parameters. In addition to the case and control male age and sperm analysis parameters, Table 1a, all the blinded test set males, Table 1b, c, age and sperm parameters were analyzed and found to also be within the mean ± SD of the case and control samples presented (Additional file 2: Figure S1). Therefore, the blinded test set of individuals was appropriate comparisons with the same clinical parameters.

Table 1 Sperm samples and clinical semen analysis.

The participant demographics and clinical information were similar between the case and control population participants. The ethnicity of all the fathers was Caucasian. No major comorbidities were observed within either the control or case populations. The date of the patient sperm collection, age of the father upon collection and age of the father at conception of child are all not statistically different and provided in Table 1. Although age can impact sperm DNA methylation, the mean age upon sperm collection, which required a 3-year collection period, for the case and control was not statistically different, as given in Table 1. In addition, no statistical difference was observed in the age of the fathers at conception of child, as given in Table 1. All the autistic children were males. Since the focus was on paternal sperm, and due to IRB restrictions, the offspring ASD spectrum severity was not considered. The human subjects approval and informed consent were obtained from all participants prior to the initiation of the study and approved by the Ethics Committee of Valencian Infertility Institute—Reproductive Medicine Associates (IVIRMA) Valencia, Spain, with code, #1311-VLC-136-FC.

Individual patient sperm samples from the collection upon enrollment were prepared for sperm analysis, and an aliquot was taken and flash frozen with liquid nitrogen and stored at − 20 °C until shipment on dry ice for the epigenetic analysis. The samples were thawed, and prior to DNA isolation, the sperm were sonicated to destroy and remove any contaminating somatic cells, as previously described [16]. Due to the sperm nuclei being resistant to sonication, any contaminating somatic cells are removed following sonication. The DNA was extracted from the sperm and then fragmented for a methylated DNA immunoprecipitation (MeDIP) procedure to obtain methylated DNA for analysis to identify differential DNA methylated regions (DMRs). The MeDIP is a genome-wide analysis examining 95% of the genome, which is comprised of low-density CpG regions in comparison with the less than 5% of the genome of high-density CpG regions such as CpG islands. The MeDIP DNA libraries were prepared for next-generation DNA sequencing by creating individual patient sequencing libraries. Samples were then sequenced for bioinformatic analysis, as described in Additional file 1 section. A comparison of the sequences between the control (non-autism children) and case (autism children) participant sperm samples identified DMRs, as shown in Fig. 1a. At a p value of p < 1e−05 there were 805 DMRs identified with the majority being a single 1-kb window with fewer (i.e., six) having multiple adjacent 1-kb windows involved. The DMRs at a number of p values are presented for p < 001 to p < 1e−07, Fig. 1a. The DMRs at EdgeR p < 1e−05 all had false discovery rate (FDR) of p < 0.05 and were used for subsequent data analysis. The p < 1e−05 was used to optimize DMR numbers and statistical considerations. A list of these DMRs with various genomic features (e.g., CpG density and chromosomal location) are presented in Additional file 3: Table S1. Observations suggest that males with autistic children have a sperm DMR signature that is distinct from males without autistic children (control).

Fig. 1
figure 1

a DMR identifications. Autism case versus control sperm DMR analysis. The number of DMRs found using different p value cutoff thresholds. The all window column shows all DMRs. The multiple window column shows the number of DMRs containing at least two adjacent significant windows and the number of DMRs with each specific number of significant windows at a p value threshold of 1e−05. b Autism case versus control patient DMR analysis. The DMR locations on the individual chromosomes are identified. All DMRs at a p value threshold of p < 1e−05 are shown with the red arrowheads and clusters of DMRs with the black boxes. c DMR CpG density in the autism case versus control patient DMRs. The number of DMRs at different CpG densities is indicated. All DMRs at a p value threshold of p < 1e−05. d Autism case versus control patient DMR lengths in kilobases. All DMRs at a p value threshold of 1e−05 are shown

The genomic features of the offspring autism susceptibility DMRs were investigated. The chromosomal locations of the DMRs at p < 1e−05 within the human genome are presented in Fig. 1b. The red arrowheads indicate the individual DMRs, and the black boxes represent clusters of DMRs. The DMRs are present on all chromosomes. The CpG density of the DMRs is generally less than 10 CpG per 100 bp with 1–3 CpG predominant for the paternal offspring autism susceptibility DMRs, as shown in Fig. 1c. The size of the DMRs was predominantly 1–3 kb for the sperm DMRs, as shown in Fig. 1d. Additional genomic features are presented in Additional file 4: Table S1. The log-fold change (LFC) in DMA methylation in Additional file 3: Table S1 demonstrated for the 805 DMR in the autism group that 303 have an increase in DNA methylation and 502 have a decrease in DNA methylation. The autism DMRs involved a 38% increase or 62% decrease in DNA methylation. Therefore, the majority of the sperm DMRs had low CpG density, termed a CpG desert, and were 1 kb in length with either an increase or decrease in DNA methylation.

The paternal offspring autism susceptibility sperm DMR-associated genes and corresponding gene functional categories were determined, as presented in Additional file 4: Table S1. The total autism 805 DMRs had 193 with no DMR gene associations (24%), and the DMRs were intergenic and not associated with genes. From the 612 DMR with gene associations (76%), there were 493 DMR that overlapped with annotated genes. There were 17 DMR in the 1–1000 bp and 62 DMR in the 1–5 kb proximal promoter regions. There were 40 DMRs in the 5–10 kb distal promoter region. Therefore, approximately 20% of the DMR are in the proximal and distal promoter region and 80% overlapping the gene annotation regions. There were 193 DMRs that were intergenic. These DMRs are intergenic and not proximal to genes, but can influence gene expression events for megabase distances through ncRNA and chromatin structure alterations, as previously described [19, 20]. Genes within 10 kb of a DMR were identified, which has previously been shown to be optimal for both proximal and distal promoter regions and epigenetic associations [21]. The functional categories corresponding to each DMR-associated gene are summarized in Fig. 2a. The signaling, transcription and metabolism functional categories are predominant. This reflects that these gene functional categories have the highest number of genes within them. A comparison of previously identified genes associated with neurodegeneration, neurodevelopment and autism with the DMR-associated genes of this study is summarized in Fig. 2b. These autism-associated genes have previously been shown to be regulated or involve genetic mutations within autism patients, and the gene symbols, descriptions and associated references are presented in Additional file 5: Table S2. The DMR-associated genes were also used in a gene pathway or gene set analysis to identify associated pathways. Interestingly, the top pathway or gene set identified was autism, and the majority of the subsequent pathways with greater than three genes were all neurodevelopmental- or neurobiology-associated pathways, as given in Table 2. All those gene sets were found to be significant, and a list of the specific DMR-associated genes is provided, as given in Table 2. As shown with all the DMR-associated genes, Additional file 5: Table S2, the associated genes in Table 2 also had approximately a 50% mixture of genes with an increase or decrease in DNA methylation. Therefore, the DMR-associated genes did correlate well with previously identified autism- and neurodevelopment-associated genes. Since the sperm DMRs will impact the embryo epigenomes and transcriptomes of all subsequent somatic cells, this dynamic cascade of developmental epigenetics needs to be considered in potential links in sperm epigenetics and potential neurological impacts on autism.

Fig. 2
figure 2

a DMR-associated gene categories. DMRs at a p value threshold p < 1e−05 are shown. b Autism case versus control DMR PCA. PCA for DMRs at p < 1e−05. The first two principal components used and samples color code index indicated. The underlying data are the RPKM read depth for all DMRs. c DMR-associated genes and autism. The paternal offspring autism-susceptible DMRs previously shown to correlate with autism and associated neurodegenerative disease are presented. DMR-associated genes from the current study were compared to genes associated with autism in the published literature using Pathway Studio software (Elsevier, Inc.). Those that were in common are depicted

Table 2 DMR-associated gene pathways or gene sets.

The final analysis examined the statistical significance and validation of the DMRs identified. A permutation analysis was performed on the DMRs to determine whether the DMRs were due to background variation in the data or randomly generated. The permutation analysis shows the number of DMRs generated from the control versus autism case comparison was significantly greater than the DMRs generated from random subset analysis, Additional file 3: Figure S2. The red line to the right indicates the comparison DMRs versus the low numbers from the random subset comparison. Another analysis involved a cross-validation of the DMRs and demonstrated approximately 80% accuracy in the confirmation of the DMRs to assess autism susceptibility [22]. A principal component analysis (PCA) of the control male sperm without an autistic child versus the male sperm with an autistic child is presented in Fig. 2c. A clear separation of the DMR principal components is seen between the groups. This demonstrates a distinction between the DMR principal components.

The final validation involved using blinded test sets of samples for analysis to identify and assess the accuracy to determine actual case or control samples. Three different molecular analyses of the original 13 cases and 13 controls were performed and combined for this test set analysis. The test set analysis involved four independent analyses that were combined for the final analysis. The first test set involved eight blinded samples selected from the main analysis samples and reanalyzed (BS1–BS8) and identified as actual case or control sample, as given in Table 1b. All were accurately identified, except one false negative which was identified as a control, but actually was a case sample. A second set of ten blinded test samples (BS9–BS18) was provided by IVI-RMA clinical collaborators and was also identified in an independent analysis as case or control, as given in Table 1c. All were accurately identified, except one false negative which was identified as a control, but was actually a case sample. Therefore, the blinded test set analysis indicated all but one in each test set were accurately identified for approximately a 90% accuracy in the analysis. Since multiple analysis was used for this blinded test set analysis, random batch effect outlier DMRs identified were removed to optimize the analysis. Although significantly more validation with larger clinical test sets is needed, the current study provides the proof of concept that epigenetic biomarkers potentially exist and may be used to diagnose that a father may potentially have a child with a susceptibility for autism.

Discussion

The frequency of autism in the population has dramatically increased over tenfold the past several decades. This increase appears to be due in part to increased diagnosis efficiency from 1975 to the early 2000s, as well as greater public awareness of the disease [3]. The more recent increase in the last couple of decades suggests environmental factors, and exposures also have a critical role in autism prevalence. Although many suggestions have been made on specific toxicants and factors being involved, more extensive analysis and better understanding of autism etiology are needed to understand this increase in autism frequency [8]. An example is the suggestion-assisted reproduction and in vitro fertilization is involved, but follow-up studies demonstrated no risk of ASD in children born after assisted reproduction [23, 24]. One factor that has been correlated with autism is paternal age [4, 9, 15] and sperm DNA methylation alterations [11, 25]. Previous studies have shown a hypermethylation of sperm DNA is associated with male infertility, abnormal sperm parameters and increasing age [13, 26, 27]. Therefore, the majority of DMR involve an increase in DNA methylation when associated with infertility or age. The current study demonstrated 60% of the DMRs have a decrease in DNA methylation, and 40% of DMRs an increase in DNA methylation, Additional file 4: Table S1. Therefore, a mixture of an increase and decrease in methylation is observed, which is distinct from the sperm hypermethylation observed in male infertility and aging [13, 26, 27]. Since all the paternal subjects were similar in age and fertile (Table 1), the current study observations appear to be distinct from infertility and aging-related DNA hypermethylation. Although some participants from both control and case populations were involved in in vitro fertilization upon sperm collection, male factor infertility was not involved. No differences in demographics or clinical variables were observed. The age upon sperm collection between the case and control was not statistically different, as given in Table 1. In addition, the age at conception of child was also not statistically different between case and control participants. The comparison was biased on age of sperm collection to control for age differences to minimize DNA methylation variation. However, similar observations were also obtained considering age of conception of child. Therefore, the current study was designed to identify sperm epigenetic alterations (i.e., biomarkers) to assess a father’s potential ability to transmit autism susceptibility to his offspring.

Altered germline epigenetics has been shown to impact offspring health later in life, and if permanently programmed, to promote the epigenetic transgenerational inheritance of disease and pathology to subsequent generations [15, 16]. Since sperm or egg epigenetics can impact the zygote epigenetics and transcription following fertilization, as well as the subsequent stem cell population in the early embryo epigenome and transcriptome, all subsequently derived somatic cells also have the potential to have an altered cell-type specific epigenomes and transcriptomes later in development [15, 28]. This molecular alteration has been shown to be associated with adult somatic cell epigenetics, transcriptomes and associated diseases [29,30,31]. The ability of an ancestral or early-life exposure to impact the germline epigenetics to subsequently impact the offspring epigenetics and susceptibility to develop pathology and disease has been established [15,16,17, 28,29,30,31] and is anticipated to be a component of autism etiology as well. The current study observations support the concept that similar events may contribute to autism etiology.

The application of a sperm molecular diagnostic is optimally used in an assisted reproduction setting. Routine semen analysis and genetic testing are used in most in vitro fertilization clinical settings. Although epigenetic analysis is not as routine, the proposal for such analysis has been made [18]. The analysis of male infertility using sperm DNA methylation alterations has been developed [32]. Epigenetic alterations (DNA methylation) in sperm have been shown to associate in fathers of families with autistic children [33]. That study used a targeted array-based approach that focused on high-density CpG islands that constitutes approximately 1% of the genome, but does demonstrate such an analysis is feasible. The current study was designed to use a genome-wide approach to identify altered DNA methylation for paternal sperm and offspring autism susceptibility.

Although genetics will be involved in autism etiology, genome-wide association studies (GWAS) have demonstrated generally less than 1% of the patients with a specific disease, such as neurodegenerative disease that has a correlated genetic mutation [34]. ASD is similar to only a few percent correlation with associated genetic mutations [35]. An additional molecular mechanism to consider for ASD disease etiology involves epigenetics. The current study uses a more epigenome-wide association study approach to investigate sperm DNA methylation in fathers with or without autistic children. A procedure to assess DNA methylation alterations in low-density CpG regions, that constitute over 95% of the human genome, was used in comparison with the high-density CpG procedures previously used. A significant signature of differential DNA methylation regions (DMRs) was identified comparing the sperm from fathers with or without autistic children. The genomic features of the DMRs were identified and demonstrated generally 1-kb lengths and low-density CpG regions. The DMR-associated genes were identified, and a number of previously identified autism-linked genes were present (Fig. 2b, Additional file 5: Table S2 and Table 2). In regard to the autism sperm DMR biomarkers, a separation in a principal component analysis (PCA) was observed. In addition to this validation, the permutation and cross-validation analyses help demonstrate the robustness and sensitivity of the analysis. The validation studies with blinded sample sets accurately identified the majority of case and control samples, but potential false-negative identification of case samples was observed. The observations demonstrate the paternal sperm epigenetic analysis is potentially effective at identifying offspring susceptibility for autism, but the current analysis needs to be improved with expanded clinical trials.

Although an epigenetic signature was identified for paternal transmission of susceptibility of autism children, which was identified and statistically significant, a limitation of the current study is the low number of samples used for the analysis. Although epigenetic alterations occur at a significantly higher frequency than genetics, expanded clinical trials are required with increased numbers, greater ethnic diversity, and more thorough assessment of the impacts of paternal age. The impacts of these variables need to be elucidated to improve and expand the accuracy of the analysis. The expanded clinical trial with greater numbers and diverse subpopulations is essential to develop a useful diagnostic. However, the current study does provide the proof of concept; such a diagnostic can be developed.

Applications of the paternal offspring autism susceptibility biomarker/diagnostic will potentially improve the health care for ASD patients. This would allow IVF patients to assess risk and determine management procedures. Importantly, this would allow clinicians to plan the offspring’s clinical management options more efficiently. Potential preventative treatments could be considered to reduce the severity of the autism spectrum disorder. The availability of the assay could also be used in a research setting to facilitate the identification of environmental factors potentially involved in the ASD etiology. Therefore, potential therapeutic and preventative options not previously considered could be taken.

The current study identified a genome-wide signature of DNA methylation sites that are associated with the paternal transmission of offspring autism susceptibility. Although a large clinical trial is needed to further validate the biomarkers and potential diagnostic, the current study provides the proof of concept for the assay and biomarkers. Therefore, the identification of offspring susceptibility can be assessed, allowing better clinical management of ASD. The potential for therapy options can be expanded to improve health care for ASD. Such epigenetic biomarkers are anticipated to exist for many disease and pathology conditions, which will facilitate the future preventative medicine strategies for health care. In addition, the current study suggests epigenetic inheritance may play a role in ASD etiology and explain the paternal transmission prevalence of the disease.

Methods summary

Clinical sample collection

A single-center (IVIRMA Valencia, Spain) prospective and open clinical study was performed. The participant approval and informed consent were obtained from all participants prior to the clinical sample collection. The study protocols were approved by the Institutional Review Board Ethics Committee of Valencian Infertility Institute—Reproductive Medicine Associates (IVIRMA) Valencia, Spain, with code, #1311-VLC-136-FC. All research was performed in accordance with relevant guidelines/regulations. The study was not designed for, nor did the IRB involve, the ability to correlate autism child clinical information to be correlated. The semen was analyzed as described in Additional file 1. Samples were immersed in liquid nitrogen and then stored at −20 °C prior to analysis.

Epigenetic analysis, statistics and bioinformatics

Somatic cell contamination was removed by sonication, and the sperm DNA was isolated as previously described [16]. Methylated DNA immunoprecipitation (MeDIP), followed by next-generation sequencing (MeDIP-Seq), was performed. MeDIP-Seq, sequencing libraries, next-generation sequencing, and bioinformatics analysis were performed as described [16] and are found in Additional file 1. The statistical analysis and validation protocols were performed as previously described [16] and are found in Additional file 1. All molecular data has been deposited into the public database at NCBI (GEO # GSE157417), and R code computational tools are available at GitHub (https://github.com/skinnerlab/MeDIP-seq) and www.skinner.wsu.edu.

Availability of data and materials

All molecular data have been deposited into the public database at NCBI (GEO # GSE157417), and R code computational tools are available at GitHub (https://github.com/skinnerlab/MeDIP-seq) and www.skinner.wsu.edu.

Abbreviations

ASD:

Autism spectrum disorder

IVF:

In vitro fertilization

DMRs:

Differential DNA methylation regions

GWAS:

Genome-wide association studies

MeDIP:

Methylated DNA immunoprecipitation

MeDIP-Seq:

Methylated DNA immunoprecipitation followed by next-generation sequencing

FDR:

False discovery rate

LFC:

Log-fold change

PCA:

Principal component analysis

BS:

Blinded samples

References

  1. King BH, Navot N, Bernier R, Webb SJ. Update on diagnostic classification in autism. Curr Opin Psych. 2014;27(2):105–9.

    Article  Google Scholar 

  2. Kumar S, Reynolds K, Ji Y, Gu R, Rai S, Zhou CJ. Impaired neurodevelopmental pathways in autism spectrum disorder: a review of signaling mechanisms and crosstalk. J Neurodev Disord. 2019;11(1):10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Weintraub K. The prevalence puzzle: autism counts. Nature. 2011;479(7371):22–4.

    Article  CAS  PubMed  Google Scholar 

  4. Kimura R, Yoshizaki K, Osumi N. Risk of neurodevelopmental disease by paternal aging: a possible influence of epigenetic alteration in sperm. Adv Exp Med Biol. 2018;1012:75–81.

    Article  CAS  PubMed  Google Scholar 

  5. Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol Autism. 2017;8:21.

    Article  CAS  Google Scholar 

  6. Breuss MW, Antaki D, George RD, Kleiber M, James KN, Ball LL, et al. Autism risk in offspring can be assessed through quantification of male sperm mosaicism. Nat Med. 2020;26(1):143–50.

    Article  CAS  PubMed  Google Scholar 

  7. Massrali A, Brunel H, Hannon E, Wong C, Baron-Cohen S, et al. Integrated genetic and methylomic analyses identify shared biology between autism and autistic traits. Mol Autism. 2019;10:31.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Herbert MR. Contributions of the environment and environmentally vulnerable physiology to autism spectrum disorders. Curr Opin Neurol. 2010;23(2):103–10.

    Article  PubMed  Google Scholar 

  9. Yassin W, Kojima M, Owada K, Kuwabara H, Gonoi W, Aoki Y, et al. Paternal age contribution to brain white matter aberrations in autism spectrum disorder. Psychiatry Clin Neurosci. 2019;73(10):649–59.

    Article  CAS  PubMed  Google Scholar 

  10. Frans EM, Sandin S, Reichenberg A, Langstrom N, Lichtenstein P, McGrath JJ, et al. Autism risk across generations: a population-based study of advancing grandpaternal and paternal age. JAMA Psych. 2013;70(5):516–21.

    Article  Google Scholar 

  11. Atsem S, Reichenbach J, Potabattula R, Dittrich M, Nava C, Depienne C, et al. Paternal age effects on sperm FOXK1 and KCNA7 methylation and transmission into the next generation. Hum Mol Genet. 2016;25(22):4996–5005.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Milekic MH, Xin Y, O’Donnell A, Kumar KK, Bradley-Moore M, Malaspina D, et al. Age-related sperm DNA methylation changes are transmitted to offspring and associated with abnormal behavior and dysregulated gene expression. Mol Psych. 2015;20(8):995–1001.

    Article  CAS  Google Scholar 

  13. Jenkins TG, Aston KI, Pflueger C, Cairns BR, Carrell DT. Age-associated sperm DNA methylation alterations: possible implications in offspring disease susceptibility. PLoS Genet. 2014;10(7):e1004458.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Rumbold AR, Sevoyan A, Oswald TK, Fernandez RC, Davies MJ, Moore VM. Impact of male factor infertility on offspring health and development. Fertil Steril. 2019;111(6):1047–53.

    Article  PubMed  Google Scholar 

  15. Nilsson E, Sadler-Riggleman I, Skinner MK. Environmentally induced epigenetic transgenerational inheritance of disease. Environ Epigenet. 2018;4(2):1–13.

    Article  CAS  Google Scholar 

  16. King SE, McBirney M, Beck D, Sadler-Riggleman I, Nilsson E, Skinner MK. Sperm epimutation biomarkers of obesity and pathologies following DDT induced epigenetic transgenerational inheritance of disease. Environ Epigenet. 2019;5(2):1–15.

    Article  CAS  Google Scholar 

  17. Anway MD, Cupp AS, Uzumcu M, Skinner MK. Epigenetic transgenerational actions of endocrine disruptors and male fertility. Science. 2005;308(5727):1466–9.

    Article  CAS  PubMed  Google Scholar 

  18. Simon L, Emery BR, Carrell DT. Review: diagnosis and impact of sperm DNA alterations in assisted reproduction. Best Pract Res Clin Obstet Gynaecol. 2017;44:38–56.

    Article  PubMed  Google Scholar 

  19. Haque MM, Nilsson EE, Holder LB, Skinner MK. Genomic Clustering of differential DNA methylated regions (epimutations) associated with the epigenetic transgenerational inheritance of disease and phenotypic variation. BMC Genom. 2016;17(418):1–13.

    Google Scholar 

  20. Skinner MK, Manikkam M, Haque MM, Zhang B, Savenkova M. Epigenetic transgenerational inheritance of somatic transcriptomes and epigenetic control regions. Genome Biol. 2012;13(10):R91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Lehman DA, Patterson B, Johnston LA, Balzer T, Britton JS, Saint R, et al. Cis-regulatory elements of the mitotic regulator, string/Cdc25. Development. 1999;126(9):1793–803.

    CAS  PubMed  Google Scholar 

  22. Hastie T, Tibshirani R, Friedman J. The elements of statistical learning: data mining, inference, and prediction. 2nd ed. New York: Springer; 2009.

    Book  Google Scholar 

  23. Hvidtjorn D, Grove J, Schendel D, Schieve LA, Svaerke C, Ernst E, et al. Risk of autism spectrum disorders in children born after assisted conception: a population-based follow-up study. J Epidemiol Commun Health. 2011;65(6):497–502.

    Article  CAS  Google Scholar 

  24. Sandin S, Nygren KG, Iliadou A, Hultman CM, Reichenberg A. Autism and mental retardation among offspring born after in vitro fertilization. JAMA. 2013;310(1):75–84.

    Article  CAS  PubMed  Google Scholar 

  25. Siu MT, Weksberg R. Epigenetics of autism spectrum disorder. Adv Exp Med Biol. 2017;978:63–90.

    Article  CAS  PubMed  Google Scholar 

  26. Laqqan M, Solomayer EF, Hammadeh M. Aberrations in sperm DNA methylation patterns are associated with abnormalities in semen parameters of subfertile males. Reprod Biol. 2017;17(3):246–51.

    Article  PubMed  Google Scholar 

  27. Alkhaled Y, Laqqan M, Tierling S, Lo Porto C, Hammadeh ME. DNA methylation level of spermatozoa from subfertile and proven fertile and its relation to standard sperm parameters. Andrologia. 2018;50(6):e13011.

    Article  CAS  PubMed  Google Scholar 

  28. Jirtle RL, Skinner MK. Environmental epigenomics and disease susceptibility. Nat Rev Genet. 2007;8(4):253–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Guerrero-Bosagna C, Savenkova M, Haque MM, Nilsson E, Skinner MK. Environmentally induced epigenetic transgenerational inheritance of altered sertoli cell transcriptome and epigenome: molecular etiology of male infertility. PLoS ONE. 2013;8(3):1–12.

    Article  CAS  Google Scholar 

  30. Nilsson E, Klukovich R, Sadler-Riggleman I, Beck D, Xie Y, Yan W, et al. Environmental toxicant induced epigenetic transgenerational inheritance of ovarian pathology and granulosa cell epigenome and transcriptome alterations: ancestral origins of polycystic ovarian syndrome and primary ovarian insufiency. Epigenetics. 2018;13(8):875–95.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Klukovich R, Nilsson E, Sadler-Riggleman I, Beck D, Xie Y, Yan W, et al. Environmental toxicant induced epigenetic transgenerational inheritance of prostate pathology and stromal-epithelial cell epigenome and transcriptome alterations: ancestral origins of prostate disease. Sci Rep. 2019;9(2209):1–17.

    CAS  Google Scholar 

  32. Aston KI, Uren PJ, Jenkins TG, Horsager A, Cairns BR, Smith AD, et al. Aberrant sperm DNA methylation predicts male fertility status and embryo quality. Fertil Steril. 2015;104(6):1388–97.

    Article  CAS  PubMed  Google Scholar 

  33. Feinberg JI, Bakulski KM, Jaffe AE, Tryggvadottir R, Brown SC, Goldman LR, et al. Paternal sperm DNA methylation associated with early signs of autism risk in an autism-enriched cohort. Int J Epidemiol. 2015;44(4):1199–210.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Visscher PM, Brown MA, McCarthy MI, Yang J. Five years of GWAS discovery. Am J Hum Genet. 2012;90(1):7–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kim YS, Leventhal BL. Genetic epidemiology and insights into interactive genetic and environmental effects in autism spectrum disorders. Biol Psychiat. 2015;77(1):66–74.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge Ms. Jayleana Barton and for technical assistance and all technicians from IVIRMA Andrology Laboratory for assistance in the clinical sample collection. We acknowledge Ms. Amanda Quilty for editing and Ms. Heather Johnson for assistance in preparation of the manuscript. We thank the Genomics Core laboratory at WSU Spokane for sequencing data. We acknowledge Epigenesys Inc., Pullman, WA, where the initial molecular analysis was performed. This study was supported by Epigenesys Inc., Pullman, WA, USA, and John Templeton Foundation (50183 and 61174) (https://templeton.org/) grants to MKS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Funding

This study was supported by Epigenesys Inc., Pullman, WA, USA, and John Templeton Foundation (50183and61174) (https://templeton.org/) grants to MKS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

NG performed clinical and sample collection oversight, data analysis and editing manuscript. FC was involved in patient’s recruitment and editing manuscript. RRE contributed to patient’s recruitment, clinical and sample collection oversight, data analysis and editing manuscript. CS conceived, initiated sample collection and edited the manuscript. ISR was involved in molecular analysis, data analysis and editing manuscript. DB performed bioinformatics, data analysis and editing manuscript. EN done sample processing, data analysis and editing manuscript. MBM was involved in sample processing, molecular analysis, data analysis and editing manuscript. MKS conceived, analyzed the data, performed funding acquisition and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Michael K. Skinner.

Ethics declarations

Ethics approval and consent to participate

The participant approval and informed consent were obtained from all participants prior to the clinical sample collection. The IRB study was approved by the Ethics Committee of Valencian Infertility Institute—Reproductive Medicine Associates (IVIRMA) Valencia, Spain, with code, #1311-VLC-136-FC. All research was performed in accordance with relevant guidelines/regulations.

Consent for publication

Not applicable.

Competing interests

The authors report no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplemental methods.

Additional file 2: Figure S1

. Clinical group statistic comparison. The various sperm/semen characteristics in the Study case and control group were compared with the Blind group. The n-value, mean, standard deviation, and standard error mean are presented. The blind groups are within the mean ± SD of the case and control study.

Additional file 3: Figure S2

. Permutation analysis. The number of DMR for autism case versus control patient comparison for all permutation analyses. The vertical red line shows the number of DMR found in the original analysis. All DMRs are defined using an edgeR p value threshold of p < 1e−05.

Additional file 4: Table S1

. DMR lists at p < 1e−05 with presentation of name, chromosomal location, DMR start and stop nucleotide number for chromosome, length (bp), number of 1 kb significant windows, minimum p value, CpG number and density, DMR maximum, log-fold change (maxLFC) (+ increase DNA methylation and − decrease DNA methylation), and gene association within 10 kb and gene functional categories.

Additional file 5: Table S2

. Gene and protein regulators of autism. The DMR-associated autism-related genes with symbol, description, and relevant references PubMed ID (PMID) numbers.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Garrido, N., Cruz, F., Egea, R.R. et al. Sperm DNA methylation epimutation biomarker for paternal offspring autism susceptibility. Clin Epigenet 13, 6 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s13148-020-00995-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13148-020-00995-2

Keywords