Introduction
A longstanding mystery surrounds the epidemiologic observation that human chorionic gonadotropin (hCG) – a placental glycoprotein hormone – is higher in maternal circulation in the case of a female than a male fetus.Reference Bremme, Lagerström, Andersson, Johansson and Eneroth 1 – Reference Yaron, Lehavi and Orr-Urtreger 5 This has been proposed as an explanation as to why women carrying a female fetus are more likely to suffer from hyperemesis gravidarum (i.e. nausea) than those carrying a male fetus;Reference Kauppila, Huhtaniemi and Ylikorkala 6 but otherwise has not been elucidated in terms of clinical significance. The stability of placental sex differences within and across populations is also unknown. We present here a population-level analysis of sexual dimorphic expression of hCG.
Placental hormones have been identified as relevant to sex-specific fetal development. hCG can bind the luteinizing hormone/hCG receptor (LHCGR) and stimulate steroidogenesis in the fetal testis in the first trimester of pregnancy.Reference Clements, Reyes, Winter and Faiman 2 , Reference Huhtaniemi, Korenbrot and Jaffe 7 This interaction is essential to normal male development in that defects in hCG-LHCGR binding, due to mutations in the LHCGR gene, are associated with hypogonadism and pseudohermaphroditism.Reference Scott, Mason and Sharpe 8 It is theorized that a male-specific negative feedback mechanism between the fetus and the placenta explains lower circulating levels in the male v. the female.Reference Clements, Reyes, Winter and Faiman 2 There is little known on the relationship of placental hCG and female-specific development. In a comparison study of hCG in fetal tissues, the highest concentrations were detected in the fetal ovary – ~five-fold higher than in the testes.Reference Huhtaniemi, Korenbrot and Jaffe 9 Hyperglycosylated hCG can bind and activate the TGF-β receptor.Reference Berndt, Blacher and Munaut 10 TGF-β plays a role in fetal ovary development and possibly in fetal origins of polycystic ovarian syndrome.Reference Hatzirodos, Bayne and Irving-Rodgers 11
The hCG observation poses an important epidemiologic question, as to whether fetal sex differences in placental function are primarily a genetic/epigenetic phenomenon (i.e. the presence/absence of the Y-chromosome and X-inactivation), or might female/male differences in hCG and other proteins also be indicators of non-genetic influences on the fetal environment such as nutrition, maternal health, maternal stress or chemical exposures. Circulating hCG levels can change in response to chemical exposures,Reference Ridano, Racca and Flores-Martín 12 – Reference Fowler, Bhattacharya, Gromoll, Monteiro and O’Shaughnessy 14 maternal stress,Reference Catalano, Saxton and Bruckner 15 maternal nutrition,Reference Lewis, Demmelmair and Gaillard 16 exogenous estrogen exposureReference Jukic, Weinberg, Baird and Wilcox 17 and other factors that can modulate the endocrine milieu of pregnancy.
hCG was put into widespread use over three decades based on its predictive power in estimating risk of fetal defects such as Down’s syndrome.Reference Brock and Sutcliffe 18 It later was observed that it can be useful in predictions of non-genetic, adverse pregnancy outcomes such as preterm labor,Reference Jelliffe-Pawlowski, Shaw and Currier 19 preeclampsiaReference Taché, Baer and Currier 20 and placental abruption.Reference Blumenfeld, Baer and Druzin 21 Lower second trimester hCG was correlated with increased risk of cryptorchidism in male infants as compared with normal controls.Reference Chedane, Puissant, Weil, Rouleau and Coutant 22 The ability of hCG to predict non-genetic fetal disorders has become a new possibility given that it is measured at the population level, and data has been collected for more than a decade. First trimester hCGβ first went into widespread use as a prenatal screening analyte in the 1990s.
The aim of the current analysis was to robustly characterize sex differences in hCG within a well enumerated, racially/ethnically and economically diverse US population. The second aim was to explore the stability of the hCG sex difference within subgroups of the population that might represent genetic as well as non-genetic differences related to social disadvantage, lifestyle, maternal stress, health, nutrition and chemical exposures.
Method
Genetic Disease Screening Program (GDSP)
Study subjects were residents of the State of California who: (a) underwent routine prenatal serum screening in the first and/or second trimesters of pregnancy from April 1, 2009 through December 2012, (b) were carrying singleton pregnancies, (c) delivered their babies in California and (d) participated in the California Newborn Screening Program. For this analysis, the sex of the baby was abstracted from the newborn screening database and merged with the prenatal screening database. The first trimester blood draw was between 10 and 13 weeks 6 days gestation (median time of prenatal screen was 12.0 weeks) and the second trimester blood draw was between 15 and 20 weeks gestation (median time of prenatal screen was 17.0 weeks). Blood samples were shipped to one of the six regional laboratories, where they were analyzed for a panel of proteins using automated methods. The Genetic Disease Laboratory reviewed results before transmission to the GDSP database. Total hCG was measured by the AutoDELFIA® time-resolved fluoroimmunoassay (PerkinElmer, Waltham, MA, USA). All data were entered directly into a centralized database along with information collected on the sample requisition form related to maternal demographics, maternal weight and health. Analyte-specific multiple of the medians (MoM) are calculated centrally, and referred to here as the clinical MoM. The Committee for the Protection of Human Subjects of the California Department of Public Health approved the protocol for this study. Data were abstracted in October 2013. All data were de-identified by the California Genetic Disease Screening Program before transfer to the collaborating institution.
Gestational age multiple of the median (ga-MoM)
The clinical MoM’s are normalized values of the serum analytes that were designed for the calculation of risks of Down’s Syndrome and other types of birth defects.Reference Wald, Cuckle and Brock 23 , Reference Wald, Rodeck and Hackshaw 24 Epidemiologic analyses have used the clinical MoM more than the analyte values. The concept behind the MoM is to compare each person’s lab value to the expected value for that individual for that gestational day. Information used in the normalization process includes dating method (nuchal translucency, ultrasound, last menstrual period, physical examination), maternal race/ethnicity, weight of the mother, insulin-dependent diabetes, smoking and ovum donor status. This information is collected through a one-page requisition form that the physician completes at the time of the screening. The result is a unitless ratio with a median of 1.0. The database is monitored periodically for analytical drift in the assays and shifts over time in gestational age and other subgroup medians. As needed, the medians are revised in order to maintain the population MoM median value at 1. For the purposes of this analysis, we calculated a ga-MoM that is only adjusted for gestational age and not the other factors. This allows us to model the effects of the other factors on the MoM without the risk of overadjustment. To calculate the ga-MoM, we fit a model of the log median analyte, calculated for each day of gestation, regressed on the day of gestation using a smoothing function. For first trimester hCG, a linear term for gestational day was included in the model. For second trimester hCG, linear and quadratic terms for gestational age were included. An exponential decay to a non-zero constant was tried for hCG but the quadratic model fit the data better. We output the predicted smoothed gestational age medians for the 3-year time period covered in our dataset. Each person’s analyte value was divided by the smoothed gestational age median to generate the ga-MoM.
Statistical analysis
The associations of maternal and fetal factors (fetal sex, maternal race/ethnicity, weight, age, smoking, diabetes and ovum donor status) with log analyte and ga-MoM values were calculated using linear regression models. All estimates were adjusted for dating method, fetal sex, maternal race/ethnicity, maternal age, the inverse of maternal weight in pounds, year of estimated due date, month of screening, smoking, diabetes and ovum donor status. The goal of the modeling approach was to describe, using the information available, the variability in the analyte and in the ga-MoM values. Race groups were generated according to the criteria used by the GDSP. Maternal age categories are according to clinical cut-offs. Maternal weight categories are based on quartiles calculated in this population. We output predicted means for each subgroup and predicted values for each subject on the log scale and then back transformed them to the linear scale. We used the predicted values per subject in order to calculate standard errors in the original units. We calculated fold differences for each category from the referent group on the linear scale. To test the interaction of these factors with fetal sex, we used the same strategy as above and included in the model an interaction term for fetal sex by maternal factor. We output predicted means for females and males, and calculated the female/male ratio for each subgroup. P-values are reported for the main effects as well as for the interactions. Mixed effects models were used to estimate the intraclass correlation coefficient (ICC) for first and second trimester hCG. To visually inspect the change in the fetal sex ratio over time in gestation, we calculated a ratio of median female to median male analyte values for each day of prenatal screening. All analyses were carried out in SAS 9.2 (Cary, NC, USA) and plots were generated in R (R Foundation, www.r-project.org).
Results
The sample presented here was restricted to those pregnancies that resulted in a livebirth; yet, the clinical MoM’s are based on all pregnancies in California. The median values for the MoM’s were equal to 1.00, indicating that our sample was not significantly different from the underlying population, even though it excluded pregnancies that ended in miscarriage, termination or otherwise were not captured in the neonatal screening program. Subject characteristics were highly consistent within the two screening intervals (Table 1). More women had second than first trimester screening, and 38% of women had screening in both trimesters. The ICC for the hCG ga-MoM was 0.74 (95% CI 0.74, 0.74), indicating fairly high reproducibility of hCG within a pregnant woman over time.
Sample characteristics of participants in Genetic Disease Screening Program of the California Department of Public Health
In the first row for first and second trimesters in Table 2, we present the estimate for the overall sex ratio for each analyte. This ratio represented the overall, or genetically dominant, fetal sex difference. The ratios below it are the subgroup fluctuations from the genetic baseline. Associations were essentially the same when calculated using the analyte values adjusted for gestational age (Table 2) and the ga-MoM values (data not shown). Demographic differences by race, weight, smoking, insulin-dependent diabetes and ovum donor status were consistent with what has been reported previously.Reference Wald, Rodeck and Hackshaw 24 – Reference Lambert-Messerlian, Dugoff and Vidaver 28
Predicted mean values of hCG (IU/l), by maternal factors, predicted mean values by placental-fetal sex, and female/male ratios in hCG levels and their fluctuations across subgroups
hCG, human chorionic gonadotropin.
a All P-values for differences with referent group are <0.0001.
Fluctuations in female/male ratios by hormone
hCG was 11% higher in female pregnancies in the first trimester and 8% higher in the second trimester (Table 2, Fig. 1). In first trimester, the greatest female/male ratio was in White women (13%) v. Asian women (9%). In second trimester, the female/male ratio was highest in Black women (10%) and lowest in Asian women (3%) (Fig. 2). All estimates were adjusted for gestational age at the time of the blood draw, maternal weight, age, year and month of blood draw, diabetes and smoking status. The female/male ratio increased with increasing maternal age in the first trimester (9–12%), but not in the second trimester (Fig. 2). Conversely, maternal weight had a stronger effect on the ratio in the second trimester (5–9%) as compared with the first (Fig. 2). We did not detect significant fluctuations in the ratio by smoking, diabetes or ovum donor even though the overall effects of these exposures on circulating hCG were large, ranging from 16% lower for first trimester insulin-dependent diabetes to 39% higher for second trimester ovum donor status.
Placental-fetal sex differences in circulating hCG in first and second trimester pregnancies as depicted by (a) hCG serum concentration (IU/l) and (b) normalized hCG. The dots and lines are median and 25th/75th percentile values for each gestational day of prenatal screening. Female values are denoted in gray, and male values are denoted in black. hCG, human chorionic gonadotropin; MoM, multiple of the medians.
Predicted means (± standard errors) in female and male placental hCG multiples of the median, normalized for ga-MoM, across categories of maternal race, weight and age in the first and second trimesters. All of these plots represent significant fluctuations in the female-to-male ratio across categories (P<0.01). hCG, human chorionic gonadotropin; ga-MoM, gestational age multiple of the median.
To assess the change in the female/male ratio over time in pregnancy, we fit a model to estimate the effect of gestational age on the female/male ratio in analyte values adjusting for other factors that may also be changing over time. The female/male ratio in hCG differed significantly with gestational day. The hCG ratio decreased from 17 to 2% in the first trimester, and increased from 2 to 19% in the second trimester (P<0.0001, Fig. 3).
Variability in female/male ratio in hCG by gestational age. Black dots represent the ratio calculated for each gestational day for the median female hCG to the median male hCG. hCG, human chorionic gonadotropin.
Discussion
We report fetal sex differences in hCG measured in early pregnancy, and the degree to which those differences fluctuate by maternal factors and by gestational age. The female/male ratio in hCG was the most variable across racial subgroups, with the smallest differences in Asian women and women younger than 20 years, and the maximum differences in African-American women. We did not detect significant fluctuations in female/male ratios by smoking, ovum donor status or insulin-dependent diabetes, which may indicate a lack of susceptibility of sex-specific hCG regulation to these types of adversity. It may also be due to the extremely low prevalence of these three exposures in the dataset (⩽1%). These types of clinical data are widely available and offer a novel framework to study maternal–placental–fetal interactions that may be responsive to genetic, physiologic and environmental influences and relevant to obstetric complications and fetal outcomes.
Approximately 75–80% of pregnancies in California during this time period are represented in our dataset, giving us confidence that we are reporting unconfounded differences that have a basis in sexually dimorphic placental function. The female/male ratios were similar with or without adjustment for maternal factors; however, they did change considerably when we calculated group-level means v. individual-level predicted values, which were then summarized as means. Proper adjustment for individual-level gestational age is critical in making valid inferences on dynamic relationships such as these in early pregnancy. It is possible that previous studies did not detect sex differences in hCG in early pregnancy because they did not incorporate individual-level information on gestational age, but instead conducted tests of group means.Reference Steier, Myking and Bergsjo 4
Race was a major source of variability in hCG as well as in the fluctuations in the female/male ratios. This may relate to polymorphic genes that encode synthesis and/or metabolism of hCG that may also cluster by race or ethnic groups; although none have been identified to our knowledge. In the context of CA, we can speculate on explanations for higher mean levels of hCG in Asian women v. Hispanic women, after controlling for other factors, that might include: (a) higher proportions of farm workers among Hispanic women or members of families of farm workers with higher than average pesticide exposure; and (b) higher proportions of Asians living in coastal cities and higher fish consumption as compared with Hispanics. Hispanic women had 45% higher odds of neural tube defects (NTDs) than White women, which may also indicate different experiences in pregnancy with respect to environmental factors and their effects on placental–fetal hormones.Reference Feuchtbaum, Carter, Dowray, Currier and Lorey 29 Risk of NTDs, anencephaly in particular, is higher for female fetuses.Reference Stumpf, Cranford and Elias 30 Similar to US trends, African-American women in CA have the highest rates of preterm birth and infant mortality which are two outcomes that also may have origins in early maternal–placental–fetal interactions and are known to have a sex bias.Reference Sabol, de Sam Lazaro and Salati 31 – Reference Cooperstock and Campbell 33 Another difference that could contribute to differences in hCG regulation by race might be maternal and/or paternal nativity or time since initial migration into California.Reference Becerra, von Ehrenstein and Heck 34 – Reference Ahmed, Quinn and Caan 36 This cannot be directly studied in this database, but should be considered in future analyses.
Maternal weight is a known contributor to the endocrine milieu of pregnancy,Reference Lewis, Demmelmair and Gaillard 16 and maternal age is tightly correlated with ovarian function.Reference Donnenfeld, Icke, Pargas and Dowman 37 Both weight and age were correlated with fluctuations in the hCG sex difference. This could be a function of cumulative environmental exposures, chemical body burden, physiologic changes or a combination of all three. Information on body mass index and gestational weight gain was not available.
In order to alter circulating hCG levels consistently, the non-genetic or genetic factor would need to work at the level of methylation, transcription of the genes [chorionic gonadotropin α (CGA), chorionic gonadotropin β (CGB, CGB1, CGB2, CGB5, CGB7)] that encode the protein, the translation of the mRNA into protein, the packaging and secretion of the protein, glycosylation (addition of carbohydrate structures), the association or coupling of the α and β subunits of the protein, and/or the metabolism and excretion of the protein from maternal circulation.Reference Hussa 38 It is likely that the types of variability that we report here reflect changes in different combinations of these mechanisms, some of which are sexually dimorphic, under the control of polymorphic genes that vary by race, and also susceptible to factors in the maternal environment.
We show here that the fetal sex differences persist in the clinical MoM after normalization for other factors. This difference has been reported across time and populations, but has never incited clinical efforts to determine and correct for fetal sex in screening protocols largely because of the difficulty in obtaining an accurate fetal sex prediction at these early time points. On the other hand, lack of adjustment for fetal sex in epidemiologic analyses may result in bias or possibly missing a sex-specific association or dose–response relationship.
The gestational age trends in the female/male ratio may tell us something about the mechanism driving the difference. The downward trend in the hCG fetal sex ratio followed by the upward trend in the second trimester suggests that there may be a switch to a different placental secretion pattern or hormonal function in the second trimester. Given that placental growth, invasion and expansion are the dominant forces in the second trimester, and male placentas are generally larger at birth than female placentas, we would predict that male hCG would be higher in the second trimester. However, that is not the case in these data at 20 weeks, and may offer support for the presence of a negative feedback mechanism in the male, or a role for hCG in ‘non-canonical’ placental functions.Reference Clements, Reyes, Winter and Faiman 2
The clinical currency in these data is the MoM. This is the value recorded in prenatal charts and used by the obstetrician to identify and properly manage a high-risk pregnancy. It is useful clinically because it ‘removes’ the subclinical sources of variability from the analyte value and yields a normalized value that can more easily be interpreted at the level of the individual. However, for the purposes of epidemiologic analyses in which we would like to study the effects of maternal race/ethnicity, weight or body mass index, smoking, and health, or other variables correlated with these factors, the clinical MoM is not useful. We propose here to use the ga-MoM, which removes the gestational age effects on the analyte yet, preserves the types of variability that we would otherwise study. In this dataset with extremely stable group mean values, the analyte and the ga-MoM essentially yielded the same results. However, in a smaller dataset the differences between the two would be greater. For the sake of capturing important gestational age variability, the analyte values would be preferable over the MoM. Precaution has to be taken in correctly modeling the gestational age dependence as it is rarely linear and mis-specification can bias estimation of associations.
In conclusion, we demonstrate that fetal sex differences in circulating hCG in early pregnancy were associated with maternal race, weight, age which may be reflective of environmental and physiologic differences in the fetal environment. These types of analyses at the population level potentially allow us to understand non-canonical placental functions, including more direct observations of fetal origins of health and disease. In the case of this analysis, changes in the size of the expected female/male ratio were characterized as potential indicators of normal v. abnormal placental regulation of sex-specific fetal development. We are making a recommendation to researchers in the field of fetal origins to consider using these data to understand placental responses to maternal environmental factors and placental contributions to sex-specific fetal outcomes.
Acknowledgments
The authors acknowledge the mentorship of Ilpo Huhtaniemi and Ulf-Håkan Stenman on hCG and Susan Fisher’s mentorship in placental biology. The authors thank Nerissa Wu for her contributions.
Financial Support
This work has been supported by the National Institute of Environmental Health Sciences (J.J.A., grant numbers 5K99ES017780-02, R00ES017780-06); and the Mt. Zion Health Fund of the Jewish Community Endowment Fund granted through the UCSF National Center of Excellence in Women’s Health (J.J.A.).
Conflicts of Interest.
None.
Ethnical Standards
The authors assert that all procedures contributing to this work comply with the ethnical standards of the Helsinki Declaration of 1975, as revised in 2008, and have been approved by institutional committees (Committee for the Protection of Human Subjects, State of California; University of Pittsburgh Institutional Review Board).