- Open Access
Association between fetal sex and maternal plasma microRNA responses to prenatal alcohol exposure: evidence from a birth outcome-stratified cohort
Biology of Sex Differences volume 11, Article number: 51 (2020)
Most persons with fetal alcohol spectrum disorders (FASDs) remain undiagnosed or are diagnosed in later life. To address the need for earlier diagnosis, we previously assessed miRNAs in the blood plasma of pregnant women who were classified as unexposed to alcohol (UE), heavily exposed with affected infants (HEa), or heavily exposed with apparently unaffected infants (HEua). We reported that maternal miRNAs predicted FASD-related growth and psychomotor deficits in infants. Here, we assessed whether fetal sex influenced alterations in maternal circulating miRNAs following prenatal alcohol exposure (PAE). To overcome the loss of statistical power due to disaggregating maternal samples by fetal sex, we adapted a strategy of iterative bootstrap resampling with replacement to assess the stability of statistical parameter estimates. Bootstrap estimates of parametric and effect size tests identified male and female fetal sex-associated maternal miRNA responses to PAE that were not observed in the aggregated sample. Additionally, we observed, in HEa mothers of female, but not male fetuses, a network of co-secreted miRNAs whose expression was linked to miRNAs encoded on the X-chromosome. Interestingly, the number of significant miRNA correlations for the HEua group mothers with female fetuses was intermediate between HEa and UE mothers at mid-pregnancy, but more similar to UE mothers by the end of pregnancy. Collectively, these data show that fetal sex predicts maternal circulating miRNA adaptations, a critical consideration when adopting maternal miRNAs as diagnostic biomarkers. Moreover, a maternal co-secretion network, predominantly in pregnancies with female fetuses, emerged as an index of risk for adverse birth outcomes due to PAE.
Prenatal alcohol exposure (PAE) is common. For instance, in a recent state-wide study in Texas, we found that 8.4% of assessed newborn infants were positive for a blood alcohol metabolite, indicative of at least 3rd trimester PAE . Moreover, a combination of factors, including high rates of untreated alcohol use disorders  and unplanned pregnancies , means that PAE is difficult to prevent. Not surprisingly, fetal alcohol spectrum disorders (FASDs), a cluster of physical and neurobehavioral anomalies that can result from PAE , are common . In the USA, a recent prospective case ascertainment study showed that between 1.1 and 9% of school-aged children met the diagnostic criteria for FASD . This study also revealed a critical need for better diagnostics, as almost all the identified FASD cases were never previously diagnosed . Moreover, vulnerable groups of children, such as those in state foster care systems, are at significant risk for remaining undiagnosed [1, 18]. The inadequacy of diagnosis results in missed opportunities for early intervention to minimize the effects of PAE and, consequently, increases the risk for secondary scholastic and other disabilities .
In attempting to address the critical need for better diagnostic tools, we and others previously reported that maternal plasma miRNAs in both sheep  and humans [5, 29] were altered due to PAE. In our study, in a prospective cohort of pregnant women in Ukraine , we identified maternal miRNAs that could be used to classify heavy alcohol-exposed pregnancies that resulted in the birth of infants with neurodevelopmental and growth deficits (HEa) as distinct from unexposed pregnancies (UE) or heavy alcohol-exposed pregnancies where the infants were apparently unaffected at birth (HEua). In that study, fetal sex contributed to the accuracy of a machine learning-based prediction model for classifying affected pregnancies , raising the intriguing possibility that fetal sex also contributes to the effects of PAE on maternal miRNA profiles.
Sex differences in miRNA secretion emerge as early as the blastocyst stage of embryogenesis and influence maternal endometrial function , and sex differences can contribute to maternal physiology into the fetal period [2, 27]. Moreover, fetal sex may also contribute to pregnancy outcomes, including increased risk for pregnancy-related hypertension and preeclampsia [58, 61, 67]. Inherent transcriptomic differences between male and female fetal tissues like the placenta may further drive differences in the fetal response to alterations in the maternal environment, such as maternal exposure to stress (reviewed in ). Some studies in animal models and human populations have also documented sex differences in the pathophysiology of FASDs, and suggest that PAE-dependent sex differences can appear during the fetal period itself [42, 49, 51, 64].
Collectively, the above studies support the hypothesis that fetal sex is one determinant of the maternal response to alcohol exposure. Moreover, if maternal plasma miRNAs are to be used as predictive biomarkers for the effects of PAE, the association between fetal sex and alterations in maternal miRNAs needs to be assessed. However, small sample sizes in typical FASD studies in human populations are a significant barrier since the disaggregation of data by fetal sex results in decreased sample size and, therefore, decreased statistical power. To overcome this limitation, we employed a statistical strategy of iterative bootstrap resampling with replacement  to simulate characteristics of the population of PAE-exposed pregnancies disaggregated by fetal sex. This methodology has recently been applied to successfully estimate measurement uncertainty in blood tests  and in imaging data  and is a useful tool for assessing highly disaggregated data for the purposes of building prediction models , but to our knowledge, has not been previously applied to studies on FASDs. We report that this strategy uncovered evidence for a stable association between fetal sex and changes in maternal miRNA profiles following PAE. Moreover, these analyses uncovered potential novel regulatory mechanisms for coordinated miRNA secretion into maternal plasma. These data suggest that fetal sex-associated maternal miRNAs are likely to constitute a highly regulated and perhaps adaptive response to PAE that can be leveraged to strengthen the biomarker potential of maternal miRNAs to predict PAE outcomes.
All study protocols were approved by the Institutional Review Boards at the University of California San Diego and Texas A&M University and for data collection at Lviv National Medical University in Ukraine. The study sample was drawn from a larger prospective cohort recruited at perinatal clinics in two regions of Ukraine, in the Rivne and Khmelnytsky provinces, by OMNI-Net teams. Recruitment methods and cohort characteristics have been described elsewhere [5, 17, 20]. The study sample included 93 pregnant women, for whom the plasma that met the purity criteria was available (for sample exclusion criteria, see ). These included second- and third-trimester plasma samples from 68 women that were profiled previously , as well as 25 third-trimester samples from newly assessed subjects. Samples were classified as UE, HEua, and HEa based on the criteria described in Balaraman et al. . Briefly, HEa met the criteria for moderate to heavy drinking, and the child had at least two characteristic alcohol-related craniofacial features (short palpebral fissures, smooth philtrum, and thin vermilion border of the upper lip), growth deficiency, and/or neurobehavioral impairment; HEua met the criteria for moderate to heavy drinking, and the child had no neurobehavioral impairment and no FASD related features; and UE did not meet the criteria for moderate to heavy drinking. Cohort demographic and other characteristics are outlined in Tables 1 and 2.
Selection of miRNAs for analyses
Plasma samples were processed and assessed for 752 miRNAs using the Exiqon/Qiagen Human panel I+II, V5, miRCURY LNA miRNA miRNome PCR Panel, as described previously, including the assessment of absorbance at 414 nm and mRNA expression of Slc4a1/Band3 to detect hemolysis . We selected miRNAs which were expressed in at least 80% of the samples in all the exposure groups collectively, 153 miRNAs met this expression criteria. For these 153 miRNAs, samples with undetected levels of a miRNA were assigned a cycle threshold (Ct) value of “45” as we previously reported . We then normalized the Ct of each expressed miRNA to the average Ct of all expressed miRNAs in that sample, a normalization strategy that has been shown, in large-scale miRNA expression profiling, to outperform alternate normalization strategies utilizing a set of miRNAs or other small RNA controls .
All statistical analyses were implemented in “R” (V3.6.1, R Foundation for Statistical Computing, Vienna, Austria). All “R” code that was developed or adapted for use in these analyses has been deposited at GitHub (https://github.com/nihalasalem/Ukraine_Analysis). For each of the 153 miRNAs, maternal samples for each exposure group (HEa, HEua, and UE) were resampled with replacement for 2000 iterations, consistent with published recommendations [22, 66], and our simulations of statistical parameter stability using an iteration range from 50 to 3000 iterations (Additional file 1). The number of subjects in each resampled group was equal to the number of subjects in each of the HEa, HEua, and UE exposure groups. In each iteration, we reported the results of a statistical test, either ANOVA or ANCOVA p values (R Base Package), Hedges’ g effect size (effsize R Package), or Pearson’s correlation coefficient (corrplot R package). We then repeated these bootstrap resampling procedures by infant sex, e.g., female HEa, HEua, and UE pregnancies. In each case, the distribution of estimates resulting from iterative resampling is hypothesized to approximate the distribution of the population .
Analysis of variance
At each time point, 2nd or 3rd trimester, we resampled each exposure group (both infant sexes combined) and conducted a one-way ANOVA on each of the miRNAs in the resampled groups. We reported the p value for each miRNA and repeated the resampling process for 2000 times (Fig. 1a). For each miRNA, we calculated the proportion of the iterations in which the ANOVA resulted in a significant p value (< 0.05). In a parallel analysis, we re-analyzed the data with bootstrap resampling of an ANCOVA model, with variables, gestational age of blood draw at mid-pregnancy (GAbd1) and late pregnancy (GAbd2) as covariates, to assess the contribution of the gestational age at blood draw (Table 1) to the stability of parametric estimates derived from the ANOVA model since GAbd2 was significantly different across the groups. We then repeated the ANOVA analyses, segregating pregnant mothers into pregnancies that resulted in the birth of male infants or female infants. We then identified miRNAs which exhibited an increased proportion of significant iterations when pregnancies were segregated by infant sex, compared to non-sex-segregated analysis, and labeled these as “likely alcohol-sensitive, fetal sex-specific” maternal miRNAs.
To investigate the directionality of expression (increased vs. decreased) and group-specific changes of the secreted miRNAs, we performed effect size calculations (with Hedges’ g correction for small sample sizes ) on the resampled groups comparing each alcohol exposure group to the unexposed group for each of the 153 miRNAs. The median effect size and 95% confidence interval estimates (computed from the 2.5 and 97.5 percentiles of the resulting effect sizes from 2000 bootstrap iterations) resulting from bootstrap resampling analysis were determined for the combined and sex-segregated groups.
We calculated Pearson’s product-moment correlation between each expressed miRNA and all other expressed miRNAs and constructed a correlation matrix for each group (significant correlations, p < 0.05). We compared the number of significant correlations between the three exposure groups in the second and third trimester separately. This was followed by a re-analysis where maternal samples were segregated by infant sex. All correlation plots were constructed using the “corrplot” package in “R” . To assess the contribution of the gestational age at blood draw (Table 1) on the correlation coefficients, we calculated the partial correlations (controlling for the gestational age at blood draw variables GAbd1 and GAbd2) and calculated R2 between the correlation and the partial correlation coefficients in each exposure group.
To assess potential infant sex-dependent maternal miRNA correlations and control for multiple hypothesis testing , we performed bootstrap resampling analyses of Pearson’s correlations. Each exposure group was resampled with replacement for 2000 iterations, where the number of subjects in each resampled group was equal to the number of subjects in each exposure group (HEa, HEua, and UE). In each iteration, we reported the number of significant correlations, and the resulting data from each group were represented in a frequency distribution histogram. We further tested the null hypothesis that the mean number of significant correlations was not different between the groups. We counted, in the alcohol-exposed groups (HEa and HEua), the number of significant correlations in each bootstrap iteration which exceeded the mean number of significant correlations in the bootstrap iterations of the non-alcohol exposed groups (UE) and divided this count by the number of bootstrap iterations. The resulted ratio was subtracted from 1 to compute the p value. If the p value was less than 0.05, the null hypothesis was rejected.
Correlations by chromosomal location
For each exposure group, we identified the chromosomal locations of the pairs of miRNAs that were significantly correlated. We next tabulated the frequency of significant correlations between each pair of chromosomes and then calculated the percent change in the number of significant correlations for each pair of chromosomes between the alcohol-exposed and unexposed groups, to derive a measure of inter-chromosome correlations in maternal miRNA expression patterns.
Identifying presumptive alcohol-sensitive fetal sex-specific maternal miRNAs
The sample group of 93 pregnant women included in this study differed significantly on the key variable of alcohol consumption (AAD and AADXP, Table 1). A greater proportion of mothers in the HEa and HEua groups also reported current or past smoking, though among self-identified smokers, there were no statistically significant differences in the number of cigarettes consumed per day. Moreover, while there were overall no statistically significant differences between the groups in socioeconomic status, a higher proportion of UE group women reported at least some college-level education (for details about sample demographics, see Table 1).
To identify maternal miRNAs that likely show a fetal sex-dependent response to PAE, we identified the miRNAs which were classified as significant (ANOVA, p < 0.05) in more than 50% of the resampling iterations following bootstrap analyses. Overall, sex segregation resulted in a decreased proportion of significant p values for most miRNAs, a predicted outcome due to the loss of statistical power with decreased sample size. However, a fraction of miRNAs, between 4 and 7%, exhibited an increased probability of achieving significance when data were segregated by sex (Fig. 1b). These miRNAs exhibited increased proportion of p values that exceeded the significance criterion, compared to the combined analyses, and, moreover, exceeded a probability of 0.5 for achieving significance (Fig. 1c, see Additional file 2 for “likely alcohol-sensitive, fetal sex-specific” miRNAs in each group). To assess the potential effects of gestational age at blood draw at mid- and late pregnancy (variables GAbd1 and GAbd2, Table 1) on the stability of parameter estimates, we additionally combined ANCOVA (with GAbd1 and GAbd2 as covariates) with bootstrap resampling analysis. We then compared, for each miRNA, the proportion of iterations that were statistically significant (p < 0.05) with ANCOVA vs. with ANOVA. This analysis showed a concordance (R2s between 0.76 and 0.95, see Additional file 3) between ANOVA bootstrap- and ANCOVA bootstrap-derived estimates of parameter stability indicating that GAbd1 and GAbd2 did not significantly influence outcomes
To assess the specific contribution of infant PAE outcomes to maternal miRNAs, we conducted a post hoc bootstrap resampling analysis, to compare the effect size (Hedges’ g) of changes in maternal miRNA expression in the HEa and HEua groups compared to the UE group (HEa vs. UE, HEua vs. UE). We assessed the stability of group differences first in the composite data and second after disaggregating data by infant sex. A number of sex-segregated miRNAs exhibited median effect sizes that were outside the range of the 95% confidence estimate for the composite sample, but the majority of PAE effects on maternal miRNAs were fetal sex-independent (Fig. 2). miRNAs whose 95% confidence estimates did not include zero were classified as exhibiting a “significant effect.” This analysis identified a number of maternal miRNAs in both HEa and HEua groups, where no effect of PAE (relative to the UE group) was observed in the composite sample, but where there was a significant effect of PAE when the samples were disaggregated by sex (Fig. 3, Additional file 4). For the lists of significant “alcohol-sensitive, fetal sex-specific” miRNAs in each group, see Additional file 5.
To determine if likely sex-specific miRNAs following PAE were due to normally occurring developmental sexual dimorphisms (i.e., differences between control group (UE) mothers who gave birth to male infants vs. those who gave birth to female infants), we segregated the UE population by infant sex, iteratively resampled each group, calculated the effect size due to fetal sex, and after 2000 iterations calculated the median effect size and 95% confidence interval. In the second trimester, 7 miRNAs were significantly decreased, while 4 miRNAs were significantly increased in the plasma of UE mothers who gave birth to male infants compared to those who gave birth to female infants. In the third trimester, 7 miRNAs were significantly decreased while 3 miRNAs were significantly increased in the plasma of UE mothers who gave birth to male infants compared to those who gave birth to female infants. Sex-specific miRNAs in the control, UE, group, identified in the second trimester, did not overlap with those identified in the third trimester. Overall, the fetal sex-specific miRNAs at baseline contributed minimally (0–25%) to the groups of identified alcohol-sensitive fetal sex-specific maternal miRNAs (Table 3).
Ethanol exposure results in increased correlated expression of secreted miRNAs in the second and third trimesters in the HEa group
The factors that control miRNA secretion into biofluids are largely unknown. However, evidence for correlated expression of miRNAs would lend support to a hypothesis that expression, or secretion, of circulating miRNAs are co-regulated by a set of common causal mechanisms. To assess the likelihood that maternal miRNAs are co-regulated and/or co-secreted in response to factors in the maternal environment, such as alcohol, we computed the correlation matrices for all expressed miRNAs in each exposure group at mid- and late gestation.
PAE increased significant miRNA cross-correlations in the HEa group in both the second and third trimester (Table 4, Fig. 4). To test the stability of those differences in miRNA correlations, we performed bootstrap resampling analyses with replacement within each group and calculated the number of significant correlations in each iteration (data represented as a bootstrap resampling distribution, Fig. 4b). We tested the null hypothesis that the bootstrap resampling-derived mean number of significant correlations was not significantly different across the exposure groups. The bootstrap analyses indicated a statistically significant overall increase in the number of significant correlations in the HEa group compared to the unexposed groups (UE) in the second and third trimesters.
Prenatal alcohol’s effect on maternal secreted miRNA cross-correlations is dependent on the fetal sex and birth outcome
We next assessed maternal miRNAs’ cross-correlations segregated by infant sex (Table 4). The numbers of significantly cross-correlated miRNAs in the plasma from mothers who gave birth to male infants were not altered by PAE in the second trimester. In contrast, maternal miRNAs from births that resulted in female infants exhibited an increased number of significant cross-correlations with PAE in the second trimester. However, in the third trimester, infant outcome (whether the infant was born affected) was the main determinant of increased correlations, such that HEa pregnancies had an increased number of significant correlations for both male and female fetuses. To test the possibility that gestational age at blood draw (variables GAbd1 and GAbd2, Table 1) affected the correlation patterns, we performed partial correlation analysis, controlling for GA, between the maternal miRNAs in each of the exposure groups segregated by infant sex. We observed concordance between correlation and partial correlation estimates (R2 ranged between 0.856 and 0.995, see Additional file 6) suggesting that the variables GAbd1 and GAbd2 did not influence patterns of correlation between miRNAs.
We used bootstrap resampling analyses to determine the stability of the increased cross-correlations. We tested the null hypothesis that the resampling mean of the number of significant correlations in each of the alcohol-exposed groups (HEa or HEua) was equal to the resampling-mean of the non-alcohol-exposed groups (UE). These data show that in the second trimester, PAE increased the coordinated expression of maternal plasma miRNAs in pregnancies with female fetuses, but not those with male fetuses (Fig. 5). In the 2nd trimester, for mothers of female fetuses, the bootstrap distribution of miRNA cross-correlations in HEua mothers was more similar to the distribution observed with HEa group mothers but shifted to become more similar to UE group mothers by mid-third trimester. Maternal miRNAs in pregnancies with male fetuses exhibited an increased number of significant correlations only in the HEa group, and only in the third trimester, when compared to the UE group.
Chromosomal location of correlated miRNA expression
The increased number of miRNAs exhibiting correlated expression in pregnancies with female fetuses suggests that either fetal sex chromosomes or fetal endocrine differences (e.g., ) may influence maternal plasma miRNA. To assess the potential contribution of chromosomal differences, we examined the cross-correlation of miRNA expression by chromosomal location. We identified the most cross-correlated chromosomes by ranking the percent enrichment of significant correlations between chromosomal pairs in the exposure groups compared to the unexposed groups (e.g., HEua vs. UE and HEa vs. UE). We identified the 10 pairs of chromosomes with the highest enrichment for significant miRNA correlations in each exposure group comparison for pregnancies with fetuses of each sex (chromosomal pairs from each group are listed in Additional file 7). In HEa mothers with female fetuses in the second trimester, X-chromosome-linked miRNAs exhibited a preferential increase in the correlated expression with miRNAs encoded on chromosomes 6 and 21 compared to the UE group mothers. In HEua group mothers with female fetuses, X-chromosome-encoded miRNAs showed an increased correlated expression with miRNAs encoded on chromosomes 4, 6, 19, and 21 compared to UE mothers. In contrast, this enrichment for X-chromosome cross-correlation was not seen in the third trimester in either HEa vs. UE or HEua vs. UE miRNAs from pregnancies with female fetuses (Fig. 6).
The genetic sex of a fetus influences important maternal physiological processes like immune tolerance, inflammation, and activation of the hypothalamic-pituitary-adrenal axis [24, 27, 30, 34, 44, 57]. Fetal sex differences in secreted placental factors may also directly influence placental vasculature and the growth and attachment of the placenta , and contribute to the differences in maternal uterine artery blood flow and to maternal vascular adaptations to pathological processes like preeclampsia . Fetal sex can also contribute to the differences in placental function and efficiency in response to maternal disease states like obesity  and asthma . Consistent with these observations, previous studies have also identified fetal sex as an important contributory factor to maternal miRNA profiles associated with other gestational pathologies that result in birthweight alterations . However, as is common with small sample size studies, in the above report, none of the identified sex differences survived correction for multiple comparisons, suggesting a need for alternate statistical approaches to disaggregating samples. Here, we utilize an innovative strategy, bootstrap resampling with replacement, as an approach to disaggregating data from small samples, while approximating a population distribution , to estimate the association between infant sex and the maternal plasma miRNA response to PAE.
Not surprisingly, ANOVA-bootstrap analyses show that most maternal miRNAs that were PAE-responsive were not influenced by fetal sex. For these miRNAs, when the data were examined in aggregate, statistical significance was observed for up to 90% of the resampling trials. This finding suggests a relatively high degree of homogeneity of samples within each exposure group. However, disaggregating data by fetal sex resulted in a decrease in sample size for each group and, consequently, a loss of statistical power. For these non-fetal sex-dimorphic miRNAs, we observed a reduction in the number of resampling iterations which resulted in statistical significance when data were disaggregated by sex, as would be predicted due to decreased sample size.
Our analyses also show that, for some miRNAs, the genetic sex of the fetus does influence maternal plasma miRNA responses to PAE. For these maternal miRNAs, disaggregating the data by fetal sex resulted in an elevation in the proportion of iterations reaching statistical significance, identifying these maternal miRNAs as “likely fetal sex-specific.” However, a majority of the presumptive fetal sexually dimorphic miRNAs reached statistical significance for less than 80% of the resampling iterations, suggesting the possibility that additional sources of variance aside from fetal sex may contribute to the maternal miRNA profiles. Further analyses are needed to define additional contributory factors that influence both basal and PAE-responsive maternal miRNAs. Subsequent post hoc analyses of effect size also indicated that both HEa and HEua pregnancies exhibited fetal sex-specific miRNAs that were significantly different from UE controls and that these fetal sex-specific patterns were observed in both second and third trimesters. These data emphasize the broader importance of assessing the contribution of fetal sex to maternal biomarkers, as part of the strategy for assay development to predict disease outcomes.
We found that PAE increased the number of significant miRNA cross-correlations, an effect which was modified by both infant sex and birth outcome. In the aggregate data (not accounting for fetal sex), we found that HEa group maternal samples showed an increased number of significant miRNA cross-correlations compared to the HEua and UE groups in both second and third trimesters. Bootstrap resampling analysis, segregating the data by infant sex, showed that the increased number of significant correlations in the HEa group in the second trimester was driven mainly by the presence of a female fetus, whereas in the third trimester, both female fetuses and, to a lesser though still statistically significant extent, male fetuses contributed to significant increases in cross-correlations in the HEa group compared to the UE group. Interestingly, the plasma miRNA cross-correlations in pregnant HEua mothers were intermediate between the HEa group at mid-pregnancy but congruent with the UE group by the end of pregnancy. It is unclear whether there is a causal link between correlated maternal miRNAs and infant birth outcomes. However, this finding raises the intriguing possibility that decreases in the correlated expression of plasma miRNAs may reflect biological mechanisms associated with resiliency, i.e., compensatory biological processes to remediate negative growth and neurobehavioral outcomes from maternal alcohol exposure. Conversely, an increased number of plasma miRNA cross-correlations through the end of pregnancy may be associated with adverse infant birth outcomes.
Research groups have explored the relationships between the expression of miRNAs and their target mRNAs, in the context of alcohol exposures [47, 56, 62] and under basal conditions . However, the correlated expression of miRNAs remains poorly investigated and, when observed, has been hypothesized to be due mainly to the proximity of their encoding locations within chromosomes [7, 19]. Such co-secretion from the miRNA cluster within the14q32 chromosome region has been documented to serve as a biomarker for heart disease  and, from the chromosome 19 microRNA cluster, to mediate trophoblast resistance to viral infection . To our knowledge, ours is the first to study to present evidence for correlated expression of miRNAs that are encoded on different chromosomes. For instance, in the 2nd trimester, in HEa and HEua mothers with female fetuses, the expression of X-chromosome encoded miRNAs were significantly correlated with the expression of miRNAs encoded on chromosomes 21 and 6 compared to the UE group.
The correlated behavior of X-chromosome encoded miRNAs was observed specifically in HEa and HEua mothers with female fetuses. However, this finding does not necessarily imply that pregnancies with female fetuses are selectively vulnerable to PAE, since we also observed PAE-linked changes in maternal miRNA profiles in pregnancies with both male and female fetuses following bootstrap-based assessments of effect size. A number of studies have shown that male fetal and postnatal development can also be uniquely sensitive to altered maternal environments due to drug and toxin exposures [6, 23, 25, 55, 65]. However, our analysis does suggest that the toxic effects of PAE may be mediated by different mechanisms in male and female fetuses.
The increased number of miRNA cross-correlations also does not imply that the miRNAs themselves are a causal link between PAE and FASD outcomes. However, the factors that contribute to the correlated miRNA expression may also contribute to the adverse birth outcomes. Firstly, the correlated expression in maternal plasma may be due to the co-transcription of miRNAs in one or a few tissues, such as the placenta, or coordinated secretion of miRNAs from multiple maternal and/or fetal tissues. In the first instance, the increased correlation may be facilitated by the alcohol-driven reorganization of chromatin structure. Chromatin is thought to be structured into topologically associated domains which act as gene regulatory “hotspots” , where eukaryotic genes located on separate chromosomes can physically interact  and adaptively interact in response to exogenous cues such has hormonal treatments , resulting in correlated gene expression changes. Such transcription factories may emerge in a limited set of tissues. For instance, we found evidence in a previous study that a number of miRNAs elevated in maternal circulation in response to alcohol exposure were preferentially transcribed in the placenta . However, correlated patterns of maternal miRNA may also emerge due to stimulated co-secretion from a number of different tissues. Circulating cytokines constitute a candidate class of molecules that may coordinate the secretion of miRNAs into the plasma. For instance, in the same cohort of pregnant women assessed in this study, another group observed that prenatal alcohol exposure also results in elevated levels of circulating cytokines in pregnant women who subsequently gave birth to affected infants . In experimental studies of liver function, cytokines have been shown to induce the release of extracellular miRNAs in exosomes . Consequently, cytokines, perhaps even fetal sex-dimorphic maternal cytokine profiles , and extracellular miRNAs may constitute a coordinated signal for future developmental disability due to prenatal alcohol exposure.
While this study does provide preliminary support for the inclusion of fetal sex in predictive models based on extracellular miRNA profiling in mothers, there are several caveats that need to be considered. The assumption behind bootstrap resampling, that the iterative process results in an approximation of the characteristics of a population, may not be generalizable. This method simulates a population to estimate the stability of a statistic, and therefore, the results may not be generalizable to other populations not reflected in the initial sample. Study results may be true for pregnant women from the selected regions in western Ukraine but may not hold true for pregnant women from other geographic locales or from diverse ethnic groups. Variation in the composition of additional risk factors like education, socioeconomic status, nutrition, and poly-substance use may also limit the generalizability of the specific miRNAs identified in the current study population. For instance, in our study population, a higher proportion of HEa and HEua mothers reported smoking during pregnancy compared to UE mothers which is consistent with the evidence for a strong association between tobacco and alcohol consumption behavior , though among the smokers in each group, the cigarette consumption levels were similar. The co-use of other psychoactive agents with alcohol may also modify the maternal miRNA profiles. A higher proportion of UE mothers also attained a college education and though not significant overall, UE mothers also tended towards higher socioeconomic status. These factors could be protective and further contribute to changes in maternal miRNA levels. In addition, we classified infants as HEa and HEua based on physical features, growth, and neurodevelopmental assessments through the first year of life. As some alcohol-related neurodevelopmental deficits are not easily detectable until later in childhood, some infants who were truly affected may have been misclassified as unaffected. We would expect this misclassification to result in the attenuation of some of the identified associations. Nevertheless, these limitations do not minimize the need to account for the role of significant variables like fetal sex in the emergence of patterns of biomarkers that predict infant disability.
Perspectives and significance
This study shows that statistical techniques like bootstrap resampling with replacement can be used successfully, to disaggregate data from small sample FASD studies to generate novel hypotheses. This study opens the door to examining the contributions of other significant maternal health variables to maternal miRNA profiles to make predictive models for the behavior of these biomarkers. These analyses also uncover potentially novel biology of miRNA co-synthesis or co-secretion as an adaptive response to the maternal environment, a phenomenon that predicts diminished fetal resiliency, and needs further investigation. It will also be important to understand why the sex of the fetus results in a change in a maternal miRNA response to PAE and if there is a common sex-dependent behavior in maternal miRNAs responses to other adverse environmental perturbations as well. Our previous studies showed that maternal miRNAs can be biologically active and control the growth of important fetal tissues like the placenta . Consequently, it will be important to ascertain what, if any, biological processes are influenced by such differences. In conclusion, the statistical approaches outlined in this study provide a means of generating hypotheses about the influence of biological variables in studies that were not initially designed to investigate such variables.
Availability of data and materials
Data on human subjects is deposited at CIFASD.org, in accordance with the NIH data repository guidelines, and is available to outside investigators following an application for data use (https://cifasd.org/data-sharing/). All “R” code that was developed or adapted for use in these analyses has been deposited at GitHub (https://github.com/nihalasalem/Ukraine_Analysis).
Analysis of variance
Fetal alcohol spectrum disorders
Heavily exposed with affected infants
Heavily exposed with apparently unaffected infants
Mature miRNA miRBase sequence entry unique accession number
Prenatal alcohol exposure
Unexposed to alcohol
Bakhireva LN, Garrison L, Shrestha S, Sharkis J, Miranda R, Rogers K. Challenges of diagnosing fetal alcohol spectrum disorders in foster and adopted children. Alcohol. 2017a;67:37–43.
Bakhireva LN, Schatz M, Jones KL, Tucker CM, Slymen DJ, Klonoff-Cohen HS, Gresham L, Johnson D, Chambers CD, Group, O.C.R. Fetal sex and maternal asthma control in pregnancy. J Asthma. 2008;45:403–7.
Bakhireva LN, Sharkis J, Shrestha S, Miranda-Sohrabji TJ, Williams S, Miranda RC. Prevalence of prenatal alcohol exposure in the State of Texas as assessed by phosphatidylethanol in newborn dried blood spot specimens. Alcohol Clin Exp Res. 2017b;41:1004–11.
Balaraman S, Lunde ER, Sawant O, Cudd TA, Washburn SE, Miranda RC. Maternal and neonatal plasma microRNA biomarkers for fetal alcohol exposure in an ovine model. Alcohol Clin Exp Res. 2014;38:1390–400.
Balaraman S, Schafer JJ, Tseng AM, Wertelecki W, Yevtushok L, Zymak-Zakutnya N, Chambers CD, Miranda RC. Plasma miRNA profiles in pregnant women predict infant outcomes following prenatal alcohol exposure. PLoS One. 2016;11:e0165081.
Bara A, Manduca A, Bernabeu A, Borsoi M, Serviado M, Lassalle O, Murphy M, Wager-Miller J, Mackie K, Pelissier-Alicot A-L, et al. Sex-dependent effects of in utero cannabinoid exposure on cortical function. eLife. 2018;7:e36234.
Baskerville S, Bartel DP. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA. 2005;11:241–7.
Bayer A, Lennemann NJ, Ouyang Y, Sadovsky E, Sheridan MA, Roberts RM, Coyne CB, Sadovsky Y. Chromosome 19 microRNAs exert antiviral activity independent from type III interferon signaling. Placenta. 2018;61:33–8.
Belyaeva A, Venkatachalapathy S, Nagarajan M, Shivashankar G, Uhler C. Network analysis identifies chromosome intermingling regions as regulatory hotspots for transcription. Proc National Acad Sci. 2017;114:13714–9.
Benes V, Collier P, Kordes C, Stolte J, Rausch T, Muckentaler MU, Häussinger D, Castoldi M. Identification of cytokine-induced modulation of microRNA expression and secretion as measured by a novel microRNA specific qPCR assay. Scientific Reports. 2015;5:11590.
Bland JM, Altman DG. Statistics notes: bootstrap resampling methods. BMJ. 2015;350:h2622.
Bobo JK, Husten C. Sociocultural influences on smoking and drinking. Alcohol Res Health. 2000;24:225–32.
Bodnar TS, Raineki C, Wertelecki W, Yevtushok L, Plotka L, Zymak-Zakutnya N, Honerkamp-Smith G, Wells A, Rolland M, Woodward TS, et al. Altered maternal immune networks are associated with adverse child neurodevelopment: impact of alcohol consumption during pregnancy. Brain Behavior Immunity. 2018;73:205–15.
Broere-Brown ZA, Schalekamp-Timmermans S, Hofman A, Jaddoe V, Steegers E. Fetal sex dependency of maternal vascular adaptation to pregnancy: a prospective population-based cohort study. BJOG. 2016;123:1087–95.
Bronson SL, Bale TL. The placenta as a mediator of stress effects on neurodevelopmental reprogramming. Neuropsychopharmacology. 2016;41:207–18.
Brown ZA, Schalekamp-Timmermans S, Tiemeier HW, Hofman A, Jaddoe VW, Steegers EA. Fetal sex specific differences in human placentation: a prospective cohort study. Placenta. 2014;35:359–64.
Chambers CD, Yevtushok L, Zymak-Zakutnya N, Korzhynskyy Y, Ostapchuk L, Akhmedzhanova D, Chan PH, Xu R, Wertelecki W. Prevalence and predictors of maternal alcohol consumption in 2 regions of Ukraine. Alcohol Clin Exp Res. 2014;38:1012–9.
Chasnoff IJ, Wells AM, King L. Misdiagnosis and missed diagnoses in foster and adopted children with prenatal alcohol exposure. Pediatrics. 2015;135:264–70.
Chaulk SG, Ebhardt HA, Fahlman RP. Correlations of microRNA: microRNA expression patterns reveal insights into microRNA clusters and global microRNA expression patterns. Molecular BioSystems. 2016;12:110–9.
Coles CD, Kable JA, Keen CL, Jones KL, Wertelecki W, Granovska IV, Pashtepa AO, Chambers CD, Cifasd. Dose and timing of prenatal alcohol exposure and maternal nutritional supplements: developmental effects on 6-month-old infants. Matern Child Health J. 2015;19:2605–14.
Cui M, Xu L, Wang H, Ju S, Xu S, Jing R. Combining Nordtest method and bootstrap resampling for measurement uncertainty estimation of hematology analytes in a medical laboratory. Clin Biochem. 2017;50:1067–72.
Davidson R, MacKinnon JG. Bootstrap tests: how many bootstraps? Econometric Reviews. 2000;19:55–68.
Dietrich KN, Krafft KM, Bornschein RL, Hammond PB, Berger O, Succop PA, Bier M. Low-level fetal lead exposure effect on neurobehavioral development in early infancy. Pediatrics. 1987;80:721–30.
DiPietro JA, Costigan KA, Kivlighan KT, Chen P, Laudenslager ML. Maternal salivary cortisol differs by fetal sex during the second half of pregnancy. Psychoneuroendocrinology. 2011;36:588–91.
DiPietro JA, Voegtline KM. The gestational foundation of sex differences in development and vulnerability. Neuroscience. 2017;342:4–20.
Efron B. Bootstrap methods: another look at the jackknife. In: Kotz S, Johnson NL, editors. Breakthroughs in statistics: methodology and distribution. New York: Springer New York; 1992. p. 569–93.
Enninga EA, Nevala WK, Creedon DJ, Markovic SN, Holtan SG. Fetal sex-based differences in maternal hormones, angiogenic factors, and immune mediators during pregnancy and the postpartum period. Am J Reprod Immunol. 2015;73:251–62.
Finer LB, Zolna MR. Declines in unintended pregnancy in the United States, 2008-2011. New England J Med. 2016;374:843–52.
Gardiner AS, Gutierrez HL, Luo L, Davies S, Savage DD, Bakhireva LN, Perrone-Bizzozero NI. Alcohol use during pregnancy is associated with specific alterations in microRNA levels in maternal serum. Alcohol Clin Exp Res. 2016;40:826–37.
Giesbrecht GF, Campbell T, Letourneau N, Team APS. Sexually dimorphic adaptations in basal maternal stress physiology during pregnancy and implications for fetal development. Psychoneuroendocrinol. 2015;56:168–78.
Grant BF, Goldstein RB, Saha TD, Chou SP, Jung J, Zhang H, Pickering RP, Ruan WJ, Smith SM, Huang B, et al. Epidemiology of DSM-5 alcohol use disorder: results from the National Epidemiologic Survey on Alcohol and Related Conditions III. JAMA Psychiatry. 2015;72:757–66.
Gross N, Kropp J, Khatib H. Sexual dimorphism of miRNAs secreted by bovine in vitro-produced embryos. Front Genetics. 2017;8:39.
Hedges LV. Distribution theory for Glass’s estimator of effect size and related estimators. J Educ Stat. 1981;6:107–28.
Kieffer TEC, Laskewitz A, Faas MM, Scherjon SA, Erwich J, Gordijn SJ, Prins JR. Lower FOXP3 mRNA expression in first-trimester decidual tissue from uncomplicated term pregnancies with a male fetus. J Immunol Res. 2018;2018:1950879.
Lange S, Probst C, Gmel G, Rehm J, Burd L, Popova S. Global prevalence of fetal alcohol spectrum disorder among children and youth: a systematic review and meta-analysis. JAMA Pediatr. 2017;171:948–56.
Le Dily F, Baù D, Pohl A, Vicent GP, Serra F, Soronellas D, Castellano G, Wright RHG, Ballare C, Filion G, et al. Distinct structural transitions of chromatin topological domains correlate with coordinated hormone-induced gene regulation. Genes Development. 2014;28:2151–62.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433:769–73.
Mando C, Calabrese S, Mazzocco MI, Novielli C, Anelli GM, Antonazzo P, Cetin I. Sex specific adaptations in placental biometry of overweight and obese women. Placenta. 2016;38:1–7.
Markiewicz PJ, Reader AJ, Matthews JC. Assessment of bootstrap resampling performance for PET data. Phys Med Biol. 2015;60:279–99.
Martínez-Micaelo N, Beltrán-Debón R, Baiges I, Faiges M, Alegret JM. Specific circulating microRNA signature of bicuspid aortic valve disease. Journal of Translational Medicine. 2017;15:76.
May PA, Chambers CD, Kalberg WO, Zellner J, Feldman H, Buckley D, Kopald D, Hasken JM, Xu R, Honerkamp-Smith G, et al. Prevalence of fetal alcohol spectrum disorders in 4 US communities. JAMA. 2018;319:474–82.
May PA, Tabachnick B, Hasken JM, Marais A-S, de Vries MM, Barnard R, Joubert B, Cloete M, Botha I, Kalberg WO. Who is most affected by prenatal alcohol exposure: boys or girls? Drug and alcohol dependence. 2017;177:258–67.
Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J. A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009;10:R64.
Mitchell AM, Palettas M, Christian LM. Fetal sex is associated with maternal stimulated cytokine production, but not serum cytokine levels, in human pregnancy. Brain Behavior Immunity. 2017;60:32–7.
Moore EM, Riley EP. What happens when children with fetal alcohol spectrum disorders become adults? Curr Dev Disord Rep. 2015;2:219–27.
Murphy VE, Gibson PG, Giles WB, Zakar T, Smith R, Bisits AM, Kessell CG, Clifton VL. Maternal asthma is associated with reduced female fetal growth. Am J Respiratory Critical Care Med. 2003;168:1317–23.
Nunez YO, Truitt JM, Gorini G, Ponomareva ON, Blednov YA, Harris RA, Mayfield RD. Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence. BMC Genomics. 2013;14:725.
O’Shaughnessy PJ, Antignac JP, Le Bizec B, Morvan ML, Svechnikov K, Soder O, Savchuk I, Monteiro A, Soffientini U, Johnston ZC, et al. Alternative (backdoor) androgen production and masculinization in the human fetus. PLoS Biol. 2019;17:e3000002.
Paolozza A, Munn R, Munoz DP, Reynolds JN. Eye movements reveal sexually dimorphic deficits in children with fetal alcohol spectrum disorder. Front Neurosci. 2015;9:76.
Pei X, Sze NN, Wong SC, Yao D. Bootstrap resampling approach to disaggregate analysis of road crashes in Hong Kong. Accid Anal Prev. 2016;95:512–20.
Raineki C, Chew L, Mok P, Ellis L, Weinberg J. Short- and long-term effects of stress during adolescence on emotionality and HPA function of animals exposed to alcohol prenatally. Psychoneuroendocrinol. 2016;74:13–23.
Riley EP, Infante MA, Warren KR. Fetal alcohol spectrum disorders: an overview. Neuropsychol Review. 2011;21:73–80.
Rodosthenous RS, Burris HH, Sanders AP, Just AC, Dereix AE, Svensson K, Solano M, Tellez-Rojo MM, Wright RO, Baccarelli AA. Second trimester extracellular microRNAs in maternal blood and fetal growth: an exploratory study. Epigenetics. 2017;12:804–10.
Romano JP, Wolf M. Exact and approximate stepdown methods for multiple hypothesis testing. Journal of the American Statistical Association. 2005;100:94–108.
Sandman CA, Glynn LM, Davis EP. Is there a viability-vulnerability tradeoff? Sex differences in fetal programming. J Psychosom Res. 2013;75:327–35.
Sathyan P, Golden HB, Miranda RC. Competing interactions between micro-RNAs determine neural progenitor survival and proliferation after ethanol exposure: evidence from an ex vivo model of the fetal cerebral cortical neuroepithelium. J Neurosci. 2007;27:8546–57.
Scott NM, Hodyl NA, Murphy VE, Osei-Kumah A, Wyper H, Hodgson DM, Smith R, Clifton VL. Placental cytokine expression covaries with maternal asthma severity and fetal sex. J Immunol. 2009;182:1411–20.
Shiozaki A, Matsuda Y, Satoh S, Saito S. Impact of fetal sex in pregnancy-induced hypertension and preeclampsia in Japan. J Reprod Immunol. 2011;89:133–9.
Wei T, Simko V. R package “corrplot”: visualization of a correlation matrix; 2017.
Spilianakis CG, Lalioti MD, Town T, Lee GR, Flavell RA. Interchromosomal associations between alternatively expressed loci. Nature. 2005;435:637–45.
Steier JA, Myking OL, Bergsjo PB. Correlation between fetal sex and human chorionic gonadotropin in peripheral maternal blood and amniotic fluid in second and third trimester normal pregnancies. Acta Obstet Gynecol Scand. 1999;78:367–71.
Tsai PC, Bake S, Balaraman S, Rawlings J, Holgate RR, Dubois D, Miranda RC. MiR-153 targets the nuclear factor-1 family and protects against teratogenic effects of ethanol exposure in fetal neural stem cells. Biology open. 2014;3:741–58.
Tseng AM, Mahnke AH, Wells AB, Salem NA, Allan AM, Roberts VH, Newman N, Walter NA, Kroenke CD, Grant KA, et al. Maternal circulating miRNAs that predict infant FASD outcomes influence placental maturation. Life Science Alliance. 2019;2:e201800252.
Varlinskaya EI, Mooney SM. Acute exposure to ethanol on gestational day 15 affects social motivation of female offspring. Behav Brain Res. 2014;261:106–9.
Wang J-Q, Hu Y-B, Gao H, Sheng J, Huang K, Zhang Y-W, Mao L-J, Zhou S-S, Cai X-X, Zhang L-J, et al. Sex-specific difference in placental inflammatory transcriptional biomarkers of maternal phthalate exposure: a prospective cohort study. J Exposure Sci Environm Epidemiol. 2020;30:835–44.
Wilcox RR. Fundamentals of modern statistical methods: substantially improving power and accuracy. New York: Springer; 2010.
Zheng Q, Deng Y, Zhong S, Shi Y. Human chorionic gonadotropin, fetal sex and risk of hypertensive disorders of pregnancy: a nested case-control study. Pregnancy Hypertens. 2016;6:17–21.
This research was supported by grants from the NIH, U01AA014835, and the Office of Dietary Supplements (Christina D. Chambers), R01AA024659 (Rajesh C. Miranda), F31AA026505 (Alexander M. Tseng), and F99NS113423 (Nihal A. Salem).
Ethics approval and consent to participate
All study protocols were approved by the Institutional Review Boards at the University of California San Diego and Texas A&M University and for data collection at Lviv National Medical University in Ukraine.
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Selection of the number of bootstrap resampling iterations using second trimester UE group. (a) The frequency distributions of the number of significant miRNA cross-correlations in each of the shown number of iterations. (b) Mean of the number of significant miRNA cross-correlations in each of the shown number of iterations. (c) Mean of the number of significant correlations in each run of 20 simulations using the shown number of resampling iterations.
Likely alcohol-sensitive, fetal sex-specific miRNAs (ANOVA-Bootstrap). Table includes MIMAT numbers, chromosomal location and names of the likely alcohol-sensitive, fetal sex-specific miRNAs for each sex in HEa and HEua in each of the second and third trimesters.
Concordance between ANCOVA-bootstrap & ANOVA bootstrap. X-Y scatterplot showing the relation between the proportion of ANCOVA significant iterations (x-axis) and the proportion of ANOVA significant iterations (y-axis)
Bootstrap-resampling to determine effect size and confidence interval estimates for the effects of fetal sex on maternal HEua miRNA expression. Median bootstrap effect sizes (solid dot) and 95% confidence interval (error bars) for miRNAs which have a non-zero containing confidence interval during sex-segregated resampling. Second and third trimester, HEua vs. UE effect sizes are shown. Dashed line indicates zero effect size.
Significant alcohol-sensitive, fetal sex-specific miRNAs (Effect-Size Bootstrap). Table includes sex-segregated miRNAs that exhibited median effect sizes that were outside the range of the 95% confidence estimate for the composite sample. The table shows the median effect size for each of the composite and sex-segregated analyses.
Concordance between partial correlation analyses correcting for gestational age (GA) at blood draw and full correlation analyses estimates. Figures showing full correlation coefficient between each miRNA pair on X-axis, and partial correlation coefficient correcting for GA on Y-axis, blue line shows the linear regression between the two coefficients. R2 between the two coefficients are shown. Second trimester groups are shown in (a), third trimester groups are shown in (b).
Table showing the 10 chromosome pairs with the greatest increase in significant cross correlations in either the HEa or HEua groups, relative to the UE group, for pregnancies with male and female fetuses. Percentage of chromosome pair change (%change)= [((no. of significant cross-correlations between miRNAs on a chromosome pair in alcohol exposed group) – (no. of significant cross-correlations between miRNAs on the same chromosome pair in the control, UE, group) * 100)/(no. of significant cross-correlations between miRNAs on the same chromosome pair in the control group)].
About this article
Cite this article
Salem, N.A., Mahnke, A.H., Wells, A.B. et al. Association between fetal sex and maternal plasma microRNA responses to prenatal alcohol exposure: evidence from a birth outcome-stratified cohort. Biol Sex Differ 11, 51 (2020). https://doi.org/10.1186/s13293-020-00327-2
- Bootstrap resampling
- Sex as a biological variable
- Extracellular miRNAs
- Fetal alcohol spectrum disorders
- Maternal miRNA co-secretion