Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits

Proteome-wide Mendelian randomization identifies potential therapeutic targets for nonalcoholic fatty liver diseases

Nonalcoholic fatty liver disease (NAFLD) is the predominant cause of liver pathology. Current evidence highlights plasma proteins as potential therapeutic targets. However, their mechanistic roles in NAFLD remain unclear. This study investigated the involvement of specific plasma proteins and intermediate risk factors in NAFLD progression. Two-sample Mendelian randomization (MR) analysis was conducted to examine the association between plasma proteins and NAFLD. Colocalization analysis determined the shared causal variants between the identified proteins and NAFLD. The MR analysis was applied separately to proteins, risk factors, and NAFLD. Mediator shares were computed by detecting the correlations among these elements. Phenome-wide association studies (phewas) were utilized to assess the safety implications of targeting these proteins. Among 1,834 cis-protein quantitative trait loci (cis-pQTLs), after-FDR correction revealed correlations between the plasma levels of four gene-predicted proteins (CSPG3, CILP2, Apo-E, and GCKR) and NAFLD. Colocalization analysis indicated shared causal variants for CSPG3 and GCKR in NAFLD (posterior probability > 0.8). Out of the 22 risk factors screened for MR analysis, only 8 showed associations with NAFLD (p ≤ 0.05), while 4 linked to CSPG3 and GCKR. The mediator shares for these associations were calculated separately. Additionally, reverse MR analysis was performed on the pQTLs, risk factors, and NAFLD, which exhibited a causal relationship with forward MR analysis. Finally, phewas summarized the potential side effects of associated-targeting proteins, including CSPG3 and GCKR. Our research emphasized the potential therapeutic targets for NAFLD and provided modifiable risk factors for preventing NAFLD.

Multi-ancestry Genome-Wide Association Study of Early Childhood Caries

Early childhood caries (ECC) is the most common non-communicable childhood disease. It is an important health problem with known environmental and social/behavioral influences that lacks evidence for specific associated genetic risk loci. To address this knowledge gap, we conducted a genome-wide association study of ECC in a multi-ancestry population of U.S. preschool-age children (n=6,103) participating in a community-based epidemiologic study of early childhood oral health. Calibrated examiners used ICDAS criteria to measure ECC with the primary trait using the dmfs index with decay classified as macroscopic enamel loss (ICDAS ≥3). We estimated heritability, concordance rates, and conducted genome-wide association analyses to estimate overall genetic effects; the effects stratified by sex, household water fluoride, and dietary sugar; and leveraged the combined gene/gene-environment effects using the 2-degree-of-freedom (2df) joint test. The common genetic variants explained 24% of the phenotypic variance (heritability) of the primary ECC trait and the concordance rate was higher with a higher degree of relatedness. We identified 21 novel non-overlapping genome-wide significant loci for ECC. Two loci, namely RP11-856F16 . 2 (rs74606067) and SLC41A3 (rs71327750) showed evidence of association with dental caries in external cohorts, namely the GLIDE consortium adult cohort (n=∼487,000) and the GLIDE pediatric cohort (n=19,000), respectively. The gene-based tests identified TAAR6 as a genome-wide significant gene. Implicated genes have relevant biological functions including roles in tooth development and taste. These novel associations expand the genomics knowledge base for this common childhood disease and underscore the importance of accounting for sex and pertinent environmental exposures in genetic investigations of oral health.

The regulation of methylation on the Z chromosome and the identification of multiple novel Male Hyper-Methylated regions in the chicken

DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. To investigate how variation in DNA methylation is regulated on the Z chromosome we utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals. The relatively large numbers of individuals allowed us to identify additional genomic MHM regions on the Z chromosome that were significantly differentially methylated between the sexes. These regions appear to down-regulate local gene expression in males, but not remove it entirely (unlike the lncRNAs identified in the initial MHM regions). These MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. In addition, quantitative trait loci (QTL) that regulate variation in methylation on the Z chromosome, and those loci that regulate methylation on the autosomes that derive from the Z chromosome were mapped. Trans-effect hotspots were also identified that were based on the autosomes but affected the Z, and also one that was based on the Z chromosome but that affected both autosomal and sex chromosome DNA methylation regulation. We show that both cis and trans loci that originate from the Z chromosome never exhibit an interaction with sex, whereas trans loci originating from the autosomes but affecting the Z chromosome always display such an interaction. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes.

Sex-dependent placental methylation quantitative trait loci provide insight into the prenatal origins of childhood onset traits and conditions

Molecular quantitative trait loci (QTLs) allow us to understand the biology captured in genome-wide association studies (GWASs). The placenta regulates fetal development and shows sex differences in DNA methylation. We therefore hypothesized that placental methylation QTL (mQTL) explain variation in genetic risk for childhood onset traits, and that effects differ by sex. We analyzed 411 term placentas from two studies and found 49,252 methylation (CpG) sites with mQTL and 2,489 CpG sites with sex-dependent mQTL. All mQTL were enriched in regions that typically affect gene expression in prenatal tissues. All mQTL were also enriched in GWAS results for growth- and immune-related traits, but male- and female-specific mQTL were more enriched than cross-sex mQTL. mQTL colocalized with trait loci at 777 CpG sites, with 216 (28%) specific to males or females. Overall, mQTL specific to male and female placenta capture otherwise overlooked variation in childhood traits.

Socio-cultural practices may have affected sex differences in stature in Early Neolithic Europe

The rules and structure of human culture impact health as much as genetics or environment. To study these relationships, we combine ancient DNA (n = 230), skeletal metrics (n = 391), palaeopathology (n = 606) and dietary stable isotopes (n = 873) to analyse stature variation in Early Neolithic Europeans from North Central, South Central, Balkan and Mediterranean regions. In North Central Europe, stable isotopes and linear enamel hypoplasias indicate high environmental stress across sexes, but female stature is low, despite polygenic scores identical to males, and suggests that cultural factors preferentially supported male recovery from stress. In Mediterranean populations, sexual dimorphism is reduced, indicating male vulnerability to stress and no strong cultural preference for males. Our analysis indicates that biological effects of sex-specific inequities can be linked to cultural influences at least as early as 7,000 yr ago, and culture, more than environment or genetics, drove height disparities in Early Neolithic Europe.

MRSamePopTest: introducing a simple falsification test for the two-sample mendelian randomisation 'same population' assumption

Two-sample MR is an increasingly popular method for strengthening causal inference in epidemiological studies. For the effect estimates to be meaningful, variant-exposure and variant-outcome associations must come from comparable populations. A recent systematic review of two-sample MR studies found that, if assessed at all, MR studies evaluated this assumption by checking that the genetic association studies had similar demographics. However, it is unclear if this is sufficient because less easily accessible factors may also be important. Here we propose an easy-to-implement falsification test. Since recent theoretical developments in causal inference suggest that a causal effect estimate can generalise from one study to another if there is exchangeability of effect modifiers, we suggest testing the homogeneity of variant-phenotype associations for a phenotype which has been measured in both genetic association studies as a method of exploring the 'same-population' test. This test could be used to facilitate designing MR studies with diverse populations. We developed a simple R package to facilitate the implementation of our proposed test. We hope that this research note will result in increased attention to the same-population assumption, and the development of better sensitivity analyses.

The role of gonadal hormones and sex chromosomes in sex-dependent effects of early nutrition on metabolic health

Early-life conditions such as prenatal nutrition can have long-term effects on metabolic health, and these effects may differ between males and females. Understanding the biological mechanisms underlying sex differences in the response to early-life environment will improve interventions, but few such mechanisms have been identified, and there is no overall framework for understanding sex differences. Biological sex differences may be due to chromosomal sex, gonadal sex, or interactions between the two. This review describes approaches to distinguish between the roles of chromosomal and gonadal sex, and summarizes findings regarding sex differences in metabolism. The Four Core Genotypes (FCG) mouse model allows dissociation of the sex chromosome genotype from gonadal type, whereas the XY* mouse model can be used to distinguish effects of X chromosome dosage vs the presence of the Y chromosome. Gonadectomy can be used to distinguish between organizational (permanent) and activational (reversible) effects of sex hormones. Baseline sex differences in a variety of metabolic traits are influenced by both activational and organizational effects of gonadal hormones, as well as sex chromosome complement. Thus far, these approaches have not been widely applied to examine sex-dependent effects of prenatal conditions, although a number of studies have found activational effects of estradiol to be protective against the development of hypertension following early-life adversity. Genes that escape X chromosome inactivation (XCI), such as Kdm5c, contribute to baseline sex-differences in metabolism, while Ogt, another XCI escapee, leads to sex-dependent responses to prenatal maternal stress. Genome-wide approaches to the study of sex differences include mapping genetic loci influencing metabolic traits in a sex-dependent manner. Seeking enrichment for binding sites of hormone receptors among genes showing sexually-dimorphic expression can elucidate the relative roles of hormones. Using the approaches described herein to identify mechanisms underlying sex-dependent effects of early nutrition on metabolic health may enable the identification of fundamental mechanisms and potential interventions.

Sex differences in risk factors for coronary heart disease events: a prospective cohort study in Iran

We investigated sex-specific associations and their differences between major cardiovascular risk factors and the risk of incident coronary heart disease (CHD) and hard CHD (defined as nonfatal myocardial infarction and CHD death). A total of 7518 (3377 men) participants from the Tehran Lipid and Glucose Study were included. Cox models were used to estimate the hazard ratios (HRs) and women-to-men ratios of HRs for CHD events associated with each risk factor. During 20 years of follow-up (1999-2018), 1068 (631 men) and 345 (238 men) new cases of CHD and hard CHD, respectively, were documented. In total population, the incidence rates per 1000 person-years were 9.5 (9.0-10.1) and 2.9 (2.6-3.2) for CHD and hard CHD, respectively. Hypertension, diabetes, pre-diabetes, and a high waist-to-hip ratio (WHR) were associated with a greater HR of hard CHD in women than men; the women-to-men HRs were 2.85 [1.36-5.98], 1.92 [1.11-3.31], 2.04 [1.09-3.80] and 1.42 [1.10-1.82], respectively. Diabetes was associated with a higher HR of CHD in women than men (ratio of HRs 1.49 (1.10-2.01). In conclusion, we found that hypertension, diabetes, pre-diabetes, and high WHR conferred a greater excess risk of CHD events in women than in men, suggesting that Iranian women may require greater attention for the prevention of CHD events.

The evolution of small molecule enzyme activators

There is a myriad of enzymes within the body responsible for maintaining homeostasis by providing the means to convert substrates to products as and when required. Physiological enzymes are tightly controlled by many signaling pathways and their products subsequently control other pathways. Traditionally, most drug discovery efforts focus on identifying enzyme inhibitors, due to upregulation being prevalent in many diseases and the existence of endogenous substrates that can be modified to afford inhibitor compounds. As enzyme downregulation and reduction of endogenous activators are observed in multiple diseases, the identification of small molecules with the ability to activate enzymes has recently entered the medicinal chemistry toolbox to afford chemical probes and potential therapeutics as an alternative means to intervene in diseases. In this review we highlight the progress made in the identification and advancement of non-kinase enzyme activators and their potential in treating various disease states.

Sex- and reproductive status-specific relationships between body composition and non-alcoholic fatty liver disease

BACKGROUND

Sex and reproductive status differences exist in both non-alcoholic fatty liver disease (NAFLD) and body composition. Our purpose was to investigate the relationship between body composition and the severity of liver steatosis and fibrosis in NAFLD in different sex and reproductive status populations.

METHODS

This cross-sectional study included 880 patients (355 men, 417 pre-menopausal women, 108 post-menopausal women). Liver steatosis and fibrosis and body composition data were measured using FibroScan and a bioelectrical impedance body composition analyzer (BIA), respectively, and the following parameters were obtained: liver stiffness measurement (LSM), controlled attenuation parameter (CAP), waist circumference (WC), body mass index (BMI), percent body fat (PBF), visceral fat area (VFA), appendicular skeletal muscle mass (ASM), appendicular skeletal muscle mass index (ASMI), fat mass (FM), fat free mass (FFM), and FFM to FM ratio (FFM/FM). Multiple ordinal logistic regression (MOLR) was used to analyze the independent correlation between body composition indicators and liver steatosis grade and fibrosis stage in different sex and menopausal status populations.

RESULTS

Men had higher WC, ASM, ASMI, FFM, and FFM/FM than pre- or post-menopausal women, while pre-menopausal women had higher PBF, VFA, and FM than the other two groups (p < 0.001). Besides, men had greater CAP and LSM values (p < 0.001). For MOLR, after adjusting for confounding factors, WC (OR, 1.07; 95% CI, 1.02-1.12; P = 0.011) and FFM/FM (OR, 0.52; 95% CI, 0.31-0.89; P = 0.017) in men and visceral obesity (OR, 4.16; 95% CI, 1.09-15.90; P = 0.037) in post-menopausal women were independently associated with liver steatosis grade. WC and visceral obesity were independently associated with liver fibrosis stage in men (OR, 1.05; 95% CI, 1.01-1.09, P = 0.013; OR, 3.92; 95% CI, 1.97-7.81; P < 0.001, respectively).

CONCLUSIONS

Increased WC and low FFM/FM in men and visceral obesity in post-menopausal women were independent correlates of more severe liver steatosis. In addition, increased WC and visceral obesity were independent correlates of worse liver fibrosis in men. These data support the sex- and reproductive status-specific management of NAFLD.

Genetics and epigenetics in the obesity phenotyping scenario

Obesity is a common complex trait that elevates the risk for various diseases, including type 2 diabetes and cardiovascular disease. A combination of environmental and genetic factors influences the pathogenesis of obesity. Advances in genomic technologies have driven the identification of multiple genetic loci associated with this disease, ranging from studying severe onset cases to investigating common multifactorial polygenic forms. Additionally, findings from epigenetic analyses of modifications to the genome that do not involve changes to the underlying DNA sequence have emerged as key signatures in the development of obesity. Such modifications can mediate the effects of environmental factors, including diet and lifestyle, on gene expression and clinical presentation. This review outlines what is known about the genetic and epigenetic contributors to obesity susceptibility, along with the albeit limited therapeutic options currently available. Furthermore, we delineate the potential mechanisms of actions through which epigenetic changes can mediate environmental influences and the related opportunities they present for future interventions in the management of obesity.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

A roadmap for sex- and gender-disaggregated health research

Sex and gender are fundamental aspects of health and wellbeing. Yet many research studies fail to consider sex or gender differences, and even when they do this is often limited to merely cataloguing such differences in the makeup of study populations. The evidence on sex and gender differences is thus incomplete in most areas of medicine. This article presents a roadmap for the systematic conduct of sex- and gender-disaggregated health research. We distinguish three phases: the exploration of sex and gender differences in disease risk, presentation, diagnosis, treatment, and outcomes; explaining any found differences by revealing the underlying mechanisms; and translation of the implications of such differences to policy and practice. For each phase, we provide critical methodological considerations and practical examples are provided, taken primarily from the field of cardiovascular disease. We also discuss key overarching themes and terminology that are at the essence of any study evaluating the relevance of sex and gender in health. Here, we limit ourselves to binary sex and gender in order to produce a coherent, succinct narrative. Further disaggregation by sex and gender separately and which recognises intersex, non-binary, and gender-diverse identities, as well as other aspects of intersectionality, can build on this basic minimum level of disaggregation. We envision that uptake of this roadmap, together with wider policy and educational activities, will aid researchers to systematically explore and explain relevant sex and gender differences in health and will aid educators, clinicians, and policymakers to translate the outcomes of research in the most effective and meaningful way, for the benefit of all.

Mendelian randomization analysis revealed potential metabolic causal factors for breast cancer

Observational studies showed that metabolic phenotypes were associated with the risk of developing breast cancer (BC). However, those results are inconsistent regarding the magnitude of the association, particularly by subtypes of breast cancer. Furthermore, the mechanisms of the association remain unclear. We performed two-sample Mendelian randomization (MR) analyses to evaluate the causal effect of metabolic risk factors on breast cancer in the European population. Assessed individually using MR, body mass index (BMI) (odds ratio [OR] 0.94, 95% Confidence interval [CI] 0.90-0.98, P = 0.007), high-density lipoprotein cholesterol (HDL-C) (OR 1.10, 95% CI 1.07-1.13, P = 6.10 × 10-11) and triglycerides (TG) (OR 0.92, 95% CI 0.90-0.96, P = 1.58 × 10-6) were causally related to breast cancer risk. In multivariable MR, only HDL-C (OR 1.08; 95% CI 1.02-1.14; P = 0.02) retained a robust effect, suggesting that the genetic association between BMI, HDL-C and TG with breast cancer risk in univariable analysis was explained via HDL-C. These findings suggest a possible causal role of HDL-C in breast cancer etiology.

Large-scale exome sequence analysis identifies sex- and age-specific determinants of obesity

Obesity contributes substantially to the global burden of disease and has a significant heritable component. Recent large-scale exome sequencing studies identified several genes in which rare, protein-coding variants have large effects on adult body mass index (BMI). Here we extended such work by performing sex-stratified associations in the UK Biobank study (N∼420,000). We identified genes in which rare heterozygous loss-of-function increases adult BMI in women (DIDO1, PTPRG, and SLC12A5) and in men (SLTM), with effect sizes up to ∼8 kg/m2. This is complemented by analyses implicating rare variants in OBSCN and MADD for recalled childhood adiposity. The known functions of these genes, as well as findings of common variant genome-wide pathway enrichment analyses, suggest a role for neuron death, apoptosis, and DNA damage response mechanisms in the susceptibility to obesity across the life-course. These findings highlight the importance of considering sex-specific and life-course effects in the genetic regulation of obesity.

Sex-specific Mendelian randomisation to assess the causality of sex differences in the effects of risk factors and treatment: spotlight on hypertension

Hypertension is a key modifiable risk factor for cardiovascular disease. Several observational studies have found a stronger association of blood pressure and cardiovascular disease risk in women compared to men. Since observational studies can be affected by sex-specific residual confounding and reverse causation, it remains unclear whether these differences reflect actual differential effects. Other study designs are needed to uncover the causality of sex differences in the strength of risk factor and treatment effects. Mendelian randomisation (MR) uses genetic variants as instrumental variables to provide evidence about putative causal relations between risk factors and outcomes. By exploiting the random allocation of genes at gamete forming, MR is unaffected by confounding and results in more reliable causal effect estimates. In this review, we discuss why and how sex-specific MR and cis-MR could be used to study sex differences in risk factor and drug target effects. Sex-specific MR can be helpful to strengthen causal inferences in the field of sex differences, where it is often challenging to distinguish nature from nurture. The challenge of sex-specific (drug target) MR lays in leveraging robust genetic instruments from sex-specific GWAS studies which are not commonly available. Knowledge on sex-specific causal effects of hypertension, or other risk factors, could improve clinical practice and health policies by tailoring interventions based on personalised risk. Drug target MR can help to determine the anticipated on-target effects of a drug compound and to identify targets to pursue in drug development.

The Casual Association Inference for the Chain of Falls Risk Factors-Falls-Falls Outcomes: A Mendelian Randomization Study

Previous associations have been observed not only between risk factors and falls but also between falls and their clinical outcomes based on some cross-sectional designs, but their causal associations were still largely unclear. We performed Mendelian randomization (MR), multivariate Mendelian randomization (MVMR), and mediation analyses to explore the effects of falls. Our study data are mainly based on White European individuals (40-69 years) downloaded from the UK Biobank. MR analyses showed that osteoporosis (p = 0.006), BMI (p = 0.003), sleeplessness (p < 0.001), rheumatoid arthritis (p = 0.001), waist circumference (p < 0.001), and hip circumference (p < 0.001) have causal effects on falls. In addition, for every one standard deviation increase in fall risk, the risk of fracture increased by 1.148 (p < 0.001), the risk of stroke increased by 2.908 (p = 0.003), and a 1.016-fold risk increase in epilepsy (p = 0.009). The MVMR found that sleeplessness is an important risk factor for falls. Finally, our mediation analyses estimated the mediation effects of falls on the hip circumference and fracture (p < 0.001), waist circumference and epilepsy (p < 0.001), and sleeplessness and fracture (p = 0.005). Our study inferred the causal effects between risk factors and falls, falls, and outcomes, and also constructed three causal chains from risk factors → falls → falls outcomes.

Sex differences of the shared genetic landscapes between type 2 diabetes and peripheral artery disease in East Asians and Europeans

Type 2 diabetes (T2D) is a critical risk factor for peripheral artery disease (PAD). However, the sex differences in genetic basis, causality, and underlying mechanisms of the two diseases are still unclear. Using sex-stratified and ethnic-based GWAS summary, we explored the genetic correlation and causal relationship between T2D and PAD in both ethnicities and sexes by linkage disequilibrium score regression, LAVA and six Mendelian Randomization approaches. We observed stronger genetic correlations between T2D and PAD in females than males in East Asians and Europeans. East Asian females exhibit higher causal effects of T2D on PAD than males. The gene-level analysis found KCNJ11 and ANK1 genes associated with the cross-trait of T2D and PAD in both sexes. Our study provides genetic evidence for the sex difference of genetic correlations and causal relationships between PAD and T2D, indicating the importance of using sex-specific strategies for monitoring PAD in T2D patients.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Quality control and analytic best practices for testing genetic models of sex differences in large populations

Phenotypic sex-based differences exist for many complex traits. In other cases, phenotypes may be similar, but underlying biology may vary. Thus, sex-aware genetic analyses are becoming increasingly important for understanding the mechanisms driving these differences. To this end, we provide a guide outlining the current best practices for testing various models of sex-dependent genetic effects in complex traits and disease conditions, noting that this is an evolving field. Insights from sex-aware analyses will not only teach us about the biology of complex traits but also aid in achieving the goals of precision medicine and health equity for all.

Recent progress in understanding the genomic architecture of sexual conflict

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.

A non-coding variant linked to metabolic obesity with normal weight affects actin remodelling in subcutaneous adipocytes

Recent large-scale genomic association studies found evidence for a genetic link between increased risk of type 2 diabetes and decreased risk for adiposity-related traits, reminiscent of metabolically obese normal weight (MONW) association signatures. However, the target genes and cellular mechanisms driving such MONW associations remain to be identified. Here, we systematically identify the cellular programmes of one of the top-scoring MONW risk loci, the 2q24.3 risk locus, in subcutaneous adipocytes. We identify a causal genetic variant, rs6712203, an intronic single-nucleotide polymorphism in the COBLL1 gene, which changes the conserved transcription factor motif of POU domain, class 2, transcription factor 2, and leads to differential COBLL1 gene expression by altering the enhancer activity at the locus in subcutaneous adipocytes. We then establish the cellular programme under the genetic control of the 2q24.3 MONW risk locus and the effector gene COBLL1, which is characterized by impaired actin cytoskeleton remodelling in differentiating subcutaneous adipocytes and subsequent failure of these cells to accumulate lipids and develop into metabolically active and insulin-sensitive adipocytes. Finally, we show that perturbations of the effector gene Cobll1 in a mouse model result in organismal phenotypes matching the MONW association signature, including decreased subcutaneous body fat mass and body weight along with impaired glucose tolerance. Taken together, our results provide a mechanistic link between the genetic risk for insulin resistance and low adiposity, providing a potential therapeutic hypothesis and a framework for future identification of causal relationships between genome associations and cellular programmes in other disorders.

A Gene-Acculturation Study of Obesity Among US Hispanic/Latinos: The Hispanic Community Health Study/Study of Latinos

OBJECTIVE

In the United States, Hispanic/Latino adults face a high burden of obesity; yet, not all individuals are equally affected, partly due in part to this ethnic group's marked sociocultural diversity. We sought to analyze the modification of body mass index (BMI) genetic effects in Hispanic/Latino adults by their level of acculturation, a complex biosocial phenomenon that remains understudied.

METHODS

Among 11,747 Hispanic/Latinos adults in the Hispanic Community Health Study/Study of Latinos aged 18 to 76 years from four urban communities (2008-2011), we a) tested our hypothesis that the effect of a genetic risk score (GRS) for increased BMI may be exacerbated by higher levels of acculturation and b) examined if GRS acculturation interactions varied by gender or Hispanic/Latino background group. All genetic modeling controlled for relatedness, age, gender, principal components of ancestry, center, and complex study design within a generalized estimated equation framework.

RESULTS

We observed a GRS increase of 0.34 kg/m 2 per risk allele in weighted mean BMI. The estimated main effect of GRS on BMI varied both across acculturation level and across gender. The difference between high and low acculturation ranged from 0.03 to 0.23 kg/m 2 per risk allele, but varied across acculturation measure and gender.

CONCLUSIONS

These results suggest the presence of effect modification by acculturation, with stronger effects on BMI among highly acculturated individuals and female immigrants. Future studies of obesity in the Hispanic/Latino community should account for sociocultural environments and consider their intersection with gender to better target obesity interventions.

Sex-Specific Features of the Correlation between GWAS-Noticeable Polymorphisms and Hypertension in Europeans of Russia

The aim of the study was directed at studying the sex-specific features of the correlation between genome-wide association studies (GWAS)-noticeable polymorphisms and hypertension (HTN). In two groups of European subjects of Russia (n = 1405 in total), such as men (n = 821 in total: n = 564 HTN, n = 257 control) and women (n = 584 in total: n = 375 HTN, n = 209 control), the distribution of ten specially selected polymorphisms (they have confirmed associations of GWAS level with blood pressure (BP) parameters and/or HTN in Europeans) has been considered. The list of studied loci was as follows: (PLCE1) rs932764 A > G, (AC026703.1) rs1173771 G > A, (CERS5) rs7302981 G > A, (HFE) rs1799945 C > G, (OBFC1) rs4387287 C > A, (BAG6) rs805303 G > A, (RGL3) rs167479 T > G, (ARHGAP42) rs633185 C > G, (TBX2) rs8068318 T > C, and (ATP2B1) rs2681472 A > G. The contribution of individual loci and their inter-locus interactions to the HTN susceptibility with bioinformatic interpretation of associative links was evaluated separately in men's and women's cohorts. The men-women differences in involvement in the disease of the BP/HTN-associated GWAS SNPs were detected. Among women, the HTN risk has been associated with HFE rs1799945 C > G (genotype GG was risky; OR_GG = 11.15 p_permGG = 0.014) and inter-locus interactions of all 10 examined SNPs as part of 26 intergenic interactions models. In men, the polymorphism BAG6 rs805303 G > A (genotype AA was protective; OR_AA = 0.30 p_permAA = 0.0008) and inter-SNPs interactions of eight loci in only seven models have been founded as HTN-correlated. HTN-linked loci and strongly linked SNPs were characterized by pronounced polyvector functionality in both men and women, but at the same time, signaling pathways of HTN-linked genes/SNPs in women and men were similar and were represented mainly by immune mechanisms. As a result, the present study has demonstrated a more pronounced contribution of BP/HTN-associated GWAS SNPs to the HTN susceptibility (due to weightier intergenic interactions) in European women than in men.

Tree-based QTL mapping with expected local genetic relatedness matrices

Understanding the genetic basis of complex phenotypes is a central pursuit of genetics. Genome-wide Association Studies (GWAS) are a powerful way to find genetic loci associated with phenotypes. GWAS are widely and successfully used, but they face challenges related to the fact that variants are tested for association with a phenotype independently, whereas in reality variants at different sites are correlated because of their shared evolutionary history. One way to model this shared history is through the ancestral recombination graph (ARG), which encodes a series of local coalescent trees. Recent computational and methodological breakthroughs have made it feasible to estimate approximate ARGs from large-scale samples. Here, we explore the potential of an ARG-based approach to quantitative-trait locus (QTL) mapping, echoing existing variance-components approaches. We propose a framework that relies on the conditional expectation of a local genetic relatedness matrix given the ARG (local eGRM). Simulations show that our method is especially beneficial for finding QTLs in the presence of allelic heterogeneity. By framing QTL mapping in terms of the estimated ARG, we can also facilitate the detection of QTLs in understudied populations. We use local eGRM to identify a large-effect BMI locus, the CREBRF gene, in a sample of Native Hawaiians in which it was not previously detectable by GWAS because of a lack of population-specific imputation resources. Our investigations can provide intuition about the benefits of using estimated ARGs in population- and statistical-genetic methods in general.

Genome-wide association study of obstructive sleep apnoea in the Million Veteran Program uncovers genetic heterogeneity by sex

BACKGROUND

Genome-wide association studies (GWAS) for obstructive sleep apnoea (OSA) are limited due to the underdiagnosis of OSA, leading to misclassification of OSA, which consequently reduces statistical power. We performed a GWAS of OSA in the Million Veteran Program (MVP) of the U.S. Department of Veterans Affairs (VA) healthcare system, where OSA prevalence is close to its true population prevalence.

METHODS

We performed GWAS of 568,576 MVP participants, stratified by biological sex and by harmonized race/ethnicity and genetic ancestry (HARE) groups of White, Black, Hispanic, and Asian individuals. We considered both BMI adjusted (BMI-adj) and unadjusted (BMI-unadj) models. We replicated associations in independent datasets, and analysed the heterogeneity of OSA genetic associations across HARE and sex groups. We finally performed a larger meta-analysis GWAS of MVP, FinnGen, and the MGB Biobank, totalling 916,696 individuals.

FINDINGS

MVP participants are 91% male. OSA prevalence is 21%. In MVP there were 18 and 6 genome-wide significant loci in BMI-unadj and BMI-adj analyses, respectively, corresponding to 21 association regions. Of these, 17 were not previously reported in association with OSA, and 13 replicated in FinnGen (False Discovery Rate p-value < 0.05). There were widespread significant differences in genetic effects between men and women, but less so across HARE groups. Meta-analysis of MVP, FinnGen, and MGB biobank revealed 17 additional, previously unreported, genome-wide significant regions.

INTERPRETATION

Sex differences in genetic associations with OSA are widespread, likely associated with multiple OSA risk factors. OSA shares genetic underpinnings with several sleep phenotypes, suggesting shared aetiology and causal pathways.

FUNDING

Described in acknowledgements.

Single Nucleotide Polymorphisms of the RAC1 Gene as Novel Susceptibility Markers for Neuropathy and Microvascular Complications in Type 2 Diabetes

Single nucleotide polymorphisms (SNP) in the RAC1 (Rac family small GTPase 1) gene have recently been linked to type 2 diabetes (T2D) and hyperglycemia due to their contribution to impaired redox homeostasis. The present study was designed to determine whether the common SNPs of the RAC1 gene are associated with diabetic complications such as neuropathy (DN), retinopathy (DR), nephropathy, angiopathy of the lower extremities (DA), and diabetic foot syndrome. A total of 1470 DNA samples from T2D patients were genotyped for six common SNPs by the MassArray Analyzer-4 system. The genotype rs7784465-T/C of RAC1 was associated with an increased risk of DR (p = 0.016) and DA (p = 0.03) in males, as well as with DR in females (p = 0.01). Furthermore, the SNP rs836478 showed an association with DR (p = 0.005) and DN (p = 0.025) in males, whereas the SNP rs10238136 was associated with DA in females (p = 0.002). In total, three RAC1 haplotypes showed significant associations (FDR < 0.05) with T2D complications in a sex-specific manner. The study's findings demonstrate, for the first time, that the RAC1 gene's polymorphisms represent novel and sex-specific markers of neuropathy and microvascular complications in type 2 diabetes, and that the gene could be a new target for the pharmacological inhibition of oxidative stress as a means of preventing diabetic complications.

Gene-Nutrient Interactions in Obesity: COBLL1 Genetic Variants Interact with Dietary Fat Intake to Modulate the Incidence of Obesity

The COBLL1 gene is associated with leptin, a hormone important for appetite and weight maintenance. Dietary fat is a significant factor in obesity. This study aimed to determine the association between COBLL1 gene, dietary fat, and incidence of obesity. Data from the Korean Genome and Epidemiology Study were used, and 3055 Korean adults aged ≥ 40 years were included. Obesity was defined as a body mass index ≥ 25 kg/m². Patients with obesity at baseline were excluded. The effects of the COBLL1 rs6717858 genotypes and dietary fat on incidence of obesity were evaluated using multivariable Cox proportional hazard models. During an average follow-up period of 9.2 years, 627 obesity cases were documented. In men, the hazard ratio (HR) for obesity was higher in CT, CC carriers (minor allele carriers) in the highest tertile of dietary fat intake than for men with TT carriers in the lowest tertile of dietary fat intake (Model 1: HR: 1.66, 95% confidence interval [CI]: 1.07-2.58; Model 2: HR: 1.63, 95% CI: 1.04-2.56). In women, the HR for obesity was higher in TT carriers in the highest tertile of dietary fat intake than for women with TT carriers in the lowest tertile of dietary fat intake (Model 1: HR: 1.49, 95% CI: 1.08-2.06; Model 2: HR: 1.53, 95% CI: 1.10-2.13). COBLL1 genetic variants and dietary fat intake had different sex-dependent effects in obesity. These results imply that a low-fat diet may protect against the effects of COBLL1 genetic variants on future obesity risk.

A case control study on the relationship between occupational stress and genetic polymorphism and dyslipidemia in coal miners

Dyslipidemia is one of the known risk factors for cardiovascular disease, and its prevalence is increasing worldwide. At present, the study of dyslipidemia has gradually shifted from simple environmental or genetic factors to environment-gene interactions. In order to further explore the etiology and mechanism of dyslipidemia, we used occupational stress(OS) and LYPLAL1, APOC3 and SOD2 gene as research variables to explore their association with dyslipidemia.Here we used a case-control study to include Han workers from a coal mining enterprise in China to determine the association between study variables and dyslipidemia. Monofactor analysis showed that smoking, drinking, physical activity level, DASH diet score, sleep quality, BMI, hypertension, hyperuricemia, shift work, OS were significantly different between the two groups (P < 0.05). In the APOC3 rs2854116 dominant model, patients with CT/CC genotype had a higher risk of dyslipidemia than those with TT genotype. In SOD2 rs4880 recessive model, patients with GG genotype had a lower risk of dyslipidemia than those with AA/AG genotype, and the difference was statistically significant. We found that rs12137855 and OS, rs2854116 and OS, rs4880 and OS had joint effects, but no interaction based on the multiplication and addition model was found (P_interaction > 0.05). GMDR model showed that the rs12137855-rs2854116-rs4880-OS four-factor model had the highest cross-validation consistency and training-validation accuracy (P < 0.05), suggesting that there was a high-order interaction between them associated with dyslipidemia. We found that dyslipidemia in coal miners was related to OS and genetic factors. Through this study, we revealed the dual regulation of environmental factors and genetic factors on dyslipidemia. At the same time, this study provides clues for understanding the etiology and mechanism of dyslipidemia.

Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Prenatal social support in low-risk pregnancy shapes placental epigenome

BACKGROUND

Poor social support during pregnancy has been linked to inflammation and adverse pregnancy and childhood health outcomes. Placental epigenetic alterations may underlie these links but are still unknown in humans.

METHODS

In a cohort of low-risk pregnant women (n = 301) from diverse ethnic backgrounds, social support was measured using the ENRICHD Social Support Inventory (ESSI) during the first trimester. Placental samples collected at delivery were analyzed for DNA methylation and gene expression using Illumina 450K Beadchip Array and RNA-seq, respectively. We examined association between maternal prenatal social support and DNA methylation in placenta. Associated cytosine-(phosphate)-guanine sites (CpGs) were further assessed for correlation with nearby gene expression in placenta.

RESULTS

The mean age (SD) of the women was 27.7 (5.3) years. The median (interquartile range) of ESSI scores was 24 (22-25). Prenatal social support was significantly associated with methylation level at seven CpGs (P_FDR < 0.05). The methylation levels at two of the seven CpGs correlated with placental expression of VGF and ILVBL (P_FDR < 0.05), genes known to be involved in neurodevelopment and energy metabolism. The genes annotated with the top 100 CpGs were enriched for pathways related to fetal growth, coagulation system, energy metabolism, and neurodevelopment. Sex-stratified analysis identified additional significant associations at nine CpGs in male-bearing pregnancies and 35 CpGs in female-bearing pregnancies.

CONCLUSIONS

The findings suggest that prenatal social support is linked to placental DNA methylation changes in a low-stress setting, including fetal sex-dependent epigenetic changes. Given the relevance of some of these changes in fetal neurodevelopmental outcomes, the findings signal important methylation targets for future research on molecular mechanisms of effect of the broader social environment on pregnancy and fetal outcomes.

TRIAL REGISTRATION

NCT00912132 ( ClinicalTrials.gov ).

X chromosome associations with chronic obstructive pulmonary disease and related phenotypes: an X chromosome-wide association study

BACKGROUND

The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations.

METHODS

Using X chromosome data from three COPD-enriched cohorts of adult smokers, we performed X chromosome specific quality control, imputation, and testing for association with COPD case-control status, lung function, and quantitative emphysema. Analyses were performed among all subjects, then stratified by sex, and subsequently combined in meta-analyses.

RESULTS

Among 10,193 subjects of non-Hispanic white or European ancestry, a variant near TMSB4X, rs5979771, reached genome-wide significance for association with lung function measured by FEV₁/FVC ([Formula: see text] 0.020, SE 0.004, p 4.97 × 10^-08), with suggestive evidence of association with FEV₁ ([Formula: see text] 0.092, SE 0.018, p 3.40 × 10^-07). Sex-stratified analyses revealed X chromosome variants that were differentially trending in one sex, with significantly different effect sizes or directions.

CONCLUSIONS

This investigation identified loci influencing lung function, COPD, and emphysema in a comprehensive genetic association meta-analysis of X chromosome genetic markers from multiple COPD-related datasets. Sex differences play an important role in the pathobiology of complex lung disease, including X chromosome variants that demonstrate differential effects by sex and variants that may be relevant through escape from X chromosome inactivation. Comprehensive interrogation of the X chromosome to better understand genetic control of COPD and lung function is important to further understanding of disease pathology. Trial registration Genetic Epidemiology of COPD Study (COPDGene) is registered at ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008). Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints Study (ECLIPSE), GlaxoSmithKline study code SCO104960, is registered at ClinicalTrials.gov, NCT00292552 (Active since February 16, 2006). Genetics of COPD in Norway Study (GenKOLS) holds GlaxoSmithKline study code RES11080, Genetics of Chronic Obstructive Lung Disease.

Fibrillin-1 and asprosin, novel players in metabolic syndrome

Fibrillin-1 is a major component of the extracellular microfibrils, where it interacts with other extracellular matrix proteins to provide elasticity to connective tissues, and regulates the bioavailability of TGFβ family members. A peptide consisting of the C-terminal 140 amino acids of fibrillin-1 has recently been identified as a glucogenic hormone, secreted from adipose tissue during fasting and targeting the liver to release glucose. This fragment, called asprosin, also signals in the hypothalamus to stimulate appetite. Asprosin levels are correlated with many of the pathologies indicative of metabolic syndrome, including insulin resistance and obesity. Previous studies and reviews have addressed the therapeutic potential of asprosin as a target in obesity, diabetes and related conditions without considering mechanisms underlying the relationship between generation of asprosin and expression of the much larger fibrillin-1 protein. Profibrillin-1 undergoes obligatory cleavage at the cell surface as part of its assembly into microfibrils, producing the asprosin peptide as well as mature fibrillin-1. Patterns of FBN1 mRNA expression are inconsistent with the necessity for regulated release of asprosin. The asprosin peptide may be protected from degradation in adipose tissue. We present evidence for an alternative possibility, that asprosin mRNA is generated independently from an internal promoter within the 3' end of the FBN1 gene, which would allow for regulation independent of fibrillin-synthesis and is more economical of cellular resources. The discovery of asprosin opened exciting possibilities for treatment of metabolic syndrome related conditions, but there is much to be understood before such therapies could be introduced into the clinic.

Strategies to identify causal common genetic variants and corresponding effector genes for paediatric obesity

BACKGROUND

Childhood obesity rates are on the rise, but there are currently no effective therapies available to slow or halt their progression. Although environmental and lifestyle factors have been implicated in its pathogenesis, childhood obesity is considered a complex disorder with a clear genetic component. Intense genome-wide association study (GWAS) efforts through large-scale collaborations have enabled the discovery of genetic loci robustly associated with childhood obesity beyond the classic FTO locus. That said, GWAS itself does not pinpoint the actual underlying causal effector genes, but rather just yields association signals in the genome.

OBJECTIVE

This review aims to outline what has been elucidated thus far on the genetic aetiology of commong childhood obesity and to describe strategies to identify and validate both causal common genetic variants and their corresponding effector genes.

RESULTS

Relevant cell types for molecular studies can be identified by gene set enrichment analysis and considering known biology of obesity-related physiological processes. Putatively causal single nucleotide polymorphisms (SNPs) can be identified by several methods including statistical fine mapping and 'assay for transposase accessible chromatin sequencing' (ATAC-seq). Variant to gene mapping can then nominate effector genes likely regulated by cis-regulatory elements harbouring putatively causal SNPs. A SNP's cis-regulatory activity can be functionally validated by several in vitro methods including luciferase assay and CRISPR approaches. These CRISPR approaches can also be used to investigate how dysregulatn of effector genes may confer obesity risk.

CONCLUSION

Uncovering the causative genes related to GWAS signals and elucidating their functional contributions to paediatric obesity with these strategies will deepen our understanding of this disease and serve better treatment outcomes.

Exploring the Genetic Association between Obesity and Serum Lipid Levels Using Bivariate Methods

It is crucial to understand the genetic mechanisms and biological pathways underlying the relationship between obesity and serum lipid levels. Structural equation models (SEMs) were constructed to calculate heritability for body mass index (BMI), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and the genetic connections between BMI and the four classes of lipids using 1197 pairs of twins from the Chinese National Twin Registry (CNTR). Bivariate genomewide association studies (GWAS) were performed to identify genetic variants associated with BMI and lipids using the records of 457 individuals, and the results were further validated in 289 individuals. The genetic background affecting BMI may differ by gender, and the heritability of males and females was 71% (95% CI [.66, .75]) and 39% (95% CI [.15, .71]) respectively. BMI was positively correlated with TC, TG and LDL-C in phenotypic and genetic correlation, while negatively correlated with HDL-C. There were gender differences in the correlation between BMI and lipids. Bivariate GWAS analysis and validation stage found 7 genes (LOC105378740, LINC02506, CSMD1, MELK, FAM81A, ERAL1 and MIR144) that were possibly related to BMI and lipid levels. The significant biological pathways were the regulation of cholesterol reverse transport and the regulation of high-density lipoprotein particle clearance (p < .001). BMI and blood lipid levels were affected by genetic factors, and they were genetically correlated. There might be gender differences in their genetic correlation. Bivariate GWAS analysis found MIR144 gene and its related biological pathways may influence obesity and lipid levels.

Insight into genetic, biological, and environmental determinants of sexual-dimorphism in type 2 diabetes and glucose-related traits

There is growing evidence that sex and gender differences play an important role in risk and pathophysiology of type 2 diabetes (T2D). Men develop T2D earlier than women, even though there is more obesity in young women than men. This difference in T2D prevalence is attenuated after the menopause. However, not all women are equally protected against T2D before the menopause, and gestational diabetes represents an important risk factor for future T2D. Biological mechanisms underlying sex and gender differences on T2D physiopathology are not yet fully understood. Sex hormones affect behavior and biological changes, and can have implications on lifestyle; thus, both sex-specific environmental and biological risk factors interact within a complex network to explain the differences in T2D risk and physiopathology in men and women. In addition, lifetime hormone fluctuations and body changes due to reproductive factors are generally more dramatic in women than men (ovarian cycle, pregnancy, and menopause). Progress in genetic studies and rodent models have significantly advanced our understanding of the biological pathways involved in the physiopathology of T2D. However, evidence of the sex-specific effects on genetic factors involved in T2D is still limited, and this gap of knowledge is even more important when investigating sex-specific differences during the life course. In this narrative review, we will focus on the current state of knowledge on the sex-specific effects of genetic factors associated with T2D over a lifetime, as well as the biological effects of these different hormonal stages on T2D risk. We will also discuss how biological insights from rodent models complement the genetic insights into the sex-dimorphism effects on T2D. Finally, we will suggest future directions to cover the knowledge gaps.

The effect of heteroscedasticity on the prediction efficiency of genome-wide polygenic score for body mass index

Globally, more than 1.9 billion adults are overweight. Thus, obesity is a serious public health issue. Moreover, obesity is a major risk factor for diabetes mellitus, coronary heart disease, and cardiovascular disease. Recently, GWAS examining obesity and body mass index (BMI) have increasingly unveiled many aspects of the genetic architecture of obesity and BMI. Information on genome-wide genetic variants has been used to estimate the genome-wide polygenic score (GPS) for a personalized prediction of obesity. However, the prediction power of GPS is affected by various factors, including the unequal variance in the distribution of a phenotype, known as heteroscedasticity. Here, we calculated a GPS for BMI using LDpred2, which was based on the BMI GWAS summary statistics from a European meta-analysis. Then, we tested the GPS in 354,761 European samples from the UK Biobank and found an effective prediction power of the GPS on BMI. To study a change in the variance of BMI, we investigated the heteroscedasticity of BMI across the GPS via graphical and statistical methods. We also studied the homoscedastic samples for BMI compared to the heteroscedastic sample, randomly selecting samples with various standard deviations of BMI residuals. Further, we examined the effect of the genetic interaction of GPS with environment (GPS×E) on the heteroscedasticity of BMI. We observed the changing variance (i.e., heteroscedasticity) of BMI along the GPS. The heteroscedasticity of BMI was confirmed by both the Breusch-Pagan test and the Score test. Compared to the heteroscedastic sample, the homoscedastic samples from small standard deviation of BMI residuals showed a decreased heteroscedasticity and an improved prediction accuracy, suggesting a quantitatively negative correlation between the phenotypic heteroscedasticity and the prediction accuracy of GPS. To further test the effects of the GPS×E on heteroscedasticity, first we tested the genetic interactions of the GPS with 21 environments and found 8 significant GPS×E interactions on BMI. However, the heteroscedasticity of BMI was not ameliorated after adjusting for the GPS×E interactions. Taken together, our findings suggest that the heteroscedasticity of BMI exists along the GPS and is not affected by the GPS×E interaction.

Polygenic signals of sex differences in selection in humans from the UK Biobank

Sex differences in the fitness effects of genetic variants can influence the rate of adaptation and the maintenance of genetic variation. For example, "sexually antagonistic" (SA) variants, which are beneficial for one sex and harmful for the other, can both constrain adaptation and increase genetic variability for fitness components such as survival, fertility, and disease susceptibility. However, detecting variants with sex-differential fitness effects is difficult, requiring genome sequences and fitness measurements from large numbers of individuals. Here, we develop new theory for studying sex-differential selection across a complete life cycle and test our models with genotypic and reproductive success data from approximately 250,000 UK Biobank individuals. We uncover polygenic signals of sex-differential selection affecting survival, reproductive success, and overall fitness, with signals of sex-differential reproductive selection reflecting a combination of SA polymorphisms and sexually concordant polymorphisms in which the strength of selection differs between the sexes. Moreover, these signals hold up to rigorous controls that minimise the contributions of potential confounders, including sequence mapping errors, population structure, and ascertainment bias. Functional analyses reveal that sex-differentiated sites are enriched in phenotype-altering genomic regions, including coding regions and loci affecting a range of quantitative traits. Population genetic analyses show that sex-differentiated sites exhibit evolutionary histories dominated by genetic drift and/or transient balancing selection, but not long-term balancing selection, which is consistent with theoretical predictions of effectively weak SA balancing selection in historically small populations. Overall, our results are consistent with polygenic sex-differential-including SA-selection in humans. Evidence for sex-differential selection is particularly strong for variants affecting reproductive success, in which the potential contributions of nonrandom sampling to signals of sex differentiation can be excluded.

Sex Differences in Adiposity and Cardiovascular Diseases

Body fat distribution is a well-established predictor of adverse medical outcomes, independent of overall adiposity. Studying body fat distribution sheds insights into the causes of obesity and provides valuable information about the development of various comorbidities. Compared to total adiposity, body fat distribution is more closely associated with risks of cardiovascular diseases. The present review specifically focuses on the sexual dimorphism in body fat distribution, the biological clues, as well as the genetic traits that are distinct from overall obesity. Understanding the sex determinations on body fat distribution and adiposity will aid in the improvement of the prevention and treatment of cardiovascular diseases (CVD).

Limited evidence for blood eQTLs in human sexual dimorphism

BACKGROUND

The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits.

METHODS

To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs.

RESULTS

Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection.

CONCLUSIONS

Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs.

Adipose Tissue Dysfunction and Obesity-Related Male Hypogonadism

Obesity is a chronic illness associated with several metabolic derangements and comorbidities (i.e., insulin resistance, leptin resistance, diabetes, etc.) and often leads to impaired testicular function and male subfertility. Several mechanisms may indeed negatively affect the hypothalamic–pituitary–gonadal health, such as higher testosterone conversion to estradiol by aromatase activity in the adipose tissue, increased ROS production, and the release of several endocrine molecules affecting the hypothalamus–pituitary–testis axis by both direct and indirect mechanisms. In addition, androgen deficiency could further accelerate adipose tissue expansion and therefore exacerbate obesity, which in turn enhances hypogonadism, thus inducing a vicious cycle. Based on these considerations, we propose an overview on the relationship of adipose tissue dysfunction and male hypogonadism, highlighting the main biological pathways involved and the current therapeutic options to counteract this condition.

Exploring and mitigating potential bias when genetic instrumental variables are associated with multiple non-exposure traits in Mendelian randomization

With the increasing size and number of genome-wide association studies, individual single nucleotide polymorphisms are increasingly found to associate with multiple traits. Many different mechanisms could result in proposed genetic IVs for an exposure of interest being associated with multiple non-exposure traits, some of which could bias MR results. We describe and illustrate, through causal diagrams, a range of scenarios that could result in proposed IVs being related to non-exposure traits in MR studies. These associations could occur due to five scenarios: (i) confounding, (ii) vertical pleiotropy, (iii) horizontal pleiotropy, (iv) reverse causation and (v) selection bias. For each of these scenarios we outline steps that could be taken to explore the underlying mechanism and mitigate any resulting bias in the MR estimation. We recommend MR studies explore possible IV-non-exposure associations across a wider range of traits than is usually the case. We highlight the pros and cons of relying on sensitivity analyses without considering particular pleiotropic paths versus systematically exploring and controlling for potential pleiotropic or other biasing paths via known traits. We apply our recommendations to an illustrative example of the effect of maternal insomnia on offspring birthweight in UK Biobank.

Inherited basis of visceral, abdominal subcutaneous and gluteofemoral fat depots

For any given level of overall adiposity, individuals vary considerably in fat distribution. The inherited basis of fat distribution in the general population is not fully understood. Here, we study up to 38,965 UK Biobank participants with MRI-derived visceral (VAT), abdominal subcutaneous (ASAT), and gluteofemoral (GFAT) adipose tissue volumes. Because these fat depot volumes are highly correlated with BMI, we additionally study six local adiposity traits: VAT adjusted for BMI and height (VATadj), ASATadj, GFATadj, VAT/ASAT, VAT/GFAT, and ASAT/GFAT. We identify 250 independent common variants (39 newly-identified) associated with at least one trait, with many associations more pronounced in female participants. Rare variant association studies extend prior evidence for PDE3B as an important modulator of fat distribution. Local adiposity traits (1) highlight depot-specific genetic architecture and (2) enable construction of depot-specific polygenic scores that have divergent associations with type 2 diabetes and coronary artery disease. These results - using MRI-derived, BMI-independent measures of local adiposity - confirm fat distribution as a highly heritable trait with important implications for cardiometabolic health outcomes.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Causality of anthropometric markers associated with polycystic ovarian syndrome: Findings of a Mendelian randomization study

Introduction

Using body mass index (BMI) as a proxy, previous Mendelian randomization (MR) studies found total causal effects of general obesity on polycystic ovarian syndrome (PCOS). Hitherto, total and direct causal effects of general- and central obesity on PCOS have not been comprehensively analyzed.

Objectives

To investigate the causality of central- and general obesity on PCOS using surrogate anthropometric markers.

Methods

Summary GWAS data of female-only, large-sample cohorts of European ancestry were retrieved for anthropometric markers of central obesity (waist circumference (WC), hip circumference (HC), waist-to-hip ratio (WHR)) and general obesity (BMI and its constituent variables–weight and height), from the IEU Open GWAS Project. As the outcome, we acquired summary data from a large-sample GWAS (118870 samples; 642 cases and 118228 controls) within the FinnGen cohort. Total causal effects were assessed via univariable two-sample Mendelian randomization (2SMR). Genetic architectures underlying causal associations were explored. Direct causal effects were analyzed by multivariable MR modelling.

Results

Instrumental variables demonstrated no weak instrument bias (F > 10). Four anthropometric exposures, namely, weight (2.69–77.05), BMI (OR: 2.90–4.06), WC (OR: 6.22–20.27), and HC (OR: 6.22–20.27) demonstrated total causal effects as per univariable 2SMR models. We uncovered shared and non-shared genetic architectures underlying causal associations. Direct causal effects of WC and HC on PCOS were revealed by two multivariable MR models containing exclusively the anthropometric markers of central obesity. Other multivariable MR models containing anthropometric markers of both central- and general obesity showed no direct causal effects on PCOS.

Conclusions

Both and general- and central obesity yield total causal effects on PCOS. Findings also indicated potential direct causal effects of normal weight-central obesity and more complex causal mechanisms when both central- and general obesity are present. Results underscore the importance of addressing both central- and general obesity for optimizing PCOS care.

A Personalized Diet Approach Study: Interaction between PPAR-γ Pro12Ala and Dietary Insulin Indices on Metabolic Markers in Diabetic Patients

BACKGROUND

The objectives were to investigate the effect of the interaction between peroxisome proliferator-activated receptor gamma (PPAR-γ) Pro12Ala polymorphisms and dietary insulin load and insulin index (DIL and DII) on Cardio-metabolic Markers among diabetic patients.

METHODS

This cross-sectional study was conducted on 393 diabetic patients. Food-frequency questionnaire (FFQ) was used for DIL and DII calculation. PPAR-γ Pro12Ala was genotyped by the PCR-RFLP method. Biochemical markers including TC, LDL, HDL, TG, SOD, CRP, TAC, PTX3, PGF2α. IL18, leptin and ghrelin were measured by standard protocol.

RESULT

Risk-allele carriers (CG, GG) had higher obesity indices WC (P _interaction =0.04), BMI (P _interaction =0.006) and, WC (P _interaction =0.04) compared with individuals with the CC genotype when they consumed a diet with higher DIL and DII respectively. Besides, carriers of the G allele who were in the highest tertile of DIL, had lower HDL (P _interaction =0.04) and higher PGF2α (P _interaction =0.03) and PTX3 (P _interaction =0.03). Moreover, the highest tertile of the DII, showed an increase in IL18 (P _interaction =0.01) and lower SOD (P _interaction =0.03) for risk allele carriers compared to those with CC homozygotes.

CONCLUSION

We revealed PPAR-γ Pro12Ala polymorphism was able to intensify the effect of DIL and DII on CVD risk factors; risk-allele carriers who consumed a diet with high DIL and DII score have more likely to be obese and have higher inflammatory markers. Also, protective factor against CVD risk factors were reduced significantly in this group compared to CC homozygotes. This article is protected by copyright. All rights reserved.

Identification of genetic effects underlying type 2 diabetes in South Asian and European populations

South Asians are at high risk of developing type 2 diabetes (T2D). We carried out a genome-wide association meta-analysis with South Asian T2D cases (n = 16,677) and controls (n = 33,856), followed by combined analyses with Europeans (neff = 231,420). We identify 21 novel genetic loci for significant association with T2D (P = 4.7 × 10-8 to 5.2 × 10-12), to the best of our knowledge at the point of analysis. The loci are enriched for regulatory features, including DNA methylation and gene expression in relevant tissues, and highlight CHMP4B, PDHB, LRIG1 and other genes linked to adiposity and glucose metabolism. A polygenic risk score based on South Asian-derived summary statistics shows ~4-fold higher risk for T2D between the top and bottom quartile. Our results provide further insights into the genetic mechanisms underlying T2D, and highlight the opportunities for discovery from joint analysis of data from across ancestral populations.

Identification of genetic loci simultaneously associated with multiple cardiometabolic traits

BACKGROUND AND AIMS

Cardiometabolic disorders (CMD) arise from a constellation of features such as increased adiposity, hyperlipidemia, hypertension and compromised glucose control. Many genetic loci have shown associations with individual CMD-related traits, but no investigations have focused on simultaneously identifying loci showing associations across all domains. We therefore sought to identify loci associated with risk across seven continuous CMD-related traits.

METHODS AND RESULTS

We conducted separate genome-wide association studies (GWAS) for systolic and diastolic blood pressure (SBP/DBP), hemoglobin A1c (HbA1c), low- and high- density lipoprotein cholesterol (LDL-C/HDL-C), waist-to-hip-ratio (WHR), and triglycerides (TGs) in the UK Biobank (N = 356,574-456,823). Multiple loci reached genome-wide levels of significance (N = 145-333) for each trait, but only four loci (in/near VEGFA, GRB14-COBLL1, KLF14, and RGS19-OPRL1) were associated with risk across all seven traits (P < 5 × 10^-8). We sought replication of these four loci in an independent set of seven trait-specific GWAS meta-analyses. GRB14-COBLL1 showed the most consistent replication, revealing nominally significant associations (P < 0.05) with all traits except DBP.

CONCLUSIONS

Our analyses suggest that very few loci are associated in the same direction of risk with traits representing the full spectrum of CMD features. We identified four such loci, and an understanding of the pathways between these loci and CMD risk may eventually identify factors that can be used to identify pathologic disturbances that represent broadly beneficial therapeutic targets.

Genetic pleiotropy between pulmonary function and age-related traits: The Long Life Family Study

BACKGROUND

Pulmonary function (PF) progressively declines with aging. Forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) are predictors of morbidity of pulmonary and cardiovascular diseases and all-cause mortality. In addition, reduced PF is associated with elevated chronic low-grade systemic inflammation, glucose metabolism, body fatness, and low muscle strength. It may suggest pleiotropic genetic effects between PF with these age-related factors.

METHODS

We evaluated whether FEV1 and FVC share common pleiotropic genetic effects factors with interleukin-6, high-sensitivity C-reactive protein, body mass index, muscle (grip) strength, plasma glucose, and glycosylated hemoglobin in 3,888 individuals (age range: 26-106). We employed sex-combined and sex-specific correlated meta-analyses to test whether combining genome-wide association p-values from two or more traits enhances the ability to detect variants sharing effects on these correlated traits.

RESULTS

We identified 32 loci for PF, including 29 novel pleiotropic loci associated with pulmonary function and (i) body fatness (CYP2U1/SGMS2), (ii) glucose metabolism (CBWD1/DOCK8 and MMUT/CENPQ), (iii) inflammatory markers (GLRA3/HPGD, TRIM9, CALN1, CTNNB1/ZNF621, GATA5/SLCO4A1/NTSR1, and NPVF/C7orf31/CYCS), and (iv) muscle strength (MAL2, AC008825.1/LINC02103, AL136418.1).

CONCLUSIONS

The identified genes/loci for PF and age-related traits suggest their underlying shared genetic effects, which can explain part of their phenotypic correlations. Integration of gene expression and genomic annotation data shows enrichment of our genetic variants in lung, blood, adipose, pancreas, and muscles, among others. Our findings highlight the critical roles of identified gene/locus in systemic inflammation, glucose metabolism, strength performance, PF, and pulmonary disease, which are involved in accelerated biological aging.

Genetic associations with carotid intima-media thickness link to atherosclerosis with sex-specific effects in sub-Saharan Africans

Atherosclerosis precedes the onset of clinical manifestations of cardiovascular diseases (CVDs). We used carotid intima-media thickness (cIMT) to investigate genetic susceptibility to atherosclerosis in 7894 unrelated adults (3963 women, 3931 men; 40 to 60 years) resident in four sub-Saharan African countries. cIMT was measured by ultrasound and genotyping was performed on the H3Africa SNP Array. Two new African-specific genome-wide significant loci for mean-max cIMT, SIRPA (p = 4.7E-08), and FBXL17 (p = 2.5E-08), were identified. Sex-stratified analysis revealed associations with one male-specific locus, SNX29 (p = 6.3E-09), and two female-specific loci, LARP6 (p = 2.4E-09) and PROK1 (p = 1.0E-08). We replicate previous cIMT associations with different lead SNPs in linkage disequilibrium with SNPs primarily identified in European populations. Our study find significant enrichment for genes involved in oestrogen response from female-specific signals. The genes identified show biological relevance to atherosclerosis and/or CVDs, sex-differences and transferability of signals from non-African studies.