- Open Access
A profile and review of findings from the Early Markers for Autism study: unique contributions from a population-based case–control study in California
Molecular Autism volume 12, Article number: 24 (2021)
The Early Markers for Autism (EMA) study is a population-based case–control study designed to learn more about early biologic processes involved in ASD.
Participants were drawn from Southern California births from 2000 to 2003 with archived prenatal and neonatal screening specimens. Across two phases, children with ASD (n = 629) and intellectual disability without ASD (ID, n = 230) were ascertained from the California Department of Developmental Services (DDS), with diagnoses confirmed according to DSM-IV-TR criteria based on expert clinical review of abstracted records. General population controls (GP, n = 599) were randomly sampled from birth certificate files and matched to ASD cases by sex, birth month and year after excluding individuals with DDS records. EMA has published over 20 papers examining immune markers, endogenous hormones, environmental chemicals, and genetic factors in association with ASD and ID. This review summarizes the results across these studies, as well as the EMA study design and future directions.
EMA enabled several key contributions to the literature, including the examination of biomarker levels in biospecimens prospectively collected during critical windows of neurodevelopment. Key findings from EMA include demonstration of elevated cytokine and chemokine levels in maternal mid-pregnancy serum samples in association with ASD, as well as aberrations in other immune marker levels; suggestions of increased odds of ASD with prenatal exposure to certain endocrine disrupting chemicals, though not in mixture analyses; and demonstration of maternal and fetal genetic influence on prenatal chemical, and maternal and neonatal immune marker and vitamin D levels. We also observed an overall lack of association with ASD and measured maternal and neonatal vitamin D, mercury, and brain-derived neurotrophic factor (BDNF) levels.
Covariate and outcome data were limited to information in Vital Statistics and DDS records. As a study based in Southern California, generalizability for certain environmental exposures may be reduced.
Results across EMA studies support the importance of the prenatal and neonatal periods in ASD etiology, and provide evidence for the role of the maternal immune response during pregnancy. Future directions for EMA, and the field of ASD in general, include interrogation of mechanistic pathways and examination of combined effects of exposures.
Motivation for the EMA study
Autism spectrum disorder (ASD) is a complex neurodevelopmental condition characterized by deficits in social communication and the presence of restricted behaviors or interests . Decades of research have firmly established the importance of genetics in the etiology of ASD [2, 3], although that role is a complex one that is still not fully understood. There is also evidence of substantial involvement of environmental factors in ASD etiology, particularly those acting during the prenatal and early neonatal periods as critical windows of neurodevelopment and susceptibility to environmental influence. Though the pathways to ASD are not fully understood, a growing body of evidence supports immune system dysregulation in risk of ASD [4,5,6,7]. Furthermore, fetal exposure to chemicals may also impact neurodevelopment, perhaps through immune pathways [8,9,10]. However, the mechanisms by which environmental agents interact with the developing immune and nervous systems to cause neurodevelopmental disorders in genetically susceptible individuals are poorly understood. Motivated by these gaps in our knowledge, and the desire to learn more about early biologic processes involved in ASD, the Early Markers for Autism (EMA) study was created in 2004. This review describes the motivation for the EMA study, the methods of this novel population-based case control project, and summarizes major findings published (or accepted) as of October 2020 that serve as unique contributions to the field of ASD epidemiology.
In the early part of the decade from 2000 to 2010, few studies had directly investigated the role of prenatal immune profiles (e.g., cytokine and chemokine profiles, antibodies to fetal brain proteins, antibodies to specific infectious agents) on neurodevelopmental outcomes in humans [11, 12]. However, findings of immune conditions in the family members of affected individuals, and aberrations in immune markers in both cases and mothers were motivating mechanistic work in this area [13,14,15,16,17,18,19,20,21,22,23,24,25,26]. Also during this time frame, no ASD studies had been conducted measuring biomarker levels of environmental contaminants in samples collected during the prenatal or early neonatal susceptibility periods. Research in related fields suggested developmental concerns for a number of environmental chemicals, including polybrominated diphenyl ethers (PBDEs), polychlorinated biphenyls (PCBs), and certain pesticides , but work examining these exposures using prenatal biomarkers for associations with neurodevelopmental outcomes more specifically was truly in its infancy. There was also a gap in knowledge regarding the potential mechanistic pathways of such exposures. When created in 2004, EMA was the first large, population-based, nested case–control study of ASD to utilize archived prenatal and newborn specimens mother-baby pairs. Key features of the EMA study are provided in Table 1.
Goals and opportunities
The EMA Study sought to advance our knowledge of the effects of immunologic, genetic, and environmental factors, individually and in combination, on ASD etiology, utilizing biomarkers of exposures and examining immune and endocrine pathways. With a large, essentially population-based sample of children with confirmed ASD, children with other developmental disabilities, as well as unaffected children, we had an unprecedented opportunity to explore the interplay between genetic susceptibility and markers of immune function and environmental exposures that contribute to ASD risk. Furthermore, because banked biospecimens were being used, the design enabled us to prospectively investigate measured levels of these factors during critical periods of fetal brain development in a novel way. We focused on exposures and markers that were known or suspected to play a role in early brain development. These included antibodies to viral infections, cytokines, chemokines, antibodies to fetal brain proteins, neurotrophins, hormones, metals, and endocrine disrupting chemicals (EDCs) that could be measured in these archived maternal blood (sera and cell pellets) collected in mid-pregnancy and neonatal blood spots collected in the first few days of life. We anticipated that findings from this study would define areas for further investigations of physiologic mechanism, provide a better understanding of the underlying biology involved in ASD etiology, and provide normative data on unaffected controls that may be useful as benchmarks for future studies utilizing prenatal and neonatal specimens.
The EMA study framework is depicted in Fig. 1. The EMA study population was drawn from a 2000–2003 cohort of births in Southern California in which mother’s prenatal screening specimens were archived as part of Project Baby’s Breath (M Kharrazi and GN DeLorenze, Co-PIs) . Women who participated in the prenatal expanded alphafetoprotein screening program (XAFP),  offered to all pregnant women in California in the second trimester of pregnancy, during 1999–2002 in Orange, Imperial, and San Diego Counties were eligible for participation in Project Baby’s Breath; in addition to this participation in screening, for inclusion in EMA, women who had a live born infant for whom a California birth certificate match could be found (approximately 90% of screening participants), and for whom both a maternal prenatal screening blood specimen and a newborn screening blood specimen were available in the State archives, were eligible. EMA was conducted in two phases; the first phase was drawn from women who were pregnant in Orange County, California, and who delivered a live born infant from July 2000 to September 2001. The second phase was drawn from women who were pregnant in Orange, San Diego, or Imperial Counties in Southern California and delivered a live born infant from January 2000 to June 2003, after excluding those in the first phase, thereby serving to expand the sample size, geographic region, and birth year range of the first phase. This geographic region included urban, suburban and rural areas in the southern part of California closest to the Mexico border with diverse geographies and race/ethnicities. Study activities were approved by the Institutional Review Board at Kaiser Permanente and the State Committee for the Protection of Human Subjects. Informed consent was not required for the EMA study, given that consent forms for screening programs were distributed at the time of collection, which stipulated that specimens and results from prenatal and neonatal testing could be used for research purposes given IRB approval.
Representativeness: Approximately 70% of pregnant women participated in the XAFP in CA during this time period [29, 30]. Screened women were more likely to be under age 35, but were similar to the total population of women delivering in study counties for these birth years with respect to maternal education, race/ethnicity, parity, and health insurance status (Additional file 1: Table S1). Basic characteristics of the EMA study population, by phase and study group, are provided in Table 2.
Group ascertainment: Three groups of children born to women in the cohort were identified: children with ASD, children with intellectual disability without ASD (ID; formerly classified as mental retardation or developmental delay in some EMA publications), and general population (GP) controls. Individuals with ID of known etiology (e.g., Down’s syndrome) were excluded. Further definition of these groups is provided below in diagnostic verification. Children with ASD or ID were ascertained from the California Department of Developmental Services (DDS), which operates a system of 21 Regional Centers (RC) that coordinate services for persons with ASD and other developmental disabilities. Services are provided without regard to citizenship or financial status; thus, the system is widely utilized across socio-economic levels and racial/ethnic groups. DDS maintains a centralized electronic database of diagnostic and identifying information on all clients, derived from information provided by the RCs on a standardized, statewide reporting form (Client Development Evaluation Report). Referrals to the RCs come from pediatricians and other clinical providers, the education system, friends, and family members. DDS and the RCs are mandated to provide services to individuals with autistic disorder and children with other Pervasive Developmental Disorders (PDDs) who have ID (IQ < 70) or are substantially impaired. Children with milder forms of ASD were less likely to be captured in these records. GP controls were randomly sampled from the birth certificate files, frequency matched to ASD cases by sex, birth month and birth year, after excluding all past or current DDS/RC clients. Controls were originally sampled to ensure a 2:1 control-ASD case ratio after retrieval of maternal and neonatal specimens. Among those with live births in the study counties and birth years, ~ 50% across all study groups linked to stored specimens. (This proportion is lower than the overall participation rate in the XAFP due to variability in participation rates over time and region, and ability to link to stored specimens). At the time of case and control ascertainment, children were approximately 3–4 years old for phase 1 and 4.5–9 years for phase 2.
Diagnostic verification and case classifications
Following a protocol initially developed by the CDC, trained medical record abstractors reviewed and abstracted detailed diagnostic and clinical data from RC records for all children identified as ASD and ID for whom archived blood specimens were found. Final case status was determined by expert clinical review of the abstracted information, according to DSM-IV-TR criteria. The overall EMA participant selection process is shown in Fig. 1. Most of those ascertained as ASD cases were confirmed (72% for phase 1 and 98% for phase 2), while 22% (phase 1) and 38% (phase 2) of those originally identified as ID were reclassified as ASD cases according to expert clinician review of records. After this expert review, final analytic samples consisted of 84 children with ASD and 49 with ID for phase 1 and 545 children with ASD and 181 children with ID for phase 2. Certain analyses across phases also considered ASD subtypes, including early onset or regressive (phase 1), and presence or absence of co-occurring ID (phase 1 and 2). Determination of ID (and co-occurring ID within the presence of ASD) was based on composite scores < 70 on all standardized cognitive and functional tests recorded in RC records (with no ID defined according to all or some scores > 70, and unknown ID defined according to the absence of standardized scores in chart).
Data were obtained through California Vital Statistics, the GDSP, and the DDS. Parental demographics, including parental age, dates of birth, education level, race, ethnicity, country of origin, and insurance status at delivery were obtained through birth certificates. Additional maternal information, including weight, height, and pregnancy-related characteristics were obtained at time of prenatal sample collection. Infant death and neonatal information, including birthweight and gestational age, was obtained through Vital Statistics. While covariates included in adjusted models varied across studies, study matching factors, maternal age, weight at prenatal specimen collection, and race/ethnicity were typically adjusted for in analyses; additional adjustment factors were considered as appropriate for analyses (see also Table 4).
Retrieval of biospecimens
Maternal mid-pregnancy (range 15–20 weeks gestation) blood specimens were obtained from Project Baby’s Breath and dried newborn bloodspots were obtained from the CA Genetic Disease Screening program (GDSP) of the California Biobank Program (https://www.cdph.ca.gov/Programs/CFH/DGDS/Pages/cbp/default.aspx). All biospecimens were retrieved from -20 C freezer archives, barcoded with a unique study identification number, and shipped to the participating laboratories. For all subjects, one whole dried bloodspot, 1–2 ml (mean 1.8 ml, SD 0.3 ml) of maternal serum, and 1-2 ml of maternal blood cell pellet were retrieved.
Laboratory and exposure analyses
Biomarkers measured in specimens from EMA participants are shown in Table 3 (see also Additional file 1: Appendix). Laboratory procedures for analytes have been previously published. We refer the reader to individual papers for more detailed information on methods for a given biomarker; those examined in multiple EMA studies are briefly described here. Cytokines and chemokines were measured in maternal serum and newborn samples using commercially available multiplex bead-based kits, with assays carried out following manufacturer protocols [31, 32]; TSH was measured in regional laboratories as part of routine screening in maternal serum using an I-labeled immunoradio-metric assay and newborn bloodspots using solid-phase, time-resolved sandwich fluoroimmunoassay [33, 34]. Vitamin D was measured as total 25-hydroxyvitamin D (sum of 25(OH)D2 and 25(OH)D3, hereinafter referred to as 25(OH)D) by a sensitive isotope dilution liquid chromatography-tandem mass spectrometry method (LC/MS/MS) in both maternal serum and dried newborn blood spots (converted to serum equivalents) [35,36,37]. Endocrine disrupting chemicals (PCBs, OCPs, PBDEs, and PFASs) were measured in maternal serum samples using gas chromatography isotope dilution high resolution mass spectrometry, [38,39,40].
Air pollution in EMA was assigned as averaged monthly exposure for the month prior to maternal serum collection based on the Environmental Protection Agency’s Air Quality System data (http://www.epa.gov/ttn/airs/airsaqs) for nitrogen dioxide (NO2), ozone, and particulate matter less than 10 and 2.5 microns in diameter (PM10, PM2.5) using the maternal residential address reported during the prenatal screening visit.
DNA was extracted from a subset of 790 maternal cell pellets and 764 neonatal bloodspots and 675,000 SNPs were genotyped across the genome using standard protocols and quality assessment, as previously described . The genotyped samples included 366 ASD mother–child pairs and 369 control pairs [42, 43]. The ID cases and ASD cases identified through reclassification of those originally enrolled as ID were not included in the genome-wide genotyping due to budget limitations, but were used for analyzing specific candidate genetic factors .
Given that EMA is comprised of mother–child pairs of cases and unrelated controls, we were able to take advantage of some features of family-based analyses. We considered models where the maternal genetic make-up could act as a determinant of the fetal environment, impacting the early risk for ASD. We tested for maternal main effects influencing risk for ASD, as well as maternal-child genetic combinations or interactions that might impact risk . With genotyping data, we also estimated copy number variants (CNVs) and tested for maternal effects that might be mediated through CNVs . Combining EMA genetic and biomarker data we conducted G × E analyses. Previous studies using genome-wide association study (GWAS) data to identify G × E interactions encountered statistical power issues compared with the detection of main effects [46, 47]. We hypothesized that the most likely genetic candidates for G × E would be the variants associated with the E biomarker itself. Our strategy was thus to perform a quantitative trait locus (QTL) analysis to identify candidate genetic factors that contribute to the inter-individual variation of the E biomarker and then test for interaction with case–control status among those top candidates [42, 43]. Given the variety of genetic ancestries in the EMA study population, methods appropriate for mixed-ancestry data were used to avoid ancestry-related false positives, such as adjustment for genome-wide principal components . In addition, because mothers and their affected children are genetically correlated with each other, we performed separate GWAS to analyze distinct maternal and neonatal effects. When we observed both maternal and fetal genetic contribution for the same biomarker, we identified the offspring alleles not inherited from the maternal lineage (fetal-specific alleles) and the maternal alleles not transmitted to the offspring (maternal-specific alleles) in order to distinguish whether both maternal and fetal genomes independently affected a specific biomarker .
We also applied the principles of Mendelian randomization to assess the possible causal relationship of a risk trait with ASD outcome. If a trait (i.e. immune biomarker or vitamin D) is causally linked to ASD, the genetic variants influencing the trait should also influence ASD risk. We thus tested the association between the mapped genetic determinants of the trait in the EMA dataset and ASD outcome [42, 43, 49].
Major EMA findings
A summary of EMA findings for associations with ASD from publications through October 2020 is provided in Table 4. Brief descriptions of key findings within each of the major target areas of interest are further provided below.
EMA has contributed 7 studies focused on the role of immune markers in ASD, many of these among the first to examine such markers in prospectively collected prenatal and neonatal biospecimens. Summaries of adjusted associations between ASD and immune markers in maternal and neonatal specimens are provided in Fig. 2a, b, respectively. Immune markers measured in EMA, by phase, are shown in these figures; those measured in phase 2 (which expanded phase 1 measures) are also listed in the Appendix. Results generally supported associations in the same direction across Phases (with the exception of neonatal eotaxin and IL-10), although the samller sample size in Phase 1 yielded less precise estimates. Across these studies, there was a general pattern of moderate increases in ASD risk with inflammatory cytokines; additional patterns are highlighted further below.
One of EMA’s first publications demonstrated a relationship between maternal antibodies to fetal brain proteins measured in mid-pregnancy samples and ASD , a result consistent with earlier work suggesting similar findings in maternal samples collected after the birth of the child, as well as in animal models [51,52,53]. Another investigation in phase 1 suggested lower levels of the immunoglobulin IgG to Toxoplasma gondii in both maternal and neonatal samples were associated with increased risk of ASD, but no differences were noted in the ID group . In addition, work conducted in the phase 2 sample suggested lower levels of C-reactive protein (CRP, acute phase protein produced in response to inflammation) in maternal samples of ASD cases, or significantly reduced risk of ASD with higher levels of CRP . Together, these results could suggest an impaired response to immune activation in mothers of ASD cases.
Examination of cytokine profiles from phase 1  demonstrated differences potentially suggestive of an asthma/allergy phenotype in mothers of ASD children, with elevated concentrations of IFN-γ, IL-4 and IL-5. A more inflammatory phenotype was found in mothers of ID children, characterized by elevated concentrations of IL-2, IL-4, and IL-6. This work was expanded to examine a larger panel of cytokines and chemokines levels in a greater number of maternal samples from phase 2, which demonstrated elevated maternal levels of several proinflammatory cytokines and chemokines in association with ASD with ID . This differentiated the ASD with ID group from all other groups including ASD without co-occurring ID, ID without ASD, and GP controls. When neonatal samples were similarly examined, higher levels of inflammatory cytokines were found in ASD cases relative to GP controls, and other inflammatory T cell cytokines and certain chemokines were elevated in ASD relative to ID .
Results across these studies of immune markers therefore suggest immune dysregulation, in some cases consistent with pro-inflammatory processes, during the prenatal and neonatal periods may be involved in the development of ASD. In addition, a number of these studies suggested these relationships may differ based on the presence or absence of co-occurring ID (Additional file 1: Table S2).
Other endogenous biomarkers
Thyroid hormones: Thyroid hormones play a critical role in brain development in regulating neurogenesis, dendrite proliferation, synaptogenesis, and myelination [57,58,59]. Furthermore, hypothyroidism early in life can lead to intellectual disability, and maternal thyroid conditions have been linked with ASD in prior work [60,61,62,63,64]. In order to address the role of thyroid hormones in ASD, thyroid stimulating hormone (TSH) was examined in maternal samples in phase 1, and neonatal bloodspot TSH in phases 1 and 2 [33, 34]. In these studies we observed inverse relationships between ASD and maternal TSH, as well as early onset ASD specifically . Relationships between ASD and neonatal TSH were not significant, though there was a suggestive inverse relationship with ASD with developmental regression in the phase 2 sample .
Vitamin D: Vitamin D is a steroid molecule whose primary source is sunlight. Neonatal levels depend on maternal levels, and deficiency has been linked with a range of maternal and child health outcomes , including neurodevelopmental outcomes such as cognitive function, language and motor development [66, 67]. Vitamin D plays a well-documented role in immune functioning, including the ability to influence growth and differentiation of innate and acquired immune cells, influence cytokine production, act as an anti-inflammatory agent, and regulate placental inflammation during pregnancy [68, 69]. Vitamin D also plays a role in protecting against oxidative stress, and can act as a neurosteroid in neurodevelopment , providing a basis for potential associations with ASD. Thus, for both maternal serum and neonatal vitamin D levels, we hypothesized potential protective effects with higher 25(OH)D levels, or increased risk with vitamin D deficiency . However, primary analyses considering vitamin D deficiency (defined as < 50 nmol/L total 25(OH)D; 14% of newborns and 9% of mothers), insufficiency (defined as 50–74 nmol/L; 26% newborns and 26% of mothers), or according to continuous 25(OH)D levels revealed no significant associations overall at either time point [37, 71]. We also did not observe different associations with ID, or ASD with or without comorbid ID (Additional file 1: Table S2). However, for both maternal and neonatal 25(OH)D levels, there were interactions with child sex and race/ethnicity, such that males and non-Hispanic whites had lower odds of ASD with higher vitamin D levels. In addition, a non-linear, inverted j-shape pattern was observed overall, and among males, between maternal levels and ASD, with the peak around 100 nmol/L 25(OH)D. This pattern was not observed for ID. Of note, 25(OH)D levels were relatively high in EMA (median ~ 85 nmol/L), as may be expected for a study based in relatively lower latitudes (in an area with higher sunlight) .
EDCs: Endocrine disrupting chemicals have been hypothesized to influence risk of ASD given evidence for associations with broader adverse neurodevelopmental outcomes. Potential links between these factors are supported by evidence of EDC disruption of thyroid-hormones [73, 74], influence on immune markers, or direct influence on developing neurons [75,76,77,78]. Three published analyses have examined classes of organohalogen endocrine disrupting chemicals in association with ASD and ID in EMA [38,39,40], and one study has examined these associations considering all chemicals simultaneously . Results from individual chemical analyses suggested significant associations with 2 PCB congeners, and a general trend of increases in odds of ASD with levels in the highest quartiles of most PCBs, but no association between ASD and the two organochlorine pesticides (OCP) detected (p,p-D’DE and trans-nonachlor) . No sex differences or differences by ID comorbidity were noted, and the general pattern of associations with PCBs was similar when examining ID as the outcome rather than ASD, though a significant, non-monotonic association was noted between the OCP p’p’-DDE and ID. In contrast, analyses of PBDEs demonstrated null or inverse associations with ASD and ID . When stratified by sex, associations with higher levels of PBDEs were below the null for boys but above for girls, suggesting potential sexual dimorphism (Additional file 1: Table S2). Analyses investigating more mechanistically-defined subgroupings of PCBs and PBDEs (such as dioxin-like or cytokine-P450 inducing PCBs, or highly brominated PBDEs), did not yield differing findings, and thus results did not directly support a specific pathway. For PFASs, associations were generally null, though inverse associations were found for the highest quartile of perfluorooctanoate (PFOA) perfluorooctane sulfonate (PFOS) and ASD . Results were similar for ID and ASD with and without ID, and stratified by sex. Assessing all these chemicals together in single Bayesian model for associations with ASD, no significant associations were found . ORs for the 25 included chemicals were all close to the null and attenuated when compared to analyses based on individual chemicals (Additional file 1: Fig. S1).
Air pollution: Air pollution may impact neurodevelopment via inflammatory or oxidative stress pathways, influences on methylation, or via direct effects of specific components on the developing brain [80, 81]. We had two goals in assessing air pollution in EMA: 1) to examine the relationship between maternal cytokine and chemokine levels measured during mid-pregnancy and previous-month (short term) air pollution exposure; and 2) to determine if immune markers measured mid-pregnancy mediate an association between mid-pregnancy air pollutant exposure and child ASD and ID. We found that previous month air pollution exposure and mid-gestational maternal cytokine and chemokine levels were significantly, though weakly, correlated . Some differences by ASD with and without ID, and for ID only, were noted for certain immune marker and pollutant combinations, though overall conclusions were similar across groups. While mid-pregnancy air pollution was not associated with any neurodevelopmental outcome in EMA, we did find that IL-6 remained associated with ASD with ID even after adjusting for air pollution exposure. These results further support a role of maternal immune activation in risk for neurodevelopmental outcomes, and suggest that prenatal air pollution exposure may result in small effects on the immune system with biological relevance.
Other environmental toxicants: One study, conducted within the first phase of EMA, examined the association between total mercury, a known neurotoxicant, and ASD . While mean mercury levels were higher in the maternal samples for ASD cases than for GP controls in the unadjusted analysis, differences did not persist after adjustment for covariates, and no case–control differences were noted in neonatal samples. Thus, mercury was not examined in phase 2.
Genetic factors are known to play a strong role in ASD, with heritability estimates ranging from 50–95%, and evidence of contributions of both rare and common genetic variants [84,85,86]. Genetic studies in EMA have investigated several hypotheses. First, we applied a novel approach considering maternal genetic effects both alone and in combination with offspring genetic effects. We identified suggestive transgenerational interactions that cannot be explained by maternal or child main effects only ; Fig. 3a). We also examined maternal autosomal burden of large, rare copy number variants (CNVs) for transgenerational effects contributing to ASD in children, addressing the hypothesis of a maternal-specific effect (during gestation, rather than interaction between maternal and offspring genotype). We observed that the mothers of children with ASD showed a higher burden of deletions and duplications (> 1 kb) compared with (1) control mothers of unaffected children, (2) both the ASD or control children’s CNV burden, and (3) the enrichment observed in females in the general population ; Fig. 3b). Thus, a maternal-specific genetic contribution of excess CNVs in proband was demonstrated.
Second, in order to ultimately assess potential G × E effects, we first applied SNP-based heritability and genome-wide association studies to several E biomarkers. Our biomarker results (unrelated to ASD risk) suggested that a number of environmental prenatal and postnatal biomarkers were highly heritable and/or strongly associated with specific loci [42,43,44, 49] (Fig. 3c). Results suggested that maternal PBDEs (5/7), maternal (2/22) and neonatal (8/42) cytokines/chemokines and vitamin D circulating during pregnancy and/or at birth are not solely determined by external stimuli (e.g. pollution or infections), but show significant contribution of specific maternal genetic loci relevant for xenobiotic and lipid metabolism , immune function , and vitamin D metabolism , respectively.
Third, we examined the relationships between biomarkers and genetic markers across mothers and neonates, independent of ASD outcome. We found evidence that maternal genetics was associated with neonatal immune biomarkers (3/42), and that fetal genetics had a strong association with mid-gestational maternal immune biomarkers (2/22) . In most cases, we observed independent genetic contribution from the maternal non-transmitted genome and fetal non-inherited genomes and/or a different set of regulatory loci independently associated with maternal or fetal genetic influences. An example is the maternal soluble IL-2 receptor (sIL2-Ra) regulated by maternal genetics through the cytokine-encoding locus (IL2RA) and by fetal genetics through a non-cytokine encoding locus (ADCYAP1). The nearby ADCYAP1 gene exerts a transcriptional activating role not previously related to immune function . Additionally, suggestive cross-associated loci influencing maternal and neonatal vitamin D were observed near immune genes. Relevant to pregnancy outcomes, we observed two distinct suggestive interactions between a maternal genetic variant, neonatal vitamin D and ID outcome, and between a second independent maternal genetic variant, neonatal IL-8 and ASD outcome [43, 44].
When we assessed whether the genetic factors controlling the measured biomarkers were associated with ASD outcome, we found none that were statistically convincing (although our power for case–control status was low due to the sample size) [42,43,44, 49]. However, these studies help to elucidate how genetics strongly influences maternal physiology during gestation with intriguing potential consequences for consideration of the environmental response and its complex interplay with genetics.
The EMA study has reported results for biological markers in association with ASD in over 20 publications to date. By capitalizing on archived biospecimens and widely-available diagnostic outcome data in the state of California, the EMA design enabled a number of “study firsts,” providing novelty in prospective biomarker assessments in ASD research. EMA has allowed researchers to address questions that had previously relied on exposure as estimated through retrospective reporting, biosamples collected after suspected critical windows, and/or had been conducted in small samples. Across work from the EMA study, results support the indication that the prenatal period is one of particular relevance and importance to ASD etiology. EMA has also facilitated comparisons between ASD and ID, specifically suggesting differences in certain immune marker relationships by ASD with and without co-occurring ID, and overall similarity between ASD and ID for associations with EDCs. Several EMA studies have also addressed sex differences in ASD, with differences suggested for vitamin D and PBDEs in particular. Below we further discuss key contributions from the EMA study, as well as future directions for EMA and related work in the field.
Some of the strongest findings from EMA point to the involvement of maternal and newborn immune dysregulation in ASD. Immunologic profiles of mothers of children with ASD, specifically ASD with ID, were marked by elevated levels of several inflammatory cytokines/chemokines that would normally be downregulated during mid-gestation. Neonates with ASD with co-occurring ID were also more likely than ID-only and GP groups to have elevated levels of some inflammatory cytokines/chemokines indicative of heightened immune activation. The findings in EMA demonstrating an association with immune markers and developmental outcome are consistent with a large and growing literature pointing to a role of prenatal and early life immune activation in ASD [87,88,89,90]. A primary function of cytokines is to act as growth factors, prime the immune system, and modulate immune activation in response to infection. However, a growing literature supports their key role in mediating signals between the immune and nervous systems, influencing neuronal response and subsequent behaviors as well as brain growth early in development. These maternal and neonatal immune findings in EMA point to promising avenues for understanding the biological mechanisms underlying specific sub-phenotypes of ASD, possible therapeutic targets, and/or early diagnostic biomarkers.
In addition to roles in brain development, cytokines also have been shown to alter neurite outgrowth and synaptic connectivity in cultured rodent neurons . Translational extensions of this work exposing human neuronal cell lines to the 5 cytokines elevated in the EMA ASD with ID group have found that this inflammatory cytokine mix induced apoptosis, altered cytokine receptor expression, and promoted apoptosis . These results indicated that concentration-dependent and interactive cytokine effects from clinically relevant mixtures were shown for the first time to affect developing human neurons, and specifically, to affect endpoints with relevance to ASD . Of future interest is the concerted effect of environmental exposures, such as EDCs, and inflammatory molecules on neuron development.
Only a few other studies have examined prenatal levels of the EDC chemical classes examined in EMA [92,93,94,95], some with similar findings to ours, and others with differences. Overall, in our study we found evidence for increased risk of ASD with higher levels of several PCBs, suggestions of sexually dimorphic effects with PBDEs, and overall null findings for OC pesticides and PFASs. We also observed similarity in risk patterns with these EDCs for ID, as well as ASD with and without comorbid ID. Recent work in the Maternal-Infant Research on Environmental Chemicals (MIREC) study in Canada supported EMA findings for PCBs in association with quantitatively-assessed ASD traits according to the Social Responsiveness Scale . Another analysis in the Finnish Prenatal Study of Autism reported no associations between prenatal exposure to PCBs, but did suggest an association with the organochlorine insecticide p’p-DDE . In EMA, p’p-DDE was not associated with ASD, but was associated with ID in a non-monotonic fashion. Mixture analyses examining these classes of EDCs together in our dataset, which have not yet been replicated in other studies, did not reveal any significant associations. The discrepancies in findings according to these different modeling strategies highlight the need for additional work in this area, and suggest either chance findings due to failure to account for correlation across chemicals in the individual chemical class analyses (though findings in those studies did persist in co-exposure models adjusting for one other ‘representative’ chemical from one other class), or perhaps challenges in identifying modest associations in models that necessarily provide conservative estimates when a large number of exposures are examined in a single framework. Because individuals are exposed to hundreds of chemicals, advancement of our understanding of developmental consequences of such exposure may necessitate further assessment of these combined effects.
EMA has also contributed evidence to our understanding of endogenous factors that possibly influence ASD development. Our findings with respect to maternal TSH in mid-pregnancy are generally consistent with literature linking ASD to disrupted thyroid regulation and thyroid conditions in mothers before and during pregnancy [63, 64, 97,98,99]. However, in the few studies that have directly measured thyroid hormone levels during the prenatal and neonatal periods, the results for associations with neonatal levels are neither consistent nor supportive of strong overall associations [100,101,102,103]. Continued work in this area, including studies of environmental and therapeutic influences on the maternal thyroid axis, could provide important insights.
Vitamin D deficiency has been hypothesized as a risk factor for ASD, but early studies were small and often based on estimated exposure (rather than measured levels), or on levels measured outside of relevant developmental windows [104,105,106,107,108]. While results for associations with vitamin D have not been uniform in the broader literature, there is somewhat greater consistency for inverse associations from lower sunlight regions . A recent meta-analysis on neurodevelopmental effects of prenatal vitamin D  suggested a protective association (pooled odds ratio of 0.42, 95% CI 0.25, 0.71, for highest vs. lowest prenatal vitamin D). However, only five studies were included, and not all had clinical ASD as the outcome. A few subsequent studies, two of which are also from high sunlight regions, have findings similar to ours, with no overall association between vitamin D deficiency or levels and ASD, but suggestions of a relationship in some sub-groups (such as whites only) [111,112,113]. Associations with vitamin D should be explored further, as it represents a risk factor with the potential for modification.
Genome-wide genetic data for mothers and their children allowed for several unique mother–child analyses in EMA such as detecting transgenerational effects, GxE interactions, and cross-genetic influence. These studies strengthen growing evidence that common variation in particular genes, including SNPs in or near autism-related genes, and maternal burden of CNVs may increase risk of ASD. Our findings at some loci were replicated or are consistent with other work . These include findings in non-pregnant individuals, including relationships between vitamin D and the GC gene ; HTR7 variation ; and ACKR1; as well as relationships between chemokines and ACKR2 and ACKR4 , and IL-2 and the gene IL2RA . However, results for additional maternal and fetal loci regulating neonatal immune biomarkers and vitamin D do not overlap with those from some large studies in European adults [115, 116]. This could be due to unique signals detected in neonates, before environmental variation is large, and in mothers indicating distinct pregnancy biology, or could be due to inclusion of substantial Hispanic and Asian representation in our study but few other GWAS. Finally, when analyzing the cross-genetic contribution to immune biomarkers and vitamin D, the significant and suggestive associated loci mapped to genes that were not already associated with the analyzed immune biomarkers or vitamin D. These results may suggest that the fetal influence on maternal physiology occurs via the placenta (which shares fetal genetics) and distinct regulatory or indirect mechanisms such as signaling at the placental interface and through a different set of genes than currently implicated in inflammation or beyond the classic regulators of the vitamin D pathway.
Our G × E approach has also provided novel, early insights into both maternal and fetal genetic determinants of maternal and neonatal immune regulation, as well as the metabolism of environmental chemicals and vitamin D, separate from pregnancy outcomes. Our findings suggest that studies focusing on risk factors for neurodevelopmental disorders should investigate G × E and cross-genetic interaction involved in early risk for diseases.
Challenges, limitations, and special considerations in the EMA study
Despite the unique strengths and contributions of the EMA study, certain limitations and considerations for placing findings from this work into the larger context of our understanding of ASD should be noted. The age at ascertainment of Phase 1 cases was on average younger than that for Phase 2, presenting the possibility that Phase 1 cases could represent greater severity. EMA relied on existing diagnostic records from regional centers of the DDS; while these records were reviewed by expert clinicians for consistency with ASD diagnosis, that determination was reliant on available data. Available information from a study that examined the completeness of DDS diagnoses suggests that DDS captures approximately 80% of ASD cases in the state, with milder cases more likely to be missed [117, 118]. We cannot rule out potential misclassification in controls, though any such misclassification would not be expected to influence findings here given the low prevalence of ASD in the general population. Likewise, information on comorbid ID status was reliant on DDS records, and it is possible that some ID cases, or comorbid ID status, may have gone undiagnosed. We do not have data to quantify the extent to which this may occur. We have not been able to examine associations with other comorbidities highly prevalent in ASD, such as anxiety and ADHD, as may be useful to address questions related to associations with other phenotypic subgroups and comorbidites. We had a relatively small number of female ASD cases, limiting power to address sex differences in some analyses. We also lacked in-depth clinical data across the child’s lifespan, preventing the study of developmental trajectories. Certain subgroup analyses (e.g., by sex) may have had more limited statistical power and require replication in larger samples. Many factors examined in phase 1 were not examined in phase 2 participants, precluding ability to pool data across phases. Prenatal confounder information was limited to birth certificate and prenatal screening data, and information on some potential confounders and modifiers, such as parental psychiatric history, or maternal diet, was not available. The EMA design did not include any direct contact with the participants; thus, additional research measures (such as parent-report questionnaires) that may be conducted in ASD case–control studies are not present in our database. These have been the trade-offs enabling a large population-based study with diagnostic information and biospecimens collected during critical developmental windows in mother–child pairs.
As a genetic study, the EMA design and dataset had several relatively unusual features. Though we were able to utilize mother–child pairs to assess some family-based analyses, we did not have paternal genetic contribution to maximize the understanding of which effects might be specific to pregnancy. On the other hand, the availability of unrelated controls allowed for exploration of G × E effects that are challenging in family-based designs. Although the genetic heterogeneity of the dataset created the risk of population stratification, we were able to adjust the analyses for genetic ancestry using the principal components derived from genome-wide maternal and fetal SNPs to account for population stratification, facilitating increased generalizability to the broader US population. The study contributed novel information regarding the relationship between genetic variation and biomarker levels, supporting strong patterns of genetic influence in such levels.
Additional considerations for EMA’s environmental exposure findings, and for placing this study’s results within the larger context of related work, include its geographic location and setting. Our overall null findings for certain factors may have been related to reduced exposure variability; specifically, vitamin D deficiency was low, while air pollution exposure was relatively high, for the Southern California counties included in EMA. EMA’s decision to match cases and controls on month and year of birth (to address concerns related to seasonality in infections, given a key focus of the study was on the immune system’s role in ASD) may also impact exposures that are seasonally associated and reduce ability to detect such associations. For exposures with decreasing (e.g., PBDEs, OC pesticides) or increasing (e.g., PFASs) use, the time period of EMA births (early 2000′s) must also be considered, as varying exposure levels may impact the ability to detect associations, and their observed effect sizes.
Finally, in placing EMA environmental findings in context, combined effects of exposures must also be considered. Environmental exposures vary across populations; these include not only air pollution levels and EDCs, but also other environmental factors not measured here but likely to either influence neurodevelopmentally-relevant pathways directly, or interact with these exposures, such as maternal diet. Such variation means that even in the face of true causal associations, results across studies of environmental risk factors may not be consistent. Future replication of findings should consider such population-level variability and seek to address exposures from a comprehensive standpoint. In compiling a dataset with over 70 biomarkers, many of which are measured in both maternal and child samples, as well as genotypes, the EMA study has been an important step in ASD research towards this overarching aim of considering multiple exposures. Future studies in the -omics and consortium era are likely to build on these findings.
Future directions and next steps for EMA
The next phases of the EMA study focus on tying together the above exposures to consider pathways and mechanisms, and additional extensions of our combined exposures work. Animal model evidence, and some emerging human studies, suggest that both thyroid hormones and vitamin D may lie along the hypothesized pathway between the environmental chemicals considered in EMA and neurodevelopmental outcomes such as ASD and ID. Emerging work has demonstrated air pollution levels may lead to lower vitamin D levels . Though evidence to date has not provided strong and consistent support for direct associations between thyroid hormones and ASD, thyroid hormones are known to play an important role in neurodevelopment, and links between thyroid hormone levels and EDC exposure have been demonstrated . Ongoing analyses in EMA will examine the potential interaction or mediating effects of vitamin D, as well as TSH, on the associations between chemical exposures and ASD. To our knowledge, EMA will be the first study to examine such pathways in pregnancy for certain classes of chemicals.
In addressing genetic contributions, EMA will also be able to apply our novel G × E and maternal–fetal genetics approaches to additional pre- and post-natal risk factors, including maternal antibodies to fetal brain, to potentially identify novel genetic involvement. In addition to further G × E considerations, EMA has the capacity to examine novel E × E interactions, including not just extensions of existing EDC mixture analyses that may incorporate air pollutant exposures, but also similar strategies to address combined effects of biomarkers available and previously examined only as single exposures (such as, for example, cytokines).
The EMA study capitalized on existing resources to address what had at one time seemed to be unaddressable questions; that is, well-powered analyses based on biomarker levels from prospectively collected biospecimens from critical periods during neurodevelopment for a rare outcome. Its success highlights the value and utility of resources like the California Biobank Program for archived biosamples and outcome registries like DDS for research purposes. Key findings from EMA include supporting the role of prenatal immune markers in ASD; suggesting potential relationships with levels of certain environmental chemicals; and providing evidence for genetic involvement in regulating or influencing levels of immune markers, endocrine disrupting chemicals, and vitamin D during pregnancy and in the neonatal period. The EMA design has created rich opportunity for analyses of multiple exposures, intermediate biomarkers, and outcomes, and can serve as a model for other studies.
Availability of data and materials
The datasets used and/or analyzed as part of EMA studies are available from the authors with permission from the California Committee for the Protection of Human Subjects and the California Biobank Program. Genetic data from prenatal and neonatal specimens are not publicly available.
Autism spectrum disorder
Copy number variants
Department of Developmental Services
Endocrine disrupting chemicals
Early Markers for Autism
- G × E:
Gene by environment
Genetic Disease Screening Program
Genome-wide association study
Polybrominated diphenyl ethers
Thyroid stimulating hormone
American Psychiatric Association, ed. Diagnostic and statistical manual of mental disorders (5th edition). 2013, American Psychiatric Association: Washington, DC.
De Rubeis S, Buxbaum JD. Genetics and genomics of autism spectrum disorder: embracing complexity. Hum Mol Genet. 2015;24(R1):R24-31.
Lyall K, et al. The changing epidemiology of autism spectrum disorders. Annu Rev Public Health. 2017a;38:81–102.
Bilbo SD, Schwarz JM. The immune system and developmental programming of brain and behavior. Front Neuroendocrinol. 2012;33(3):267–86.
Goines P, Van de Water J. The immune system’s role in the biology of autism. Curr Opin Neurol. 2010;23(2):111–7.
Nardone S, Elliott E. The interaction between the immune system and epigenetics in the etiology of autism spectrum disorders. Front Neurosci. 2016;10:329.
Patterson PH. Maternal infection and immune involvement in autism. Trends Mol Med. 2011;17(7):389–94.
Kim KH, et al. Diverse influences of androgen-disrupting chemicals on immune responses mounted by macrophages. Inflammation. 2014;37(3):649–56.
Kortenkamp A, Ten, . years of mixing cocktails: a review of combination effects of endocrine-disrupting chemicals. Environ Health Perspect. 2007;115(Suppl 1):98–105.
Weiss B. The intersection of neurotoxicology and endocrine disruption. Neurotoxicology. 2012;33(6):1410–9.
Ashwood P, Van de Water J. A review of autism and the immune response. Clin Dev Immunol. 2004a;11(2):165–74.
Ashwood P, Van de Water J. Is autism an autoimmune disease? Autoimmun Rev. 2004b;3(7–8):557–62.
Comi AM, et al. Familial clustering of autoimmune disorders and evaluation of medical risk factors in autism. J Child Neurol. 1999;14(6):388–94.
Sweeten TL, et al. Increased prevalence of familial autoimmunity in probands with pervasive developmental disorders. Pediatrics. 2003;112(5):e420.
Ashwood P, et al. Intestinal lymphocyte populations in children with regressive autism: evidence for extensive mucosal immunopathology. J Clin Immunol. 2003;23(6):504–17.
Ashwood P, et al. Spontaneous mucosal lymphocyte cytokine profiles in children with autism and gastrointestinal symptoms: mucosal immune activation and reduced counter regulatory interleukin-10. J Clin Immunol. 2004;24(6):664–73.
Ashwood P, et al. Decreased transforming growth factor beta1 in autism: a potential link between immune dysregulation and impairment in clinical behavioral outcomes. J Neuroimmunol. 2008a;204(1–2):149–53.
Ashwood P, et al. Brief report: plasma leptin levels are elevated in autism: association with early onset phenotype? J Autism Dev Disord. 2008b;38(1):169–75.
Ashwood P, et al. Preliminary evidence of the in vitro effects of BDE-47 on innate immune responses in children with autism spectrum disorders. J Neuroimmunol. 2009;208(1–2):130–5.
Jyonouchi H, Sun S, Le H. Proinflammatory and regulatory cytokine production associated with innate and adaptive immune responses in children with autism spectrum disorders and developmental regression. J Neuroimmunol. 2001;120(1–2):170–9.
Jyonouchi H, Sun S, Itokazu N. Innate immunity associated with inflammatory responses and cytokine production against common dietary proteins in patients with autism spectrum disorder. Neuropsychobiology. 2002;46(2):76–84.
Jyonouchi H, et al. Dysregulated innate immune responses in young children with autism spectrum disorders: their relationship to gastrointestinal symptoms and dietary intervention. Neuropsychobiology. 2005;51(2):77–85.
Enstrom A, et al. Increased IgG4 levels in children with autism disorder. Brain Behav Immun. 2009a;23(3):389–95.
Enstrom A, et al. Peripheral blood leukocyte production of bdnf following mitogen stimulation in early onset and regressive autism. Am J Biochem Biotechnol. 2008;4(2):121–9.
Heuer L, et al. Reduced levels of immunoglobulin in children with autism correlates with behavioral symptoms. Autism Res. 2008;1(5):275–83.
Enstrom AM, et al. Altered gene expression and function of peripheral blood natural killer cells in children with autism. Brain Behav Immun. 2009b;23(1):124–33.
Ribas-Fito N, et al. Breastfeeding, exposure to organochlorine compounds, and neurodevelopment in infants. Pediatrics. 2003;111(5 Pt 1):e580–5.
Steinmaus C, et al. Thyroid hormones and moderate exposure to perchlorate during pregnancy in women in Southern California. Environ Health Perspect. 2016;124(6):861–7.
Cunningham GC, Tompkinison DG. Cost and effectiveness of the California triple marker prenatal screening program. Genet Med. 1999;1(5):199–206.
Currier R, et al. Integrated and first trimester prenatal screening in California: program implementation and patient choice for follow-up services. Prenat Diagn. 2012;32(11):1077–83.
Goines PE, et al. Increased midgestational IFN-gamma, IL-4 and IL-5 in women bearing a child with autism: a case-control study. Mol Autism. 2011;2:13.
Zerbo O, et al. Neonatal cytokines and chemokines and risk of Autism Spectrum Disorder: the Early Markers for Autism (EMA) study: a case-control study. J Neuroinflammation. 2014;11:113–113.
Ames JL, et al. Neonatal thyroid stimulating hormone and subsequent diagnosis of autism spectrum disorders and intellectual disability. Autism Res. 2020;13(3):444–55.
Yau VM, et al. Prenatal and neonatal thyroid stimulating hormone levels and autism spectrum disorders. J Autism Dev Disord. 2015;45(3):719–30.
Kvaskoff D, et al. minimizing matrix effects for the accurate quantification of 25-hydroxyvitamin D metabolites in dried blood spots by LC-MS/MS. Clin Chem. 2016;62(4):639–46.
Eyles D, et al. A sensitive LC/MS/MS assay of 25OH vitamin D3 and 25OH vitamin D2 in dried blood spots. Clin Chim Acta. 2009;403(1–2):145–51.
Windham GC, et al. Newborn vitamin D levels in relation to autism spectrum disorders and intellectual disability: a case-control study in california. Autism Res. 2019;12(6):989–98.
Lyall K, et al. Prenatal maternal serum concentrations of per- and polyfluoroalkyl substances in association with autism spectrum disorder and intellectual disability. Environ Health Perspect. 2018;126(1):017001.
Lyall K, et al. Prenatal serum concentrations of brominated flame retardants and autism spectrum disorder and intellectual disability in the early markers of autism study: a population-based case-control study in California. Environ Health Perspect. 2017b;125(8):087023.
Lyall K, et al. Polychlorinated biphenyl and organochlorine pesticide concentrations in maternal mid-pregnancy serum samples: association with autism spectrum disorder and intellectual disability. Environ Health Perspect. 2017c;125(3):474–80.
Tsang KM, et al. A genome-wide survey of transgenerational genetic effects in autism. PLoS ONE. 2013;8(10):e76978.
Traglia M, et al. Independent maternal and fetal genetic effects on midgestational circulating levels of environmental pollutants. G3 (Bethesda). 2017;7(4):1287–99.
Traglia M, et al. Cross-genetic determination of maternal and neonatal immune mediators during pregnancy. Genome Med. 2018;10(1):67.
Traglia M, et al. Genetic contributions to maternal and neonatal vitamin D levels. Genetics. 2020;214(4):1091–102.
Desachy G, et al. Increased female autosomal burden of rare copy number variants in human populations and in autism families. Mol Psychiatry. 2015;20(2):170–5.
Dahl A, et al. A robust method uncovers significant context-specific heritability in diverse complex traits. Am J Hum Genet. 2020;106(1):71–91.
Frost HR, et al. Identifying significant gene-environment interactions using a combination of screening testing and hierarchical false discovery rate control. Genet Epidemiol. 2016;40(7):544–57.
Novembre J, et al. Genes mirror geography within Europe. Nature. 2008;456(7218):98–101.
Zerbo O, et al. Maternal mid-pregnancy C-reactive protein and risk of autism spectrum disorders: the early markers for autism study. Transl Psychiatry. 2016;6:e783.
Croen LA, et al. Maternal mid-pregnancy autoantibodies to fetal brain protein: the early markers for autism study. Biol Psychiatry. 2008;64(7):583–8.
Braunschweig D, et al. Autism: maternally derived antibodies specific for fetal brain proteins. Neurotoxicology. 2008;29(2):226–31.
Zimmerman AW, et al. Maternal antibrain antibodies in autism. Brain Behav Immun. 2007;21(3):351–7.
Dalton P, et al. Maternal neuronal antibodies associated with autism and a language disorder. Ann Neurol. 2003;53(4):533–7.
Grether JK, et al. Prenatal and newborn immunoglobulin levels from mother-child pairs and risk of autism spectrum disorders. Front Neurosci. 2016;10:218.
Jones KL, et al. Autism with intellectual disability is associated with increased levels of maternal cytokines and chemokines during gestation. Mol Psychiatry. 2017;22(2):273–9.
Heuer LS, et al. An exploratory examination of neonatal cytokines and chemokines as predictors of autism risk: the early markers for autism study. Biol Psychiatry. 2019;86(4):255–64.
Morreale de Escobar G, Obregon MJ, Escobar del Rey F. Role of thyroid hormone during early brain development. Eur J Endocrinol. 2004;151 Suppl 3:U25–37.
Gillberg IC, Gillberg C, Kopp S. Hypothyroidism and autism spectrum disorders. J Child Psychol Psychiatry. 1992;33(3):531–42.
Roman GC. Autism: transient in utero hypothyroxinemia related to maternal flavonoid ingestion during pregnancy and to other environmental antithyroid agents. J Neurol Sci. 2007;262(1–2):15–26.
Fetene DM, Betts KS, Alati R. Mechanisms in endocrinology: Maternal thyroid dysfunction during pregnancy and behavioural and psychiatric disorders of children: a systematic review. Eur J Endocrinol. 2017;177(5):R261–73.
Thompson W, et al. Maternal thyroid hormone insufficiency during pregnancy and risk of neurodevelopmental disorders in offspring: a systematic review and meta-analysis. Clin Endocrinol (Oxf). 2018;88(4):575–84.
Andersen SL, et al. Attention deficit hyperactivity disorder and autism spectrum disorder in children born to mothers with thyroid dysfunction: a Danish nationwide cohort study. BJOG. 2014;121(11):1365–74.
Getahun D, et al. Association between maternal hypothyroidism and autism spectrum disorders in children. Pediatr Res. 2018;83(3):580–8.
Roman GC, et al. Association of gestational maternal hypothyroxinemia and increased autism risk. Ann Neurol. 2013;74(5):733–42.
Schwalfenberg GK. A review of the critical role of vitamin D in the functioning of the immune system and the clinical implications of vitamin D deficiency. Mol Nutr Food Res. 2011;55(1):96–108.
Janbek J, Specht IO, Heitmann BL. Associations between vitamin D status in pregnancy and offspring neurodevelopment: a systematic literature review. Nutr Rev. 2019;77(5):330–49.
Veena SR, et al. Association between maternal nutritional status in pregnancy and offspring cognitive function during childhood and adolescence; a systematic review. BMC Pregnancy Childbirth. 2016;16:220.
Veldman CM, Cantorna MT, DeLuca HF. Expression of 1,25-dihydroxyvitamin D(3) receptor in the immune system. Arch Biochem Biophys. 2000;374(2):334–8.
Walker VP, et al. Cord blood vitamin D status impacts innate immune responses. J Clin Endocrinol Metab. 2011;96(6):1835–43.
Eyles D, Burne T, McGrath J. Vitamin D in fetal brain development. Semin Cell Dev Biol. 2011;22(6):629–36.
Windham GC, et al. Maternal vitamin D levels during pregnancy in association with autism spectrum disorders or intellectual disability in offspring; exploring non-linear patterns and demographic sub-groups. Autism Res. 2020;13(12):2216–29.
Holick MF. Vitamin D deficiency. N Engl J Med. 2007;357(3):266–81.
Lema SC, et al. Dietary exposure to 2,2’,4,4’-tetrabromodiphenyl ether (PBDE-47) alters thyroid status and thyroid hormone-regulated gene transcription in the pituitary and brain. Environ Health Perspect. 2008;116(12):1694–9.
Ghassabian A, Trasande L. Disruption in thyroid signaling pathway: a mechanism for the effect of endocrine-disrupting chemicals on child neurodevelopment. Front Endocrinol (Lausanne). 2018;9:204.
Hertz-Picciotto I, et al. Prenatal exposures to persistent and non-persistent organic compounds and effects on immune system development. Basic Clin Pharmacol Toxicol. 2008;102(2):146–54.
Kimura-Kuroda J, Nagata I, Kuroda Y. Disrupting effects of hydroxy-polychlorinated biphenyl (PCB) congeners on neuronal development of cerebellar Purkinje cells: a possible causal factor for developmental brain disorders? Chemosphere. 2007;67(9):S412–20.
Chevrier J, et al. Prenatal exposure to polybrominated diphenyl ether flame retardants and neonatal thyroid-stimulating hormone levels in the CHAMACOS study. Am J Epidemiol. 2011;174(10):1166–74.
Liu C, et al. PCB153 disrupts thyroid hormone homeostasis by affecting its biosynthesis, biotransformation, feedback regulation, and metabolism. Horm Metab Res. 2012;44(9):662–9.
Hamra GB, et al. Prenatal exposure to endocrine-disrupting chemicals in relation to autism spectrum disorder and intellectual disability. Epidemiology. 2019;30(3):418–26.
Genc S, et al. The adverse effects of air pollution on the nervous system. J Toxicol. 2012;2012:782462.
Costa LG, et al. Neurotoxicants are in the air: convergence of human, animal, and in vitro studies on the effects of air pollution on the brain. Biomed Res Int. 2014;2014:736385.
Volk HE, et al. Maternal immune response and air pollution exposure during pregnancy: insights from the early markers for autism (EMA) Study. J Neurodevelop Disord. 2020;12(1):42. https://doi.org/10.1186/s11689-020-09343-0.
Yau VM, et al. Prenatal and neonatal peripheral blood mercury levels and autism spectrum disorders. Environ Res. 2014;133:294–303.
Bai D, et al. Association of genetic and environmental factors with autism in a 5-country cohort. JAMA Psychiatry. 2019;76(10):1035–43.
Hallmayer J, et al. Genetic heritability and shared environmental factors among twin pairs with autism. Arch Gen Psychiatry. 2011;68(11):1095–102.
Sandin S, et al. The familial risk of autism. JAMA. 2014;311(17):1770–7.
Thom RP, et al. Beyond the brain: A multi-system inflammatory subtype of autism spectrum disorder. Psychopharmacology (Berl). 2019;236(10):3045–61. https://doi.org/10.1007/s00213-019-05280-6.
Guma E, Plitman E, Chakravarty MM. The role of maternal immune activation in altering the neurodevelopmental trajectories of offspring: A translational review of neuroimaging studies with implications for autism spectrum disorder and schizophrenia. Neurosci Biobehav Rev. 2019;104:141–57.
Matelski L, et al. Effects of cytokines on nuclear factor-kappa B, cell viability, and synaptic connectivity in a human neuronal cell line. Mol Psychiatry. 2021;26(3):875–87. https://doi.org/10.1038/s41380-020-0647-2.
Estes ML, et al. Baseline immunoreactivity before pregnancy and poly(I:C) dose combine to dictate susceptibility and resilience of offspring to maternal immune activation. Brain Behav Immun. 2020;88:619–30. https://doi.org/10.1016/j.bbi.2020.04.061.
Li L, et al. Endogenous interferon gamma directly regulates neural precursors in the non-inflammatory brain. J Neurosci. 2010;30(27):9038–50.
Braun JM, et al. Gestational Exposure to Endocrine-Disrupting Chemicals and Reciprocal Social, Repetitive, and Stereotypic Behaviors in 4- and 5-Year-Old Children: The HOME Study. Environ Health Perspect. 2014;122(5):513–20.
Brown, A.S., et al. Association of maternal insecticide levels with autism in offspring from a National Birth Cohort. Am J Psychiatry, 2018: p. appiajp201817101129.
Cheslack-Postava K, et al. Maternal serum persistent organic pollutants in the Finnish Prenatal Study of Autism: a pilot study. Neurotoxicol Teratol. 2013;38:1–5.
Liew Z, et al. Attention deficit/hyperactivity disorder and childhood autism in association with prenatal exposure to perfluoroalkyl substances: a nested case-control study in the Danish National Birth Cohort. Environ Health Perspect. 2015;123(4):367–73.
Bernardo BA, et al. Assessing the relation between plasma PCB concentrations and elevated autistic behaviours using bayesian predictive odds ratios. Int J Environ Res Public Health. 2019;16(3):457.
Rotem RS, et al. Maternal thyroid disorders and risk of autism spectrum disorder in progeny. Epidemiology. 2020;31(3):409–17.
Andersen SL, et al. Maternal thyroid function in early pregnancy and child neurodevelopmental disorders: a Danish Nationwide Case-Cohort Study. Thyroid. 2018;28(4):537–46.
Ge GM, et al. Maternal thyroid dysfunction during pregnancy and the risk of adverse outcomes in the offspring: a systematic review and meta-analysis. J Clin Endocrinol Metab. 2020;105(12):dgaa555.
Hoshiko S, et al. Are thyroid hormone concentrations at birth associated with subsequent autism diagnosis? Autism Res. 2011;4(6):456–63.
Korzeniewski SJ, et al. Association between transient hypothyroxinaemia of prematurity and adult autism spectrum disorder in a low-birthweight cohort: an exploratory study. Paediatr Perinat Epidemiol. 2013;27(2):182–7.
Lyall K, et al. Neonatal thyroid hormone levels in association with autism spectrum disorder. Autism Res. 2017d;10(4):585–92.
Soldin OP, et al. Lack of a relation between human neonatal thyroxine and pediatric neurobehavioral disorders. Thyroid. 2003;13(2):193–8.
Fernell E, et al. Serum levels of 25-hydroxyvitamin D in mothers of Swedish and of Somali origin who have children with and without autism. Acta Paediatr. 2010;99(5):743–7.
Humble MB, Gustafsson S, Bejerot S. Low serum levels of 25-hydroxyvitamin D (25-OHD) among psychiatric out-patients in Sweden: relations with season, age, ethnic origin and psychiatric diagnosis. J Steroid Biochem Mol Biol. 2010;121(1–2):467–70.
Meguid NA, et al. Reduced serum levels of 25-hydroxy and 1,25-dihydroxy vitamin D in Egyptian children with autism. J Altern Complement Med. 2010;16(6):641–5.
Molloy CA, et al. Plasma 25(OH)D concentration in children with autism spectrum disorder. Dev Med Child Neurol. 2010;52(10):969–71.
Mostafa GA, Al-Ayadhi LY. Reduced serum concentrations of 25-hydroxy vitamin D in children with autism: relation to autoimmunity. J Neuroinflammation. 2012;9:201.
Zhong C, et al. Maternal dietary factors and the risk of autism spectrum disorders: a systematic review of existing evidence. Autism Res. 2020;13(10):1634–58. https://doi.org/10.1002/aur.2402.
Garcia-Serna AM, Morales E. Neurodevelopmental effects of prenatal vitamin D in humans: systematic review and meta-analysis. Mol Psychiatry. 2020;25(10):2468–81. https://doi.org/10.1038/s41380-019-0357-9.
Lee BK, et al. Developmental vitamin D and autism spectrum disorders: findings from the Stockholm Youth Cohort. Mol Psychiatry. 2019. https://doi.org/10.1038/s41380-019-0578-y.
Schmidt RJ, et al. Neonatal vitamin D status in relation to autism spectrum disorder and developmental delay in the CHARGE case-control study. Autism Res. 2019;12(6):976–88.
López-Vicente M, et al. Maternal circulating Vitamin D(3) levels during pregnancy and behaviour across childhood. Sci Rep. 2019;9(1):14792.
Connolly S, et al. A genome-wide investigation into parent-of-origin effects in autism spectrum disorder identifies previously associated genes including SHANK3. Eur J Hum Genet. 2017;25(2):234–9.
Ahola-Olli AV, et al. Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet. 2017;100(1):40–50.
Jiang X, et al. Genome-wide association study in 79,366 European-ancestry individuals informs the genetic architecture of 25-hydroxyvitamin D levels. Nat Commun. 2018;9(1):260.
Croen LA, et al. The changing prevalence of autism in California. J Autism Dev Disord. 2002;32(3):207–15.
Windham GC, et al. Birth prevalence of autism spectrum disorders in the San Francisco Bay area by demographic and ascertainment source characteristics. J Autism Dev Disord. 2011;41(10):1362–72.
Mousavi SE, et al. Air pollution, environmental chemicals, and smoking may trigger vitamin D deficiency: evidence and potential mechanisms. Environ Int. 2019;122:67–90.
Boas M, Feldt-Rasmussen U, Main KM. Thyroid effects of endocrine disrupting chemicals. Mol Cell Endocrinol. 2012;355(2):240–8.
Project Baby’s Breath was funded by the California Tobacco-Related Disease Research Program, Grant # 8RT-0115; EMA was funded by NIH Grant R01-ES016669.
Ethics approval and consent to participate
Study activities were approved by the Institutional Review Board at Kaiser Permanente and the State Committee for the Protection of Human Subjects. Informed consent was not required for the EMA study; consent forms for screening programs were distributed at the time of collection, which stipulated that specimens and results from prenatal and neonatal testing could be used for research purposes given IRB approval.
Consent for publication
Not applicable; no individual data is presented here.
The authors have no financial conflicts of interest to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lyall, K., Ames, J.L., Pearl, M. et al. A profile and review of findings from the Early Markers for Autism study: unique contributions from a population-based case–control study in California. Molecular Autism 12, 24 (2021). https://doi.org/10.1186/s13229-021-00429-7
- Risk factors
- Immune response
- Early Markers for Autism