Sex-biased genetic effects on gene regulation in humans

Segregation GWAS to linearize a non-additive locus with incomplete penetrance: an example of horn status in sheep

BACKGROUND

The objective of this study was to introduce a genome-wide association study (GWAS) in conjunction with segregation analysis on monogenic categorical traits. Genotype probabilities calculated from phenotypes, mode of inheritance and pedigree information, are expressed as the expected allele count (EAC) (range 0 to 2), and are inherited additively, by definition, unlike the original phenotypes, which are non-additive and could be of incomplete penetrance. The EAC are regressed on the single nucleotide polymorphism (SNP) genotypes, similar to an additive GWAS. In this study, horn phenotypes in Merino sheep are used to illustrate the advantages of using the segregation GWAS, a trait believed to be monogenic, affected by dominance, sex-dependent expression and likely affected by incomplete penetrance. We also used simulation to investigate whether incomplete penetrance can cause prediction errors in Merino sheep for horn status.

RESULTS

Estimated penetrance values differed between the sexes, where males showed almost complete penetrance, especially for horned and polled phenotypes, while females had low penetrance values for the horned status. This suggests that females homozygous for the 'horned allele' have a horned phenotype in only 22% of the cases while 78% will be knobbed or have scurs. The GWAS using EAC on 4001 animals and 510,174 SNP genotypes from the Illumina Ovine high-density (600k) chip gave a stronger association compared to using actual phenotypes. The correlation between the EAC and the allele count of the SNP with the highest -log10(p-value) was 0.73 in males and 0.67 in females. Simulations using penetrance values found by the segregation analyses resulted in higher correlations between the EAC and the causative mutation (0.95 for males and 0.89 for females, respectively), suggesting that the most predictive SNP is not in full LD with the causative mutation.

CONCLUSIONS

Our results show clear differences in penetrance values between males and female Merino sheep for horn status. Segregation analysis for a trait with mutually exclusive phenotypes, non-additive inheritance, and/or incomplete penetrance can lead to considerably more power in a GWAS because the linearized genotype probabilities are additive and can accommodate incomplete penetrance. This method can be extended to any monogenic controlled categorical trait of which the phenotypes are mutually exclusive.

Dynamic stress- and inflammatory-based regulation of psychiatric risk loci in human neurons

The prenatal environment can alter neurodevelopmental and clinical trajectories, markedly increasing risk for psychiatric disorders in childhood and adolescence. To understand if and how fetal exposures to stress and inflammation exacerbate manifestation of genetic risk for complex brain disorders, we report a large-scale context-dependent massively parallel reporter assay (MPRA) in human neurons designed to catalogue genotype x environment (GxE) interactions. Across 240 genome-wide association study (GWAS) loci linked to ten brain traits/disorders, the impact of hydrocortisone, interleukin 6, and interferon alpha on transcriptional activity is empirically evaluated in human induced pluripotent stem cell (hiPSC)-derived glutamatergic neurons. Of ~3,500 candidate regulatory risk elements (CREs), 11% of variants are active at baseline, whereas cue-specific CRE regulatory activity range from a high of 23% (hydrocortisone) to a low of 6% (IL-6). Cue-specific regulatory activity is driven, at least in part, by differences in transcription factor binding activity, the gene targets of which show unique enrichments for brain disorders as well as co-morbid metabolic and immune syndromes. The dynamic nature of genetic regulation informs the influence of environmental factors, reveals a mechanism underlying pleiotropy and variable penetrance, and identifies specific risk variants that confer greater disorder susceptibility after exposure to stress or inflammation. Understanding neurodevelopmental GxE interactions will inform mental health trajectories and uncover novel targets for therapeutic intervention.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Extensive Alternative Splicing Patterns in Systemic Lupus Erythematosus Highlight Sexual Differences

Substantial evidence highlights divergences in immune responses between men and women. Women are more susceptible to autoimmunity, whereas men suffer from the more severe presentation of autoimmune disorders. The molecular mechanism of this sexual dimorphism remains elusive. Herein, we conducted a comprehensive analysis of sex differences in whole-blood gene expression focusing on alternative splicing (AS) events in systemic lupus erythematosus (SLE), which is a prototype sex-biased disease. This study included 79 SLE patients with active disease and 58 matched healthy controls who underwent whole-blood RNA sequencing. Sex differences in splicing events were widespread, existent in both SLE and a healthy state. However, we observed distinct gene sets and molecular pathways targeted by sex-dependent AS in SLE patients as compared to healthy subjects, as well as a notable sex dissimilarity in intron retention events. Sexually differential spliced genes specific to SLE patients were enriched for dynamic cellular processes including chromatin remodeling, stress and inflammatory responses. Remarkably, the extent of sexual differences in AS in the SLE patients and healthy individuals exceeded those in gene expression. Overall, this study reveals an unprecedent variation in sex-dependent splicing events in SLE and the healthy state, with potential implications for understanding the molecular basis of sexual dimorphism in autoimmunity.

Genetics and epigenetics of primary Sjögren syndrome: implications for future therapies

In primary Sjögren syndrome (pSS), chronic inflammation of exocrine glands results in tissue destruction and sicca symptoms, primarily of the mouth and eyes. Fatigue, arthralgia and myalgia are also common symptoms, whereas extraglandular manifestations that involve the respiratory, nervous and vascular systems occur in a subset of patients. The disease predominantly affects women, with an estimated female to male ratio of 14 to 1. The aetiology of pSS, however, remains incompletely understood, and effective treatment is lacking. Large-scale genetic and epigenetic investigations have revealed associations between pSS and genes in both innate and adaptive immune pathways. The genetic variants mediate context-dependent effects, and both sex and environmental factors can influence the outcome. As such, genetic and epigenetic studies can provide insight into the dysregulated molecular mechanisms, which in turn might reveal new therapeutic possibilities. This Review discusses the genetic and epigenetic features that have been robustly connected with pSS, putting them into the context of cellular function, carrier sex and environmental challenges. In all, the observations point to several novel opportunities for early detection, treatment development and the pathway towards personalized medicine.

Sex chromosome mechanisms in cardiac development and disease

Many human diseases, including cardiovascular disease, show differences between men and women in pathology and treatment outcomes. In the case of cardiac disease, sex differences are exemplified by differences in the frequency of specific types of congenital and adult-onset heart disease. Clinical studies have suggested that gonadal hormones are a factor in sex bias. However, recent research has shown that gene and protein networks under non-hormonal control also account for cardiac sex differences. In this review, we describe the sex chromosome pathways that lead to sex differences in the development and function of the heart and highlight how these findings affect future care and treatment of cardiac disease.

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.

A systematic review of analytical methods used in genetic association analysis of the X-chromosome

Genetic association studies have been very successful at elucidating the genetic background of many complex diseases/traits. However, the X-chromosome is often neglected in these studies because of technical difficulties and the fact that most tools only utilize genetic data from autosomes. In this review, we aim to provide an overview of different practical approaches that are followed to incorporate the X-chromosome in association analysis, such as Genome-Wide Association Studies and Expression Quantitative Trait Loci Analysis. In general, the choice of which test statistics is most appropriate will depend on three main criteria: (1) the underlying X-inactivation model, (2) if Hardy-Weinberg equilibrium holds and sex-specific allele frequencies are expected and (3) whether adjustment for confounding variables is required. All in all, it is recommended that a combination of different association tests should be used for the analysis of X-chromosome.

Characterising sex differences of autosomal DNA methylation in whole blood using the Illumina EPIC array

BACKGROUND

Sex differences are known to play a role in disease aetiology, progression and outcome. Previous studies have revealed autosomal epigenetic differences between males and females in some tissues, including differences in DNA methylation patterns. Here, we report for the first time an analysis of autosomal sex differences in DNAme using the Illumina EPIC array in human whole blood by performing a discovery (n = 1171) and validation (n = 2471) analysis.

RESULTS

We identified and validated 396 sex-associated differentially methylated CpG sites (saDMPs) with the majority found to be female-biased CpGs (74%). These saDMP's are enriched in CpG islands and CpG shores and located preferentially at 5'UTRs, 3'UTRs and enhancers. Additionally, we identified 266 significant sex-associated differentially methylated regions overlapping genes, which have previously been shown to exhibit epigenetic sex differences, and novel genes. Transcription factor binding site enrichment revealed enrichment of transcription factors related to critical developmental processes and sex determination such as SRY and ESR1.

CONCLUSION

Our study reports a reliable catalogue of sex-associated CpG sites and elucidates several characteristics of these sites using large-scale discovery and validation data sets. This resource will benefit future studies aiming to investigate sex specific epigenetic signatures and further our understanding of the role of DNA methylation in sex differences in human whole blood.

Sex and gender in interstitial lung diseases

Sex and gender differences influence key domains of research, lung health, healthcare access and healthcare delivery. In interstitial lung diseases (ILDs), mouse models of pulmonary fibrosis are clearly influenced by sex hormones. Additionally, short telomeres, a biomarker of telomere regulation gene mutations, are impacted by sex, while heritability unexplained by genetic variation may be attributable to gendered environmental factors that drive epigenetic control. Diseases like idiopathic pulmonary fibrosis, hypersensitivity pneumonitis, occupational ILDs, connective tissue-associated ILDs and lymphangioleiomyomatosis have different prevalence and prognosis between men and women. These differences arise from a complex interplay between biological sex and sociocultural gender influencing genetics, epigenomic modifiers, hormones, immune function, response to treatment and interaction with healthcare systems. Much work remains to be done to systematically integrate sex and gender analysis into relevant domains of science and clinical care in ILD, from strategic considerations for establishing research priorities to guidelines for establishing best clinical practices. Accounting for sex and gender in ILD is essential to the practice of individualised, patient-centred medicine.

Sex differences in the genetic regulation of the blood transcriptome response to glucocorticoid receptor activation

Substantial sex differences have been reported in the physiological response to stress at multiple levels, including the release of the stress hormone, cortisol. Here, we explore the genomic variants in 93 females and 196 males regulating the initial transcriptional response to cortisol via glucocorticoid receptor (GR) activation. Gene expression levels in peripheral blood were obtained before and after GR-stimulation with the selective GR agonist dexamethasone to identify differential expression following GR-activation. Sex stratified analyses revealed that while the transcripts responsive to GR-stimulation were mostly overlapping between males and females, the quantitative trait loci (eQTLs) regulation differential transcription to GR-stimulation was distinct. Sex-stratified eQTL SNPs (eSNPs) were located in different functional genomic elements and sex-stratified transcripts were enriched within postmortem brain transcriptional profiles associated with Major Depressive Disorder (MDD) specifically in males and females in the cingulate cortex. Female eSNPs were enriched among SNPs linked to MDD in genome-wide association studies. Finally, transcriptional sensitive genetic profile scores derived from sex-stratified eSNPS regulating differential transcription to GR-stimulation were predictive of depression status and depressive symptoms in a sex-concordant manner in a child and adolescent cohort (n = 584). These results suggest the potential of eQTLs regulating differential transcription to GR-stimulation as biomarkers of sex-specific biological risk for stress-related psychiatric disorders.

Role of Sex on the Genetic Susceptibility to Childhood Asthma in Latinos and African Americans

Asthma is a respiratory disease whose prevalence changes throughout the lifespan and differs by sex, being more prevalent in males during childhood and in females after puberty. In this study, we assessed the influence of sex on the genetic susceptibility to childhood asthma in admixed populations. Sex-interaction and sex-stratified genome-wide association studies (GWAS) were performed in 4291 Latinos and 1730 African Americans separately, and results were meta-analyzed. Genome-wide (p ≤ 9.35 × 10-8) and suggestive (p ≤ 1.87 × 10-6) population-specific significance thresholds were calculated based on 1000 Genomes Project data. Additionally, protein quantitative trait locus (pQTL) information was gathered for the suggestively associated variants, and enrichment analyses of the proteins identified were carried out. Four independent loci showed interaction with sex at a suggestive level. The stratified GWAS highlighted the 17q12-21 asthma locus as a contributor to asthma susceptibility in both sexes but reached genome-wide significance only in females (p-females < 9.2 × 10-8; p-males < 1.25 × 10-2). Conversely, genetic variants upstream of ligand-dependent nuclear receptor corepressor-like gene (LCORL), previously involved in height determination and spermatogenesis, were associated with asthma only in males (minimum p = 5.31 × 10-8 for rs4593128). Enrichment analyses revealed an overrepresentation of processes related to the immune system and highlighted differences between sexes. In conclusion, we identified sex-specific polymorphisms that could contribute to the differences in the prevalence of childhood asthma between males and females.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians

Background

Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors).

Methods

In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors.

Results

The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors.

Conclusions

Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01900-y.

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.

Sexual dimorphism in cancer: insights from transcriptional signatures in kidney tissue and renal cell carcinoma

Sexual dimorphism in cancer incidence and outcome is widespread. Understanding the underlying mechanisms is fundamental to improve cancer prevention and clinical management. Sex disparities are particularly striking in kidney cancer: across diverse populations, men consistently show unexplained 2-fold increased incidence and worse prognosis. We have characterized genome-wide expression and regulatory networks of 609 renal tumors and 256 non-tumor renal tissues. Normal kidney displayed sex-specific transcriptional signatures, including higher expression of X-linked tumor suppressor genes in women. Sex-dependent genotype-phenotype associations unraveled women-specific immune regulation. Sex differences were markedly expanded in tumors, with male-biased expression of key genes implicated in metabolism, non-malignant diseases with male predominance and carcinogenesis, including markers of tumor infiltrating leukocytes. Analysis of sex-dependent RCC progression and survival uncovered prognostic markers involved in immune response and oxygen homeostasis. In summary, human kidney tissues display remarkable sexual dimorphism at the molecular level. Sex-specific transcriptional signatures further shape renal cancer, with relevance for clinical management.

Sex and Gender Omic Biomarkers in Men and Women With COPD: Considerations for Precision Medicine

Sex and gender differences in lung health and disease are imperative to consider and study if precision pulmonary medicine is to be achieved. The development of reliable COPD biomarkers has been elusive, and the translation of biomarkers to clinical care has been limited. Useful and effective biomarkers must be developed with attention to clinical heterogeneity of COPD; inherent heterogeneity exists related to grouping women and men together in the studies of COPD. Considering sex and gender differences and influences related to -omics may represent progress in susceptibility, diagnostic, prognostic, and therapeutic biomarker development and clinical innovation to improve the lung health of men and women.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes

BACKGROUND

HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk.

METHODS

One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed.

RESULTS

A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1.

CONCLUSIONS

The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

Sex-specific differences in peripheral blood leukocyte transcriptional response to LPS are enriched for HLA region and X chromosome genes

Sex-specific differences in prevalence are well documented for many common, complex diseases, especially for immune-mediated diseases, yet the precise mechanisms through which factors associated with biological sex exert their effects throughout life are not well understood. We interrogated sex-specific transcriptional responses of peripheral blood leukocytes (PBLs) to innate immune stimulation by lipopolysaccharide (LPS) in 46 male and 66 female members of the Hutterite community, who practice a communal lifestyle. We identified 1217 autosomal and 54 X-linked genes with sex-specific responses to LPS, as well as 71 autosomal and one X-linked sex-specific expression quantitative trait loci (eQTLs). Despite a similar proportion of the 15 HLA genes responding to LPS compared to all expressed autosomal genes, there was a significant over-representation of genes with sex by treatment interactions among HLA genes. We also observed an enrichment of sex-specific differentially expressed genes in response to LPS for X-linked genes compared to autosomal genes, suggesting that HLA and X-linked genes may disproportionately contribute to sex disparities in risk for immune-mediated diseases.

Sex-dimorphic genetic effects and novel loci for fasting glucose and insulin variability

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

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.

Understanding immune variation for improved translational medicine

The goal of translational medicine is to use an improved understanding of human biology to develop new clinical approaches. Immune responses are highly variable from one person to another, with this variability strongly impacting clinical outcome. Variable immunity can determine differential risks for infection, for development of autoimmunity, and for response to therapeutic interventions. Therefore, a better understanding of the causes of such differences has huge potential to improve patient management through precision medicine strategies. Variability in immunity is determined by intrinsic (e.g. age, sex), extrinsic (e.g. environment, diet), and genetic factors. There is a growing consensus that genetics factors account for 20-40% of immune variability between individuals. The remaining unexplained variability is likely due to direct environmental influences, as well as specific gene-environmental interactions, which are more challenging to quantify and study. However, population based cohort studies with systems immunology approaches are now providing new understanding into these associations.

Genetics and Epigenetics of Sex Bias: Insights from Human Cancer and Autoimmunity

High-throughput sequencing and genome-wide association studies have revealed a sex bias in human diseases. The underlying molecular mechanisms remain, however, unknown. Here, we cover recent advances in cancer and autoimmunity focusing on intrinsic genetic and epigenetic differences underlying sex biases in human disease. These studies reveal a central role of genome regulatory mechanisms including genome repair, chromosome folding, and epigenetic regulation in dictating the sex bias. These highlight the importance of considering sex as a variable in both basic science and clinical investigations. Understanding the molecular mechanisms underlying sex bias in human diseases will be instrumental in making a first step forwards into the era of personalized medicine.

Clinical Advances in Sex- and Gender-Informed Medicine to Improve the Health of All: A Review

IMPORTANCE

Biological sex and sociocultural gender represent major sources of diversity among patients, and recent research has shown the association of sex and gender with health. A growing body of literature describes widespread associations of sex and gender with cells, organs, and the manner in which individual patients interact with health care systems. Sex- and gender-informed medicine is a young paradigm of clinical practice and medical research founded on this literature that considers the association of sex and gender with each element of the disease process from risk, to presentation, to response to therapy.

OBSERVATIONS

Characteristics that underlie sex and gender involve both endogenous and exogenous factors that change throughout the life course. This review details clinical examples with broad applicability that highlight sex and gender differences in the key domains of genetics, epigenomic modifiers, hormonal milieu, immune function, neurocognitive aging process, vascular health, response to therapeutics, and interaction with health care systems. These domains interact with one another in multidimensional associations, contributing to the diversity of the sex and gender spectra. Novel research has identified differences of clinical relevance with the potential to improve care for all patients.

CONCLUSIONS AND RELEVANCE

Clinicians should consider incorporating sex and gender in their decision-making to practice precision medicine that integrates fundamental components of patient individuality. Recognizing the biological and environmental factors that affect the disease course is imperative to optimizing care for each patient. Research highlights the myriad ways sex and gender play a role in health and disease. However, these clinically relevant insights have yet to be systematically incorporated into care. The framework described in this review serves as a guide to help clinicians consider sex and gender as they practice precision medicine.

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.

A New Liver Expression Quantitative Trait Locus Map From 1,183 Individuals Provides Evidence for Novel Expression Quantitative Trait Loci of Drug Response, Metabolic, and Sex-Biased Phenotypes

Expression quantitative trait locus (eQTL) studies in human liver are crucial for elucidating how genetic variation influences variability in disease risk and therapeutic outcomes and may help guide strategies to obtain maximal efficacy and safety of clinical interventions. Associations between expression microarray and genome-wide genotype data from four human liver eQTL studies (n = 1,183) were analyzed. More than 2.3 million cis-eQTLs for 15,668 genes were identified. When eQTLs were filtered against a list of 1,496 drug response genes, 187,829 cis-eQTLs for 1,191 genes were identified. Additionally, 1,683 sex-biased cis-eQTLs were identified, as well as 49 and 73 cis-eQTLs that colocalized with genome-wide association study signals for blood metabolite or lipid levels, respectively. Translational relevance of these results is evidenced by linking DPYD eQTLs to differences in safety of chemotherapy, linking the sex-biased regulation of PCSK9 expression to anti-lipid therapy, and identifying the G-protein coupled receptor GPR180 as a novel drug target for hypertriglyceridemia.

Distinct molecular etiologies of male and female hepatocellular carcinoma

BACKGROUND

Sex-differences in cancer occurrence and mortality are evident across tumor types; men exhibit higher rates of incidence and often poorer responses to treatment. Targeted approaches to the treatment of tumors that account for these sex-differences require the characterization and understanding of the fundamental biological mechanisms that differentiate them. Hepatocellular Carcinoma (HCC) is the second leading cause of cancer death worldwide, with the incidence rapidly rising. HCC exhibits a male-bias in occurrence and mortality, but previous studies have failed to explore the sex-specific dysregulation of gene expression in HCC.

METHODS

Here, we characterize the sex-shared and sex-specific regulatory changes in HCC tumors in the TCGA LIHC cohort using combined and sex-stratified differential expression and eQTL analyses.

RESULTS

By using a sex-specific differential expression analysis of tumor and tumor-adjacent samples, we uncovered etiologically relevant genes and pathways differentiating male and female HCC. While both sexes exhibited activation of pathways related to apoptosis and cell cycle, males and females differed in the activation of several signaling pathways, with females showing PPAR pathway enrichment while males showed PI3K, PI3K/AKT, FGFR, EGFR, NGF, GF1R, Rap1, DAP12, and IL-2 signaling pathway enrichment. Using eQTL analyses, we discovered germline variants with differential effects on tumor gene expression between the sexes. 24.3% of the discovered eQTLs exhibit differential effects between the sexes, illustrating the substantial role of sex in modifying the effects of eQTLs in HCC. The genes that showed sex-specific dysregulation in tumors and those that harbored a sex-specific eQTL converge in clinically relevant pathways, suggesting that the molecular etiologies of male and female HCC are partially driven by differential genetic effects on gene expression.

CONCLUSIONS

Sex-stratified analyses detect sex-specific molecular etiologies of HCC. Overall, our results provide new insight into the role of inherited genetic regulation of transcription in modulating sex-differences in HCC etiology and provide a framework for future studies on sex-biased cancers.

Putative contributions of the sex chromosome proteins SOX3 and SRY to neurodevelopmental disorders

The male-biased prevalence of certain neurodevelopmental disorders and the sex-biased outcomes associated with stress exposure during gestation have been previously described. Here, we hypothesized that genes distinctively targeted by only one or both homologous proteins highly conserved across therian mammals, SOX3 and SRY, could induce sexual adaptive changes that result in a differential risk for neurodevelopmental disorders. ChIP-seq/chip data showed that SOX3/SRY gene targets were expressed in different brain cell types in mice. We used orthologous human genes in rodent genomes to extend the number of SOX3/SRY set (1,721). These genes were later found to be enriched in five modules of coexpressed genes during the early and mid-gestation periods (FDR < 0.05), independent of sexual hormones. Genes with differential expression (24, p < 0.0001) and methylation (40, p < 0.047) between sexes were overrepresented in this set. Exclusive SOX3 or SRY target genes were more associated with the late gestational and postnatal periods. Using autism as a model sex-biased disorder, the SOX3/SRY set was enriched in autism gene databases (FDR ≤ 0.05), and there were more de novo variations from the male autism spectrum disorder (ASD) samples under the SRY peaks compared to the random peaks (p < 0.024). The comparison of coexpressed networks of SOX3/SRY target genes between male autism and control samples revealed low preservation in gene modules related to stress response (99 genes) and neurogenesis (78 genes). This study provides evidence that while SOX3 is a regulatory mechanism for both sexes, the male-exclusive SRY also plays a role in gene regulation, suggesting a potential mechanism for sex bias in ASD.

The role of sex in the genomics of human complex traits

Nearly all human complex traits and disease phenotypes exhibit some degree of sex differences, including differences in prevalence, age of onset, severity or disease progression. Until recently, the underlying genetic mechanisms of such sex differences have been largely unexplored. Advances in genomic technologies and analytical approaches are now enabling a deeper investigation into the effect of sex on human health traits. In this Review, we discuss recent insights into the genetic models and mechanisms that lead to sex differences in complex traits. This knowledge is critical for developing deeper insight into the fundamental biology of sex differences and disease processes, thus facilitating precision medicine.

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.

Sex-Interacting mRNA- and miRNA-eQTLs and Their Implications in Gene Expression Regulation and Disease

Despite sex being an important epidemiological and physiological factor, not much is known about how sex works to interact with genotypes to result in different phenotypes. Both messenger RNA (mRNA) and microRNA (miRNA) may be differentially expressed between the sexes in different physiological conditions, and both may be differentially regulated between males and females. Using whole transcriptome data on lymphoblastoid cell lines from 338 samples of European origin, we tried to uncover genes differentially expressed between the two sexes and sex-interacting expression quantitative trait loci (ss-eQTLs). Two miRNAs were found to be differentially expressed between the two sexes, both of which were found to be functionally implicated in breast cancer. Using two stage linear regression analysis, 21 mRNA ss-eQTL and 3 miRNA ss-eQTLs were discovered. We replicated two of the mRNA ss-eQTLs (p < 0.1) using a separate dataset of gene expression data derived from monocytes. Three mRNA ss-eQTLs are in high linkage disequilibrium with variants also found to be associated with sexually dimorphic traits. Taken together, we believe the ss-eQTLs presented will assist researchers in uncovering the basis of sex-biased gene expression regulation, and ultimately help us understand the genetic basis of differences in phenotypes between sexes.

Gene-by-Sex Interactions in Mitochondrial Functions and Cardio-Metabolic Traits

We studied sex differences in over 50 cardio-metabolic traits in a panel of 100 diverse inbred strains of mice. The results clearly showed that the effects of sex on both clinical phenotypes and gene expression depend on the genetic background. In support of this, genetic loci associated with the traits frequently showed sex specificity. For example, Lyplal1, a gene implicated in human obesity, was shown to underlie a sex-specific locus for diet-induced obesity. Global gene expression analyses of tissues across the panel implicated adipose tissue "beiging" and mitochondrial functions in the sex differences. Isolated mitochondria showed gene-by-sex interactions in oxidative functions, such that some strains (C57BL/6J) showed similar function between sexes, whereas others (DBA/2J and A/J) showed increased function in females. Reduced adipose mitochondrial function in males as compared to females was associated with increased susceptibility to obesity and insulin resistance. Gonadectomy studies indicated that gonadal hormones acting in a tissue-specific manner were responsible in part for the sex differences.

Identification and analysis of the human sex-biased genes

Tremendous differences between human sexes are universally observed. Therefore, identifying and analyzing the sex-biased genes are becoming basically important for uncovering the mystery of sex differences and personalized medicine. Here, we presented a computational method to identify sex-biased genes from public gene expression databases. We obtained 1407 female-biased genes (FGs) and 1096 male-biased genes (MGs) across 14 different tissues. Bioinformatics analysis revealed that compared with MGs, FGs have higher evolutionary rate, higher single-nucleotide polymorphism density, less homologous gene numbers and smaller phyletic age. FGs have lower expression level, higher tissue specificity and later expressed stage in body development. Moreover, FGs are highly involved in immune-related functions, whereas MGs are more enriched in metabolic process. In addition, cellular network analysis revealed that MGs have higher degree, more cellular activating signaling and tend to be located in cellular inner space, whereas FGs have lower degree, more cellular repressing signaling and tend to be located in cellular outer space. Finally, the identified sex-biased genes and the discovered biological insights together can be a valuable resource helpful for investigating sex-biased physiology and medicine, for example sex-biased disease diagnosis and therapy, which represents one important aspect of personalized and precision medicine.

Estrogen-Dependent Upregulation of Adcyap1r1 Expression in Nucleus Accumbens Is Associated With Genetic Predisposition of Sex-Specific QTL for Alcohol Consumption on Rat Chromosome 4

Humans show sex differences related to alcohol use disorders (AUD). Animal model research has the potential to provide important insight into how sex differences affect alcohol consumption, particularly because female animals frequently drink more than males. In previous work, inbred strains of the selectively bred alcohol-preferring (P) and non-preferring (NP) rat lines revealed a highly significant quantitative trait locus (QTL) on rat chromosome 4, with a logarithm of the odds score of 9.2 for alcohol consumption. Recently, interval-specific congenic strains (ISCS) were developed by backcrossing the congenic P.NP line to inbred P (iP) rats to further refine the chromosome 4 QTL region. Two ISCS sub-strains, ISCS-A and ISCS-B, were obtained with a narrowed QTL, where the smallest region of overlap consisted of 8.9 Mb in ISCS-B. Interestingly, we found that females from both ISCS lines consumed significantly less alcohol than female iP controls (p < 0.05), while no differences in alcohol consumption were observed between male ISCS and iP controls. RNA-sequencing was performed on the nucleus accumbens of alcohol-naïve female ISCS-B and iP rats, which revealed differentially expressed genes (DEG) with greater than 2-fold change and that were functionally relevant to behavior. These DEGs included down-regulation of Oxt, Asb4, Gabre, Gabrq, Chat, Slc5a7, Slc18a8, Slc10a4, and Ngfr, and up-regulation of Ttr, Msln, Mpzl2, Wnt6, Slc17a7, Aldh1a2, and Gstm2. Pathway analysis identified significant alterations in gene networks controlling nervous system development and function, as well as cell signaling, GABA and serotonin receptor signaling and G-protein coupled receptor signaling. In addition, β-estradiol was identified as the most significant upstream regulator. The expression levels of estrogen-responsive genes that mapped to the QTL interval and have been previously associated with alcohol consumption were measured using RT-qPCR. We found that expression of the Adcyap1r1 gene, encoding the pituitary adenylate cyclase-activating polypeptide type 1 (PAC₁) receptor, was upregulated in female ISCS-B compared to female iP controls, while no differences were exhibited in males. In addition, sequence variants in the Adcyap1r1 promoter region showed a differential response to estrogen stimulation in vitro. These findings demonstrate that rat chromosome 4 QTL contains genetic variants that respond to estrogen and are associated with female alcohol consumption.

Thrombin-Activatable Fibrinolysis Inhibitor Polymorphisms and Cerebral Venous Thrombosis in Mexican Mestizo Patients

Thrombin-activatable fibrinolysis inhibitor (TAFI) gene polymorphisms have been proposed as a predisposing factor for cerebral venous thrombosis (CVT). We analyzed the association between CVT and TAFI single-nucleotide polymorphisms (rs3742264, rs2146881, and rs1926447) compared to healthy controls. Mexico Mestizo confirmed cases with CVT and age- and sex-matched controls with no history of venous thrombotic events were recruited from July 2006 to July 2015. Demographic, clinical, and imaging information was included in the analysis. Genotyping single-nucleotide polymorphisms were performed by allele-specific polymerase chain reaction. Allelic univariate analysis, haplotype association, and Hardy-Weinberg equilibrium were assessed. A total of 113 CVT cases (94 females [83.2%]; median age 35 years [interquartile range 27-43 years]) and 134 age- and sex-matched controls were included. The main risk factors for CVT were pregnancy/puerperium (30.9%), oral contraceptive use (19.5%), and hereditary thrombophilia (7.1%). We found no significant association for heterozygous and homozygous models for rs3742264 ( P = .30 and P = .69, respectively), rs2146881 ( P = .90 and P = .17, respectively), or rs1926447 ( P = .40 and P = .52, respectively) compared to controls; these findings were consistent in subgroup and haplotype analyses. In conclusion, TAFI rs3742264, rs2146881, and rs1926447 polymorphisms do not increase the risk of CVT in comparison to healthy controls.

Familial Analysis of Epistatic and Sex-Dependent Association of Genes of the Renin-Angiotensin-Aldosterone System and Blood Pressure

BACKGROUND

Renin-angiotensin-aldosterone system genes have been inconsistently associated with blood pressure, possibly because of unrecognized influences of sex-dependent genetic effects or gene-gene interactions (epistasis).

METHODS AND RESULTS

We tested association of systolic blood pressure with single-nucleotide polymorphisms (SNPs) at renin (REN), angiotensinogen (AGT), angiotensin-converting enzyme (ACE), angiotensin II type 1 receptor (AGTR1), and aldosterone synthase (CYP11B2), including sex-SNP or SNP-SNP interactions. Eighty-eight tagSNPs were tested in 2872 white individuals in 809 pedigrees from the Victorian Family Heart Study using variance components models. Three SNPs (rs8075924 and rs4277404 at ACE and rs12721297 at AGTR1) were individually associated with lower systolic blood pressure with significant (P<0.00076) effect sizes ≈1.7 to 2.5 mm Hg. Sex-specific associations were seen for 3 SNPs in men (rs2468523 and rs2478544 at AGT and rs11658531 at ACE) and 1 SNP in women (rs12451328 at ACE). SNP-SNP interaction was suggested (P<0.005) for 14 SNP pairs, none of which had shown individual association with systolic blood pressure. Four SNP pairs were at the same gene (2 for REN, 1 for AGT, and 1 for AGTR1). The SNP rs3097 at CYP11B2 was represented in 5 separate pairs.

CONCLUSIONS

SNPs at key renin-angiotensin-aldosterone system genes associate with systolic blood pressure individually in both sexes, individually in one sex only and only when combined with another SNP. Analyses that incorporate sex-dependent and epistatic effects could reconcile past inconsistencies and account for some of the missing heritability of blood pressure and are generally relevant to SNP association studies for any phenotype.

Distinctive roles of age, sex, and genetics in shaping transcriptional variation of human immune responses to microbial challenges

Identifying the drivers of the interindividual diversity of the human immune system is crucial to understand their consequences on immune-mediated diseases. By examining the transcriptional responses of 1,000 individuals to various microbial challenges, we show that age and sex influence the expression of many immune-related genes, but their effects are overall moderate, whereas genetic factors affect a smaller gene set but with a stronger effect. We identify numerous genetic variants that affect transcriptional variation on infection, many of which are associated with autoimmune or inflammatory disorders. These results enable additional exploration of the role of regulatory variants in the pathogenesis of immune-related diseases and improve our understanding of the respective effects of age, sex, and genetics on immune response variation.

The contribution of host genetic and nongenetic factors to immunological differences in humans remains largely undefined. Here, we generated bacterial-, fungal-, and viral-induced immune transcriptional profiles in an age- and sex-balanced cohort of 1,000 healthy individuals and searched for the determinants of immune response variation. We found that age and sex affected the transcriptional response of most immune-related genes, with age effects being more stimulus-specific relative to sex effects, which were largely shared across conditions. Although specific cell populations mediated the effects of age and sex on gene expression, including CD8⁺ T cells for age and CD4⁺ T cells and monocytes for sex, we detected a direct effect of these intrinsic factors for the majority of immune genes. The mapping of expression quantitative trait loci (eQTLs) revealed that genetic factors had a stronger effect on immune gene regulation than age and sex, yet they affected a smaller number of genes. Importantly, we identified numerous genetic variants that manifested their regulatory effects exclusively on immune stimulation, including a Candida albicans-specific master regulator at the CR1 locus. These response eQTLs were enriched in disease-associated variants, particularly for autoimmune and inflammatory disorders, indicating that differences in disease risk may result from regulatory variants exerting their effects only in the presence of immune stress. Together, this study quantifies the respective effects of age, sex, genetics, and cellular heterogeneity on the interindividual variability of immune responses and constitutes a valuable resource for further exploration in the context of different infection risks or disease outcomes.

Regulatory and evolutionary signatures of sex-biased genes on both the X chromosome and the autosomes

Background

Sex is an important but understudied factor in the genetics of human diseases. Analyses using a combination of gene expression data, ENCODE data, and evolutionary data of sex-biased gene expression in human tissues can give insight into the regulatory and evolutionary forces acting on sex-biased genes.

Methods

In this study, we analyzed the differentially expressed genes between males and females. On the X chromosome, we used a novel method and investigated the status of genes that escape X-chromosome inactivation (escape genes), taking into account the clonality of lymphoblastoid cell lines (LCLs). To investigate the regulation of sex-biased differentially expressed genes (sDEG), we conducted pathway and transcription factor enrichment analyses on the sDEGs, as well as analyses on the genomic distribution of sDEGs. Evolutionary analyses were also conducted on both sDEGs and escape genes.

Results

Genome-wide, we characterized differential gene expression between sexes in 462 RNA-seq samples and identified 587 sex-biased genes, or 3.2% of the genes surveyed. On the X chromosome, sDEGs were distributed in evolutionary strata in a similar pattern as escape genes. We found a trend of negative correlation between the gene expression breadth and nonsynonymous over synonymous mutation (dN/dS) ratios, showing a possible pleiotropic constraint on evolution of genes. Genome-wide, nine transcription factors were found enriched in binding to the regions surrounding the transcription start sites of female-biased genes. Many pathways and protein domains were enriched in sex-biased genes, some of which hint at sex-biased physiological processes.

Conclusions

These findings lend insight into the regulatory and evolutionary forces shaping sex-biased gene expression and their involvement in the physiological and pathological processes in human health and diseases.

Electronic supplementary material

The online version of this article (10.1186/s13293-017-0156-4) contains supplementary material, which is available to authorized users.

Mutation spectrum of RB1 mutations in retinoblastoma cases from Singapore with implications for genetic management and counselling

Retinoblastoma (RB) is a rare childhood malignant disorder caused by the biallelic inactivation of RB1 gene. Early diagnosis and identification of carriers of heritable RB1 mutations can improve disease outcome and management. In this study, mutational analysis was conducted on fifty-nine matched tumor and peripheral blood samples from 18 bilateral and 41 unilateral unrelated RB cases by a combinatorial approach of Multiplex Ligation-dependent Probe Amplification (MLPA) assay, deletion screening, direct sequencing, copy number gene dosage analysis and methylation assays. Screening of both blood and tumor samples yielded a mutation detection rate of 94.9% (56/59) while only 42.4% (25/59) of mutations were detected if blood samples alone were analyzed. Biallelic mutations were observed in 43/59 (72.9%) of tumors screened. There were 3 cases (5.1%) in which no mutations could be detected and germline mutations were detected in 19.5% (8/41) of unilateral cases. A total of 61 point mutations were identified, of which 10 were novel. There was a high incidence of previously reported recurrent mutations, occurring at 38.98% (23/59) of all cases. Of interest were three cases of mosaic RB1 mutations detected in the blood from patients with unilateral retinoblastoma. Additionally, two germline mutations previously reported to be associated with low-penetrance phenotypes: missense-c.1981C>T and splice variant-c.607+1G>T, were observed in a bilateral and a unilateral proband, respectively. These findings have implications for genetic counselling and risk prediction for the affected families. This is the first published report on the spectrum of mutations in RB patients from Singapore and shows that further improved mutation screening strategies are required in order to provide a definitive molecular diagnosis for every case of RB. Our findings also underscore the importance of genetic testing in supporting individualized disease management plans for patients and asymptomatic family members carrying low-penetrance, germline mosaicism or heritable unilateral mutational phenotypes.

Sex-specific effect of CPB2 Ala147Thr but not Thr325Ile variants on the risk of venous thrombosis: A comprehensive meta-analysis

BACKGROUND

Thrombin activatable fibrinolysis inhibitor (TAFI), encoded by the Carboxypeptidase B2 gene (CPB2), is an inhibitor of fibrinolysis and plays a role in the pathogenesis of venous thrombosis. Experimental findings support a functional role of genetic variants in CPB2, while epidemiological studies have been unable to confirm associations with risk of venous thrombosis. Sex-specific effects could underlie the observed inconsistent associations between CPB2 genetic variants and venous thrombosis.

METHODS

A comprehensive literature search was conducted for associations between Ala147Thr and Thr325Ile variants with venous thrombosis. Authors were contacted to provide sex-specific genotype counts from their studies. Combined and sex-specific random effects meta-analyses were used to estimate a pooled effect estimate for primary and secondary genetic models.

RESULTS

A total of 17 studies met the inclusion criteria. A sex-specific meta-analysis applying a dominant model supported a protective effect of Ala147Thr on venous thrombosis in females (OR = 0.81, 95%CI: 0.68,0.97; p = 0.018), but not in males (OR = 1.06, 95%CI:0.96-1.16; p = 0.263). The Thr325Ile did not show a sex-specific effect but showed variation in allele frequencies by geographic region. A subgroup analysis of studies in European countries showed decreased risk, with a recessive model (OR = 0.83, 95%CI:0.71-0.97, p = 0.021) for venous thrombosis.

CONCLUSIONS

A comprehensive literature review, including unpublished data, provided greater statistical power for the analyses and decreased the likelihood of publication bias influencing the results. Sex-specific analyses explained apparent discrepancies across genetic studies of Ala147Thr and venous thrombosis. While, careful selection of genetic models based on population genetics, evolutionary and biological knowledge can increase power by decreasing the need to adjust for testing multiple models.

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.

The rs61764370 Functional Variant in the KRAS Oncogene is Associated with Chronic Myeloid Leukemia Risk in Women

BACKGROUND

Chronic myeloid leukemia (CML) is one of the most frequent hematopoietic malignancies in the elderly population; however, knowledge is limited regarding the genetic factors associated with increased risk for CML. Polymorphisms affecting microRNA (miRNA) biogenesis or mRNA:miRNA interactions are important risk factors in the development of different types of cancer. Thus, we carried out a case-control study to test the association with CML susceptibility of gene variants located in the miRNA machinery genes AGO1 (rs636832) and GEMIN4 (rs2740348), as well as in the miRNA binding sites of the genes BRCA1 (rs799917) and KRAS (rs61764370).

MATERIALS AND METHODS

We determined the genotype of 781 Mexican-Mestizo individuals (469 healthy subjects and 312 CML cases) for the four polymorphisms using TaqMan probes to test the association with CML susceptibility.

RESULTS

We found a borderline association of the minor homozygote genotype of the KRAS_rs61764370 polymorphism with an increased risk for CML susceptibility (P = 0.06). After gender stratification, this association was significant only for women (odds ratio [OR] = 13.41, P = 0.04). The distribution of the allelic and genotypic frequencies of the four studied SNPs was neither associated with advanced phases of CML nor treatment response.

CONCLUSIONS

To the best of our knowledge, this study is the first to show a significant association of the KRAS_rs61764370 SNP with CML. To further determine such an association of with CML susceptibility, our results must be replicated in different ethnic groups.

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

BACKGROUND

Despite their nearly identical genomes, males and females differ in risk, incidence, prevalence, severity and age-at-onset of many diseases. Sexual dimorphism is also seen in human autosomal gene expression, and has largely been explored by examining the contribution of genotype-by-sex interactions to variation in gene expression.

RESULTS

In this study, we use data from a mixture of pedigree and unrelated individuals with verified European ancestry to investigate the sex-specific genetic architecture of gene expression measured in whole blood across n=1048 males and n=1005 females by treating gene expression intensities in the sexes as two distinct traits and estimating the genetic correlation (r _G) between them. These correlations measure the similarity of the combined additive genetic effects of all single-nucleotide polymorphisms across the autosomal chromosomes, and thus the level of common genetic control of gene expression across the sexes. Genetic correlations are estimated across the sexes for the expression levels of 12,528 autosomal gene expression probes using bivariate GREML, and tested for differences in autosomal genetic control of gene expression across the sexes. Overall, no deviation of the distribution of test statistics is observed from that expected under the null hypothesis of a common autosomal genetic architecture for gene expression across the sexes.

CONCLUSIONS

These results suggest that males and females share the same common genetic control of gene expression.

A childhood acute lymphoblastic leukemia genome-wide association study identifies novel sex-specific risk variants

Childhood acute lymphoblastic leukemia (ALL) occurs more frequently in males. Reasons behind sex differences in childhood ALL risk are unknown. In the present genome-wide association study (GWAS), we explored the genetic basis of sex differences by comparing genotype frequencies between male and female cases in a case-only study to assess effect-modification by sex.The case-only design included 236 incident cases of childhood ALL consecutively recruited at the Texas Children's Cancer Center in Houston, Texas from 2007 to 2012. All cases were non-Hispanic whites, aged 1 to 10 years, and diagnosed with confirmed B-cell precursor ALL. Genotyping was performed using the Illumina HumanCoreExome BeadChip on the Illumina Infinium platform. Besides the top 100 statistically most significant results, results were also analyzed by the top 100 highest effect size with a nominal statistical significance (P <0.05).The statistically most significant sex-specific association (P = 4 × 10) was with the single nucleotide polymorphism (SNP) rs4813720 (RASSF2), an expression quantitative trait locus (eQTL) for RASSF2 in peripheral blood. rs4813720 is also a strong methylation QTL (meQTL) for a CpG site (cg22485289) within RASSF2 in pregnancy, at birth, childhood, and adolescence. cg22485289 is one of the hypomethylated CpG sites in ALL compared with pre-B cells. Two missense SNPs, rs12722042 and 12722039, in the HLA-DQA1 gene yielded the highest effect sizes (odds ratio [OR] ∼ 14; P <0.01) for sex-specific results. The HLA-DQA1 SNPs belong to DQA1*01 and confirmed the previously reported male-specific association with DQA1*01. This finding supports the proposed infection-related etiology in childhood ALL risk for males. Further analyses revealed that most SNPs (either direct effect or through linkage disequilibrium) were within active enhancers or active promoter regions and had regulatory effects on gene expression levels.Cumulative data suggested that RASSF2 rs4813720, which correlates with increased RASSF2 expression, may counteract the suppressor effect of estrogen-regulated miR-17-92 on RASSF2 resulting in protection in males. Given the amount of sex hormone-related mechanisms suggested by our findings, future studies should examine prenatal or early postnatal programming by sex hormones when hormone levels show a large variation.

Functional genomics bridges the gap between quantitative genetics and molecular biology

Deep characterization of molecular function of genetic variants in the human genome is becoming increasingly important for understanding genetic associations to disease and for learning to read the regulatory code of the genome. In this paper, I discuss how recent advances in both quantitative genetics and molecular biology have contributed to understanding functional effects of genetic variants, lessons learned from eQTL studies, and future challenges in this field.

BAFF Expression is Modulated by Female Hormones in Human Immune Cells

Among several autoimmune diseases, one of the main risk factors is the female gender, and much consideration has been given to the involvement of female hormones in their etiology. B-cell activating factor (BAFF) is a key factor in survival and maturation of B cells and is overexpressed in several autoimmune patients although the mechanism behind this feature is unclear. In murine models, BAFF expression could be upregulated by exogenous estrogen treatment in splenocytes; however, no evidence of this relationship was available in humans. Here, leukocytes from healthy male and female individuals were collected and cultivated in the presence or absence of estrogen or progesterone. BAFF gene expression was accessed by quantitative PCR and compared between treated and untreated group of cells. In the presence of estrogen, BAFF expression was upregulated by more than 5 times in both genders. When exposed to progesterone, the female-originated cells showed increased expression, while the cells of male origin a significant downregulation of BAFF. Our results suggest that female hormones can modulate the expression of BAFF, a key cytokine in autoimmune pathways, in human immune cells. These data might contribute to the understanding of the etiology as well as the gender bias featured by several autoimmune disorders.

Impact of the X Chromosome and sex on regulatory variation

The X Chromosome, with its unique mode of inheritance, contributes to differences between the sexes at a molecular level, including sex-specific gene expression and sex-specific impact of genetic variation. Improving our understanding of these differences offers to elucidate the molecular mechanisms underlying sex-specific traits and diseases. However, to date, most studies have either ignored the X Chromosome or had insufficient power to test for the sex-specific impact of genetic variation. By analyzing whole blood transcriptomes of 922 individuals, we have conducted the first large-scale, genome-wide analysis of the impact of both sex and genetic variation on patterns of gene expression, including comparison between the X Chromosome and autosomes. We identified a depletion of expression quantitative trait loci (eQTL) on the X Chromosome, especially among genes under high selective constraint. In contrast, we discovered an enrichment of sex-specific regulatory variants on the X Chromosome. To resolve the molecular mechanisms underlying such effects, we generated chromatin accessibility data through ATAC-sequencing to connect sex-specific chromatin accessibility to sex-specific patterns of expression and regulatory variation. As sex-specific regulatory variants discovered in our study can inform sex differences in heritable disease prevalence, we integrated our data with genome-wide association study data for multiple immune traits identifying several traits with significant sex biases in genetic susceptibilities. Together, our study provides genome-wide insight into how genetic variation, the X Chromosome, and sex shape human gene regulation and disease.

Imputing Gene Expression in Uncollected Tissues Within and Beyond GTEx

Gene expression and its regulation can vary substantially across tissue types. In order to generate knowledge about gene expression in human tissues, the Genotype-Tissue Expression (GTEx) program has collected transcriptome data in a wide variety of tissue types from post-mortem donors. However, many tissue types are difficult to access and are not collected in every GTEx individual. Furthermore, in non-GTEx studies, the accessibility of certain tissue types greatly limits the feasibility and scale of studies of multi-tissue expression. In this work, we developed multi-tissue imputation methods to impute gene expression in uncollected or inaccessible tissues. Via simulation studies, we showed that the proposed methods outperform existing imputation methods in multi-tissue expression imputation and that incorporating imputed expression data can improve power to detect phenotype-expression correlations. By analyzing data from nine selected tissue types in the GTEx pilot project, we demonstrated that harnessing expression quantitative trait loci (eQTLs) and tissue-tissue expression-level correlations can aid imputation of transcriptome data from uncollected GTEx tissues. More importantly, we showed that by using GTEx data as a reference, one can impute expression levels in inaccessible tissues in non-GTEx expression studies.

Non-negligible Occurrence of Errors in Gender Description in Public Data Sets

Due to advances in omics technologies, numerous genome-wide studies on human samples have been published, and most of the omics data with the associated clinical information are available in public repositories, such as Gene Expression Omnibus and ArrayExpress. While analyzing several public datasets, we observed that errors in gender information occur quite often in public datasets. When we analyzed the gender description and the methylation patterns of gender-specific probes (glucose-6-phosphate dehydrogenase [G6PD], ephrin-B1 [EFNB1], and testis specific protein, Y-linked 2 [TSPY2]) in 5,611 samples produced using Infinium 450K HumanMethylation arrays, we found that 19 samples from 7 datasets were erroneously described. We also analyzed 1,819 samples produced using the Affymetrix U133Plus2 array using several gender-specific genes (X (inactive)-specific transcript [XIST], eukaryotic translation initiation factor 1A, Y-linked [EIF1AY], and DEAD [Asp-Glu-Ala-Asp] box polypeptide 3, Y-linked [DDDX3Y]) and found that 40 samples from 3 datasets were erroneously described. We suggest that the users of public datasets should not expect that the data are error-free and, whenever possible, that they should check the consistency of the data.