Environmental and genetic contributors to salivary testosterone levels in infants

Transient activation of the hypothalamic–pituitary–gonadal axis in early infancy plays an important role in male genital development and sexual differentiation of the brain, but factors contributing to individual variation in testosterone levels during this period are poorly understood. We measured salivary testosterone levels in 222 infants (119 males, 103 females, 108 singletons, 114 twins) between 2.70 and 4.80 months of age. We tested 16 major demographic and medical history variables for effects on inter-individual variation in salivary testosterone. Using the subset of twins, we estimated genetic and environmental contributions to salivary testosterone levels. Finally, we tested single nucleotide polymorphisms (SNPs) within ±5 kb of genes involved in testosterone synthesis, transport, signaling, and metabolism for associations with salivary testosterone using univariate tests and random forest (RF) analysis. We report an association between 5 min APGAR scores and salivary testosterone levels in males. Twin modeling indicated that individual variability in testosterone levels was primarily explained by environmental factors. Regarding genetic variation, univariate tests did not reveal any variants significantly associated with salivary testosterone after adjusting for false discovery rate. The top hit in males was rs10923844, an SNP of unknown function located downstream of HSD3B1 and HSD3B2. The top hits in females were two SNPs located upstream of ESR1 (rs3407085 and rs2295190). RF analysis, which reflects joint and conditional effects of multiple variants, indicated that genes involved in regulation of reproductive function, particularly LHCGR, are related to salivary testosterone levels in male infants, as are genes involved in cholesterol production, transport, and removal, while genes involved in estrogen signaling are related to salivary testosterone levels in female infants.

Introduction

Transient activation of the hypothalamic–pituitary–gonadal (HPG) axis in the early post-natal period results in elevated levels of gonadotropins and testosterone in human males, a phenomenon known as the “minipuberty” or “neonatal surge” (1–5). The consequences of this activation are not fully understood, but it likely plays an important role in genital development (6–9) and has been linked to future fertility (10). There is also an increasing body of evidence that this is a critical period for the development of sexually dimorphic behavior and psychopathology (11–16). Understanding the causes of individual variation in testosterone levels during the minipuberty is thus of considerable theoretical and clinical interest.

Levels of pubertal and post-pubertal testosterone in males are determined, in part, genetically, with heritability estimates around 40–70% (17–23) for both salivary and plasma measures. Studies in females are less consistent with some reporting minimal heritable variation (23) and others reporting heritabilities similar to those seen in males (17, 20, 22). Recently, genome-wide association studies (GWAS) have been used to identify single nucleotide polymorphisms (SNPs) associated with serum levels of testosterone in adults. Ohlsson et al. (24) identified two loci that met genome-wide significance for serum testosterone levels in men of European ancestry, one in the gene for sex hormone-binding globulin (SHBG, lead SNPs rs12150660 and rs6258) and one near family with sequence similarity 9, member B (FAM9B) on Xp22 (lead SNP rs5934505). Jin et al. (25) confirmed both loci in an independent study, although the lead SNP in the SHBG locus differed (rs727428). They also reported a new hit at rs10822184, which is located in receptor accessory protein 3 (REEP3). Chen et al. (26) reported a significant association between rs2075230 in the SHBG locus and serum testosterone in a sample of Chinese men. A GWAS of 1600 post-menopausal women failed to find any genome-wide significant associations with testosterone (27).

There is evidence to suggest that distinct genetic mechanisms influence testosterone levels across developmental time (17), so it is theoretically possible that the relative importance of genetic factors and the specific variants involved may be different in the minipuberty. Indeed, a study of salivary testosterone in twins between 4 and 8 months of age reported that variation was completely explained by common and unique environmental factors (28). This result is extremely intriguing, but the age range studied is slightly beyond the minipuberty as classically defined. Most studies suggest that the post-natal surge in testosterone peaks in males between 1 and 3 months of age and declines significantly by 6 months of age (4, 5, 29, 30). The gonads are the main source of testosterone in males during the minipuberty, but after 6 months the adrenal becomes the primary source of testosterone (31) as a consequence of the regression and degeneration of neonatal Leydig cells (32). Complementary studies at earlier ages are therefore needed as are studies looking at the impact of specific environmental factors.

The current study takes advantage of an ongoing GWAS of infant brain development being carried out at the University of North Carolina (UNC) at Chapel Hill. The study includes both twins and singletons. A subset of participating children donated saliva samples for the assessment of salivary testosterone around 3 months of age [see Ref. (33)]. We used this sample to estimate genetic and environmental contributions to salivary testosterone levels during the minipuberty. We also examined the impact of 16 major demographic and medical history variables on inter-individual variation in salivary testosterone. In addition, we tested SNPs within ±5 kb of genes involved in testosterone synthesis, transport, signaling, and metabolism for their relationship with salivary testosterone. We also tested SNPs in REEP3 and the Xp22 loci identified in GWAS studies of serum testosterone in adult males to see if these were predictive in infants. Finally, to address the issue of potentially small effect sizes for individual SNPs, as well as non-additive effects, we used random forest (RF) methods to identify combinations of SNPs contributing to variation in salivary testosterone in infancy.

Materials and Methods

Subjects

Two-hundred twenty-two infants (119 males, 103 females) between 2.70 and 4.80 months of age (mean age 3.34 months post-date of birth ± 0.36) with high quality genetic information and saliva samples suitable for hormone assay are included in this analysis. The sample includes 108 singletons and 114 twins (27 same-sex DZ pairs, 27 same-sex MZ pairs, 2 opposite-sex DZ pairs, and 2 unpaired twins). All children were participating in prospective longitudinal studies of early brain development for which the senior author is a co-investigator. Mothers were recruited during the second trimester of pregnancy from the outpatient obstetrics and gynecology clinics at UNC hospitals. Exclusion criteria at enrollment were the presence of abnormalities on fetal ultrasound or major medical illness in the mother. Demographic variables (maternal age, paternal age, maternal education, paternal education, maternal ethnicity, paternal ethnicity, maternal psychiatric history, paternal psychiatric history, and total household income) were collected via maternal report at the time of enrollment (see Supplementary Material). For the purpose of the current study, maternal psychiatric history and paternal psychiatric history were treated as binary variables. Individuals were counted as positive for psychiatric history if they had received any psychiatric diagnosis. Medical history variables (birth weight, gestational age at birth, 5 min APGAR scores, stay in neonatal intensive care unit over 24 h, gestation number, and delivery method) were collected from maternity and pediatric medical records shortly after birth (see Supplementary Material). Maternal smoking during pregnancy was collected via maternal report at two timepoints during pregnancy and shortly after birth. Demographic and medical history data are summarized in Table 1. Experiments were undertaken with the understanding and written consent of each subject’s mother or father, with the approval of the Institutional Review Board of the UNC School of Medicine.

TABLE 1

Table 1. Demographic and medical history information.

Salivary Assays

Saliva samples were collected during visits to participants’ homes. Visits were scheduled for 9:00 a.m. Mean ± SD for collection start time was 9:14 a.m. ± 16 min. Parents were advised not to feed their children for at least 15 min prior to collection. One milliliter of passive drool was collected from each participant using a suction catheter (Centurion Healthcare Products, Howell, MO, USA). Mean ± SD for duration of collection was 4 ± 2 min. All samples were frozen within 5 h and stored in a −80 °C freezer. Mean ± SD for time until storage was 60 ± 31 min.

Salivary testosterone levels were measured by enzyme immunoassay using a commercially available kit (Salimetrics, State College, PA, USA). The intra-assay precision for samples with low testosterone levels (mean 18.12 pg/ml) is 6.7%; for high testosterone levels (mean 188.83 pg/ml) it is 2.5%. Inter-assay precision for samples with low testosterone levels (mean 19.6 pg/ml) is 14.05%; for high testosterone levels (mean 199.08 pg/ml) it is 5.6%. Percent recovery for this assay varies from 92 to 111.4%. The minimal concentration of testosterone that can be distinguished from 0 is <1.0 pg/ml. Only one sample had a concentration below the detection limit. This sample was coded as 0 pg/ml in subsequent analyses. We also evaluated all samples for blood contamination using the Salimetrics Salivary Blood Contamination Enzyme Immunoassay kit, which quantitatively measures transferrin, a large protein, which is present in abundance in blood, but that is normally present in only trace amounts in saliva. Intra-assay precision for samples with high (3.88 mg/dl) transferrin levels is 10.2%, for samples with low (0.42 mg/dl) transferrin levels it is 4.9%. Inter-assay precision is 7.1% for low (1.02 mg/dl) and 7.2% for high (4.93 mg/dl) transferrin levels. Percent recovery varies from 91.9 to 101.5%. The minimal concentration of transferrin that can be distinguished from zero is 0.08 mg/dl.

Genotyping

DNA was extracted from buccal cells using standard methods as described in the Puregene^® DNA Purification Kit (Gentra Systems) or using phenol/chloroform. After extraction, samples were stored in a −80 °C freezer until analysis. Genotyping was carried out at the Bionomics Research and Technology Center at Rutgers (Piscataway, NJ, USA) using Affymetrix Axiom Genome-Wide LAT and Exome arrays. Samples were randomized across 96-well plates. Each plate contained a common control sample. Genotype calling was performed with the Affymetrix Genotyping Console. Rigorous quality control procedures were carried out prior to analysis. In brief, we excluded samples with low DishQC (<0.82 for LAT array and <0.79 for Exome array), low call rates (<95%), outliers for homozygosity, sex, or zygosity from genotypes inconsistent with reported phenotypes, ancestry outliers, excessive relatedness, and unexpected relatedness. We also removed individual SNPs that deviated from Hardy–Weinberg equilibrium (P_HWE < 1 × 10⁻⁸), had low call rate (<95%), high deviation of allele frequency from 1000G EUR/AFR founders, and that did not match 1000G EUR/AFR founders. Population stratification was assessed using principal component analysis (PCA) (34, 35). Population stratification is the presence of systematic differences in allele frequencies between subpopulations with different genetic ancestry and can lead to spurious results in genetic association studies. We first pruned the genotyped SNPs using the pruning tool in PLINK (36). We then created a subset of our total GWAS sample without related subjects by randomly selecting one individual from twin and sibling pairs (all non-twins/non-sibs were also included). Smartpca in EIGENSOFT was applied to compute the principal components (PCs or eigenvectors). The first three PCs explain 6.8% of the variability in genetic variables. See Figures 1 and 2 for a scree plot and 3D scatter plots of the first three PCs color-coded by maternal and paternal reported ethnicity. Strong clustering of subjects with similar reported ethnicity confirms that the first three PCs index genetic ancestry. Imputation was then performed with MaCH-Admix using 1000G reference panel (phase1_release_v3.20101123) followed by post-imputation quality control based on imputation quality score.

FIGURE 1

Figure 1. Scree plot demonstrating that the majority of genetic variation is explained by the first three PCs.

FIGURE 2

Figure 2. 3D scatter plots of the first three PCs color-coded by maternal and paternal reported ethnicity.

For the current analysis, we extracted all genotyped SNPs within ±5 kb of genes involved in testosterone synthesis, transport, signaling, and metabolism (see Table 2). We also included genotyped SNPs in/near REEP3 and the Xp22 loci flagged in GWAS studies of serum testosterone in adult males. Three of the top GWAS hits in SHBG were not genotyped by either array, but were imputed with good quality (rs727428, rs12150660, and rs2075230). The imputed SNPs were used in the current analysis. rs6258 was directly genotyped, but was excluded due to low minor-allele frequency (MAF). Top hits in REEP3 and the Xp22 loci were also imputed (rs10822184 and rs5934505). Additional imputed SNPs near the Xp22 loci were also included due to the low number of genotyped SNPs in this region. SNPs with MAF <0.05 were not included in our analyses. In total, 512 SNPs were examined in males and 473 SNPs were examined in females.

TABLE 2

Table 2. Genes probed in the current analysis.

Environmental Association Analysis

In order to determine the potential impact of major demographic and medical history variables on inter-individual variation in salivary testosterone, we used a moment-based method to select fixed effects in linear mixed effects models (37). Twins are treated as repeated measures. For fixed effects selection, we applied an adaptive Lasso penalty using the feasible generalized least squares estimator as an initial. In the model, we always include log(10) transferrin and age since DOB as predictors. We focused on selecting a list of predictors including maternal age, paternal age, maternal education, paternal education, maternal ethnicity, paternal ethnicity, maternal psychiatric history, paternal psychiatric history, total household income, maternal smoking during pregnancy, birth weight, gestational age at birth, 5 min APGAR scores, stay in neonatal intensive care unit over 24 h, gestation number, and delivery method. We used the BIC statistic to select the tuning parameter of the adaptive Lasso as in Section 4 of Ahn et al. (37). Before applying our variable selection method, we standardized all covariates and centered the response variable (testosterone). We also applied bootstrap methods 1000 times to assess the stability of our results.

After model selection, we ran a mixed effect model using the selected variables for significance testing and to estimate r² values. Mixed effect models were also run including all variables for comparison. Males and females were analyzed separately for this and all subsequent analyses. We separated the sexes because the source of testosterone is primarily gonadal in males and adrenal in females during the minipuberty (31). Studies in adults also suggest that different genetic factors influence testosterone concentrations in men and women (17). Variables were considered significant if they were selected in the Lasso model and survived Bonferroni correction in the mixed effect model including all variables.

Intraclass Correlation Coefficient and Heritability Estimation

The intraclass correlation coefficient (ICCs) and their confidence intervals were estimated separately for MZM (monozygotic male), DZM (dizygotic male), MZF (monozygotic female), and DZF (dizygotic female) using R package called ICC. The heritability was estimated using linear mixed effect model (ACE model) while restricting variance of genetic, shared environment, and random error to be >0. One-tailed t-tests were used to test whether the estimates were significantly >0 (p-values <0.05 were considered significant). For males, log(10)transferrin, age since DOB, and 5 min APGAR score were included as covariates. For females, log(10)transferrin and age since DOB were included as covariates.

Univariate Genetic Association Analysis

Males and females were analyzed separately. For males, log(10)transferrin, age since DOB, and 5 min APGAR score were included as covariates. For females, log(10)transferrin and age since DOB were included as covariates. In addition, the first three principle components derived from all genotyped SNPs were included as covariates to control for possible population stratification. Association analysis was performed using mixed effect models with likelihood ratio tests. p-values were adjusted by false discovery rate (FDR); p-values <0.05 after FDR correction were considered significant.

Univariate genetic association analysis provides a straightforward approach to identifying genetic variants associated with phenotypes of interest, in this case salivary testosterone levels. However, given the large number of tests performed, this approach is underpowered for identifying SNPs with small effect sizes. In addition, it does not account for genetic interactions, i.e., causal effects that are only observed when specific combinations of mutations and/or non-mutations are jointly present. In order to address these limitations, we also carried out RF analysis to identify combinations of SNPs contributing to variation in salivary testosterone in infancy.

To gain a better understanding of the power of our univariate genetic association analysis, we generated power plots using different SNP MAFs and effect sizes (beta).

We used the same number of subjects (and proportion of MZs, DZs, and singletons) as our real data and similar scale of effect size, mean, variance, and additive genetic and shared environment effects estimated from the real data. Ten thousand simulations were done for each scenario and the significance level was set to 0.05/500, which is similar to our data after Bonferroni correction.

Random Forests

Random forest analysis was implemented using R software. The result of RF analysis reflects joint and conditional effects of multiple variables. As variables with strong effects may mask weaker, yet important effects, we removed the effects of log(10)transferrin, age since DOB, 5 min APGAR score, and the first three principle components (for males) and log(10)transferrin, age since DOB, and the first three principle components (for females). Specifically, we (1) used the residual from the regression Testosterone ~log(10)transferrin + age since DOB + 5 min APGAR score (males only) + PC1 + PC2 + PC3 as the new response, and (2) used the residual from the regression SNPs ~log(10)transferrin + age since DOB + 5 min APGAR score (males only) + PC1 + PC2 + PC3 as the new predictor variables for the RF analysis. The number of randomly preselected variables (m_try) was set to the number of SNPs divided by 3 (p/3) as recommended by Liaw and Wiener (38). We also examined the error rate for m_try = 2p/3, 0.5p/3, and p. Error rate was highly similar across different m_try, but p/3 produced the lowest error rate. The number of trees (ntree) was set to 2000. RF analysis was run 50 times, each run using different seeds for random number generation (RNG). In each run, we recorded the 30 SNPs with the highest variable importance. After all 50 runs were completed, we identified those SNPs, which appeared in the “top 30 list” in every 1 of 50 runs. Because RF analysis is not set-up to treat twins as a repeated measure, we randomly picked one twin from each twin pair to ensure that no related subjects were included in the RF analysis.

Results

Mean ± SD concentration of salivary testosterone was 40.39 ± 13.39 pg/ml in males and 39.70 ± 16.64 pg/ml in females and did not differ significantly (p = 0.73). This is likely a consequence of differing correlations between salivary testosterone and serum testosterone in males and females. According to Salimetrics, the relationship between serum and saliva for males as determined by linear regression is y (total serum testosterone in nanogram per milliliter) = 0.2421 + 0.0496 × x (salivary testosterone in picogram per milliliter). The linear regression equation for females is y (total serum testosterone in nanogram per milliliter) = 0.1415 + 0.0055 × x (salivary testosterone in picogram per milliliter). Assuming that the relationship between total serum testosterone and salivary testosterone in infants is similar to that seen in adults; we would estimate a mean serum level of 2.25 ng/ml for male infants and 0.36 ng/ml for female infants, which is comparable to published reports on serum testosterone levels in this age range (39, 40).

Mean ± SD for transferrin was 0.844 ± 0.93 mg/dl in males and 0.91 ± 0.98 mg/dl in females. A moderate correlation between transferrin and testosterone was observed in both males and females (r = 0.28, p = 0.002 and r = 0.43, p < 0.001, respectively). Transferrin levels were higher than expected, which raises the possibility of blood contamination. We addressed this issue in two ways. (1) In all our primary analyses, we included transferrin as a covariate. Specifically, we adjusted for log(10)transferrin as the untransformed variable showed high levels of skewness and kurtosis. (2) For the environmental association analysis, univariate genetic analyses, and RF analysis, we also performed sensitivity analyses in infants with transferrin levels <0.50 mg/dl [cut-off based on the recommendation of Granger et al. (41)]. For ICC and heritability estimation, we did not perform sensitivity analyses due to insufficient sample size.

Environmental Association Analysis

The final model for males using adaptive Lasso included 5 min APGAR score in addition to the fixed variables, log(10)transferrin, and age since DOB. Bootstrapping supported the importance of 5 min APGAR score in that it was selected in 837 out of 1000 tests (see Table 3). A linear mixed effect model including these variables explained approximately 32% of the variance in salivary testosterone. Examination of a linear mixed effect model including all predictors confirmed that 5 min APGAR score was the only major demographic or medical history variable examined which significantly predicted salivary testosterone in males after Bonferroni correction (0.05/29 predictors = 0.0017) and explained approximately 15% of the variance. In males with transferrin levels <0.50 mg/dl, a linear mixed model including age since DOB and 5 min APGAR score explained 22% of the variance in salivary testosterone (see Table 4). The final model for females using adaptive Lasso only included the fixed variables, log(10)transferrin, and age since DOB. A mixed effect model including these variables explained approximately 18% of the variance in salivary testosterone. Examination of a mixed effect model including all predictors confirmed that none of the demographic or medical history variables examined significantly predicted salivary testosterone in females after Bonferroni correction (0.05/29 predictors = 0.0017). In females with transferrin levels <0.50 mg/dl, a mixed model including age since DOB explained 1% of the variance in salivary testosterone (see Table 5). Given the converging evidence for the importance of 5 min APGAR score on salivary testosterone in males, it was included as a covariate in subsequent analyses.

TABLE 3

Table 3. Bootstrapping results.

TABLE 4

Table 4. Association of salivary testosterone with demographic and medical history variables in males.

TABLE 5

Table 5. Association of salivary testosterone with demographic and medical history variables in females.

Intraclass Correlation Coefficient and Heritability Estimation

The intraclass correlations for MZ and DZ twins are shown in Figure 3 and suggest a high environmental component and a low genetic component for both sexes. The ACE model confirmed that the majority of variation in salivary testosterone was explained by shared environmental factors in both sexes (see Table 6).

FIGURE 3

Figure 3. Intraclass correlations within twin pairs are shown. MZM (monozygotic male), DZM (dizygotic male), MZF (monozygotic female), and DZF (dizygotic female).

TABLE 6

Table 6. Univariate genetic model.

Univariate Genetic Association Analysis

In males, no SNPs were significantly associated with salivary testosterone levels after adjusting for FDR in the full sample. The most significant association was for rs10923844, a variant of unknown function located downstream of HSD3B1 and HSD3B2 (adjusted p-value = 0.07, unadjusted p = 0.0002). This same SNP had an unadjusted p-value of 0.02 in the subsample of males with transferrin <0.50 mg/dl. Effect sizes for this variant were highly similar in the total male sample [adjusting for log(10)transferrin] and the subsample of males with transferrin <0.50 mg/dl (Beta = 6.26 and 6.02, respectively). Indeed, effect sizes for the full male sample are highly similar to the corresponding effect sizes in the subset of males with transferrin <0.50 mg/dl (see Figure 4). None of the previous GWAS hits for serum testosterone levels in adult males were significantly associated with salivary testosterone in infant males before or after adjusting for FDR in either the full sample or the subset of males with transferrin <0.50 mg/dl. In females, no SNPs were significantly associated with salivary testosterone levels after adjusting for FDR. The most significant associations were for two SNPs located upstream of ESR1 (rs3407085 and rs2295190) within the intron region of SYNE1 (spectrin repeat containing, nuclear envelope (1); adjusted and unadjusted p-values were 0.09 and 0.0004 for both variants (these variants are in complete linkage disequilibrium within our sample). Both variants had an unadjusted p-value of 0.04 in the subsample of females with transferrin <0.50 mg/dl. Effect sizes were highly similar for the total female sample [adjusting for log(10)transferrin] and the subsample of females with transferrin <0.50 mg/dl (Beta = −13.7 and −11.2, respectively), though we note that effect sizes for the full female sample do not correspond as well to the effect sizes in the reduced sample as compared to males (See Figure 4). Full results are available as supplemental data. Power calculations indicate that theoretically, our analysis is well-powered to identify variants with Beta values >6 in the MAF range examined (see Figure 5).

FIGURE 4

Figure 4. Correlation between effect sizes in full samples [correcting for log(10)transferrin] and subsamples with transferrin levels <0.50 mg/dl.

FIGURE 5

Figure 5. Plots demonstrating the relationship between statistical power, minor-allele frequencies, and effect size are shown.

Random Forests

Tables 7 and 8 display the top SNPs in males and females, respectively, defined as the intersection of the top 30 variables obtained from each of 50 RF runs. Also included in the table are the beta.hat, likelihood ratio, and unadjusted and FDR adjusted p-values from the univariate association tests. In males, regulators of reproductive function appear to play an important role. Specifically, multiple variants in or near LHCGR were consistently selected as high importance. Multiple SNPs in genes related to cholesterol production, transport, and removal were also selected. These genes included DHCR7, TSPO, CYP7A1, and POR. The top hit from the univariate analyses (rs10923844) was selected, as were two variants near the Xp22 loci identified in GWAS studies of serum testosterone in adult males. Analyses in the subset of males with transferrin <0.50 mg/dl supported the importance of regulators of reproductive function and cholesterol-related genes. Five SNPs were selected as high importance in both the full and reduced male sample; these were rs4952922, rs10495960, and rs2301267 (all in/near LHCGR), rs47340 (TSPO), and rs74091680 (HSD17B6). In females, genes related to estrogen signaling appear to play an important role. Specifically, variants in ESR1, CYP19A1, CYP1B1, and HSD17B8 were consistently selected as high importance in the full sample. Analyses in the subset of females with transferrin <0.50 mg/dl supported the importance of estrogen signaling, although the specific SNPs differed. Two SNPs were selected as high importance in both the full and reduced female sample; these were rs4149448 (SULT2A1) and rs4659175 (HSD3B2).

TABLE 7

Table 7. Random forest results males.

TABLE 8

Table 8. Random forest results females.

Discussion

This study provides the first detailed analysis of environmental and genetic contributors to variation in salivary testosterone during the minipuberty, a developmental period that plays a critical role in genital development and sexual differentiation of the brain. As such, it should be considered exploratory and requires replication. Twin modeling indicated that individual variability in testosterone levels in the minipuberty is primarily explained by environmental factors. In terms of specific environmental contributors, we observed a consistent and robust association between 5 min APGAR scores and salivary testosterone levels in males. In terms of specific genetic contributors, univariate tests did not reveal any variants significantly associated with salivary testosterone after adjusting for FDR. However, we note that this approach is underpowered for identifying SNPs with small effect sizes. The top hit in males was rs10923844, an SNP of unknown function located downstream of HSD3B1 and HSD3B2. The top hits in females were two SNPs located upstream of ESR1 (rs3407085 and rs2295190). RF analysis, which reflects joint and conditional effects of multiple variables, including those with small individual effect sizes, suggests that genes involved in regulation of reproductive function and cholesterol production, transport, and removal are involved in individual variation in salivary testosterone in males, while genes involved in estrogen signaling are important in females.

Major Demographic and Medical History Variables and Salivary Testosterone

Higher APGAR scores taken 5 min post-birth were associated with lower salivary testosterone levels in males during the minipuberty. APGAR scores were originally designed to quickly evaluate a newborn’s physical condition and to identify any immediate need for extra medical or emergency care. The APGAR score includes five components: heart rate, respiratory effort, muscle tone, reflex irritability, and color, each of which is given a score of 0, 1, or 2 with 2 being the best score (42). A number of factors may influence an APGAR score, including hypoxia, exposure to drugs, trauma, congenital anomalies, infections, and hypovolemia (43). While primarily conceived as a measure of short peripartum stress, some researchers have suggested that it also indexes a suboptimal fetal environment (44). Studies in rodent models provide compelling evidence that the stress response induced by physical or emotional challenges in fetal life affects later reproductive function (45–47). Regarding the HPG axis, maternal stress during pregnancy disrupts the prenatal surge of testosterone that normally occurs in the developing male rat (48, 49) and is associated with reduced testosterone levels in adult rats (50). No studies have examined stress response during labor and delivery and its impact on reproductive function. Physical stressors associated with lower APGAR scores may explain the association between 5 min APGAR scores and salivary testosterone in males observed in our study, although the direction of effect is opposite to that reported for prenatal stress in rodents. We note that while 5-min Apgar score is a valid predictor of neonatal mortality, its predictive value for other outcomes continues to be debated (43). Never-the-less, this scoring system remains the only widely used and accepted tool for assessing the vitality of newborn infants across the world. We also note that the majority of children in our sample had 5-min Apgar scores between 7 and 10, which is considered normal, with only a few falling in the intermediate range (scores between 4 and 6). Additional research is necessary to replicate these findings, determine if they extend to individuals with lower APGAR scores, and understand underlying mechanisms. None of the major demographic or maternal history variables examined in this study impacted salivary testosterone levels in females.

Heritability of Salivary Testosterone

Our results suggest that individual variability in salivary testosterone levels in the minipuberty is predominantly explained by environmental factors. Our findings are similar to those reported by Caramaschi et al. (28) in older infants, but contrast with studies in adolescents and adults, which show high heritability, especially in males (17–23). Caramaschi et al. (28) suggested that the pubertal rise in male testosterone levels has a strong genetic component and that testosterone levels are less influenced by genetic factors when male and female levels are similar and very low as is the case in later infancy. Salivary testosterone levels in our study were significantly higher than those reported by Caramaschi et al. (28): 40.39 vs. 9.60 pg/ml in boys and 39.70 vs. 8.93 pg/ml in girls, but are lower than those reported for adults: 163.81 and 60.86 pg/ml for males and females, respectively (51). Salivary testosterone did not differ significantly between males and females in our sample, but estimated serum levels exhibited a sex difference, which is comparable to published reports on serum testosterone levels in this age range (39, 40). Therefore, our results suggest that environmental factors play a stronger role than genetic factors during the minipuberty despite moderately high testosterone levels and sex differences in serum testosterone. Identifying the environmental factors involved is an important area for future research.

Genetic Factors Influencing Salivary Testosterone

Although our heritability analysis indicated that environmental risk factors play a larger role in determining individual variation in salivary testosterone during the minipuberty than genetic factors, deciphering the genetic component is still important in terms of understanding underlying biochemical pathways. We were also mindful of the fact that heritability by itself does not provide information about the genetic architecture of traits. In principle, a trait with a low heritability can have a single locus that causes variation and a trait with a high heritability can have hundreds of contributing loci (52). The presence of gene–environment correlations and interactions also introduces systematic biases in heritability estimates made under the independence assumption in twin studies, which can result in an underestimation of the genetic component (53) For all these reasons, we performed univariate genetic association analyses and RF analyses of SNPs within ±5 kb of genes involved in testosterone synthesis, transport, signaling, and metabolism for their relationship with salivary testosterone.

No SNPs were significantly associated with salivary testosterone levels after adjusting for FDR in the full male sample or the full female sample. Simulations indicated that theoretically we were well-powered to detect significant SNPs with Beta values >6 within the range of minor-allele frequencies we examined. Thus, our study suggests that common variants with large effect sizes do not play a role in individual differences in salivary testosterone in infants, at least within the genes examined. The top hit in males was rs10923844, a variant of unknown function located downstream of HSD3B1 and HSD3B2, two genes that play a critical role in steroidogenesis. This variant was also identified as an important variable in the RF analysis. Both isoforms of 3-beta-HSD convert pregnenolone to progesterone, 17α-hydroxypregnenolone to 17α-hydroxyprogesterone, dehydroepiandrosterone (DHEA) to androstenedione, and androstenediol to testosterone. Type 1 is primarily expressed in the placenta and non-steroidogenic tissues, while type 2 is primarily expressed in the adrenals and gonads (54). The top hits in females were two SNPs located upstream of ESR1 (rs3407085 and rs2295190), the gene coding for estrogen receptor alpha, one of two main types of estrogen receptor.

It is increasingly recognized that individual variation in complex phenotypes results, in part, from the joint and conditional effects of many common SNPs whose individual effect sizes are relatively small, making them difficult to identify via standard univariate analyses (55, 56). RF analysis represents a promising approach to this problem and is well-suited to high-dimension, low-sample-size data typical of genetic association studies. RF analysis suggested that regulators of reproductive function and genes related to cholesterol production, transport, and removal play an important role in salivary testosterone levels in males. Five SNPs were selected as high importance in both the full and reduced male sample and thus are particularly promising candidates for follow-up. rs4952922, rs10495960, and rs2301267 are all located in/near LHCGR, which codes for the LH/choriogonadotropin receptor. rs47340 is located downstream of TSPO, which codes for a protein that interacts with StAR (steroidogenic acute regulatory protein) to transport cholesterol into mitochondria to permit steroid synthesis [data on gene functions from www.genecards.org (57), SNP locations from dbSNP (58) and UCSC genome browser (59)]. rs74091680 is an intronic variant in HSD17B6, which codes for an enzyme involved in androgen catabolism. Specifically, the oxidoreductase activity can convert 3 alpha-adiol to dihydrotestosterone, while the epimerase activity can convert androsterone to epi-androsterone. Interestingly, two variants near the Xp22 loci identified in GWAS studies of serum testosterone in adult males (24, 25) were selected in the primary RF analysis. This suggests that some genetic factors influencing testosterone levels are active throughout the lifespan, although this finding was not recapitulated in the subsample analysis.

In females, RF analysis suggested that genes related to estrogen signaling play an important role. Specifically, variants in/near ESR1, CYP19A1, CYP1B1, and HSD17B8 were consistently selected as high importance in the full sample. ESR1 codes for estrogen receptor alpha, one of two main types of estrogen receptor. CYP19A1 codes for aromatase, the enzyme responsible for the conversion of androgens into estrogens. CYP1B1 codes for an enzyme that metabolizes multiple compounds, including 17β-estradiol. HSD17B8 inactivates estrogens and androgens, with high activity toward estrogen and low activity toward testosterone. It can also synthesize estradiol from estrone. Analyses in the subset of females with transferrin <0.50 mg/dl supported the importance of estrogen signaling, although the specific SNPs implicated differed. Two SNPs were selected as high importance in both the full and reduced female sample; these were rs4149448 (SULT2A1) and rs4659175 (HSD3B2). SULT2A1 codes for a protein, which catalyzes the sulfation of steroids and bile acids in the liver and adrenal glands, and may have a role in the inherited adrenal androgen excess in women with polycystic ovary syndrome (60). HSD3B2 codes for an enzyme with a critical role in the biosynthesis of all classes of steroid hormones.

Limitations

The primary limitation of the current study is the use of salivary testosterone levels rather than plasma testosterone levels. While plasma and serum are the traditional matrices for the determination of endocrine parameters including testosterone, saliva offers a non-invasive and stress-free alternative, which enjoys widespread acceptance and has been used for over 40 years (61). In men, salivary testosterone is strongly correlated with the free fraction of testosterone in serum (r = 0.64–0.97) (51, 62–64). In women, reported correlations are more moderate (r = 0.37–0.85) (51, 63, 65). There is no a priori reason to presume that correlations between serum and saliva are different in the neonate, although this has not been directly tested. Collection of saliva maximized parental acceptability and minimized child distress. It is unlikely that we would have achieved a reasonable sample size if we opted to use blood. Substitution of saliva assay results for serum values is known to underestimate testosterone–behavior associations, primarily in females, and this problem may also be applicable to the environmental and genetic associations, which were the focus of this study. Null results, especially in females, should be treated with caution.

We also observed higher levels of transferrin in saliva than expected and moderate correlations between transferrin and testosterone in both sexes (stronger in females), which raises the possibility of blood contamination. For comparison, Granger et al. (41), in a study of children between 6 and 13 years of age, reported a mean ± SD for transferrin of 0.37 ± 0.46 mg/dl and a correlation with testosterone of 0.058. Unfortunately, although the transferrin assay (66) has been available for over a decade, it has not been routinely used in studies of infants making it difficult to determine if the levels that we observed are typical for this age or not. Reassuringly, the mean ± SD for salivary testosterone in our study was highly similar to that reported by other studies in our age range (11) and estimated serum levels are comparable to published reports on serum testosterone in this age range (39, 40). We also addressed this issue statistically. In all our primary analyses, we included transferrin as a covariate. In addition, for the environmental association analysis, univariate genetic analyses, and RF analysis, we performed sensitivity analyses in infants with transferrin levels <0.50 mg/dl (cut-off based on the recommendation of Granger et al. (41). In general, relationships identified in the full sample were also present in the subsample with transferrin levels <0.50 mg/dl. Finally, we note that collecting samples on multiple days would have allowed a more precise evaluation of individual differences in testosterone levels, but was deemed impractical. The sample size for twin modeling was also somewhat underpowered.

Conclusion

Despite these limitations, the current study provides novel information about the environmental and genetic contributors to testosterone levels in the early post-natal period. In contrast to the strong genetic contributions observed during puberty, transient activation of the HPG axis in the early post-natal period appears to be heavily influenced by environmental factors. We identified 5 min APGAR score as a significant predictor of salivary testosterone levels in males. Further research is needed to elucidate the biological mechanisms underlying this relationship. Our study also suggests that genetic variants in regulators of reproductive function and cholesterol play an important role in salivary testosterone levels in males, while genes related to estrogen signaling play an important role in females. All results require replication, preferably with serially collected serum samples. Ideally, such a study would also include measures of prenatal testosterone in amniotic fluid or through cordocentesis, providing insight into the genetic architecture of the prenatal testosterone surge as well as the neonatal surge.

Author Contributions

Kai Xia, Yang Yu, Mihye Ahn, Hongtu Zhu, and Fei Zou made substantial contributions to the analysis and interpretation of data for this project. John H. Gilmore made substantial contributions to the acquisition of data for this project. Rebecca C. Knickmeyer made substantial contributions to the conception and design of the work; as well as the acquisition, analysis, and interpretation of data. Kai Xia, Yang Yu, Mihye Ahn, and Rebecca C. Knickmeyer drafted the manuscript. All coauthors revised the manuscript for important intellectual content, and approved the final version to be published. RK agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Institutes of Health (MH064065 and MH070890 to John H. Gilmore, MH083045 and MH092335 to Rebecca C. Knickmeyer, and RR025747, P01CA142538, MH086633, EB005149, and AG033387 to Hongtu Zhu) and by Autism Speaks. We would like to thank the participating families who made this project possible as well as the staff of the UNC Early Brain Development Program and David Barrow, Director of the UNC Cytokine Analysis Facility where enzyme immunoassays were performed.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fendo.2014.00187/abstract

References

1. Forest MG, Cathiard AM. Pattern of plasma testosterone and ⁴-androstenedione in normal newborns: evidence for testicular activity at birth. J Clin Endocrinol Metab (1975) 41:977–80. doi:10.1210/jcem-41-5-977

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

2. Winter JS, Faiman C, Hobson WC, Prasad AV, Reyes FI. Pituitary-gonadal relations in infancy. I. Patterns of serum gonadotropin concentrations from birth to four years of age in man and chimpanzee. J Clin Endocrinol Metab (1975) 40(4):545–51. doi:10.1210/jcem-40-4-545

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

3. Forest MG, de Peretti E, Bertrand J. Testicular and adrenal androgens and their binding to plasma proteins in the perinatal period: developmental patterns of plasma testosterone, 4-androstenedione, dehydroepiandrosterone and its sulfate in premature and small for date infants as compared with that of full-term infants. J Steroid Biochem (1980) 12:25–36.

Google Scholar

4. Kulle AE, Riepe FG, Melchior D, Hiort O, Holterhus PM. A novel ultrapressure liquid chromatography tandem mass spectrometry method for the simultaneous determination of androstenedione, testosterone, and dihydrotestosterone in pediatric blood samples: age- and sex-specific reference data. J Clin Endocrinol Metab (2010) 95(5):2399–409. doi:10.1210/jc.2009-1670

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

5. Kuiri-Hanninen T, Seuri R, Tyrvainen E, Turpeinen U, Hamalainen E, Stenman UH, et al. Increased activity of the hypothalamic-pituitary-testicular axis in infancy results in increased androgen action in premature boys. J Clin Endocrinol Metab (2011) 96(1):98–105. doi:10.1210/jc.2010-1359

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

6. Mann DR, Gould KG, Collins DC. Blockade of neonatal activation of the pituitary-testicular axis: effect on peripubertal luteinizing hormone and testosterone secretion and on testicular development in male monkeys. J Clin Endocrinol Metab (1989) 68:600–7. doi:10.1210/jcem-68-3-600

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

7. Brown GR, Dixson AF. Investigation of the role of postnatal testosterone in the expression of sex differences in behavior in infnats rhesus macaques (Macaca mulatta). Horm Behav (1999) 35:186–94. doi:10.1006/hbeh.1999.1512

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

8. Main KM, Schmidt IM, Skakkebaek NE. A possible role for reproductive hormones in newborn boys: progressive hypogonadism without the postnatal testosterone peak. J Clin Endocrinol Metab (2000) 85(12):4905–7. doi:10.1210/jcem.85.12.7058

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

9. Sharpe RM, Fraser HM, Brougham MF, McKinnell C, Morris KD, Kelnar CJ, et al. Role of the neonatal period of pituitary-testicular activity in germ cell proliferation and differentiation in the primate testis. Hum Reprod (2003) 18(10):2110–7. doi:10.1093/humrep/deg413

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

10. Zivkovic D, Fratric I. Disturbances of sperm maturation and minipuberty: is there a connection? Biomed Res Int (2014) 2014:912746. doi:10.1155/2014/912746

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

11. Alexander GM, Wilcox T, Farmer ME. Hormone-behavior associations in early infancy. Horm Behav (2009) 56(5):498–502. doi:10.1016/j.yhbeh.2009.08.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

12. Alexander GM, Saenz J. Postnatal testosterone levels and temperament in early infancy. Arch Sex Behav (2011) 40(6):1287–92. doi:10.1007/s10508-010-9701-5

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

13. Lamminmaki A, Hines M, Kuiri-Hanninen T, Kilpelainen L, Dunkel L, Sankilampi U. Testosterone measured in infancy predicts subsequent sex-typed behavior in boys and in girls. Horm Behav (2012) 61(4):611–6. doi:10.1016/j.yhbeh.2012.02.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

14. Saenz J, Alexander GM. Digit ratios (2D:4D), postnatal testosterone and eye contact in toddlers. Biol Psychol (2013) 94(1):106–8. doi:10.1016/j.biopsycho.2013.05.010

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

15. Saenz J, Alexander GM. Postnatal testosterone levels and disorder relevant behavior in the second year of life. Biol Psychol (2013) 94(1):152–9. doi:10.1016/j.biopsycho.2013.05.011

CrossRef Full Text | Google Scholar

16. Alexander GM. Postnatal testosterone concentrations and male social development. Front Endocrinol (2014) 5:15. doi:10.3389/fendo.2014.00015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

17. Harris JA, Vernon PA, Boomsma DI. The heritability of testosterone: a study of Dutch adolescent twins and their parents. Behav Genet (1998) 28:165–71. doi:10.1023/A:1021466929053

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

18. Sluyter F, Keijser JM, Boomsma DI, van Doornen LJ, van den Oord EJCG, Snieder H. Genetics of testosterone and the aggression-hostility-anger (AHA) syndrome: a study of middle aged twins. Twin Res (2000) 3:266–76. doi:10.1375/twin.3.4.266

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

19. Ring HZ, Lessov CN, Reed T, Marcus R, Holloway L, Swan GE, et al. Heritability of plasma sex hormones and hormone binding globulin in adult male twins. J Clin Endocrinol Metab (2005) 90(6):3653–8. doi:10.1210/jc.2004-1025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

20. Hoekstra RA, Bartels M, Boomsma DI. Heritability of testosterone levels in 12-year-old twins and its relation to pubertal development. Twin Res Hum Genet (2006) 9(4):558–65. doi:10.1375/twin.9.4.558

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

21. Kuijper EA, Lambalk CB, Boomsma DI, van der Sluis S, Blankenstein MA, de Geus EJ, et al. Heritability of reproductive hormones in adult male twins. Hum Reprod (2007) 22(8):2153–9. doi:10.1093/humrep/dem145

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

22. Koenis MM, Brouwer RM, van Baal GC, van Soelen IL, Peper JS, van Leeuwen M, et al. Longitudinal study of hormonal and physical development in young twins. J Clin Endocrinol Metab (2013) 98(3):E518–27. doi:10.1210/jc.2012-3361

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

23. Harden KP, Kretsch N, Tackett JL, Tucker-Drob EM. Genetic and environmental influences on testosterone in adolescents: evidence for sex differences. Dev Psychobiol (2014) 56(6):1278–89. doi:10.1002/dev.21207

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

24. Ohlsson C, Wallaschofski H, Lunetta KL, Stolk L, Perry JR, Koster A, et al. Genetic determinants of serum testosterone concentrations in men. PLoS Genet (2011) 7(10):e1002313. doi:10.1371/journal.pgen.1002313

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

25. Jin G, Sun J, Kim ST, Feng J, Wang Z, Tao S, et al. Genome-wide association study identifies a new locus JMJD1C at 10q21 that may influence serum androgen levels in men. Hum Mol Genet (2012) 21(23):5222–8. doi:10.1093/hmg/dds361

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

26. Chen Z, Tao S, Gao Y, Zhang J, Hu Y, Mo L, et al. Genome-wide association study of sex hormones, gonadotropins and sex hormone-binding protein in Chinese men. J Med Genet (2013) 50(12):794–801. doi:10.1136/jmedgenet-2013-101705

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

27. Prescott J, Thompson DJ, Kraft P, Chanock SJ, Audley T, Brown J, et al. Genome-wide association study of circulating estradiol, testosterone, and sex hormone-binding globulin in postmenopausal women. PLoS One (2012) 7(6):e37815. doi:10.1371/journal.pone.0037815

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

28. Caramaschi D, Booij L, Petitclerc A, Boivin M, Tremblay RE. Genetic and environmental contributions to saliva testosterone levels in male and female infant twins. Psychoneuroendocrinology (2012) 37(12):1954–9. doi:10.1016/j.psyneuen.2012.04.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

29. Forest MG, Sizonenko PC, Cathiard AM, Bertrand J. Hypophyso-gonadal function in humans during the first year of life: I. Evidence for testicular activity in early infancy. J Clin Invest (1974) 53:819–28. doi:10.1172/JCI107621

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

30. Konforte D, Shea JL, Kyriakopoulou L, Colantonio D, Cohen AH, Shaw J, et al. Complex biological pattern of fertility hormones in children and adolescents: a study of healthy children from the CALIPER cohort and establishment of pediatric reference intervals. Clin Chem (2013) 59(8):1215–27. doi:10.1373/clinchem.2013.204123

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

31. Bidlingmaier F, Dorr HG, Eisenmenger W, Kuhnle U, Knorr D. Contribution of the adrenal gland to the production of androstenedione and testosterone during the first two years of life. J Clin Endocrinol Metab (1986) 62(2):331–5. doi:10.1210/jcem-62-2-331

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

32. Prince FP. The triphasic nature of Leydig cell development in humans, and comments on nomenclature. J Endocrinol (2001) 168(2):213–6. doi:10.1677/joe.0.1680213

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

33. Knickmeyer RC, Woolson S, Hamer RM, Konneker T, Gilmore JH. 2D:4D ratios in the first 2 years of life: stability and relation to testosterone exposure and sensitivity. Horm Behav (2011) 60(3):256–63. doi:10.1016/j.yhbeh.2011.05.009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

34. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS Genet (2006) 2(12):e190. doi:10.1371/journal.pgen.0020190

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

35. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet (2006) 38(8):904–9. doi:10.1038/ng1847

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet (2007) 81(3):559–75. doi:10.1086/519795

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

37. Ahn M, Zhang HH, Lu W. Moment-based method for random effects selection in linear mixed models. Stat Sin (2012) 22(4):1539–62. doi:10.5705/ss.2011.054

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

38. Liaw A, Wiener M. Classification and regression by randomforest. R News (2002) 2(3):18–22. Available from: http://cran.r-project.org/doc/Rnews/Rnews_2002-3.pdf

Google Scholar

39. Garagorri JM, Rodriguez G, Lario-Elboj AJ, Olivares JL, Lario-Munoz A, Orden I. Reference levels for 17-hydroxyprogesterone, 11-desoxycortisol, cortisol, testosterone, dehydroepiandrosterone sulfate and androstenedione in infants from birth to six months of age. Eur J Pediatr (2007) 167(6):647–53. doi:10.1007/s00431-007-0565-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

40. Ji C, Huang XW, Yang RW, Wang X, Zhao ZY. Gonadotropins and sex hormones in healthy Chinese infants. Indian Pediatr (2008) 45(6):489–92. Available from: http://indianpediatrics.net/june2008/june-489-492.htm

Pubmed Abstract | Pubmed Full Text | Google Scholar

41. Granger DA, Cicchetti D, Rogosch FA, Hibel LC, Teisl M, Flores E. Blood contamination in children’s saliva: prevalence, stability, and impact on the measurement of salivary cortisol, testosterone, and dehydroepiandrosterone. Psychoneuroendocrinology (2007) 32(6):724–33. doi:10.1016/j.psyneuen.2007.05.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

42. Apgar V. A proposal for a new method of evaluation of the newborn infant. Curr Res Anesth Analg (1953) 32(4):260–7. doi:10.1213/00000539-195301000-00041

CrossRef Full Text | Google Scholar

43. AAP ACOG. The Apgar score. Pediatrics (2006) 117(4):1444–7. doi:10.1542/peds.2006-0325

CrossRef Full Text | Google Scholar

44. Li J, Cnattingus S, Gissler M, Vestergaard M, Obel C, Ahrensberg J, et al. The 5-minute Apgar score as a predictor of childhood cancer: a population-based cohort study in five million children. BMJ Open (2012) 2(4):e001095. doi:10.1136/bmjopen-2012-001095

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

45. Velazquez-Moctezuma J, Dominguez Salazar E, Cruz Rueda ML. The effect of prenatal stress on adult sexual behavior in rats depends on the nature of the stressor. Physiol Behav (1993) 53(3):443–8. doi:10.1016/0031-9384(93)90137-5

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

46. Wang CT, Huang RL, Tai MY, Tsai YF, Peng MT. Dopamine release in the nucleus accumbens during sexual behavior in prenatally stressed adult male rats. Neurosci Lett (1995) 200(1):29–32. doi:10.1016/0304-3940(95)12071-B

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

47. Rhees RW, Al-Saleh HN, Kinghorn EW, Fleming DE, Lephart ED. Relationship between sexual behavior and sexually dimorphic structures in the anterior hypothalamus in control and prenatally stressed male rats. Brain Res Bull (1999) 50(3):193–9. doi:10.1016/S0361-9230(99)00191-4

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

48. Ward IL, Ward OB, Hayden T, Weisz J, Orth JM. Naltrexone normalizes the suppression but not the surge of delta 5-3 beta-hydroxysteroid dehydrogenase activity in Leydig cells of stressed rat fetuses. Endocrinology (1990) 127(1):88–92. doi:10.1210/endo-127-1-88

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

49. Ward IL, Ward OB, Affuso JD, Long WD III, French JA, Hendricks SE. Fetal testosterone surge: specific modulations induced in male rats by maternal stress and/or alcohol consumption. Horm Behav (2003) 43(5):531–9. doi:10.1016/S0018-506X(03)00061-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

50. Pallares ME, Adrover E, Baier CJ, Bourguignon NS, Monteleone MC, Brocco MA, et al. Prenatal maternal restraint stress exposure alters the reproductive hormone profile and testis development of the rat male offspring. Stress (2013) 16(4):429–40. doi:10.3109/10253890.2012.761195

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

51. Shirtcliff EA, Granger DA, Likos A. Gender differences in the validity of testosterone measured in saliva by immunoassay. Horm Behav (2002) 42(1):62–9. doi:10.1006/hbeh.2002.1798

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

52. Visscher PM, Hill WG, Wray NR. Heritability in the genomics era – concepts and misconceptions. Nat Rev Genet (2008) 9(4):255–66. doi:10.1038/nrg2322

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

53. Johnson W, Penke L, Spinath FM. Heritability in the era of molecular genetics: some thoughts for understanding genetic influences on behavioral traits. Eur J Pers (2011) 25(4):254–66. doi:10.1002/per.834

CrossRef Full Text | Google Scholar

54. Miller WL, Auchus RJ. The molecular biology, biochemistry, and physiology of human steroidogenesis and its disorders. Endocr Rev (2011) 32(1):81–151. doi:10.1210/er.2010-0013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

55. Yang JA, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, et al. Common SNPs explain a large proportion of the heritability for human height. Nat Genet (2010) 42(7):565–U131. doi:10.1038/ng.608

CrossRef Full Text | Google Scholar

56. Visscher PM, Goddard ME, Derks EM, Wray NR. Evidence-based psychiatric genetics, AKA the false dichotomy between common and rare variant hypotheses. Mol Psychiatry (2012) 17(5):474–85. doi:10.1038/mp.2011.65

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

57. Safran M, Chalifa-Caspi V, Shmueli O, Olender T, Lapidot M, Rosen N, et al. Human gene-centric databases at the weizmann institute of science: genecards, UDB, CroW 21 and HORDE. Nucleic Acids Res (2003) 31(1):142–6. doi:10.1093/nar/gkg050

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

58. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res (2001) 29(1):308–11. doi:10.1093/nar/29.1.308

CrossRef Full Text | Google Scholar

59. Karolchik DM, Barber GP, Clawson H, Cline MS, Diekhans M, Dreszer TR, et al. The UCSC Genome Browser database: 2014 update. Nucelic Acids Res (2014) 42:D764–70. doi:10.1093/nar/gkt1168

CrossRef Full Text | Google Scholar

60. Goodarzi MO, Antoine HJ, Azziz R. Genes for enzymes regulating dehydroepiandrosterone sulfonation are associated with levels of dehydroepiandrosterone sulfate in polycystic ovary syndrome. J Clin Endocrinol Metab (2007) 92(7):2659–64. doi:10.1210/jc.2006-2600

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

61. Groschl M. Current status of salivary hormone analysis. Clin Chem (2008) 54(11):1759–69. doi:10.1373/clinchem.2008.108910

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

62. Baxendale PM, Reed MJ, James VHT. Testosterone in saliva of normal men and its relationship with unbound and total testosterone levels in plasma. J Endocrinol (1980) 4:46–7.

Google Scholar

63. Vittek J, L’Hommedieu DG, Gordon GG, Rappaport SC, Southren AL. Direct radioimmunoassay (RIA) of salivary testosterone: correlation with free and total serum testosterone. Life Sci (1985) 37(8):711–6. doi:10.1016/0024-3205(85)90540-5

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

64. Navarro MA, Juan L, Bonnin MR, Villabona C. Salivary testosterone: relationship to total and free testosterone in serum. Clin Chem (1986) 32(1 Pt 1):231–2.

Google Scholar

65. Magrini G, Chiodoni G, Rey F, Felber JP. Further evidence for the usefulness of the salivary testosterone radioimmunoassay in the assessment of androgenicity in man in basal and stimulated conditions. Horm Res (1986) 23(2):65–73. doi:10.1159/000180291

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

66. Schwartz EB, Granger DA. Transferrin enzyme immunoassay for quantitative monitoring of blood contamination in saliva. Clin Chem (2004) 50(3):654–6. doi:10.1373/clinchem.2003.028266

CrossRef Full Text | Google Scholar