Evaluating the genetic effects of sex hormone traits on the development of mental traits: a polygenic score analysis and genome-wide genetic interaction analysis in UK Biobank cohort

Objective: To evaluate the genetic effects of sex hormone traits on the development of mental traits in middle-aged adults. Methods: The SNPs associated with sex hormone traits were derived from a two-stage genome-wide association study (GWAS). Four sex hormone traits were selected in the current study, including sex hormone-binding globulin (SHBG), testosterone, bioavailable testosterone and estradiol. The polygenic risk score (PRS) of sex hormone traits were calculated from individual-level genotype data of the United Kingdom (UK) Biobank cohort. We then used logistic and linear regression models to assess the associations between individual PRS of sex hormone traits and the frequency of alcohol consumption, anxiety, intelligence and so on. Finally, genome-wide genetic interaction study (GWGIS) was performed to detect novel candidate genes interacting with the sex hormone traits on the development of uid intelligence and the frequency of smoking and alcohol consumption by PLINK2.0. Results: We observed positive association between SHBG and the frequency of alcohol consumption (b=0.0101, p=3.84×10 –11 ) in middle-aged males and females. In addition, estradiol was positively associated with the frequency of alcohol consumption (b=0.0128, p=1,96×10 –8 ) in middle-aged males. Moreover, SHBG was associated with the frequency of alcohol consumption (b=0.0067, p=2.60×10 –3 ) in middle-aged females. Finally, GWGIS identied one signicant loci, Tenascin R (TNR) (rs34633780, p=3.45×10 –8 ) interacting with total testosterone for uid intelligence. Conclusion: Our study results support the genetic effects of sex hormone traits on the development of intelligence and the frequency of alcohol consumption in middle-aged adults in UK. the frequency of alcohol consumption and uid intelligence in middle-aged adults. The most signicant interaction effect was observed between total testosterone and TNR for uid intelligence. Our study could provide novel insights into the impact of sex hormone traits on mental traits and highlight the importance of sex specic effects of sex hormone traits on mental traits.


Introduction
Mental disorders are highly prevalent and disabling globally, which lead to heavy burden on the health care system and society [1]. Based on the Global Burden of Disease, Injuries, and Risk Factors Study 2017 (GBD 2017), mental disorders consistently accounted for more than 14% of age-standardized years lived with disability for nearly 30 years [2]. It has been reported that 17.6% of adults suffered from a common mental disorder within the past 12 months and 29.2% across their lifetime [3]. Anxiety, alcohol use disorder and addiction are common mental disorders [3][4][5]. Interestingly, tobacco addiction is the most common co-occurring disorder among persons with serious mental illness [6]. The intelligence is linked to the risk of whole range of mental disorders [7].
Mental traits are multi-factorial, which are the result of multiple genetic and environmental factors that may interact in complicated ways to impact mental traits and disorders susceptibility. Addictions to alcohol and nicotine are heritable disorders [8]. The estimated heritability was 43% for tobacco usage and 19%~29% for alcohol dependence [8]. Addictive disorder has a strong heritable component with an estimated heritability around 30%~70% [9]. Sex differences in brain have known function to generate differences in the control of gonadotropic hormones, reproductive behavior or cognitive functions [10].
The sex differences in the dopamine response in the nucleus accumbens may affect the different vulnerability to addictive disorders in males and females [11]. Sex hormone activity strongly affects the individual's behavior and the constitution of the brain [12].
Over past decade, there are robust evidences to support the relationships between hormones and mental disorders [13][14][15]. For instance, Lenz et al. [13] indicated that exposure to sex hormone in utero and during early development would contribute to the risk of alcohol addiction later in life. Interestingly, they also observed bidirectional relationship between the sex hormone axis and alcohol drinking behavior [13]. The association of hormone with alcohol intake was also supported by experimental animal study [16], which has identi ed that estradiol can stimulate alcohol consumption and aggression in male mice. In addition, it has been found that the change in hormone levels overtime and the ratio of progesterone to estradiol were the strongest hormonal predictors of smoking behavior [15]. Sex differences in anxiety-like behavior partially were affected by aged-related testosterone decline in male rats [17]. Stanikov et al. [18] suggested that higher free testosterone imbalance may mediate depression in overweight premenopausal women. However, limited efforts have been paid to explore the interactions effects between genetic factors and sex hormone traits for mental traits from sex-speci c genetic perspective.
Polygenic risk score (PRS) is a sum of risk alleles, weighted by their effect size estimated from previous published genome-wide association study (GWAS) . Utilizing identi ed susceptibility loci, PRS analysis can evaluate the effects of susceptible loci on disease risks and explore the genetic relationships between complex diseases and traits [9]. PRS has been applied in many studies on neuropsychiatric disorders [19,20]. For instance, Jacqueline et al. [20] conducted PRS analysis to explore the possibility of overlapping genetic factors between smoking and the use of alcohol and cannabis. They found that PRS of cigarettes per day was associated with the number of glasses alcohol per week and cannabis initiation [20]. In another study, researchers found that the PRS of nicotine metabolism can predict nicotine metabolism biomarkers [21]. Additionally, Belsky et al. [22] observed that individuals with higher PRS were more likely to persist longer in smoking heavily and develop nicotine dependence more frequently.
Recently, Katherine et al. [23] conducted a two-stage GWAS in 425,097 United Kingdom (UK) biobank study participants and identi ed 2,571 genetic variant-sex hormone associations. Using the PRS of sex hormone traits as instrumental variables, we can calculate the PRS of sex hormone traits and explore the correlations between PRS of sex hormone traits and mental traits.
GWAS has great power to identify susceptibility genetic loci associated with mental disorders [24,25]. However, the signi cant loci identi ed by GWAS are usually limited and functionally independent. Genetic effects are different between individuals due to gene-environmental (G×E) interactions, which resulted from individuals responding differently to environmental stimuli depending on their genotype [26]. Identifying G×E interactions would improve risk-assessment for complex diseases and reveal underlying biological pathways [26]. The genome-wide genetic interaction study (GWGIS) can estimate the effect of G×E interactions [27]. GWGIS can investigate the genetic interaction effect at a genome-wide scale, which can improve the ability of detecting genotype-phenotype associations missed in GWAS [27,28].
In this study, we rst calculated the PRSs of sex hormone traits in middle-aged adults in UK Biobank.
Logistic and linear regression analyses were then performed to detect the associations between individual PRSs value of sex hormone traits and the phenotypic data of mental traits in UK Biobank.
Finally, GWGIS was conducted to explore novel candidate genes interacting with sex hormone traits on the development of uid intelligence and the frequency of alcohol consumption and smoking.

UK Biobank samples and mental phenotypes
The phenotypic and genotypic data of this study were derived from UK Biobank health resource under UK Biobank application 46478, which was a population-based prospective cohort study. Between 2006 to 2010, UK Biobank collected 502,656 participants aged 40 and 69 at recruitment. UK Biobank cohort has collected a rich variety of phenotypic, health-related information on each participant, including physical and biological measurements, lifestyle indicators, imaging of the body and brain and genome-wide genotyping. Longitudinal follow-up for a wide range of health-related information are provided by linking health and medical records.
Several potential measures of smoking behavior were selected to de ne the phenotype of ever smoking. The UK Biobank data eld of 20432 was described as ongoing behavioural or miscellanous addiction.
Anxiety and depression were de ned according to the previous study [29], which were based on the general anxiety disorder (GAD-7) and Patient Health Questionnaire (PHQ-9) [30,31]. Fluid intelligence score was described as a simple unweighted sum of the number of correct answers given to the 13 uid intelligence questions. The maximum number of reported past or current cigarettes (or pipes/cigars) consumed per day was used to de ne the frequency of smoking (UK Biobank data elds: 20116, 2887 and 3456). In addition, the frequency of alcohol consumption (UK Biobank data eld: 20117) was de ned as the sum of all alcoholic beverages per week. Ethical approval of UK Biobank study was granted by the National Health Service National Research Ethics Service (reference 11/NW/0382). The detailed de nition of phenotypes are shown in Additional le 1.

UK Biobank genotyping, imputation and quality control
A total of 488,377 middle-aged adults have genome-wide genotype data, which were assayed by two similar genotyping array. DNA was extracted from stored blood samples when participants visited to a UK Biobank assessment Centre. Genotyping was carried out by Affymetrix UK BiLEVE Axiom Array or the Affymetrix UK Biobank Axiom arrays (Santa Clara, CA, USA), which shared 95% of marker content. Imputation was conducted by IMPUTE4 (https://jmarchini.org/software/) to carry out in chunks of approximately 50,000 imputed markers with a 250 kb buffer region. Routine quality checks were carried throughout the process, including sample retrieval, DNA extraction and genotype calling. Statistical tests were performed to identify poor quality markers by checking for consistency of genotype calling across experimental factors, including batch effects, plate effects, departures from Hardy-Weinberg equilibrium (HWE), sex effects, array effects, and discordance across control replicates. Based on self-reported ethnicity (UK Biobank data eld: 21000), the individuals were restricted to only "white British". Finally, the participants who reported inconsistencies between self-reported gender or genetic gender, who were genotyped but not imputed, and who withdraw their consents, were excluded in the current study. Detailed description of array design, genotyping and quality control procedures can be found in the previous studies [32,33].

GWAS data of sex hormone traits
The SNPs associated with sex hormone traits were derived from a published GWAS [23]. Brie y, the published GWAS analyzed four sex hormone traits, including sex hormone-binding globulin (SHBG), testosterone, bioavailable testosterone and estradiol. Association test was conducted to account for cryptic population structure and relatedness by linear mixed models implemented in BOLT-LMM (v2.3.2). Genotypic data was derived from the 'v3' release of UK Biobank [32], which contained the full set of Haplotype Reference Consortium (HRC) and 1000 Genomes imputed variants. The SNPs with signi cant threshold of p value < 5×10 -8 were selected to calculate PRSs. Detailed description of sample characteristics, array design, quality control and statistical analysis can be found in the previous study [23].

PRS of sex hormone traits
Using the genotype data of UK Biobank cohort, PRS calculation was performed by using the PLINK's "-score" command [34]. Brie y, PRS denotes the PRS of the sex hormone traits for the ith subjects, de ned as , where n (n = 1, 2, 3, …, t) and i (i=1, 2, 3, …, k) denote the number of genetic markers and the sample size, respectively. β n is the effect parameter of risk allele of the nth signi cant SNP related to sex hormone traits, which obtained from the published study. is the dosage (0 to 2) of the risk allele of the nth SNP for the ith subject. In addition, we have excluded Linkage Disequilibrium (LD) when calculating PRSs by using the command "--indep-pairwise" implemented in PLINK with parameters window size (500kb), step size (5 variant ct) and r 2 < 0.5.

Statistical analysis
Four serum sex hormone traits, including SHBG, testosterone, bioavailable testosterone and estradiol, were analyzed both within and across sexes, with the exception of estradiol where analyses were performed only in men. Logistic regression model was performed to assess the associations between individual PRSs of sex hormone traits and ever smoking and ongoing behavioural or miscellanous addiction, respectively. 6 statistical tests (3 serum sex hormone traits × 2 mental traits) were analyzed by logistic regression in total samples and females. The signi cant correlation was identi ed at p value < 8.33 × 10 -3 (0.05/6) after Bonferroni correction. And 8 statistical tests (4 serum sex hormone traits × 2 mental traits) were analyzed by logistic regression in males. The signi cant correlation was identi ed at p value < 6.25 × 10 -3 (0.05/8) after Bonferroni correction. Correspondingly, linear regression model was conducted to evaluate the correlations between individual PRSs of sex hormone traits and anxiety score, depression score, uid intelligence score, and the frequency of alcohol consumption and smoking, respectively. 15 statistical tests (3 serum sex hormone traits × 5 mental traits) were analyzed by linear regression in total samples and females. The signi cant correlation was identi ed at p value < 3.33 × 10 -3 (0.05/15) after Bonferroni correction. And 20 statistical tests (4 serum sex hormone traits × 5 mental traits) were analyzed by linear regression in males. The signi cant correlation was identi ed at p value < 2.50 × 10 -3 (0.05/20) after Bonferroni correction. The regression analyses were conducted by R software (version 3.5.3). Sex, age, and 10 principle components of population structure were used as covariates in the regression model.

Genome-wide genetic interaction study
Based on the result of regression model, GWGIS was performed to assess the interaction effects between genetic factors and sex hormone traits for uid intelligence, and the frequency of smoking per day and alcohol consumption per week in UK Biobank cohort. The GWGIS was conducted by PLINK2.0 [34,35]. Letting D is the disease outcome variable, the penetrance models form is described as the following: where G is genetic factors and E is the environmental factors [36]. In this study, the outcome variables were uid intelligence score and the frequency of smoking per day and alcohol consumption per week. And the instrumental variables were the PRS of serum sex hormone traits. The HWE p value < 0.001 or minor allele frequencies (MAFs) < 0.01 or the SNPs with low call rates (< 0.90) were excluded in the current study for quality control. To avoid the in uence of population strati cation, cryptic relatedness and null model misspeci cations on our results, we calculated the in ation factor of GWGIS. Signi cant interactions were identi ed at p value < 5.0 × 10 -8 in this study. Rectangular Manhattan plot was generated using the "CMplot" R script (httcps://github.com/YinLiLin/R-CMplot).

Result
A total of 9 signi cant associations were identi ed in this study, and the general characteristics of the subjects are presented in Table 1.

Associations between PRSs of sex hormone traits and mental traits in total middle-aged samples
Three sex hormone traits were analyzed in total middle-aged samples, including SHBG, total testosterone, and bioavailable testosterone. We observed positive associations between SHBG and the frequency of alcohol consumption (b=0.0101, p=3.84×10 -11 ) and ever smoking (b=0.0084, p=4.27×10 -3 ). In addition, total testosterone was positively associated with the frequency of alcohol consumption (b=0.0067, p=1.59×10 -5 ). No signi cant association was found between bioavailable testosterone and mental traits in total middle-aged samples. The basic characteristics of study subjects and detailed information are presented in Additional le 2.

Associations between PRSs of sex hormone traits and mental traits in middle-aged males
Brie y, four sex hormone traits were analyzed in middle-aged males, including SHBG, total testosterone, bioavailable testosterone and estradiol. Estradiol was positively associated with the frequency of alcohol consumption (b=0.0128, p=1,96×10 -8 ). In addition, total testosterone was associated with the frequency of smoking (b=-0.0108, p=2.07×10 -6 ). SHBG was positively associated with the frequency of alcohol consumption (b=0.0090, p=8.18×10 -5 ). No signi cant association was found between bioavailable testosterone and mental traits in middle-aged males. The basic characteristics of study subjects and detailed information are presented in Additional le 3.

Associations between PRSs of sex hormone traits and mental traits in middle-aged females
Three sex hormone traits were analyzed in females, including SHBG, total testosterone, and bioavailable testosterone. SHBG was associated with the frequency of alcohol consumption (b=0.0067, p=2.60×10 -3 ). In addition, total testosterone was associated with the frequency of alcohol consumption (b=0.0102, p=4.55×10 -6 ). We also observed negative association between bioavailable testosterone and uid intelligence (b=-0.0136, p=5.74×10 -5 ). The basic characteristics of study subjects and detailed information are presented in Additional le 4.

Genome-wide genetic interaction analysis results
The Rectangular Manhattan plot is shown in Figure 1. The in ation factors was between 0.99~1.04 in this study, which indicates a low possibility of false-positive association resulting from population strati cation and null model misspeci cations. GWGIS identi ed one signi cant loci, Tenascin R (TNR) (rs34633780, p=3.45×10 -8 ) interacting with total testosterone for uid intelligence in middle-aged adults.

Discussion
Sex hormone supplementation has commonly effects on metabolic traits, sexual function and bone health [23]. Epidemiological study has indicated strong correlations between sex hormone supplements and health conditions [15]. However, the interactions effects between genetic factors and sex hormone traits for mental traits remain largely unknown now. Considering the effects of sex hormone traits on individual's behavior and the constitution of the brain exhibiting fundamental differences between males and females, our study focuses on the potential impact of sex hormone traits on mental traits from sex-speci c genetic perspective. We observed that sex hormone traits were associated with ever smoking and the frequency of alcohol consumption in middle-aged adults in UK. GWGIS further identi ed that the interaction of total testosterone with rs34633780 that mapped to the TNR gene can modulate the phenotype of uid intelligence, possibly through in uencing cognitive function.
TNR, a member of the tenascin family of neural extracellular matrix glycoproteins, primarily expressed in the central nervous system. This protein affects neural cell adhesion, neurite outgrowth and modulation of sodium channel function. It has been implicated that TNR was known to function in many neurological diseases [37,38], such as attention de cit hyperactivity disorder (ADHD) [38] and neurodegenerative disorders [37]. David et al. [39] reported that an individual with a homozygous deletion of the TNR gene was associated with intellectual disability, supporting the role of TNR in brain development and cognition in humans. Moreover, Anna et al. [40] reported a case with intellectual disability had a 6.14 Mb duplication at 1q25.1-q25.2 by utilizing array comparative genomic hybridization (array CGH). Interestingly, the TNR gene located within this region (1q25.1-q25.2). Additionally, genome-wide association analysis found a novel association of ADHD with TNR gene [38]. In another study, researchers found that TNR de ciency would cause an early onset and nonprogressive neurodevelopmental disorder [41]. To the best of our knowledge, the study involving biological processes of TNR in uid intelligence is still less. Although indirect, the "TNR-fluid intelligence" of GWGIS nding, combined with ndings reported for other psychiatric disorders [37,38,41], supporting the evidence of TNR in the aetiology of psychiatric conditions. Another signi cant nding of this study is the disclosure of the association between SHBG and the frequency of alcohol consumption in middle-aged males and females. SHBG, the major and speci c binding protein for testosterone and estradiol, is known to regulate the bioavailability of sex steroids. Meanwhile, SHBG can assess bioavailable testosterone level [42]. Lee et al. [43] have reported that serum SHBG level is an independent predictive factor for extraprostatic extension of tumor in prostate cancer patients. Higher SHBG concentrations were observed in the premenopausal women who consumed alcohol [44], which consistent with our result in females. Interestingly, Markus et al. [45] demonstrated that serum SHBG levels can be regulated by metabolic factors, including alcohol consumption and several drugs. Another study found that high concentrations of SHBG were consistently related to type II alcoholism [46]. However, Shiels et al. observed the paradoxical results. They found that decreased levels of SHBG were associated with increased alcohol consumption in United States men [47]. The inconsistent result of the association of alcohol with SHBG may be due to the samples size and different genetic background. The biological mechanism explanation for the correlation of SHBG with alcohol consumption is still unclear. No association was observed between smoking status and SHBG in this study, which was consistent with previous study nding [47]. In addition, total testosterone was positively correlated with the frequency of alcohol consumption in middle-aged adults in UK. It has been identi ed that the frequency of alcohol consumption was positively associated with the concentrations of total testosterone among adult men [47]. High concentrations of total testosterone were shown to be associated with type II alcoholism [46]. Martin et al. [48] observed that females with higher levels of testosterone were more likely to be using alcohol currently. Tina Kold et al. [49] suggested that alcohol consumption was associated with changes in testosterone levels among young men.
Estradiol played an important role in the establishment of sex differences in brain structure and function, which may act primary target for the investigation of sex-related differences in alcohol effects [50]. Our study results also supported that estradiol was positively associated with the frequency of alcohol consumption among middle-aged men. Population-based study in adolescent males found higher salivary estradiol level was associated with earlier onset and higher quantity of alcohol use [51]. Experimental animal study [16] suggested that estradiol can in uence voluntary alcohol consumption and alcohol related behaviors in male mice, including aggression and depression. Most interestingly, these effects are strongly gender dependence [16]. Besides, previous study reported that alcohol consumption can affect gonadal hormone. For instance, Sarkola et al. [52] found that estradiol levels were increased after intake of alcohol among subjects who used oral contraceptives. In addition, we observed that total testosterone was correlated with the frequency of smoking in middle-aged males.
Ponholzer et al. observed that nicotine consumption was associated with serum levels of testosterone or free testosterone [53]. Another study demonstrated that current smokers of ve or more cigarettes/day showed signi cantly higher levels of testosterone [54]. Furthermore, signi cant increases in serum testosterone level were observed in smokers group [55].
It is important to emphasize that our study has two limitations. First, all the samples were collected from UK Biobank cohort, aged between 40 and 69 at recruitment. Therefore, our ndings should be carefully interpreted when applied to others ages and different genetic background populations. Second, the SNPs associated with sex hormone traits were derived from previous GWAS. The accuracy of our regression analyses may be in uenced by the power of previous GWAS on the sex hormone traits. Further replication studies with other genetic background individuals and experimental studies are required to verify the results of this study.
In conclusion, the standardized collection of genotype and sex hormone supplementation data in UK Biobank give us an opportunity to access the interaction effect between sex hormone traits and genetic factors for mental traits. We observed correlations between sex hormone traits and the frequency of alcohol consumption and uid intelligence in middle-aged adults. The most signi cant interaction effect was observed between total testosterone and TNR for uid intelligence. Our study could provide novel insights into the impact of sex hormone traits on mental traits and highlight the importance of sex speci c effects of sex hormone traits on mental traits.

Declarations
Ethics approval and consent to participate Ethical approval of UK Biobank study was granted by the National Health Service National Research Ethics Service (reference 11/NW/0382).

Consent for publication
Not applicable.

Availability of data and material
The UKB data are available through the UK Biobank Access Management System (https://www.ukbiobank.ac.uk/). We will return the derived data elds following UKB policy; in due course, they will be available through the UK Biobank Access Management System.

Competing interests
There are no competing interests to declare.

Funding
There are no funding in this study.

Authors' contributions
Xiao Liang contributes to the design of the work, analysis of data and draft the paper; ShiQiang Cheng contributes to the part of data analysis and revises the paper; Jing Ye contributes to the part of data analysis; XiaoMeng Chu corrects the grammar issue and revises the paper; Yan Wen contributes to the part of data analysis; Li Liu revises the paper; Xin Qi drafts the Note Bioavailable testosterone (Bioavailable T); sex hormone-binding globulin (SHBG); Total testosterone (Total T).