Toward a deeper understanding of gene-by-sex interaction models.

In this issue of Cell Genomics, Zhu et al.1 propose amplification as the primary mode for gene-by-sex interactions in complex traits. Khramtsova et al. preview their modeling approach and discuss implications for the future work on the genomics of sex differences.

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.

A deep transcriptome meta-analysis reveals sex differences in multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS), a chronic auto-immune, inflammatory, and degenerative disease of the central nervous system, affects both males and females; however, females suffer from a higher risk of developing MS (2-3:1 ratio relative to males). The precise sex-based factors influencing risk of MS are currently unknown. Here, we explore the role of sex in MS to identify molecular mechanisms underlying observed MS sex differences that may guide novel therapeutic approaches tailored for males or females. METHODS: We performed a rigorous and systematic review of genome-wide transcriptome studies of MS that included patient sex data in the Gene Expression Omnibus and ArrayExpress databases following PRISMA statement guidelines. For each selected study, we analyzed differential gene expression to explore the impact of the disease in females (IDF), in males (IDM) and our main goal: the sex differential impact of the disease (SDID). Then, for each scenario (IDF, IDM and SDID) we performed 2 meta-analyses in the main tissues involved in the disease (brain and blood). Finally, we performed a gene set analysis in brain tissue, in which a higher number of genes were dysregulated, to characterize sex differences in biological pathways. RESULTS: After screening 122 publications, the systematic review provided a selection of 9 studies (5 in blood and 4 in brain tissue) with a total of 474 samples (189 females with MS and 109 control females; 82 males with MS and 94 control males). Blood and brain tissue meta-analyses identified, respectively, 1 (KIR2DL3) and 13 (ARL17B, CECR7, CEP78, IFFO2, LOC401127, NUDT18, RNF10, SLC17A5, STMP1, TRAF3IP2-AS1, UBXN2B, ZNF117, ZNF488) MS-associated genes that differed between males and females (SDID comparison). Functional analyses in the brain revealed different altered immune patterns in females and males (IDF and IDM comparisons). The pro-inflammatory environment and innate immune responses related to myeloid lineage appear to be more affected in females, while adaptive responses associated with the lymphocyte lineage in males. Additionally, females with MS displayed alterations in mitochondrial respiratory chain complexes, purine, and glutamate metabolism, while MS males displayed alterations in stress response to metal ion, amine, and amino acid transport. CONCLUSION: We found transcriptomic and functional differences between MS males and MS females (especially in the immune system), which may support the development of new sex-based research of this disease. Our study highlights the importance of understanding the role of biological sex in MS to guide a more personalized medicine.

Deciphering genetic causes for sex differences in human health through drug metabolism and transporter genes.

Sex differences have been widely observed in human health. However, little is known about the underlying mechanism behind these observed sex differences. We hypothesize that sex-differentiated genetic effects are contributors of these phenotypic differences. Focusing on a collection of drug metabolism enzymes and transporters (DMET) genes, we discover sex-differentiated genetic regulatory mechanisms between these genes and human complex traits. Here, we show that sex-differentiated genetic effects were present at genome-level and at DMET gene regions for many human complex traits. These sex-differentiated regulatory mechanisms are reflected in the levels of gene expression and endogenous serum biomarkers. Through Mendelian Randomization analysis, we identify putative sex-differentiated causal effects in each sex separately. Furthermore, we identify and validate sex differential gene expression of a subset of DMET genes in human liver samples. We observe higher protein abundance and enzyme activity of CYP1A2 in male-derived liver microsomes, which leads to higher level of an active metabolite formation of clozapine, a commonly prescribed antipsychotic drug. Taken together, our results demonstrate the presence of sex-differentiated genetic effects on DMET gene regulation, which manifest in various phenotypic traits including disease risks and drug responses.

Sex dimorphism in response to targeted therapy and immunotherapy in non-small cell lung cancer patients: a narrative review.

Background and Objective: Multiple agents have been developed for treating non-small cell lung cancer (NSCLC). However, patients' response to these therapies vary drastically, which indicates a need to tailor therapy. Sex is a readily usable clinical characteristic that has been shown to impact patients' response to drugs. The main objective of this narrative review is to summarize the current state of knowledge, compiled from meta-analyses, on sex differences in treatment efficacy for targeted therapy and immunotherapy in NSCLC. We discuss the interplay of patient characteristics, both molecular and demographic, with sex on how they impact therapeutic response. Methods: PubMed search was performed with the term "sex/gender differences" with currently FDA approved targeting therapy and immunotherapy agents in treating NSCLC. Key Content and Findings: For targeted therapy, women tend to benefit more in terms of progression-free survival upon receiving first-generation anti-epidermal growth factor receptor (EGFR) treatment than men. On the other hand, there is an ongoing debate on sex differences in response to immunotherapy. Although preliminary, whether sex differences were observed depends on treatment settings, patient characteristics, and molecular features. Importantly, incorporating sex as a biological component in the biomarker discovery seems to reveal novel insights in immunotherapy response. Conclusions: Taken together, sex differences in responding to standard care have been observed in clinical settings for NSCLC patients. A better understanding of sex-associated treatment response and the underlying biology will improve cancer prognosis and eliminate these sex differences.

Systematic evaluation of transcriptomics-based deconvolution methods and references using thousands of clinical samples.

Estimating cell type composition of blood and tissue samples is a biological challenge relevant in both laboratory studies and clinical care. In recent years, a number of computational tools have been developed to estimate cell type abundance using gene expression data. Although these tools use a variety of approaches, they all leverage expression profiles from purified cell types to evaluate the cell type composition within samples. In this study, we compare 12 cell type quantification tools and evaluate their performance while using each of 10 separate reference profiles. Specifically, we have run each tool on over 4000 samples with known cell type proportions, spanning both immune and stromal cell types. A total of 12 of these represent in vitro synthetic mixtures and 300 represent in silico synthetic mixtures prepared using single-cell data. A final 3728 clinical samples have been collected from the Framingham cohort, for which cell populations have been quantified using electrical impedance cell counting. When tools are applied to the Framingham dataset, the tool Estimating the Proportions of Immune and Cancer cells (EPIC) produces the highest correlation, whereas Gene Expression Deconvolution Interactive Tool (GEDIT) produces the lowest error. The best tool for other datasets is varied, but CIBERSORT and GEDIT most consistently produce accurate results. We find that optimal reference depends on the tool used, and report suggested references to be used with each tool. Most tools return results within minutes, but on large datasets runtimes for CIBERSORT can exceed hours or even days. We conclude that deconvolution methods are capable of returning high-quality results, but that proper reference selection is critical.

Examining Sex-Differentiated Genetic Effects Across Neuropsychiatric and Behavioral Traits.

BACKGROUND: The origin of sex differences in prevalence and presentation of neuropsychiatric and behavioral traits is largely unknown. Given established genetic contributions and correlations, we tested for a sex-differentiated genetic architecture within and between traits. METHODS: Using European ancestry genome-wide association summary statistics for 20 neuropsychiatric and behavioral traits, we tested for sex differences in single nucleotide polymorphism (SNP)-based heritability and genetic correlation (rg < 1). For each trait, we computed per-SNP z scores from sex-stratified regression coefficients and identified genes with sex-differentiated effects using a gene-based approach. We calculated correlation coefficients between z scores to test for shared sex-differentiated effects. Finally, we tested for sex differences in across-trait genetic correlations. RESULTS: We observed no consistent sex differences in SNP-based heritability. Between-sex, within-trait genetic correlations were high, although <1 for educational attainment and risk-taking behavior. We identified 4 genes with significant sex-differentiated effects across 3 traits. Several trait pairs shared sex-differentiated effects. The top genes with sex-differentiated effects were enriched for multiple gene sets, including neuron- and synapse-related sets. Most between-trait genetic correlation estimates were not significantly different between sexes, with exceptions (educational attainment and risk-taking behavior). CONCLUSIONS: Sex differences in the common autosomal genetic architecture of neuropsychiatric and behavioral phenotypes are small and polygenic and unlikely to fully account for observed sex-differentiated attributes. Larger sample sizes are needed to identify sex-differentiated effects for most traits. For well-powered studies, we identified genes with sex-differentiated effects that were enriched for neuron-related and other biological functions. This work motivates further investigation of genetic and environmental influences on sex differences.

The Evolutionary History of Common Genetic Variants Influencing Human Cortical Surface Area.

Structural brain changes along the lineage leading to modern Homo sapiens contributed to our distinctive cognitive and social abilities. However, the evolutionarily relevant molecular variants impacting key aspects of neuroanatomy are largely unknown. Here, we integrate evolutionary annotations of the genome at diverse timescales with common variant associations from large-scale neuroimaging genetic screens. We find that alleles with evidence of recent positive polygenic selection over the past 2000-3000 years are associated with increased surface area (SA) of the entire cortex, as well as specific regions, including those involved in spoken language and visual processing. Therefore, polygenic selective pressures impact the structure of specific cortical areas even over relatively recent timescales. Moreover, common sequence variation within human gained enhancers active in the prenatal cortex is associated with postnatal global SA. We show that such variation modulates the function of a regulatory element of the developmentally relevant transcription factor HEY2 in human neural progenitor cells and is associated with structural changes in the inferior frontal cortex. These results indicate that non-coding genomic regions active during prenatal cortical development are involved in the evolution of human brain structure and identify novel regulatory elements and genes impacting modern human brain structure.

Integration of genomics and transcriptomics predicts diabetic retinopathy susceptibility genes.

We determined differential gene expression in response to high glucose in lymphoblastoid cell lines derived from matched individuals with type 1 diabetes with and without retinopathy. Those genes exhibiting the largest difference in glucose response were assessed for association with diabetic retinopathy in a genome-wide association study meta-analysis. Expression quantitative trait loci (eQTLs) of the glucose response genes were tested for association with diabetic retinopathy. We detected an enrichment of the eQTLs from the glucose response genes among small association p-values and identified folliculin (FLCN) as a susceptibility gene for diabetic retinopathy. Expression of FLCN in response to glucose was greater in individuals with diabetic retinopathy. Independent cohorts of individuals with diabetes revealed an association of FLCN eQTLs with diabetic retinopathy. Mendelian randomization confirmed a direct positive effect of increased FLCN expression on retinopathy. Integrating genetic association with gene expression implicated FLCN as a disease gene for diabetic retinopathy.

One of the side effects of diabetes is loss of vision from diabetic retinopathy, which is caused by injury to the light sensing tissue in the eye, the retina. Almost all individuals with diabetes develop diabetic retinopathy to some extent, and it is the leading cause of irreversible vision loss in working-age adults in the United States. How long a person has been living with diabetes, the extent of increased blood sugars and genetics all contribute to the risk and severity of diabetic retinopathy. Unfortunately, virtually no genes associated with diabetic retinopathy have yet been identified. When a gene is activated, it produces messenger molecules known as mRNA that are used by cells as instructions to produce proteins. The analysis of mRNA molecules, as well as genes themselves, can reveal the role of certain genes in disease. The studies of all genes and their associated mRNAs are respectively called genomics and transcriptomics. Genomics reveals what genes are present, while transcriptomics shows how active genes are in different cells. Skol et al. developed methods to study genomics and transcriptomics together to help discover genes that cause diabetic retinopathy. Genes involved in how cells respond to high blood sugar were first identified using cells grown in the lab. By comparing the activity of these genes in people with and without retinopathy the study identified genes associated with an increased risk of retinopathy in diabetes. In people with retinopathy, the activity of the folliculin gene (FLCN) increased more in response to high blood sugar. This was further verified with independent groups of people and using computer models to estimate the effect of different versions of the folliculin gene. The methods used here could be applied to understand complex genetics in other diseases. The results provide new understanding of the effects of diabetes. They may also help in the development of new treatments for diabetic retinopathy, which are likely to improve on the current approach of using laser surgery or injections into the eye.

Association of CNVs with methylation variation.

Germline copy number variants (CNVs) and single-nucleotide polymorphisms (SNPs) form the basis of inter-individual genetic variation. Although the phenotypic effects of SNPs have been extensively investigated, the effects of CNVs is relatively less understood. To better characterize mechanisms by which CNVs affect cellular phenotype, we tested their association with variable CpG methylation in a genome-wide manner. Using paired CNV and methylation data from the 1000 genomes and HapMap projects, we identified genome-wide associations by methylation quantitative trait locus (mQTL) analysis. We found individual CNVs being associated with methylation of multiple CpGs and vice versa. CNV-associated methylation changes were correlated with gene expression. CNV-mQTLs were enriched for regulatory regions, transcription factor-binding sites (TFBSs), and were involved in long-range physical interactions with associated CpGs. Some CNV-mQTLs were associated with methylation of imprinted genes. Several CNV-mQTLs and/or associated genes were among those previously reported by genome-wide association studies (GWASs). We demonstrate that germline CNVs in the genome are associated with CpG methylation. Our findings suggest that structural variation together with methylation may affect cellular phenotype.

The impact of sex on gene expression across human tissues.

Many complex human phenotypes exhibit sex-differentiated characteristics. However, the molecular mechanisms underlying these differences remain largely unknown. We generated a catalog of sex differences in gene expression and in the genetic regulation of gene expression across 44 human tissue sources surveyed by the Genotype-Tissue Expression project (GTEx, v8 release). We demonstrate that sex influences gene expression levels and cellular composition of tissue samples across the human body. A total of 37% of all genes exhibit sex-biased expression in at least one tissue. We identify cis expression quantitative trait loci (eQTLs) with sex-differentiated effects and characterize their cellular origin. By integrating sex-biased eQTLs with genome-wide association study data, we identify 58 gene-trait associations that are driven by genetic regulation of gene expression in a single sex. These findings provide an extensive characterization of sex differences in the human transcriptome and its genetic regulation.

Determinants of telomere length across human tissues.

Telomere shortening is a hallmark of aging. Telomere length (TL) in blood cells has been studied extensively as a biomarker of human aging and disease; however, little is known regarding variability in TL in nonblood, disease-relevant tissue types. Here, we characterize variability in TLs from 6391 tissue samples, representing >20 tissue types and 952 individuals from the Genotype-Tissue Expression (GTEx) project. We describe differences across tissue types, positive correlation among tissue types, and associations with age and ancestry. We show that genetic variation affects TL in multiple tissue types and that TL may mediate the effect of age on gene expression. Our results provide the foundational knowledge regarding TL in healthy tissues that is needed to interpret epidemiological studies of TL and human health.

Cell type-specific genetic regulation of gene expression across human tissues.

The Genotype-Tissue Expression (GTEx) project has identified expression and splicing quantitative trait loci in cis (QTLs) for the majority of genes across a wide range of human tissues. However, the functional characterization of these QTLs has been limited by the heterogeneous cellular composition of GTEx tissue samples. We mapped interactions between computational estimates of cell type abundance and genotype to identify cell type-interaction QTLs for seven cell types and show that cell type-interaction expression QTLs (eQTLs) provide finer resolution to tissue specificity than bulk tissue cis-eQTLs. Analyses of genetic associations with 87 complex traits show a contribution from cell type-interaction QTLs and enables the discovery of hundreds of previously unidentified colocalized loci that are masked in bulk tissue.

Pharmacogenomic genotypes define genetic ancestry in patients and enable population-specific genomic implementation.

The importance of genetic ancestry characterization is increasing in genomic implementation efforts, and clinical pharmacogenomic guidelines are being published that include population-specific recommendations. Our aim was to test the ability of focused clinical pharmacogenomic SNP panels to estimate individual genetic ancestry (IGA) and implement population-specific pharmacogenomic clinical decision-support (CDS) tools. Principle components and STRUCTURE were utilized to assess differences in genetic composition and estimate IGA among 1572 individuals from 1000 Genomes, two independent cohorts of Caucasians and African Americans (AAs), plus a real-world validation population of patients undergoing pharmacogenomic genotyping. We found that clinical pharmacogenomic SNP panels accurately estimate IGA compared to genome-wide genotyping and identify AAs with ≥70 African ancestry (sensitivity >82%, specificity >80%, PPV >95%, NPV >47%). We also validated a new AA-specific warfarin dosing algorithm for patients with ≥70% African ancestry and implemented it at our institution as a novel CDS tool. Consideration of IGA to develop an institutional CDS tool was accomplished to enable population-specific pharmacogenomic guidance at the point-of-care. These capabilities were immediately applied for guidance of warfarin dosing in AAs versus Caucasians, but also provide a real-world model that can be extended to other populations and drugs as actionable genomic evidence accumulates.

5-Hydroxymethylcytosine Profiles in Circulating Cell-Free DNA Associate with Disease Burden in Children with Neuroblastoma.

PURPOSE: 5-Hydroxymethylcytosine (5-hmC) is an epigenetic marker of open chromatin and active gene expression. We profiled 5-hmC with Nano-hmC-Seal technology using 10 ng of plasma-derived cell-free DNA (cfDNA) in blood samples from patients with neuroblastoma to determine its utility as a biomarker. EXPERIMENTAL DESIGN: For the Discovery cohort, 100 5-hmC profiles were generated from 34 well children and 32 patients (27 high-risk, 2 intermediate-risk, and 3 low-risk) at various time points during the course of their disease. An independent Validation cohort encompassed 5-hmC cfDNA profiles (n = 29) generated from 21 patients (20 high-risk and 1 intermediate-risk). Metastatic burden was classified as high, moderate, low, or none per Curie metaiodobenzylguanidine scores and percentage of tumor cells in bone marrow. Genes with differential 5-hmC levels between samples according to metastatic burden were identified using DESeq2. RESULTS: Hierarchical clustering using 5-hmC levels of 347 genes identified from the Discovery cohort defined four clusters of samples that were confirmed in the Validation cohort and corresponded to high, high-moderate, moderate, and low/no metastatic burden. Samples from patients with increased metastatic burden had increased 5-hmC deposition on genes in neuronal stem cell maintenance and epigenetic regulatory pathways. Further, 5-hmC cfDNA profiles generated with 1,242 neuronal pathway genes were associated with subsequent relapse in the cluster of patients with predominantly low or no metastatic burden (sensitivity 65%, specificity 75.6%). CONCLUSIONS: cfDNA 5-hmC profiles in children with neuroblastoma correlate with metastatic burden and warrants development as a biomarker of treatment response and outcome.

ImmGen report: sexual dimorphism in the immune system transcriptome.

Sexual dimorphism in the mammalian immune system is manifested as more frequent and severe infectious diseases in males and, on the other hand, higher rates of autoimmune disease in females, yet insights underlying those differences are still lacking. Here we characterize sex differences in the immune system by RNA and ATAC sequence profiling of untreated and interferon-induced immune cell types in male and female mice. We detect very few differentially expressed genes between male and female immune cells except in macrophages from three different tissues. Accordingly, very few genomic regions display differences in accessibility between sexes. Transcriptional sexual dimorphism in macrophages is mediated by genes of innate immune pathways, and increases after interferon stimulation. Thus, the stronger immune response of females may be due to more activated innate immune pathways prior to pathogen invasion.

Author Correction: The role of sex in the genomics of human complex traits.

In Box 4 of the originally published article, the text describing the Miami plot in part c of the figure contained some minor errors. Specifically, during pre-publication revision of the article, the authors updated the illustrative Miami plot (shown in figure part c) from that of reference 80 (Randall et al. PLoS Genet. (2013)) to a more recent study in reference 82 (Winkler et al. PLoS Genet. (2015)). The box text has now been updated to reflect that change. In paragraph 2, the trait has been updated from "hip circumference adjusted for body mass index" to "waist-to-hip ratio adjusted for body mass index (under 50 years old)" and it has been clarified that female GWAS data are shown on the top half of the plot with male data at the bottom. The original two citations of reference 80 in the Box 4 text have been updated to reference 82. Finally, a typographical artefact was corrected on the Y axis of the Miami plot, whereby the labels '14' and '16' in the top half of the plot we originally both shown as '12'. None of these corrections alter the overall illustrative point that genetic architectures for traits can differ between males and females, which was a conclusion of both reference 80 and reference 82.

5-Hydroxymethylcytosine Profiles Are Prognostic of Outcome in Neuroblastoma and Reveal Transcriptional Networks That Correlate With Tumor Phenotype.

PURPOSE: Whole-genome profiles of the epigenetic modification 5-hydroxymethylcytosine (5-hmC) are robust diagnostic biomarkers in adult patients with cancer. We investigated if 5-hmC profiles would serve as novel prognostic markers in neuroblastoma, a clinically heterogeneous pediatric cancer. Because this DNA modification facilitates active gene expression, we hypothesized that 5-hmC profiles would identify transcriptomic networks driving the clinical behavior of neuroblastoma. PATIENTS AND METHODS: Nano-hmC-Seal sequencing was performed on DNA from Discovery (n = 51), Validation (n = 38), and Children's Oncology Group (n = 20) cohorts of neuroblastoma tumors. RNA was isolated from 48 tumors for RNA sequencing. Genes with differential 5-hmC or expression between clusters were identified using DESeq2. A 5-hmC model predicting outcome in high-risk patients was established using linear discriminant analysis. RESULTS: Comparison of low- versus high-risk tumors in the Discovery cohort revealed 577 genes with differential 5-hmC. Hierarchical clustering of tumors from the Discovery and Validation cohorts using these genes identified two main clusters highly associated with established prognostic markers, clinical risk group, and outcome. Genes with increased 5-hmC and expression in the favorable cluster were enriched for pathways of neuronal differentiation and KRAS activation, whereas genes involved in inflammation and the PRC2 complex were identified in the unfavorable cluster. The linear discriminant analysis model trained on high-risk Discovery cohort tumors was prognostic of outcome when applied to high-risk tumors from the Validation and Children's Oncology Group cohorts (hazard ratio, 3.8). CONCLUSION: 5-hmC profiles may be optimal DNA-based biomarkers in neuroblastoma. Analysis of transcriptional networks regulated by these epigenomic modifications may lead to a deeper understanding of drivers of neuroblastoma phenotype.

Sex differences in gene expression in response to ischemia in the human left ventricular myocardium.

Sex differences exist in the prevalence, presentation and outcomes of ischemic heart disease (IHD). Females have higher risk of heart failure post-myocardial infarction relative to males and are two to three times more likely to die after coronary artery bypass grafting surgery. We examined sex differences in human myocardial gene expression in response to ischemia. Left ventricular biopsies from 68 male/46 female patients undergoing aortic valve replacement surgery were obtained at baseline and after a median 74 min of cold cardioplegic arrest/ischemia. Transcriptomes were quantified by RNA-sequencing. Cell-type enrichment analysis was used to estimate the identity and relative proportions of different cell types in each sample. A sex-specific response to ischemia was observed for 271 genes. Notably, the expression FAM5C, PLA2G4E and CYP1A1 showed an increased expression in females compared to males due to ischemia and DIO3, MT1G and CMA1 showed a decreased expression in females compared to males due to ischemia. Functional annotation analysis revealed sex-specific modulation of the oxytocin signaling pathway and common pathway of fibrin clot formation. Expression quantitative trait locus (eQTL) analysis identified variant-by-sex interaction eQTLs, indicative of sex differences in the genotypic effects on gene expression. Cell-type enrichment analysis showed sex-bias in proportion of specific cell types. Common lymphoid progenitor cells and M2 macrophages were found to increase in female samples from pre- to post-ischemia, but no change was observed in male samples. These differences in response to myocardial ischemia provide insight into the sexual dimorphism of IHD and may aid in the development of sex-specific therapies that reduce myocardial injury.