Professor Craig Pennell

Grip strength (GS) is a proxy measure for muscular strength and a predictor for bone fracture risk among other diseases. Previous genome-wide association studies (GWASs) have been conducted in large cohorts of adults focusing on scores collected for the dominant hand, therefore increasing the likelihood of confounding effects by environmental factors. Here, we perform the first GWAS meta-analyses on maximal GS with the dominant (GSD) and non-dominant (GSND) hand in two cohorts of children (ALSPAC, N = 5450; age range = 10.65¿13.61; Raine Study, N = 1162, age range: 9.42¿12.38 years). We identified a novel significant association for GSND (rs9546244, LINC02465, p = 3.43e-08) and replicated associations previously reported in adults including with a HOXB3 gene marker that shows an expression quantitative trait locus (eQTL) effect. Despite a much smaller sample (~3%) compared with the UK Biobank we replicated correlation analyses previously reported in this much larger adult cohort, such as a negative correlation with coronary artery disease. Although the results from the polygenic risk score (PRS) analyses did not survive multiple testing correction, we observed nominally significant associations between GS and risk of overall fracture, as previously reported, as well ADHD which will require further investigations. Finally, we observed a higher SNP-heritability (24%¿41%) compared with previous studies (4%¿24%) in adults. Overall, our results suggest that cohorts of children might be better suited for genetic studies of grip strength, possibly due to the shorter exposure to confounding environmental factors compared with adults.

Background: Obesity rates have nearly tripled in the past 50 years, and by 2030 more than 1 billion individuals worldwide are projected to be obese. This creates a significant economic strain due to the associated non-communicable diseases. The root cause is an energy expenditure imbalance, owing to an interplay of lifestyle, environmental, and genetic factors. Obesity has a polygenic genetic architecture; however, single genetic variants with large effect size are etiological in a minority of cases. These variants allowed the discovery of novel genes and biology relevant to weight regulation and ultimately led to the development of novel specific treatments. Methods: We used a case-control approach to determine metabolic differences between individuals homozygous for a loss-of-function genetic variant in the small integral membrane protein 1 (SMIM1) and the general population, leveraging data from five cohorts. Metabolic characterization of SMIM1-/- individuals was performed using plasma biochemistry, calorimetric chamber, and DXA scan. Findings: We found that individuals homozygous for a loss-of-function genetic variant in SMIM1 gene, underlying the blood group Vel, display excess body weight, dyslipidemia, altered leptin to adiponectin ratio, increased liver enzymes, and lower thyroid hormone levels. This was accompanied by a reduction in resting energy expenditure. Conclusion: This research identified a novel genetic predisposition to being overweight or obese. It highlights the need to investigate the genetic causes of obesity to select the most appropriate treatment given the large cost disparity between them. Funding: This work was funded by the National Institute of Health Research, British Heart Foundation, and NHS Blood and Transplant.

This study examined associations between pregnancy and infant birth outcomes with child telomere length at age 17 years; and investigated if there are sex differences between pregnancy complications and telomere length. We utilised the population-based prospective Raine cohort study in Western Australia, Australia. 2900 pregnant women were recruited at 16¿20 weeks' gestation (Gen 1), and their children (Gen 2) were followed up over several years. Generalised linear models were used to examine relationships between pregnancy or birth outcomes (gestational diabetes, pre-eclampsia, preterm birth, low birth weight, macrosomia), and as a composite, with telomere length, measured via a DNA sample from blood at 17 years of age. Analyses were adjusted for a range of confounders. Among the 1202 included children, there were no differences in child telomere length for any of the individual maternal or birth weight pregnancy outcomes nor were there any significant interactions between each of the complications (individual or composite) and the sex of the child. However, females born from any of the 5 adverse outcomes had shorter telomeres (estimated mean (SE) = -0.159 (0.061), p = 0.010) than females born in the absence of these complications. Specifically, females born from a pre-eclamptic pregnancy had shorter telomeres than females not born from a preeclamptic pregnancy (estimated mean (SE) = -0.166 (0.072), p = 0.022). No relationships were observed in males. Further longitudinal studies are needed to understand mediating factors that are important in predicting offspring telomere length and the necessity to investigate females and males independently.

PURPOSE. Changes in refractive error during young adulthood is common yet risk factors at this age are largely unexplored. This study explored risk factors for these changes, including gene¿environmental interactions. METHODS. Spherical equivalent refraction (SER) and axial length (AL) for 624 community-based adults were measured at 20 (baseline) and 28 years old. Participants were genotyped and their polygenic scores (PGS) for refractive error calculated. Self-reported screen time (computer, television, and mobile devices) from 20 to 28 years old were collected prospectively and longitudinal trajectories were generated. Past sun exposure was quantified using conjunctival ultraviolet autofluorescence (CUVAF) area. RESULTS. Median change in SER and AL were -0.023 diopters (D)/year (interquartile range [IQR] = -0.062 to ¿0.008) and +0.01 mm/year (IQR = 0.000 to 0.026), respectively. Sex, baseline myopia, parental myopia, screen time, CUVAF, and PGS were significantly associated with myopic shift. Collectively, these factors accounted for approximately 20% of the variance in refractive error change, with screen time, CUVAF, and PGS each explaining approximately 1% of the variance. Four trajectories for total screen time were found: "consistently low" (n = 148), "consistently high" (n = 250), "consistently very high" (n = 76), and "increasing" (n = 150). Myopic shift was faster in those with "consistently high" or "consistently very high" screen time compared to "consistently-low" (P = 0.031). For each z-score increase in PGS, changes in SER and AL increased by -0.005 D/year and 0.002 mm/year (P = 0.045). Of the three types of screen time, only computer time was associated with myopic shift (P = 0.040). There was no two- or three-way interaction effect between PGS, CUVAF, or screen time (P = 0.26). CONCLUSIONS. Higher total or computer screen time, less sun exposure, and genetic predisposition are each independently associated with greater myopic shifts during young adulthood. Given that these factors explained only a small amount of the variance, there are likely other factors driving refractive error change during young adulthood.

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Aims/hypothesis: Characterisation of genetic variation that influences the response to glucose-lowering medications is instrumental to precision medicine for treatment of type 2 diabetes. The Study to Understand the Genetics of the Acute Response to Metformin and Glipizide in Humans (SUGAR-MGH) examined the acute response to metformin and glipizide in order to identify new pharmacogenetic associations for the response to common glucose-lowering medications in individuals at risk of type 2 diabetes. Methods: One thousand participants at risk for type 2 diabetes from diverse ancestries underwent sequential glipizide and metformin challenges. A genome-wide association study was performed using the Illumina Multi-Ethnic Genotyping Array. Imputation was performed with the TOPMed reference panel. Multiple linear regression using an additive model tested for association between genetic variants and primary endpoints of drug response. In a more focused analysis, we evaluated the influence of 804 unique type 2 diabetes- and glycaemic trait-associated variants on SUGAR-MGH outcomes and performed colocalisation analyses to identify shared genetic signals. Results: Five genome-wide significant variants were associated with metformin or glipizide response. The strongest association was between an African ancestry-specific variant (minor allele frequency [MAFAfr]=0.0283) at rs149403252 and lower fasting glucose at Visit 2 following metformin (p=1.9×10-9); carriers were found to have a 0.94 mmol/l larger decrease in fasting glucose. rs111770298, another African ancestry-specific variant (MAFAfr=0.0536), was associated with a reduced response to metformin (p=2.4×10-8), where carriers had a 0.29 mmol/l increase in fasting glucose compared with non-carriers, who experienced a 0.15 mmol/l decrease. This finding was validated in the Diabetes Prevention Program, where rs111770298 was associated with a worse glycaemic response to metformin: heterozygous carriers had an increase in HbA1c of 0.08% and non-carriers had an HbA1c increase of 0.01% after 1 year of treatment (p=3.3×10-3). We also identified associations between type 2 diabetes-associated variants and glycaemic response, including the type 2 diabetes-protective C allele of rs703972 near ZMIZ1 and increased levels of active glucagon-like peptide 1 (GLP-1) (p=1.6×10-5), supporting the role of alterations in incretin levels in type 2 diabetes pathophysiology. Conclusions/interpretation: We present a well-phenotyped, densely genotyped, multi-ancestry resource to study gene¿drug interactions, uncover novel variation associated with response to common glucose-lowering medications and provide insight into mechanisms of action of type 2 diabetes-related variation. Data availability: The complete summary statistics from this study are available at the Common Metabolic Diseases Knowledge Portal (https://hugeamp.org) and the GWAS Catalog (www.ebi.ac.uk/gwas/ , accession IDs: GCST90269867 to GCST90269899). Graphical Abstract: [Figure not available: see fulltext.]

Refractive error, measured here as mean spherical equivalent (SER), is a complex eye condition caused by both genetic and environmental factors. Individuals with strong positive or negative values of SER require spectacles or other approaches for vision correction. Common genetic risk factors have been identified by genome-wide association studies (GWAS), but a great part of the refractive error heritability is still missing. Some of this heritability may be explained by rare variants (minor allele frequency [MAF] = 0.01.). We performed multiple gene-based association tests of mean Spherical Equivalent with rare variants in exome array data from the Consortium for Refractive Error and Myopia (CREAM). The dataset consisted of over 27,000 total subjects from five cohorts of Indo-European and Eastern Asian ethnicity. We identified 129 unique genes associated with refractive error, many of which were replicated in multiple cohorts. Our best novel candidates included the retina expressed PDCD6IP, the circadian rhythm gene PER3, and P4HTM, which affects eye morphology. Future work will include functional studies and validation. Identification of genes contributing to refractive error and future understanding of their function may lead to better treatment and prevention of refractive errors, which themselves are important risk factors for various blinding conditions.

Reading and writing are crucial life skills but roughly one in ten children are affected by dyslexia, which can persist into adulthood. Family studies of dyslexia suggest heritability up to 70%, yet few convincing genetic markers have been found. Here we performed a genome-wide association study of 51,800 adults self-reporting a dyslexia diagnosis and 1,087,070 controls and identified 42 independent genome-wide significant loci: 15 in genes linked to cognitive ability/educational attainment, and 27 new and potentially more specific to dyslexia. We validated 23 loci (13 new) in independent cohorts of Chinese and European ancestry. Genetic etiology of dyslexia was similar between sexes, and genetic covariance with many traits was found, including ambidexterity, but not neuroanatomical measures of language-related circuitry. Dyslexia polygenic scores explained up to 6% of variance in reading traits, and might in future contribute to earlier identification and remediation of dyslexia.

Background: The onset of the term human parturition involves myometrial gene expression changes to transform the uterus from a quiescent to a contractile phenotype. It is uncertain whether the same changes occur in the uterus during preterm labor. Objective: This study aimed to compare the myometrial gene expression between term and preterm labor and to determine whether the presence of acute clinical chorioamnionitis or twin gestation affects these signatures. Study Design: Myometrial specimens were collected during cesarean delivery from the following 7 different groups of patients: term not in labor (n=31), term labor (n=13), preterm not in labor (n=21), preterm labor with acute clinical chorioamnionitis (n=6), preterm labor with no acute clinical chorioamnionitis (n=9), twin preterm not in labor (n=8), and twin preterm labor with no acute clinical chorioamnionitis (n=5). RNA was extracted, reverse transcribed and quantitative polymerase chain reactions were performed on 44 candidate genes (with evidence for differential expression in human term labor) using the Fluidigm platform. Computational analysis was performed using 2-class unpaired Wilcoxon tests and principal component analysis. Results: Computational analysis revealed that gene expression in the preterm myometrium, irrespective of whether in labor or not in labor, clustered tightly and is clearly different from the term labor and term not-in-labor groups. This was true for both singleton and twin pregnancies. Principal component analysis showed that 57% of the variation was explained by 3 principal components. These 44 genes interact in themes of prostaglandin activity and inflammatory signaling known to be important during term labor, but are not a full representation of the myometrium transcriptional activity. Conclusion: The myometrial contractions associated with preterm labor are associated with a pattern of gene expression that is distinct from term labor. Therefore, preterm labor may be initiated by a different myometrial process or processes outside the myometrium.

Background: Whilst cortisol reactivity has been associated with depression and anxiety disorders, research examining cortisol reactivity with early symptoms of these conditions in males and females is limited. Methods: At age 18, 748 males and females from Gen2 of the Raine Study were assessed for their salivary cortisol response to a psychosocial stressor using the Trier Social Stress Test (TSST). Participants later completed the Depression Anxiety Stress Scale (DASS-21) at age 20 which was used as the outcome measure in regression models. Results: We found differences in DASS-21 across TSST responder categories in females but not males. Female reactive-responders (RR) and non-responders (NR) had increased symptoms of depression and anxiety compared to anticipatory-responders (AR). AR were associated with the lowest symptomology in females. We found limited evidence for an association between salivary cortisol summary measures (CBL, CMAX, CMIN, CRANGE, AUCG and AUCR) and depression/anxiety symptoms at age 20. Conclusions: This study sheds new light on adaptive and maladaptive physiological responses to psychosocial stress in terms of depression and anxiety symptoms. These preliminary findings indicate the pattern of response to a psychosocial stressor may contribute to individual vulnerability for stress-related diseases in a sex-specific manner.

Objective: To determine whether vaginal progesterone reduces spontaneous preterm birth (sPTB) before 37 weeks in asymptomatic high-risk women with a singleton pregnancy and normal mid-gestation cervical length. Study design: Databases were searched (from inception to December 2020) with the search terms "progesterone" and "premature birth" or "preterm birth". Studies were screened and included if they assessed vaginal progesterone compared to placebo in women with normal cervical length. Data were pooled and synthesized in a meta-analysis using a random effects model. Data sources: MEDLINE and Embase databases. Study synthesis: Following PRISMA screening guidelines, data from 1127 women across three studies were available for synthesis. All studies had low risk of bias and were of high quality. The primary outcome was sPTB <37 weeks, with secondary outcomes of sPTB <34 weeks. Vaginal progesterone did not significantly reduce sPTB before 37 weeks, or before 34 weeks with a relative risk (RR) of 0.76 (95% CI 0.37¿1.55, p =.45) and 0.51 (95% CI 0.12¿2.13, p =.35), respectively. Conclusions: Vaginal progesterone does not decrease the risk of sPTB in high-risk singleton pregnancies with a normal mid-gestation cervical length.

We test whether genetic influences that explain individual differences in aggression in early life also explain individual differences across the life-course. In two cohorts from The Netherlands (N = 13,471) and Australia (N = 5628), polygenic scores (PGSs) were computed based on a genome-wide meta-analysis of childhood/adolescence aggression. In a novel analytic approach, we ran a mixed effects model for each age (Netherlands: 12¿70¿years, Australia: 16¿73¿years), with observations at the focus age weighted as 1, and decaying weights for ages further away. We call this approach a 'rolling weights' model. In The Netherlands, the estimated effect of the PGS was relatively similar from age 12 to age 41, and decreased from age 41¿70. In Australia, there was a peak in the effect of the PGS around age 40¿years. These results are a first indication from a molecular genetics perspective that genetic influences on aggressive behavior that are expressed in childhood continue to play a role later in life.

It is well established that genetics, environment, and interplay between them play a crucial role in adult disease. We aimed to evaluate the role of genetics, early life nutrition, and the interaction between them, on optimal adult health. As part of a large international consortium (n ~ 154,000), we identified 60 SNPs associated with both birthweight and adult disease. Utilising the Raine Study, we developed a birthweight polygenic score (BW-PGS) based on the 60 SNPs and examined relationships between BW-PGS and adulthood cardiovascular risk factors, specifically evaluating interactions with early life nutrition. Healthy nutrition was beneficial for all individuals; longer duration of any breastfeeding was particularly associated with lower BMI and lower Systolic Blood Pressure in those with higher BW-PGS. Optimal breastfeeding offers the greatest benefit to reduce adult obesity and hypertension in those genetically predisposed to high birthweight. This provides an example of how precision medicine in early life can improve adult health.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sexdifferential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10-8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10-6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV1) (P=3.15x10-15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV1 more in males (untransformed FEV1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein (HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10-6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Background: Early life stress exposures may cause dysregulation of the Hypothalamic Pituitary Adrenal (HPA)-axis and cortisol production, with timing and sex-specific effects. Studies examining the impact of early life stress on cortisol responses to stress have focused on severe trauma and have produced inconsistent results. The aim of this study was to investigate whether common early life stressors, experienced prenatally or throughout childhood and adolescence, play a role in the dysregulation of the HPA-axis in early adulthood. Methods: Exposures to common life stress events were examined prenatally and as longitudinal trajectories of stress exposure from birth to age 17 in males and females from Gen2 of the Raine Study. At age 18 years, 986 participants were assessed for their salivary cortisol response to a psychosocial stressor - the Trier Social Stress Test (TSST). Results: In males there was an association between high prenatal stress exposure at 18 weeks gestation and a heightened TSST response. We found evidence for sex-specific associations with increasing stress exposure during adolescence (the ascending trajectory) whereby males were more likely to be non-responders to the TSST and females were more likely to be responders. Conclusion: Our results point to sex differences in how stress exposure in-utero and exposure increasing during adolescence may affect regulation of the HPA-axis later in life. However, overall common life stress events experienced in-utero, during childhood and adolescence show limited impact on the HPA-axis stress response in early adulthood.

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Objective: Genetic factors underpin the narrow intraindividual variability of thyroid function, although precise contributions of environmental vs genetic factors remain uncert ain. We sought to clarify the heritability of thyroid function traits and thyroid peroxidase antibody (TPOAb) positiv ity and identify single nucleotide polymorphisms (SNPs) contributing to the trait variance. Methods: Heritability of thyroid-stimulating hormone (TSH), free T4 (fT4 ), free T3 (fT3) and TPOAb in a cohort of 2854 euthyroid, dizygous and monozygous twins (age range 11.9-16.9 y ears) from the Brisbane Longitudinal Twin Study (BLTS) was assessed using structural equation modelling. A geno me-wide analysis was conducted on 2832 of these individuals across 7 522 526 SNPs as well as gene-based associa tion analyses. Replication analysis of the association results was performed in the Raine Study (n = 1115) followed by meta-analysis to maximise power for discover y. Results: Heritability of thyroid function parameters in the BLTS was 70. 8% (95% CI: 66.7-74.9%) for TSH, 67.5% (59.8-75.3%) for fT4, 59.7% (54.4-65.0%) for fT3 and 48.8% (40. 6-56.9%) for TPOAb. The genome-wide association study (GWAS) in the discovery cohort identified a novel associat ion between rs2026401 upstream of NCOA3 and TPOAb. GWAS meta-analysis found associations between TPOAb and rs445219, also near NCOA3, and fT3 and rs12687280 near SERPINA7. Gene-based association analysis highlighted SERPINA7 for fT3 and NPAS3 for fT4. Conclusion: Our findings resolve former contention regarding heritability es timates of thyroid function traits and TPOAb positivity. GWAS and gene-based association analysis iden tified variants accounting for a component of this heritability 2021 European Society of Endocrinology Printed in Great Britain.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15¿17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype¿phenotype map than previously anticipated.

Epidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study, we evaluated the causal association of low birthweight with adulthood hypertension following a standard analytic protocol using the study-level data of 183,433 participants from 60 studies (CHARGE-BIG consortium), as well as that with blood pressure using publicly available summary-level genome-wide association data from EGG consortium of 153,781 participants, ICBP consortium and UK Biobank cohort together of 757,601 participants. We used seven SNPs as the instrumental variable in the study-level analysis and 47 SNPs in the summary-level analysis. In the study-level analyses, decreased birthweight was associated with a higher risk of hypertension in adults (the odds ratio per 1 standard deviation (SD) lower birthweight, 1.22; 95% CI 1.16 to 1.28), while no association was found between genetically instrumented birthweight and hypertension risk (instrumental odds ratio for causal effect per 1 SD lower birthweight, 0.97; 95% CI 0.68 to 1.41). Such results were consistent with that from the summary-level analyses, where the genetically determined low birthweight was not associated with blood pressure measurements either. One SD lower genetically determined birthweight was not associated with systolic blood pressure (ß = - 0.76, 95% CI - 2.45 to 1.08¿mmHg), 0.06¿mmHg lower diastolic blood pressure (ß = - 0.06, 95% CI - 0.93 to 0.87¿mmHg), or pulse pressure (ß = - 0.65, 95% CI - 1.38 to 0.69¿mmHg, all p > 0.05). Our findings suggest that the inverse association of birthweight with hypertension risk from observational studies was not supported by large Mendelian randomization analyses.

Humans, similar to many other species, exhibit developmental plasticity ¿ the ability to modify the structure and function of key vital organs during specific developmental windows. Adverse exposures during these critical periods may alter organ function to enhance survival; however, over many years ¿ these alterations may predispose toward chronic disease. The hypothalamic¿pituitary¿adrenal axis plays an important role in maintaining metabolic homoeostasis and coordinating the neuropsychiatric response to stress. Decades of animal and human research suggests that the long-term function of this axis may be particularly sensitive to adverse exposures during critical developmental windows. The aim of this review is to summarise the existing knowledge of hypothalamic¿pituitary¿adrenal axis programming in humans and briefly discuss its implications for long-term health.

BACKGROUND: Associations between dairy intake and body composition and cardiometabolic traits have been inconsistently observed in epidemiological studies, and the causal relationship remains ill-defined. METHODS: We performed Mendelian randomization analysis using an established genetic variant located upstream of the lactase gene (LCT-13910 C/T, rs4988235) associated with dairy intake as an instrumental variable (IV). The causal effects of dairy intake on body composition and cardiometabolic traits (lipids, glycemic traits, and inflammatory factors) were quantified by IV estimators among 182041 participants from 18 studies. RESULTS: Each 1 serving/day higher dairy intake was associated with higher lean mass [ß (SE) = 0.117 kg (0.035); P = 0.001], higher hemoglobin A1c [0.009% (0.002); P < 0.001], lower LDL [-0.014 mmol/L (0.006); P = 0.013], total cholesterol (TC) [-0.012 mmol/L (0.005); P = 0.023], and non-HDL [-0.012 mmol/L (0.005); P = 0.028]. The LCT-13910 C/T CT + TT genotype was associated with 0.214 more dairy servings/day (SE = 0.047; P < 0.001), 0.284 cm higher waist circumference (SE = 0.118; P = 0.017), 0.112 kg higher lean mass (SE = 0.027; P = 3.8 × 10-5), 0.032 mmol/L lower LDL (SE = 0.009; P = 0.001), and 0.032 mmol/L lower TC (SE = 0.010; P = 0.001). Genetically higher dairy intake was associated with increased lean mass [0.523 kg per serving/day (0.170); P = 0.002] after correction for multiple testing (0.05/18). However, we find that genetically higher dairy intake was not associated with lipids and glycemic traits. CONCLUSIONS: The present study provides evidence to support a potential causal effect of higher dairy intake on increased lean mass among adults. Our findings suggest that the observational associations of dairy intake with lipids and glycemic traits may be the result of confounding.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Objective: To compare rates of detectability of circulating Rh(D)-immunoglobulin (anti-D) at delivery with single and two-dose antenatal anti-D prophylaxis (RAADP) regimens; to compare compliance with the two regimens. Design: Open label, randomised controlled trial between May 2013 and November 2015. Setting, participants: 277 women who attended a tertiary obstetric referral hospital in Perth for antenatal care and were at least 18 years of age, less than 30 weeks pregnant and yet to receive RAADP, Rh(D)-negative (negative antibody screen), and who intended to deliver their baby at the hospital. Exclusion criteria were prior anti-D sensitisation, any contraindication of anti-D administration, and a history of isolated IgA deficiency. Interventions: One 1500¿IU anti-D dose at 28 weeks of pregnancy (single dose regimen); two doses of 625¿IU each at 28 and 34 weeks of pregnancy (two-dose regimen). Main outcome measures: The primary outcome was the proportion of women with detectable anti-D levels at delivery; the secondary outcome was compliance with the allocated RAADP regimen. Results: Circulating anti-D was detectable at delivery in a greater proportion of women in the two-dose group (111 of 129, 86%) than in the single dose group (70 of 125, 56%; P¿<¿0.001). Compliance was not significantly different between the single dose (86 of 138, 61%) and two-dose groups (70 of 139, 50%; P¿=¿0.06). Conclusions: The two-dose RAADP schedule currently recommended in Australia provides better protection against Rh(D) sensitisation than a one-dose regimen. Trial registration: Australian and New Zealand Clinical Trials Registry (ACTRN12613000661774).

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF =5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Macronutrient intake, the proportion of calories consumed from carbohydrate, fat, and protein, is an important risk factor for metabolic diseases with significant familial aggregation. Previous studies have identified two genetic loci for macronutrient intake, but incomplete coverage of genetic variation and modest sample sizes have hindered the discovery of additional loci. Here, we expanded the genetic landscape of macronutrient intake, identifying 12 suggestively significant loci (P < 1 × 10-6) associated with intake of any macronutrient in 91,114 European ancestry participants. Four loci replicated and reached genome-wide significance in a combined meta-analysis including 123,659 European descent participants, unraveling two novel loci; a common variant in RARB locus for carbohydrate intake and a rare variant in DRAM1 locus for protein intake, and corroborating earlier FGF21 and FTO findings. In additional analysis of 144,770 participants from the UK Biobank, all identified associations from the two-stage analysis were confirmed except for DRAM1. Identified loci might have implications in brain and adipose tissue biology and have clinical impact in obesity-related phenotypes. Our findings provide new insight into biological functions related to macronutrient intake.

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight¿blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (=95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Objective: To investigate prevalence rates and the risk of anteand intrapartum stillbirth in Western Australia with respect to maternal country of birth and ethnic origin. Design, setting and participants: Whole population retrospective cohort analysis of de-identified, linked routinely collected birth, perinatal and mortality data for all births to non- Indigenous women in WA during 2005e2013. Main outcome measures: Crude and adjusted odds ratios (aORs) with 95% confidence intervals were estimated by logistic regression and adjusted for confounding factors, for all stillbirths, antepartum stillbirths and intrapartum stillbirths, stratified by migrant status and ethnic background (white, Asian, Indian, African, Ma¯ori, other). Results: Women born overseas were more likely to have a stillbirth than Australian-born women (aOR, 1.26; 95% CI, 1.09e1.37). There was no significant difference for any type of stillbirth between Australian-born women of white and nonwhite backgrounds, but non-white migrant women were more likely than white migrants to have a stillbirth (OR, 1.42; 95% CI, 1.19e1.70). Compared with Australian-born women, migrants of Indian (aOR, 1.71; 95% CI, 1.17e2.47), African (aOR, 2.12; 95% CI, 1.46e3.08), and "other" ethnic origins (aOR, 1.43; 95% CI, 1.06e1.93) were more likely to have antepartum stillbirths; women of African (aOR, 5.08; 95% CI, 3.14e8.22) and "other" (aOR, 1.86; 95% CI, 1.15e3.00) background were more likely to have an intrapartum stillbirth. Conclusions: Immigrants of African or Indian background appear to be at greater risk of ante- and intrapartum stillbirth in WA. Specific strategies are needed reduce the prevalence of stillbirth in these communities.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sexspecific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N=15 661, with replication N=75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (b=0.06; P=0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

Background:An association between depression and obesity is well recognised, but longitudinal studies of depressive symptoms in adolescents as a predictor of body composition are lacking.Objective:We examined depressive symptoms at age 14, 17 and 20 years as predictors of lean, fat and bone mass at age 20 years in a birth cohort.Subjects/Methods:In 1161 participants (569 females) in the Western Australia Pregnancy Cohort (Raine) Study, depressive symptoms were assessed using the Beck Depression Inventory for Youth at age 14 and 17 years, and the Depression, Anxiety and Stress Scale 21 at age 20 years. Participants were further classified into two trajectories using latent class analysis: no/transient and persistent/recurrent depression. At age 20 years, lean body mass (LBM), fat body mass (FBM) and total body bone mass were measured by dual-energy X-ray absorptiometry.Results:In females, accounting for age and lifestyle factors, depression scores at age 14 and 20 years were positively associated with body weight, body mass index (BMI), FBM and % FBM (r=0.110-0.184, P<0.05) but negatively correlated with % LBM (r=-0.120, P<0.05) at age 20 years. Females in the persistent/recurrent depression trajectory (n=99) had significantly higher body weight (+5.1 kg), BMI (+1.8 kg m -2), FBM (+3.9 kg) and % FBM (+2.2%) and significantly lower % LBM (-2.2%) at age 20 years than those with no/transient depression (n=470; all P<0.05). In males, depression scores at age 17 and 20 years were negatively associated with LBM but not weight or BMI, and depression trajectory was not a predictor of body composition at age 20 years. Depression scores and trajectories did not predict bone mass in either males or females.Conclusions:Depressive symptoms and persistent/recurrent depression in adolescence are predictors of greater adiposity at age 20 years in females, but not males, but do not predict bone mass in either gender.

Hair cortisol concentration (HCC) is a promising measure of long-Term hypothalamus-pituitary-Adrenal (HPA) axis activity. Previous research has suggested an association between HCC and psychological variables, and initial studies of inter-individual variance in HCC have implicated genetic factors. However, whether HCC and psychological variables share genetic risk factors remains unclear. The aims of the present twin study were to: (i) assess the heritability of HCC; (ii) estimate the phenotypic and genetic correlation between HPA axis activity and the psychological variables perceived stress, depressive symptoms, and neuroticism; using formal genetic twin models and molecular genetic methods, i.e. polygenic risk scores (PRS). HCC was measured in 671 adolescents and young adults. These included 115 monozygotic and 183 dizygotic twin-pairs. For 432 subjects PRS scores for plasma cortisol, major depression, and neuroticism were calculated using data from large genome wide association studies. The twin model revealed a heritability for HCC of 72%. No significant phenotypic or genetic correlation was found between HCC and the three psychological variables of interest. PRS did not explain variance in HCC. The present data suggest that HCC is highly heritable. However, the data do not support a strong biological link between HCC and any of the investigated psychological variables.

Twin studies indicate that dyscalculia (or mathematical disability) is caused partly by a genetic component, which is yet to be understood at the molecular level. Recently, a coding variant (rs133885) in the myosin-18B gene was shown to be associated with mathematical abilities with a specific effect among children with dyslexia. This association represents one of the most significant genetic associations reported to date for mathematical abilities and the only one reaching genome-wide statistical significance. We conducted a replication study in different cohorts to assess the effect of rs133885 maths-related measures. The study was conducted primarily using the Avon Longitudinal Study of Parents and Children (ALSPAC), (N=3819). We tested additional cohorts including the York Cohort, the Specific Language Impairment Consortium (SLIC) cohort and the Raine Cohort, and stratified them for a definition of dyslexia whenever possible. We did not observe any associations between rs133885 in myosin-18B and mathematical abilities among individuals with dyslexia or in the general population. Our results suggest that the myosin-18B variant is unlikely to be a main factor contributing to mathematical abilities. We could not replicate the association of the myosin-18B gene with mathematical ability.

Background: Young adulthood is a critical life period for health and health behaviours. Related measurements collected before and after birth, and during childhood and adolescence can provide a life-course analysis of important factors that contribute to health and behaviour in young adulthood. The Western Australian Pregnancy Cohort (Raine) Study has collected a large number of such measurements during the fetal, perinatal, infancy, childhood and adolescence periods and plans to relate them to common health issues and behaviours in young adults, including spinal pain, asthma, sleep disorders, physical activity and sedentary behaviour and, work absenteeism and presenteeism. The aim of this paper is to describe the rationale, design and methods of the 22 year follow-up of the Raine Study cohort. Methods/Design: The Raine Study is a prospective cohort study. Participants still active in the cohort (n = 2,086) were contacted around the time of their 22nd birthday and invited to participate in the 22 year follow-up. Each was asked to complete a questionnaire, attend a research facility for physical assessment and an overnight sleep study, wear activity monitors for a week, and to maintain a sleep and activity diary over this week. The questionnaire was broad and included questions related to sociodemographics, medical history, quality of life, psychological factors, lifestyle factors, spinal pain, respiratory, sleep, activity and work factors. Physical assessments included anthropometry, blood pressure, back muscle endurance, tissue sensitivity, lung function, airway reactivity, allergic status, 3D facial photographs, cognitive function, and overnight polysomnography. Discussion: Describing the prevalence of these health issues and behaviours in young adulthood will enable better recognition of the issues and planning of health care resources. Providing a detailed description of the phenotype of these issues will provide valuable information to help educate health professionals of the needs of young adults. Understanding the life-course risk factors of health issues and behaviours in young adulthood will have important health planning implications, supporting the development of targeted interventions to improve current health status and reduce the onset and development of further ill-health across adulthood.

Summary: The relationships between fat mass and bone mass in young adults are unclear. In 1,183 young Australians, lean body mass had a strong positive relationship with total body bone mass in both genders. Fat mass was a positive predictor of total body bone mass in females, with weaker association in males. Introduction: Body weight and lean body mass are established as major determinants of bone mass, but the relationships between fat mass (including visceral fat) and peak bone mass in young adults are unclear. The aim of this study was to evaluate the associations between bone mass in young adults and three body composition measurements: lean body mass, fat mass and trunk-to-limb fat mass ratio (a surrogate measure of visceral fat). Methods: Study participants were 574 women and 609 men aged 19-22 years from the Raine study. Body composition, total body bone mineral content (TBBMC), bone area and areal bone mineral density (TBBMD) were measured using DXA. Results: In multivariate linear regression models with height, lean body mass, fat mass and trunk-to-limb fat mass ratio as predictor variables, lean mass was uniquely associated with the largest proportion of variance of TBBMC and TBBMD in males (semi-partial R 2 0.275 and 0.345, respectively) and TBBMC in females (semi-partial R 2 0.183). Fat mass was a more important predictor of TBBMC and TBBMD in females (semi-partial R 2 0.126 and 0.039, respectively) than males (semi-partial R 2 0.006 and 0.018, respectively). Trunk-to-limb fat mass ratio had a weak, negative association with TBBMC and bone area in both genders (semi-partial R 2 0.004 to 0.034). Conclusions: Lean body mass has strong positive relationship with total body bone mass in both genders. Fat mass may play a positive role in peak bone mass attainment in women but the association was weaker in men; different fat compartments may have different effects. © 2014 International Osteoporosis Foundation and National Osteoporosis Foundation.

Objective Hypertension in pregnancy and preeclampsia have been linked to poor outcomes in cognitive, mental and psychomotor development; however, few longitudinal studies have researched their effect on offspring motor development, particularly in late childhood and adolescence. The purpose of this study was to determine if maternal hypertensive diseases during pregnancy are a risk factor for compromised motor development at 10, 14, and 17 years. Study design Longitudinal cohort study using data from the Western Australian Pregnancy Cohort Study (Raine). Main outcome measure Offspring (n = 2868) were classified by their maternal blood pressure profiles during pregnancy: normotension (n = 2133), hypertension (n = 626) and preeclampsia (n = 109). Offspring motor development, at 10, 14, and 17 years was measured by the Neuromuscular Developmental Index (NDI) of the McCarron Assessment of Motor Development (MAND). Methods Linear mixed models were used to compare outcomes between pregnancy groups. Results Offspring from pregnancies complicated by preeclampsia had poorer motor outcomes at all ages than offspring from either normotensive mothers (p = 0.001) or those with hypertension (p = 0.002). Conclusion Hypertensive diseases during pregnancy, in particular preeclampsia, have long term and possibly permanent consequences for motor development of offspring. © 2014 International Society for the Study of Hypertension in Pregnancy Published by Elsevier B.V. All rights reserved.

Twin studies suggest that expressive vocabulary at ~24 months is modestly heritable. However, the genes influencing this early linguistic phenotype are unknown. Here we conduct a genome-wide screen and follow-up study of expressive vocabulary in toddlers of European descent from up to four studies of the EArly Genetics and Lifecourse Epidemiology consortium, analysing an early (15-18 months, 'one-word stage', NTotal=8,889) and a later (24-30 months, 'two-word stage', NTotal=10,819) phase of language acquisition. For the early phase, one single-nucleotide polymorphism (rs7642482) at 3p12.3 near ROBO2, encoding a conserved axon-binding receptor, reaches the genome-wide significance level (P=1.3×10-8) in the combined sample. This association links language-related common genetic variation in the general population to a potential autism susceptibility locus and a linkage region for dyslexia, speech-sound disorder and reading. The contribution of common genetic influences is, although modest, supported by genome-wide complex trait analysis (meta-GCTA h15-18-months2=0.13, meta-GCTA h24-30-months2 =0.14) and in concordance with additional twin analysis (5,733 pairs of European descent, h24-months2=0.20).

Objective To evaluate the risk of placenta praevia accreta following primary (first) elective or primary emergency caesarean section in a pregnancy complicated by placenta praevia. Design Retrospective matched case-control study, employing variable matching. Setting Tertiary referral centre between 1993 and 2008. Population Sixty-five cases and 102 controls were used for the analysis from a total of 82 667 births during the study period. Methods Relevant data were abstracted from clinical records. Matching of cases with controls was based on co-existing placenta praevia, number of previous caesarean sections, and age, with one or two controls per case. Results are presented as odds ratios (ORs) with 95% confidence intervals (95% CIs). Main outcome measures Placenta accreta in a pregnancy complicated by placenta praevia following a primary elective or emergency caesarean section, and morbidity associated with pregnancies complicated by placenta accreta. Results Significantly more cases than controls had an elective caesarean section for their primary caesarean delivery (46.2 versus 18.6%; P < 0.001). There were no differences between groups for previous pregnancy loss, uterine surgery, and vaginal delivery, before or after the primary caesarean section. Compared with primary emergency caesarean section, primary elective caesarean section significantly increased the risk of placenta accreta in a subsequent pregnancy in the presence of placenta praevia (OR 3.00; 95% CI 1.47-6.12; P = 0.025). Conclusions Our results suggest that women with a primary elective caesarean section without labour are more likely, compared with those undergoing primary emergency caesarean section with labour, to develop an accreta in a subsequent pregnancy with placenta praevia. © 2013 The Authors BJOG An International Journal of Obstetrics and Gynaecology © 2013 RCOG.

Dysfunctional regulation of the hypothalamic-pituitary-adrenal (HPA) axis has been proposed as an important biological mechanism underlying stress-related diseases; however, a better understanding of the interlinked neuroendocrine events driving the release of cortisol by this stress axis is essential for progress in preventing or halting irreversible development of adverse HPA-function. We aimed to investigate basal HPA-activity in a normal population in late adolescence, the time of life believed to overlap with HPA-axis maturation and establishment of a lasting set point level of HPA function. A total of 1258 participants (mean age 16.6 years) recruited from the Western Australian Pregnancy (Raine) Cohort provided fasting morning blood and saliva samples for basal HPA activity assessment. Irrespective of gender, linear regression modelling identified a positive correlation between the main components of the HPA-cascade of events, ACTH, total cortisol and free cortisol in saliva. Corticosteroid binding globulin (CBG) was inversely associated with free cortisol in saliva, an effect most clearly observed in boys. ACTH levels were lower, but cortisol levels were higher in girls than in boys. Girls may also be exposed to more bioactive cortisol, based on higher average free cortisol measured in saliva at awakening. These relatively higher female free cortisol levels were significantly reduced by oral contraceptive use, eliminating the gender specific difference in salivary cortisol. Free plasma cortisol, calculated from total circulating cortisol and CBG concentrations, was also significantly reduced in girls using oral contraceptives, possibly via an enhancing effect of oral contraceptives on blood CBG content. This study highlights a clear gender difference in HPA activity under non-stressful natural conditions. This finding may be relevant for research into sex-specific stress-related diseases with a typical onset in late adolescence. © 2012 Elsevier Ltd.