Autosomal genetic control of human gene expression does not differ across the sexes

Leveraging sex-genetic interactions to understand brain disorders: recent advances and current gaps

It is established that there are sex differences in terms of prevalence, age of onset, clinical manifestations, and response to treatment for a variety of brain disorders, including neurodevelopmental, psychiatric, and neurodegenerative disorders. Cohorts of increasing sample sizes with diverse data types collected, including genetic, transcriptomic and/or phenotypic data, are providing the building blocks to permit analytical designs to test for sex-biased genetic variant-trait associations, and for sex-biased transcriptional regulation. Such molecular assessments can contribute to our understanding of the manifested phenotypic differences between the sexes for brain disorders, offering the future possibility of delivering personalized therapy for females and males. With the intention of raising the profile of this field as a research priority, this review aims to shed light on the importance of investigating sex-genetic interactions for brain disorders, focusing on two areas: (i) variant-trait associations and (ii) transcriptomics (i.e. gene expression, transcript usage and regulation). We specifically discuss recent advances in the field, current gaps and provide considerations for future studies.

Genome-wide association study of varicose veins identifies a protective missense variant in GJD3 enriched in the Finnish population

Varicose veins is the most common manifestation of chronic venous disease that displays female-biased incidence. To identify protein-inactivating variants that could guide identification of drug target genes for varicose veins and genetic evidence for the disease prevalence difference between the sexes, we conducted a genome-wide association study of varicose veins in Finns using the FinnGen dataset with 17,027 cases and 190,028 controls. We identified 50 associated genetic loci (P < 5.0 × 10-8) of which 29 were novel including one near ERG with female-specificity (rs2836405-G, OR[95% CI] = 1.09[1.05-1.13], P = 3.1 × 10-8). These also include two X-chromosomal (ARHGAP6 and SRPX) and two autosomal novel loci (TGFB2 and GJD3) with protein-coding lead variants enriched above 56-fold in Finns over non-Finnish non-Estonian Europeans. A low-frequency missense variant in GJD3 (p.Pro59Thr) is exclusively associated with a lower risk for varicose veins (OR = 0.62 [0.55-0.70], P = 1.0 × 10-14) in a phenome-wide scan of the FinnGen data. The absence of observed pleiotropy and its membership of the connexin gene family underlines GJD3 as a potential connexin-modulating therapeutic strategy for varicose veins. Our results provide insights into varicose veins etiopathology and highlight the power of isolated populations, including Finns, to discover genetic variants that inform therapeutic development.

Effect of all-but-one conditional analysis for eQTL isolation in peripheral blood

Expression quantitative trait locus detection has become increasingly important for understanding how noncoding variants contribute to disease susceptibility and complex traits. The major challenges in expression quantitative trait locus fine-mapping and causal variant discovery relate to the impact of linkage disequilibrium on signals due to one or multiple functional variants that lie within a credible set. We perform expression quantitative trait locus fine-mapping using the all-but-one approach, conditioning each signal on all others detected in an interval, on the Consortium for the Architecture of Gene Expression cohorts of microarray-based peripheral blood gene expression in 2,138 European-ancestry human adults. We contrast these results with traditional forward stepwise conditional analysis and a Bayesian localization method. All-but-one conditioning significantly modifies effect-size estimates for 51% of 2,351 expression quantitative trait locus peaks, but only modestly affects credible set size and location. On the other hand, both conditioning approaches result in unexpectedly low overlap with Bayesian credible sets, with just 57% peak concordance and between 50% and 70% SNP sharing, leading us to caution against the assumption that any one localization method is superior to another. We also cross reference our results with ATAC-seq data, cell-type-specific expression quantitative trait locus, and activity-by-contact-enhancers, leading to the proposal of a 5-tier approach to further reduce credible set sizes and prioritize likely causal variants for all known inflammatory bowel disease risk loci active in immune cells.

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.

Sex differences in genetic architecture in the UK Biobank

Males and females present differences in complex traits and in the risk of a wide array of diseases. Genotype by sex (GxS) interactions are thought to account for some of these differences. However, the extent and basis of GxS are poorly understood. In the present study, we provide insights into both the scope and the mechanism of GxS across the genome of about 450,000 individuals of European ancestry and 530 complex traits in the UK Biobank. We found small yet widespread differences in genetic architecture across traits. We also found that, in some cases, sex-agnostic analyses may be missing trait-associated loci and looked into possible improvements in the prediction of high-level phenotypes. Finally, we studied the potential functional role of the differences observed through sex-biased gene expression and gene-level analyses. Our results suggest the need to consider sex-aware analyses for future studies to shed light onto possible sex-specific molecular mechanisms.

Genetic factors contributing to extensive variability of sex-specific hepatic gene expression in Diversity Outbred mice

Sex-specific transcription characterizes hundreds of genes in mouse liver, many implicated in sex-differential drug and lipid metabolism and disease susceptibility. While the regulation of liver sex differences by growth hormone-activated STAT5 is well established, little is known about autosomal genetic factors regulating the sex-specific liver transcriptome. Here we show, using genotyping and expression data from a large population of Diversity Outbred mice, that genetic factors work in tandem with growth hormone to control the individual variability of hundreds of sex-biased genes, including many long non-coding RNA genes. Significant associations between single nucleotide polymorphisms and sex-specific gene expression were identified as expression quantitative trait loci (eQTLs), many of which showed strong sex-dependent associations. Remarkably, autosomal genetic modifiers of sex-specific genes were found to account for more than 200 instances of gain or loss of sex-specificity across eight Diversity Outbred mouse founder strains. Sex-biased STAT5 binding sites and open chromatin regions with strain-specific variants were significantly enriched at eQTL regions regulating correspondingly sex-specific genes, supporting the proposed functional regulatory nature of the eQTL regions identified. Binding of the male-biased, growth hormone-regulated repressor BCL6 was most highly enriched at trans-eQTL regions controlling female-specific genes. Co-regulated gene clusters defined by overlapping eQTLs included sets of highly correlated genes from different chromosomes, further supporting trans-eQTL action. These findings elucidate how an unexpectedly large number of autosomal factors work in tandem with growth hormone signaling pathways to regulate the individual variability associated with sex differences in liver metabolism and disease.

Genome-Wide Sex and Gender Differences in Cancer

Despite their known importance in clinical medicine, differences based on sex and gender are among the least studied factors affecting cancer susceptibility, progression, survival, and therapeutic response. In particular, the molecular mechanisms driving sex differences are poorly understood and so most approaches to precision medicine use mutational or other genomic data to assign therapy without considering how the sex of the individual might influence therapeutic efficacy. The mandate by the National Institutes of Health that research studies include sex as a biological variable has begun to expand our understanding on its importance. Sex differences in cancer may arise due to a combination of environmental, genetic, and epigenetic factors, as well as differences in gene regulation, and expression. Extensive sex differences occur genome-wide, and ultimately influence cancer biology and outcomes. In this review, we summarize the current state of knowledge about sex-specific genetic and genome-wide influences in cancer, describe how differences in response to environmental exposures and genetic and epigenetic alterations alter the trajectory of the disease, and provide insights into the importance of integrative analyses in understanding the interplay of sex and genomics in cancer. In particular, we will explore some of the emerging analytical approaches, such as the use of network methods, that are providing a deeper understanding of the drivers of differences based on sex and gender. Better understanding these complex factors and their interactions will improve cancer prevention, treatment, and outcomes for all individuals.

Tissue-specific sex differences in human gene expression

Despite extensive sex differences in human complex traits and disease, the male and female genomes differ only in the sex chromosomes. This implies that most sex-differentiated traits are the result of differences in the expression of genes that are common to both sexes. While sex differences in gene expression have been observed in a range of different tissues, the biological mechanisms for tissue-specific sex differences (TSSDs) in gene expression are not well understood. A total of 30 640 autosomal and 1021 X-linked transcripts were tested for heterogeneity in sex difference effect sizes in n = 617 individuals across 40 tissue types in Genotype-Tissue Expression (GTEx). This identified 65 autosomal and 66 X-linked TSSD transcripts (corresponding to unique genes) at a stringent significance threshold. Results for X-linked TSSD transcripts showed mainly concordant direction of sex differences across tissues and replicate previous findings. Autosomal TSSD transcripts had mainly discordant direction of sex differences across tissues. The top cis-expression quantitative trait loci (eQTLs) across tissues for autosomal TSSD transcripts are located a similar distance away from the nearest androgen and estrogen binding motifs and the nearest enhancer, as compared to cis-eQTLs for transcripts with stable sex differences in gene expression across tissue types. Enhancer regions that overlap top cis-eQTLs for TSSD transcripts, however, were found to be more dispersed across tissues. These observations suggest that androgen and estrogen regulatory elements in a cis region may play a common role in sex differences in gene expression, but TSSD in gene expression may additionally be due to causal variants located in tissue-specific enhancer regions.

Genetic constraints on microevolutionary divergence of sex-biased gene expression

The evolution of sex-specific phenotypes is an important dimension of diversification and local adaptation. The sex-dependent regulation of gene expression is considered a key genomic mechanism facilitating sex-dependent adaptation. In many species, genes with male-biased expression evolve faster in DNA sequence and expression level than genes with female-biased or sexually monomorphic expression. While positive selection may be responsible for rapid DNA sequence evolution, why expression of male-biased genes also evolves rapidly remains unclear. Beyond sex differences in selection, some aspects of the genetic architecture of gene expression could contribute to the rapid evolution of male-biased gene expression. First, male-biased genes might simply have greater standing genetic variance than female-biased genes. Second, male-biased genes could be less constrained by pleiotropy, either within or between sexes. Here, we evaluate these alternative explanations on an intraspecific scale using a series of quantitative genetic experiments conducted on natural variation in male and female gene expression in the fly Drosophila serrata Male-biased genes had significantly higher genetic variance than female-biased genes and were generally more narrowly expressed across tissues, suggesting lower within-individual pleiotropy. However, consistent with stronger constraints due to between-sex pleiotropy, their between-sex genetic correlations, r_MF, were higher than for female-biased genes and more strongly negatively associated with sex bias. Using an extensive clinal dataset, we tested whether sex differences in gene expression divergence among populations have been shaped by pleiotropy. Here too, male-biased gene divergence was more strongly associated with between-sex pleiotropy than was female-biased gene divergence. Systematic differences in genetic variance and pleiotropy may be important factors influencing sex-specific adaptation arising through changes in gene expression.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Conserved properties of genetic architecture of renal and fat transcriptomes in rat models of insulin resistance

To define renal molecular mechanisms that are affected by permanent hyperglycaemia and might promote phenotypes relevant to diabetic nephropathy, we carried out linkage analysis of genome-wide gene transcription in the kidneys of F2 offspring from the Goto-Kakizaki (GK) rat model of type 2 diabetes and normoglycaemic Brown Norway (BN) rats. We mapped 2526 statistically significant expression quantitative trait loci (eQTLs) in the cross. More than 40% of eQTLs mapped in the close vicinity of the linked transcripts, underlying possible cis-regulatory mechanisms of gene expression. We identified eQTL hotspots on chromosomes 5 and 9 regulating the expression of 80-165 genes, sex or cross direction effects, and enriched metabolic and immunological processes by segregating GK alleles. Comparative analysis with adipose tissue eQTLs in the same cross showed that 496 eQTLs, in addition to the top enriched biological pathways, are conserved in the two tissues. Extensive similarities in eQTLs mapped in the GK rat and in the spontaneously hypertensive rat (SHR) suggest a common aetiology of disease phenotypes common to the two strains, including insulin resistance, which is a prominent pathophysiological feature in both GK rats and SHRs. Our data shed light on shared and tissue-specific molecular mechanisms that might underlie aetiological aspects of insulin resistance in the context of spontaneously occurring hyperglycaemia and hypertension.

Summary: Kidney and fat expression QTL mapping in rat models of spontaneously occurring insulin resistance associated with either diabetes or hypertension reveals conserved gene expression regulation, suggesting shared aetiology of disease phenotypes.

Investigating the relationship between iron and depression

Lower levels of circulating iron have been associated with depression. Our objective was to investigate the phenotypic and genetic relationship between measures of circulating levels of iron (serum iron, transferrin, transferrin saturation, and ferritin) and depressive symptoms. Data were available from ongoing studies at QIMR Berghofer Medical Research Institute (QIMRB), including twin adolescents (mean age 15.1 years, standard deviation (SD) 3.2 years), and twin adults (mean age 23.2 years, SD 2.2 years). In the adolescent cohort, there were 3416 participants from 1688 families. In the adult cohort there were 9035 participants from 4533 families. We estimated heritabilities of, and phenotypic and genetic correlations between, traits. We conducted analyses that linked results from published large-scale genome-wide association studies (including iron and Major Depressive Disorder) with our study samples using single SNP and multi-SNP genetic risk score analyses, and LD score regression analyses. In both cohorts, measures of iron, transferrin, transferrin saturation, and log 10 of ferritin (L10Fer) were all highly heritable, while depressive measures were moderately heritable. In adolescents, depression measures were higher in those in the middle 10th versus top 10th percentile of transferrin saturation measures (p = 0.002). Genetic profile risk scores of the iron measures were not significantly associated with depression in study participants. LD score analyses showed no significant genetic relationship between iron and depression. Genetic factors strongly influence iron measures in adolescents and adults. Using several different strategies we find no evidence for a genetic contribution to the relationship between blood measures of iron and measures of depression.

The landscape of sex-differential transcriptome and its consequent selection in human adults

BACKGROUND

The prevalence of several human morbid phenotypes is sometimes much higher than intuitively expected. This can directly arise from the presence of two sexes, male and female, in one species. Men and women have almost identical genomes but are distinctly dimorphic, with dissimilar disease susceptibilities. Sexually dimorphic traits mainly result from differential expression of genes present in both sexes. Such genes can be subject to different, and even opposing, selection constraints in the two sexes. This can impact human evolution by differential selection on mutations with dissimilar effects on the two sexes.

RESULTS

We comprehensively mapped human sex-differential genetic architecture across 53 tissues. Analyzing available RNA-sequencing data from 544 adults revealed thousands of genes differentially expressed in the reproductive tracts and tissues common to both sexes. Sex-differential genes are related to various biological systems, and suggest new insights into the pathophysiology of diverse human diseases. We also identified a significant association between sex-specific gene transcription and reduced selection efficiency and accumulation of deleterious mutations, which might affect the prevalence of different traits and diseases. Interestingly, many of the sex-specific genes that also undergo reduced selection efficiency are essential for successful reproduction in men or women. This seeming paradox might partially explain the high incidence of human infertility.

CONCLUSIONS

This work provides a comprehensive overview of the sex-differential transcriptome and its importance to human evolution and human physiology in health and in disease.