Abstract

Background

Steroidogenic Factor 1/Nuclear Receptor Subfamily 5 Group A Member 1 (SF-1/NR5A1) is critical for the development and function of sex organs, influencing steroidogenesis and reproduction. While rare deleterious NR5A1/SF-1 variants have been identified in individuals with various differences of sex development (DSD), primary ovarian insufficiency, and infertility, their impact on the general population remains unclear.

Methods

We analyzed health records and exome sequencing data from up to 420 162 individuals (227 858 women) from the UK Biobank study to assess the impact of rare (frequency < 0.1%) predicted deleterious NR5A1/SF-1 variants on age at menopause and 26 other traits.

Results

No carriers of rare protein truncating variants in NR5A1/SF-1 were identified. We found that the previously reported association of rare deleterious missense NR5A1/SF-1 variants with earlier age at menopause is driven by variants in the DNA binding domain (DBD) and ligand binding domain (LBD) (combined test: beta = −2.36 years/allele, [95% CI: 3.21, −1.51], N = 107 carriers, P = 4.6 × 10−8). Carriers also had a higher risk of adult obesity (OR = 1.061, [95% CI: 1.003, 1.104], N = 344, P = .015), particularly among women (OR = 1.095 [95% CI: 1.034, 1.163, P = 3.87 × 10−3], N = 176), but not men (OR = 1.019, [95% CI: 0.955, 1.088], P = .57, N = 168).

Conclusion

Deleterious missense variants in the DBD and LBD likely disrupt NR5A1/SF-1 function. This study broadens the relevance of deleterious NR5A1/SF-1 variants beyond rare DSDs, suggesting the need for extended phenotyping and monitoring of affected individuals.

Rare deleterious mutations in Steroidogenic Factor 1/Nuclear Receptor Subfamily 5 Group A Member 1 (SF-1/NR5A1) are prevalent among individuals with differences of sex development (DSD) presenting with a wide spectrum of phenotype ranging from mild to severe. They occur mostly in heterozygous but also in homozygous state [1-4]. But, NR5A1/SF-1 variants can also be identified in individuals without DSD, including healthy carriers, infertile men, or women with primary ovarian insufficiency (POI) [4-8]. Broader health impacts of NR5A1/SF-1 variations are possible, given the key role of this gene in adrenal and reproductive development and function [9]. So far, our knowledge regarding the effects of NR5A1/SF-1 variants beyond DSD is limited.

Early studies on SF-1 knockout (KO) mice hinted at the involvement of SF-1 in energy balance, as evidenced by the compromised structure of the ventromedial hypothalamic nucleus (VMH) [10, 11], a recognized regulator of energy homeostasis and body weight [12-15]. SF-1 neurons, which are a subset of neurons primarily located in the VMH, are regulated by various factors such as nutrients (eg, glucose) and hormones (eg, leptin and insulin) [16]. SF-1 KO mice exhibited late-onset obesity due to decreased energy expenditure rather than hyperphagia, establishing them as potential model for hypothalamic late-onset obesity [17].

Subsequent investigations [18, 19] focused on the role of leptin in SF-1 neurons and obesity. Conditional KO mice lacking leptin receptors (LepRs) in SF-1 neurons, displayed obesity, particularly under a high-fat diet condition. The idea of a distinct functional subsets of leptin-activated and leptin-inhibited SF-1 neurons has been proposed [20, 21]. Moreover, mice with genetic modifications affecting various elements in SF-1 neurons displayed a modified metabolic phenotype, primarily evident under high-fat diet conditions, rather than a standardized nutritional balanced diet [16]. These observations suggested that SF-1 neurons may play a role in metabolic adaptation, particularly during the initial phases of obesity.

Optogenetic and chemogenetic studies have further elucidated the functions of SF-1 neurons in energy balance and glucose regulation. Optogenetic stimulation showed that SF-1 neurons exert differential effects depending on the frequency of activation [22]. Low-frequency activation of SF-1 neurons in mice resulted in lower food intake and less interest in food, shown by a decrease in time the mice spent near their food. Similarly, chemogenetic activation of SF-1 neurons led to a reduction in food intake and an increase in energy expenditure and fat oxidation in fasted mice [23], while inhibition of SF-1 neurons led to increase in cumulative food intake in mice with ad libitum feeding [24]. Conversely, in another study, optogenetic activation of SF-1 neurons resulted in diabetes-range hyperglycemia [25]. Conflicting results across these studies support the idea that glucose regulation by SF-1 neurons may be mediated by distinct subsets of neurons, leading to either an increase in insulin sensitivity or hyperglycemia.

While numerous studies in mouse models establish the role of SF-1 in metabolic regulation, limited data are available for DSD individuals with rare deleterious NR5A1/SF-1 variants and their association to obesity or overweight [26-30]. In one study in the Han Chinese population, the specific NR5A1/SF-1 variant p.Gly146Ala was associated with type 2 diabetes [31]. Data on metabolic health outcomes in individuals with rare NR5A1/SF-1 variants are limited due to the rarity of the disorder and the fact that, so far, the clinical focus was primarily on DSD (eg, pubertal and fertility issues) [4]. Thus, the younger age of studied individuals may lead to the oversight and under-reporting of late-onset obesity or other metabolic health issues in these individuals.

Apart from SF-1 involvement in metabolism, spleen anomalies (eg, spleen hypoplasia) are reported in a significant number of individuals with rare NR5A1/SF-1 variants [4, 30, 32, 33], supported also by mouse model data [34].

These collective findings suggest a potential role of SF-1 in broader aspects of human health beyond reproduction. To address this question, we analyzed multiple sources of health record data and exome sequencing data of up to 420 162 individuals from the population-based UK Biobank study. We investigated the impact of rare (frequency < 0.1%) predicted deleterious NR5A1/SF-1 variants on potentially relevant traits with reproduction, puberty timing, metabolism, and spleen function.

Methods

UK Biobank Data Processing and Quality Control

To conduct rare variant burden analyses, we obtained whole exome sequencing (WES) data for 454 787 individuals provided by the UK Biobank (UKBB) study [35]. Individuals were excluded from subsequent analyses if they exhibited excessive heterozygosity, had autosomal variant missingness on genotyping arrays of ≥5%, or were not encompassed within the subset of phased samples as defined by Bycroft et al [36]. Additionally, we excluded individuals who did not have a broadly European genetic ancestry, resulting in a final WES dataset of 420 162 individuals. Additional exclusions occurred during the analysis of specific phenotypes due to incomplete data.

We used the UKBB Research Analysis Platform (RAP; https://ukbiobank.dnanexus.com/) for variant quality control and annotation. The RAP is a cloud-based computing environment and central data repository for UKBB WES and phenotypic data. Customized applets designed for the RAP facilitated additional quality control on exome sequencing data [37], expanding on the procedures outlined in Backman et al [35]. Using provided population-level Variant Call Format (VCF) files, we first split and left-corrected multi-allelic variants into separate alleles using “bcftools norm.” Next, we performed genotype-level filtering using “bcftools filter” separately for single nucleotide variants (SNVs) and insertions/deletions (InDels) using a missingness-based approach. With this approach, SNV genotypes with depth < 7 and genotype quality < 20 or InDel genotypes with a depth < 10 and genotype quality < 20 were set to missing (ie, ./.). We further tested for an expected alternate allele contribution of 50% for heterozygous SNVs using a binomial test; SNV genotypes with a binomial test P value ≤ 1 × 10−3 (indicative of a genotype calling error or somatic mutation) were set to missing. Following genotype-level filtering we recalculated the proportion of individuals with a missing genotype for each variant and filtered out all variants with a missingness value > 50%.

Following that, we annotated variants by predicting their functional consequences using ENSEMBL Variant Effect Predictor (VEP) v104 [38]. We applied the “–everything” flag and activating plugins for REVEL [39], CADD [40], and LOFTEE [41]. For each variant, we assigned priority to a specific ENSEMBL transcript based on whether the annotated transcript was protein-coding, MANE select v1.0 [42], or the VEP Canonical transcript, respectively. The consequence for each individual variant was determined by prioritizing severity according to VEP's definitions. After the annotation process, we consolidated stop gained, frameshift, splice acceptor, and splice donor variants into a unified category known as protein truncating variants (PTVs). The consequences for missense and synonymous variants were determined following the definitions provided by VEP. Subsequently, only autosomal or chrX variants within ENSEMBL protein-coding transcripts and within transcripts covered by the UKBB WES assay were retained for subsequent burden testing.

Exome-Wide Association Analyses in the UK Biobank

Here, we used the previously reported impact of rare predicted deleterious variants in the NR5A1/SF-1 gene on age at natural menopause (ANM) in this cohort [43], to optimize our definition of deleterious variant class for association testing. We then explored the effect of high-confidence PTVs and missense variants with CADD score ≥25 on other potentially relevant traits in UKBB: age at menarche in women (AAM), age at voice breaking in men (AVB), birth weight, blood pressure (diastolic and systolic), body fat percentage, body mass index (BMI), childlessness, circulating IGF-1 (insulin-like growth factor 1) levels, glucose levels, glycated hemoglobin (HbA1c), high-density lipoprotein (HDL) cholesterol, height, low-density lipoprotein (LDL) cholesterol, number of live births, obesity (BMI >30 kg/m2), red blood cell count, relative adiposity in childhood, SHBG (sex hormone binding globulin) levels, testosterone levels, triglycerides, type 2 diabetes, and waist hip ratio (Supplementary Table S1) [44].

To conduct rare variant burden tests, we used a custom implementation of BOLT-LMM v2.3.6 [45] for the RAP. BOLT-LMM requires 2 main inputs: (i) genotypes at variants with a minor allele count > 100 obtained from genotyping arrays to construct a null linear mixed effects model; and (ii) a larger set of variants collapsed on ENSEMBL transcript to perform association tests. For the former, we extracted genotyping data available on the RAP, limiting it to the same set of individuals used for rare variant association tests. For the latter, since BOLT-LMM expects imputed genotyping data rather than per-gene carrier status, we created dummy genotype files where each variant represented one gene, and individuals with a qualifying variant within that gene were coded as heterozygous, irrespective of the number of variants they carried in that gene.

We generated dummy genotype files for 3 overlapping variant classes: (i) high-confidence PTVs defined by LOFTEE; (ii) missense variants with CADD ≥ 25; and (iii) damaging variants, which included both high-confidence PTVs and missense variants with CADD ≥ 25. For each phenotype tested, BOLT-LMM was executed with default parameters, with the inclusion of the “lmmInfOnly” flag. To obtain association statistics for individual markers, we provided all 26 657 229 markers, regardless of their filtering status, as input to BOLT-LMM. All tested phenotypes were treated as continuous traits or binary traits, adjusted for age, age squared (age2), sex, the first 10 genetic principal components as calculated in Bycroft et al [36], and WES batch as a categorical covariate (50k, 200k, or 450k).

Phenotype Derivation

For the phenotype derivation and processing refer to Supplementary Table S1 [44].

SF-1 Protein Domain-Specific Analyses

To explore how rare deleterious missense variants disrupt SF-1 function, we next categorized variants by protein domain. Briefly, we included predicted deleterious missense variants with CADD score ≥ 25 with MAF < 0.1% which map to either of the 2 main functional domains of SF-1 protein (DNA binding domain [DBD] and ligand binding domain [LBD]), as defined by UniProt. We implemented a generalized linear model (GLM) using the Python package “statsmodels” [46]. For the domain-level burden test, the number of variant alleles across the 2 domains was summed up into a single score under a simple additive model. This score was used as a predictor of the traits in a three-step regression. All traits assessed with BOLT-LMM were included in the GLM burden test, and in addition, we expanded the analysis to encompass additional traits, including spleen volume, platelet count, and reticulocyte count (Supplementary Table S1) [44]. For all association tests we corrected for age, age2, the first 10 genetic principal components provided by Bycroft et al [36], and study participant WES batch as a categorical covariate.

All analyses and data manipulations were conducted in R (v4.1.2) and in the UKB Research Analysis Platform (RAP; https://ukbiobank.dnanexus.com/).

Results

Rare Deleterious NR5A1/SF-1 Variants and Earlier Age at Natural Menopause

Using WES data from 227 858 women of European genetic ancestry, we confirmed that individuals carrying rare deleterious missense variants (with CADD score ≥ 25) in the NR5A1/SF-1 gene had earlier ANM (beta = −1.76 years, [−2.50, −1.01], P = 3.80 × 10−6, N = 131) (Fig. 1, Supplementary Table S2) [44]. Notably, we did not identify any carrier of rare PTVs alleles in NR5A1/SF-1.

Exome associations between NR5A1/SF-1 gene with age at natural menopause in the UKBB. Included were NR5A1/SF-1 variants with MAF smaller than 0.1% and annotated as missense variants with a high CADD score ≥ 25, distributed across exons of NR5A1 gene. Each variant is depicted as a separate line, extending to its association P value (−log10), with the direction indicating the effect on age at natural menopause in carriers of the alternate allele. The point size corresponds to the respective variant's allele count, as indicated in the figure legend.
Figure 1.

Exome associations between NR5A1/SF-1 gene with age at natural menopause in the UKBB. Included were NR5A1/SF-1 variants with MAF smaller than 0.1% and annotated as missense variants with a high CADD score ≥ 25, distributed across exons of NR5A1 gene. Each variant is depicted as a separate line, extending to its association P value (−log10), with the direction indicating the effect on age at natural menopause in carriers of the alternate allele. The point size corresponds to the respective variant's allele count, as indicated in the figure legend.

Trait Associations Beyond Reproduction

Using the same variant mask of rare deleterious missense variants with CADD score ≥ 25 (Supplementary Table S3) [44] to explore associations with additional traits (Supplementary Table S1) [44], we observed nominal associations of the NR5A1/SF-1 gene with obesity (odds ratio [OR] = 1.041, [1, 1.078], P = .046, N = 472) (Fig. 2), HbA1c levels (beta = .54, [0.11, 0.97], P = .013, N = 451), HDL (beta = .04, [0.012, 0.073], P = 6.5 × 10−3, N = 395), circulating IGF-1 levels (beta = −.48, [−0.95, 0], P = .05, N = 415) and systolic blood pressure (beta = .46, [−1.41, 2.32], P = .03, N = 442) (Supplementary Table S2) [44]. Surprisingly, we found no association between this variant mask with other expected reproductive traits, eg, puberty timing in either sex (Supplementary Table S2) [44].

Exome associations between NR5A1/SF-1 gene with obesity in the UKBB. Included were NR5A1/SF-1 variants with MAF smaller than 0.1% and annotated as missense variants with a high CADD score ≥ 25, distributed across exons of NR5A1 gene. Each variant is depicted as a separate line, extending to its association P value (−log10), with the direction indicating the effect on obesity in carriers of the alternate allele. The point size corresponds to the respective variant's allele count, as indicated in the figure legend.
Figure 2.

Exome associations between NR5A1/SF-1 gene with obesity in the UKBB. Included were NR5A1/SF-1 variants with MAF smaller than 0.1% and annotated as missense variants with a high CADD score ≥ 25, distributed across exons of NR5A1 gene. Each variant is depicted as a separate line, extending to its association P value (−log10), with the direction indicating the effect on obesity in carriers of the alternate allele. The point size corresponds to the respective variant's allele count, as indicated in the figure legend.

Individual Variant-Level Analyses

Next, we conducted individual rare deleterious NR5A1/SF-1 variant-level association tests with ANM. We found 7 variants that showed nominally significant individual associations (P ≤ .05) with earlier ANM. Notably, one of these variants (c.43G>A; p.Val15Met, Allele count (AC) = 1) has previously been reported in rare cases of 46,XY sex reversal [47], and premature ovarian insufficiency [48] (Supplementary Table S4) [44]. A second variant (c.721C>T; Arg241Trp, AC = 2) has been also reported in a rare case of 46,XY sex reversal [49], (Supplementary Table S4) [44]. Lastly, a third variant (c.1063G>A; p.Val355Met, AC = 135) has been reported in a boy with bilateral anorchia, micropenis, and progressive testicular degeneration. His mother, who also carried this variant, had 2 spontaneous miscarriages, and had left ovariectomy and removal of the corresponding fallopian tube due to ovarian cysts at the age of 22 (Supplementary Table S4) [44, 50].

Rare Deleterious Missense Variants Within SF-1 Protein Domains

To explore how rare deleterious missense NR5A1/SF-1 variants associated with ANM disrupt its protein function, we categorized them by the SF-1 protein domains. Only 3 variants (n = 24 carriers) were found in the hinge region of the SF-1 protein, and none showed an association with ANM (all P > .05). Thus, we focused only on variants in the DBD and LBD.

Using GLM domain-level burden testing, we found that carriers of rare deleterious missense NR5A1/SF-1 variants in each domain had earlier ANM (DBD: beta = −4.03 years, [−5.71, −2.34], P = 2.80 × 10−6, N = 27; LBD: beta = −1.80, [−2.78, −0.82], P = 3.18 × 10−4, N = 80). We also tested if considering both DBD and LBD simultaneously increased our strength of association with ANM and found that the association between NR5A1/SF-1 and ANM was strengthened (beta = −2.36, [−3.20, −1.51], P = 4.6 × 10−8, N = 107). This combined analysis demonstrates a more robust and significant relationship between NR5A1/SF-1 variants and ANM compared to the domain-agnostic test, underscoring the importance of domain-specific effects in influencing NR5A1 function and menopausal timing.

We used this “domains” mask of grouped variants in DBD and LBD to search for NR5A1/SF-1 gene-level associations with other traits beyond ANM, using GLM burden testing. Associations with domain-unrestricted variants are shown for comparison in Supplementary Table S5 [44]. Our analysis revealed that rare deleterious variants in the DBD and LBD of the NR5A1/SF-1 gene were associated with a higher risk of obesity (sex combined: OR = 1.061, [1.01, 1.106], P = .015, N = 344), particularly in women (OR = 1.057, [1.034, 1.163], P = 3.87 × 10−3, N = 176), but not in men (OR = 1.019, [0.955, 1.088], P = .57, N = 168) (Supplementary Table S5) [44]. Notably, 32.3% of women who carried rare deleterious missense NR5A1/SF-1 variants located in the DBD and LBD had obesity, compared to 23% of noncarriers. Additionally, we found a positive association of these variants with higher platelet count (beta = 8.55 × 109 cells/L, [2.36 × 109, 14.73 × 109], P = 6.8 × 10−3, N = 334) in both sexes (Supplementary Table S5) [44].

Discussion

Here we present the results of an exome-wide association study assessing the collective impact of rare (frequency < 0.1%), predicted deleterious NR5A1/SF-1 variants in individuals from the population-based UKBB study. We aimed to understand the potential associations of these variants with traits putatively linked to the function of NR5A1; a gene known for its pivotal role in reproductive development and function, yet not thoroughly characterized in humans for effects on other health outcomes such as metabolism. A similar analysis reported an association of rare NR5A1/SF-1 variants with earlier age at menopause in women in the same UKBB cohort [43]. In the current study, we found that this association is driven by deleterious variants that specifically code the DBD and LBD of the protein. Through further protein domain association testing, we identified that women who carry rare deleterious missense NR5A1/SF-1 variants located in the DBD and LBD also have an elevated risk of obesity.

These consistent findings on earlier ANM align with the known role of rare deleterious NR5A1/SF-1 variants in POI in women [2, 5, 7, 51]. SF-1 regulates crucial genes of ovarian development, steroidogenesis, growth and maturation of follicles, the disruption of which can lead to ovarian dysfunction [5]. However, incomplete penetrance, together with the suggested role of oligogenic contribution [4, 52, 53] to the phenotype of individuals with NR5A1/SF-1 variants, underscores the complexity of genotype-phenotype correlations. Further investigations are needed to identify potential modifying factors that contribute to the variable clinical outcomes associated with NR5A1/SF-1 variants in the context of ovarian function. This understanding will be crucial for refining risk assessment and facilitating more accurate genetic counseling for affected individuals and their families.

The observed association of deleterious NR5A1/SF-1 variants with obesity supports the proposed involvement of SF-1 in regulating energy metabolism, beyond its role in reproductive function. Previous research suggested a role for SF-1 in energy balance, as evidenced by studies on SF-1 KO mice and the regulatory role of SF-1 neurons in the VMH [10, 11]. Subsequent investigations into the impact of SF-1 on leptin signaling and the development of obesity, especially under high-fat diet conditions, further informed on the multifaceted role of SF-1 in metabolic adaptation [16-24]. Optogenetic and chemogenetic studies provided additional insights, revealing the diverse effects of SF-1 neurons on feeding behavior, energy expenditure, and glucose regulation [22-25].

There is limited previous evidence on the possible link between rare NR5A1/SF-1 variants and obesity in humans. Previous reports are limited to individuals with DSD [26-30], revealing a critical gap in translating mouse model findings to human populations. This scarcity is exacerbated by challenges in studying metabolic health outcomes in individuals with rare NR5A1/SF-1 variants due to the few numbers [26-30]. In individuals with rare NR5A1/SF-1 variants and DSD, the clinical focus is primarily on reproductive and urogenital issues during younger ages, and the majority of affected individuals have a 46,XY karyotype [1, 4]. Women with rare NR5A1/SF-1 variants and DSD or POI are less frequently reported [3, 5, 7], which limits the understanding of the sex-specific role of SF-1. In our study the association of rare NR5A1/SF-1 variants with obesity particularly in women, could be linked to sex-specific hormonal dysregulation, potentially rendering women with these variants more susceptible to obesity compared to men. Our study not only emphasizes the genetic association between rare predicted deleterious NR5A1/SF-1 variants and obesity but also underscores the need for more comprehensive research to bridge the translational gap between mouse models and humans, younger and older age, as well as sex differences.

In this study, we also observed modestly higher platelet counts in carriers of deleterious missense NR5A1/SF-1 variants. Spleen anomalies including hypoplasia are reported in some individuals with rare NR5A1/SF-1 variants. While NR5A1 carriers had higher platelet counts, spleen volume was measured only in a subset of UKBB and only 29 carriers had this data. Thrombocytosis (high platelet count > 450 000/μL) [54] has previously been reported in individuals with DSD due to NR5A1/SF-1 variants, possibly indicating spleen dysfunction [4, 33, 55]. In combination with our findings, we suggest that these variants affect spleen function, due to the impaired activity of SF-1 during spleen development [30, 32, 33].

It is noteworthy that we did not find carriers of rare protein truncating variants in NR5A1/SF-1 in the UKBB population. This likely indicates the gene's critical importance and its conservation within human populations. While there are no known mutation hotspots in the SF-1 protein, our study revealed a concentration of rare deleterious missense NR5A1/SF-1 variants in the DBD and LBD of the SF-1 protein. This observation is consistent with findings from other studies, where the majority of reported rare deleterious variants are commonly located in these 2 domains, which constitute a significant portion of the protein [1, 4].

We acknowledge limitations of our study. Firstly, consistent with other large-scale exome sequencing studies, we are unable to validate and replicate our findings in an independent cohort. While the association with ANM was our primary hypothesis, associations with other potentially relevant traits should be considered as exploratory and should be confirmed in other large cohorts or patient series. Additionally, our reliance on in silico predictions to assess the functional implications of rare NR5A1/SF-1 variants underscores the need for future extensive experimental characterization to accurately assess their impact on protein function. Furthermore, we acknowledge the influence of “healthy volunteer bias” in the UK Biobank [56], which could attenuate disease associations.

Despite these limitations, our study proposes links between rare predicted deleterious missense NR5A1/SF-1 variants with obesity and adverse metabolic outcomes. It is one of the first studies investigating the impact of rare NR5A1/SF-1 variants in the general population. Further research is needed to identify the possible sex hormone related or independent pathways linking NR5A1/SF-1 to obesity risk. Meanwhile, our study suggests the need for clinical awareness and characterization to detect and mitigate possible adverse health metabolic outcomes in individuals with rare NR5A1/SF-1 variants.

Acknowledgments

We thank the participants and investigators in the UK Biobank study who made this work possible (Resource Application Number 26041; 9905) the UK Biobank Exome Sequencing Consortium (UKB-ESC) members AbbVie, Alnylam Pharmaceuticals, AstraZeneca, Biogen, Bristol-Myers Squibb, Pfizer, Regeneron, and Takeda for funding the generation of the data; the Regeneron Genetics Center for completing the sequencing and initial quality control of the exome sequencing data; and the AstraZeneca Centre for Genomics Research Analytics and Informatics team for processing and analysis of sequencing and phenotype data.

Funding

C.K., C.e.F: Swiss National Science Foundation (320030-197725), Swiss National Science Foundation Mobility grant (320030-197725) and Walter und Anne Marie  Boveri Stiftung, Switzerland. R.Y.J., K.A.K., E.J.G., J.R.B.P., and K.K.O. are supported by the Medical Research Council (MC_UU_00006/2).

Author Contributions

All authors read and approved the final version of the manuscript. C.K.: Data analysis, interpretation, and presentation; creation of tables and figures; manuscript writing. R.Y.J.: Discussion of data interpretation and presentation; manuscript reviewing and proofreading. K.A.K.: Data analysis; discussion of data interpretation and presentation; manuscript reviewing and proofreading. E.J.G.: Discussion of data interpretation and presentation; manuscript reviewing and proofreading. J.R.B.P.: Discussion of data interpretation and presentation; manuscript reviewing and proofreading. C.e.F.: Project idea and study design; discussion of data interpretation and presentation; Manuscript reviewing and proofreading. K.K.O.: Project idea and study design; discussion of data interpretation and presentation; manuscript reviewing and proofreading.

Disclosures

J.R.B.P. and E.J.G. are employees and shareholders of Insmed Ltd.

Data Availability

Original data generated and analyzed during this study are included in this published article or in the data repositories listed in References. UK Biobank data are available on application (https://www.ukbiobank.ac.uk/enable-your-research/register).

References

1

Fabbri-Scallet
 
H
,
de Sousa
 
LM
,
Maciel-Guerra
 
AT
,
Guerra-Júnior
 
G
,
de Mello
 
MP
.
Mutation update for the NR5A1 gene involved in DSD and infertility
.
Hum Mutat
.
2020
;
41
(
1
):
58
68
.

2

Domenice
 
S
,
Machado
 
AZ
,
Ferreira
 
FM
, et al.  
Wide spectrum of NR5A1-related phenotypes in 46,XY and 46,XX individuals
.
Birth Defects Res C Embryo Today
.
2016
;
108
(
4
):
309
320
.

3

Baetens
 
D
,
Stoop
 
H
,
Peelman
 
F
, et al.  
NR5A1 is a novel disease gene for 46,XX testicular and ovotesticular disorders of sex development
.
Genet Med
.
2017
;
19
(
4
):
367
376
.

4

Kouri
 
C
,
Sommer
 
G
,
Martinez de Lapiscina
 
I
, et al.  
Clinical and genetic characteristics of a large international cohort of individuals with rare NR5A1/SF-1 variants of sex development
.
EBioMedicine
.
2024
:
99
:
104941
.

5

Lourenço
 
D
,
Brauner
 
R
,
Lin
 
L
, et al.  
Mutations in NR5A1 associated with ovarian insufficiency
.
N Engl J Med
.
2009
;
360
(
12
):
1200
1210
.

6

Bashamboo
 
A
,
Ferraz-de-Souza
 
B
,
Lourenço
 
D
, et al.  
Human male infertility associated with mutations in NR5A1 encoding steroidogenic factor 1
.
Am J Hum Genet
.
2010
;
87
(
4
):
505
512
.

7

Camats
 
N
,
Pandey
 
AV
,
Fernández-Cancio
 
M
, et al.  
Ten novel mutations in the NR5A1 gene cause disordered sex development in 46,XY and ovarian insufficiency in 46,XX individuals
.
J Clin Endocrinol Metab
.
2012
;
97
(
7
):
E1294
E1306
.

8

Suntharalingham
 
JP
,
Buonocore
 
F
,
Duncan
 
AJ
,
Achermann
 
JC
.
DAX-1 (NR0B1) and steroidogenic factor-1 (SF-1, NR5A1) in human disease
.
Best Pract Res Clin Endocrinol Metab
.
2015
;
29
(
4
):
607
619
.

9

Schimmer
 
BP
,
White
 
PC
.
Minireview: steroidogenic factor 1: its roles in differentiation, development, and disease
.
Mol Endocrinol
.
2010
;
24
(
7
):
1322
1337
.

10

Ikeda
 
Y
,
Luo
 
X
,
Abbud
 
R
,
Nilson
 
JH
,
Parker
 
KL
.
The nuclear receptor steroidogenic factor 1 is essential for the formation of the ventromedial hypothalamic nucleus
.
Mol Endocrinol
.
1995
;
9
(
4
):
478
486
.

11

Shinoda
 
K
,
Lei
 
H
,
Yoshii
 
H
, et al.  
Developmental defects of the ventromedial hypothalamic nucleus and pituitary gonadotroph in the Ftz-F1 disrupted mice
.
Dev Dyn
.
1995
;
204
(
1
):
22
29
.

12

Kim
 
JH
,
Choi
 
JH
.
Pathophysiology and clinical characteristics of hypothalamic obesity in children and adolescents
.
Ann Pediatr Endocrinol Metab
.
2013
;
18
(
4
):
161
167
.

13

Gaur
 
A
,
Pal
 
GK
,
Ananthanarayanan
 
PH
,
Pal
 
P
.
Role of ventromedial hypothalamus in high fat diet induced obesity in male rats: association with lipid profile, thyroid profile and insulin resistance
.
Ann Neurosci
.
2014
;
21
(
3
):
104
107
.

14

Timper
 
K
,
Brüning
 
JC
.
Hypothalamic circuits regulating appetite and energy homeostasis: pathways to obesity
.
Dis Model Mech
.
2017
;
10
(
6
):
679
689
.

15

Tran
 
LT
,
Park
 
S
,
Kim
 
SK
,
Lee
 
JS
,
Kim
 
KW
,
Kwon
 
O
.
Hypothalamic control of energy expenditure and thermogenesis
.
Exp Mol Med
.
2022
;
54
(
4
):
358
369
.

16

Choi
 
YH
,
Fujikawa
 
T
,
Lee
 
J
,
Reuter
 
A
,
Kim
 
KW
.
Revisiting the ventral medial nucleus of the hypothalamus: the roles of SF-1 neurons in energy homeostasis
.
Front Neurosci
.
2013
;
7
:
71
.

17

Majdic
 
G
,
Young
 
M
,
Gomez-Sanchez
 
E
, et al.  
Knockout mice lacking steroidogenic factor 1 are a novel genetic model of hypothalamic obesity
.
Endocrinology
.
2002
;
143
(
2
):
607
614
.

18

Dhillon
 
H
,
Zigman
 
JM
,
Ye
 
C
, et al.  
Leptin directly activates SF1 neurons in the VMH, and this action by leptin is required for normal body-weight homeostasis
.
Neuron
.
2006
;
49
(
2
):
191
203
.

19

Bingham
 
NC
,
Anderson
 
KK
,
Reuter
 
AL
,
Stallings
 
NR
,
Parker
 
KL
.
Selective loss of leptin receptors in the ventromedial hypothalamic nucleus results in increased adiposity and a metabolic syndrome
.
Endocrinology
.
2008
;
149
(
5
):
2138
2148
.

20

Sohn
 
JW
,
Oh
 
Y
,
Kim
 
KW
,
Lee
 
S
,
Williams
 
KW
,
Elmquist
 
JK
.
Leptin and insulin engage specific PI3K subunits in hypothalamic SF1 neurons
.
Mol Metab
.
2016
;
5
(
8
):
669
679
.

21

Fosch
 
A
,
Zagmutt
 
S
,
Casals
 
N
,
Rodríguez-Rodríguez
 
R
.
New insights of SF1 neurons in hypothalamic regulation of obesity and diabetes
.
Int J Mol Sci
.
2021
;
22
(
12
):
6186
.

22

Viskaitis
 
P
,
Irvine
 
EE
,
Smith
 
MA
, et al.  
Modulation of SF1 neuron activity coordinately regulates both feeding behavior and associated emotional states
.
Cell Rep
.
2017
;
21
(
12
):
3559
3572
.

23

Coutinho
 
EA
,
Okamoto
 
S
,
Ishikawa
 
AW
, et al.  
Activation of SF1 neurons in the ventromedial hypothalamus by DREADD technology increases insulin sensitivity in peripheral tissues
.
Diabetes
.
2017
;
66
(
9
):
2372
2386
.

24

Bonaventura
 
J
,
Eldridge
 
MAG
,
Hu
 
F
, et al.  
High-potency ligands for DREADD imaging and activation in rodents and monkeys
.
Nat Commun
.
2019
;
10
(
1
):
4627
.

25

Meek
 
TH
,
Nelson
 
JT
,
Matsen
 
ME
, et al.  
Functional identification of a neurocircuit regulating blood glucose
.
Proc Natl Acad Sci U S A
.
2016
;
113
(
14
):
E2073
E2082
.

26

Achermann
 
JC
,
Ito
 
M
,
Ito
 
M
,
Hindmarsh
 
PC
,
Jameson
 
JL
.
A mutation in the gene encoding steroidogenic factor-1 causes XY sex reversal and adrenal failure in humans
.
Nat Genet
.
1999
;
22
(
2
):
125
126
.

27

Correa
 
RV
,
Domenice
 
S
,
Bingham
 
NC
, et al.  
A microdeletion in the ligand binding domain of human steroidogenic factor 1 causes XY sex reversal without adrenal insufficiency
.
J Clin Endocrinol Metab
.
2004
;
89
(
4
):
1767
1772
.

28

Hasegawa
 
T
,
Fukami
 
M
,
Sato
 
N
, et al.  
Testicular dysgenesis without adrenal insufficiency in a 46,XY patient with a heterozygous inactive mutation of steroidogenic factor-1
.
J Clin Endocrinol Metab
.
2004
;
89
(
12
):
5930
5935
.

29

Teoli
 
J
,
Mallet
 
D
,
Renault
 
L
, et al.  
Case report: longitudinal follow-up and testicular sperm extraction in a patient with a pathogenic NR5A1 (SF-1) frameshift variant: p.(Phe70Serfs*5)
.
Front Endocrinol (Lausanne)
.
2023
;
14
:
1171822
.

30

Colson
 
C
,
Aubry
 
E
,
Cartigny
 
M
, et al.  
SF1 and spleen development: new heterozygous mutation, literature review and consequences for NR5A1-mutated patient's management
.
Clin Genet
.
2017
;
92
(
1
):
99
103
.

31

Liu
 
W
,
Liu
 
M
,
Fan
 
W
,
Nawata
 
H
,
Yanase
 
T
.
The Gly146Ala variation in human SF-1 gene: its association with insulin resistance and type 2 diabetes in Chinese
.
Diabetes Res Clin Pract
.
2006
;
73
(
3
):
322
328
.

32

Zangen
 
D
,
Kaufman
 
Y
,
Banne
 
E
, et al.  
Testicular differentiation factor SF-1 is required for human spleen development
.
J Clin Invest
.
2014
;
124
(
5
):
2071
2075
.

33

Cools
 
M
,
Grijp
 
C
,
Neirinck
 
J
, et al.  
Spleen function is reduced in individuals with NR5A1 variants with or without a difference of sex development: a cross-sectional study
.
Eur J Endocrinol
.
2024
;
190
(
1
):
34
43
.

34

Morohashi
 
K
,
Tsuboi-Asai
 
H
,
Matsushita
 
S
, et al.  
Structural and functional abnormalities in the spleen of an mFtz-F1 gene-disrupted mouse
.
Blood
.
1999
;
93
(
5
):
1586
1594
.

35

Backman
 
JD
,
Li
 
AH
,
Marcketta
 
A
, et al.  
Exome sequencing and analysis of 454,787 UK Biobank participants
.
Nature
.
2021
;
599
(
7886
):
628
634
.

36

Bycroft
 
C
,
Freeman
 
C
,
Petkova
 
D
, et al.  
The UK Biobank resource with deep phenotyping and genomic data
.
Nature
.
2018
;
562
(
7726
):
203
209
.

37

Gardner
 
EJ
,
Kentistou
 
KA
,
Stankovic
 
S
, et al.  
Damaging missense variants in IGF1R implicate a role for IGF-1 resistance in the etiology of type 2 diabetes
.
Cell Genom
.
2022
;
2
(
12
). Doi:

38

McLaren
 
W
,
Gil
 
L
,
Hunt
 
SE
, et al.  
The ensembl variant effect predictor
.
Genome Biol
.
2016
;
17
(
1
):
122
.

39

Ioannidis
 
NM
,
Rothstein
 
JH
,
Pejaver
 
V
, et al.  
REVEL: an ensemble method for predicting the pathogenicity of rare missense variants
.
Am J Hum Genet
.
2016
;
99
(
4
):
877
885
.

40

Rentzsch
 
P
,
Witten
 
D
,
Cooper
 
GM
,
Shendure
 
J
,
Kircher
 
M
.
CADD: predicting the deleteriousness of variants throughout the human genome
.
Nucleic Acids Res
.
2019
;
47
(
D1
):
D886
Dd94
.

41

Karczewski
 
KJ
,
Francioli
 
LC
,
Tiao
 
G
, et al.  
The mutational constraint spectrum quantified from variation in 141,456 humans
.
Nature
.
2020
;
581
(
7809
):
434
443
.

42

Morales
 
J
,
Pujar
 
S
,
Loveland
 
JE
, et al.  
A joint NCBI and EMBL-EBI transcript set for clinical genomics and research
.
Nature
.
2022
;
604
(
7905
):
310
315
.

43

Shekari
 
S
,
Stankovic
 
S
,
Gardner
 
EJ
, et al.  
Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency
.
Nat Med
.
2023
;
29
(
7
):
1692
1699
.

44

Kouri
 
C
,
Jia
 
RY
,
Kentistou
 
KA
, et al.  Supplementary material of population-based study of rare coding variants in NR5A1/SF-1.
2024
. https://zenodo.org/doi/10.5281/zenodo.13739320

45

Loh
 
PR
,
Tucker
 
G
,
Bulik-Sullivan
 
BK
, et al.  
Efficient Bayesian mixed-model analysis increases association power in large cohorts
.
Nat Genet
.
2015
;
47
(
3
):
284
290
.

46

Seabold
 
S
,
Perktold
 
J.
 
Statsmodels: econometric and statistical modeling with python. SciPy Proceedings, Proceedings of the 9th Python in Science Conference, Austin, Texas; 2010.

47

Lin
 
L
,
Philibert
 
P
,
Ferraz-de-Souza
 
B
, et al.  
Heterozygous missense mutations in steroidogenic factor 1 (SF1/Ad4BP, NR5A1) are associated with 46,XY disorders of sex development with normal adrenal function
.
J Clin Endocrinol Metab
.
2007
;
92
(
3
):
991
999
.

48

Jaillard
 
S
,
Sreenivasan
 
R
,
Beaumont
 
M
, et al.  
Analysis of NR5A1 in 142 patients with premature ovarian insufficiency, diminished ovarian reserve, or unexplained infertility
.
Maturitas
.
2020
;
131
:
78
86
.

49

Hu
 
X
,
Li
 
N
,
Xu
 
Y
, et al.  
Proband-only medical exome sequencing as a cost-effective first-tier genetic diagnostic test for patients without prior molecular tests and clinical diagnosis in a developing country: the China experience
.
Genet Med
.
2018
;
20
(
9
):
1045
1053
.

50

Philibert
 
P
,
Zenaty
 
D
,
Lin
 
L
, et al.  
Mutational analysis of steroidogenic factor 1 (NR5a1) in 24 boys with bilateral anorchia: a French collaborative study
.
Hum Reprod
.
2007
;
22
(
12
):
3255
3261
.

51

Ferraz-de-Souza
 
B
,
Lin
 
L
,
Achermann
 
JC
.
Steroidogenic factor-1 (SF-1, NR5A1) and human disease
.
Mol Cell Endocrinol
.
2011
;
336
(
1-2
):
198
205
.

52

Camats
 
N
,
Flück
 
CE
,
Audí
 
L
.
Oligogenic origin of differences of sex development in humans
.
Int J Mol Sci
.
2020
;
21
(
5
):
1809
.

53

Kouri
 
C
,
Sommer
 
G
,
Flück
 
CE
.
Oligogenic causes of human differences of sex development: facing the challenge of genetic complexity
.
Horm Res Paediatr
.
2023
;
96
(
2
):
169
179
.

54

Sałacki
 
AJ
,
Wysokiński
 
A
.
Reactive thrombocytosis resulting from frequent blood donations as an extremely rare cause of ST segment elevation myocardial infarction in the case of a 19-year-old male
.
Ann Agric Environ Med
.
2018
;
25
(
4
):
602
604
.

55

Kirkineska
 
L
,
Perifanis
 
V
,
Vasiliadis
 
T
.
Functional hyposplenism
.
Hippokratia
.
2014
;
18
(
1
):
7
11
.

56

Fry
 
A
,
Littlejohns
 
TJ
,
Sudlow
 
C
, et al.  
Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population
.
Am J Epidemiol
.
2017
;
186
(
9
):
1026
1034
.

Abbreviations

     
  • ANM

    age at natural menopause

  •  
  • DBD

    DNA binding domain

  •  
  • DSD

    differences of sex development

  •  
  • GLM

    generalized linear model

  •  
  • HbA1c

    glycated hemoglobin

  •  
  • HDL

    high-density lipoprotein

  •  
  • IGF-1

    insulin-like growth factor 1

  •  
  • KO

    knockout

  •  
  • LBD

    ligand binding domain

  •  
  • NR5A1

    Nuclear Receptor Subfamily 5 Group A Member 1

  •  
  • OR

    odds ratio

  •  
  • PTV

    protein truncating variant

  •  
  • RAP

    UKBB Research Analysis Platform

  •  
  • SF-1

    steroidogenic factor 1

  •  
  • SNV

    single nucleotide variation

  •  
  • UKBB

    UK Biobank

  •  
  • VEP

    Variant Effect Predictor

  •  
  • VMH

    ventromedial hypothalamic nucleus

  •  
  • WES

    whole exome sequencing

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. See the journal About page for additional terms.