Sex differences in the human peripheral blood transcriptome

Combined transcriptomics and proteomics unveil the impact of vitamin C in modulating specific protein abundance in the mouse liver

BACKGROUND

Vitamin C (ascorbate) is a water-soluble antioxidant and an important cofactor for various biosynthetic and regulatory enzymes. Mice can synthesize vitamin C thanks to the key enzyme gulonolactone oxidase (Gulo) unlike humans. In the current investigation, we used Gulo-/- mice, which cannot synthesize their own ascorbate to determine the impact of this vitamin on both the transcriptomics and proteomics profiles in the whole liver. The study included Gulo-/- mouse groups treated with either sub-optimal or optimal ascorbate concentrations in drinking water. Liver tissues of females and males were collected at the age of four months and divided for transcriptomics and proteomics analysis. Immunoblotting, quantitative RT-PCR, and polysome profiling experiments were also conducted to complement our combined omics studies.

RESULTS

Principal component analyses revealed distinctive differences in the mRNA and protein profiles as a function of sex between all the mouse cohorts. Despite such sexual dimorphism, Spearman analyses of transcriptomics data from females and males revealed correlations of hepatic ascorbate levels with transcripts encoding a wide array of biological processes involved in glucose and lipid metabolisms as well as in the acute-phase immune response. Moreover, integration of the proteomics data showed that ascorbate modulates the abundance of various enzymes involved in lipid, xenobiotic, organic acid, acetyl-CoA, and steroid metabolism mainly at the transcriptional level, especially in females. However, several proteins of the mitochondrial complex III significantly correlated with ascorbate concentrations in both males and females unlike their corresponding transcripts. Finally, poly(ribo)some profiling did not reveal significant enrichment difference for these mitochondrial complex III mRNAs between Gulo-/- mice treated with sub-optimal and optimal ascorbate levels.

CONCLUSIONS

Thus, the abundance of several subunits of the mitochondrial complex III are regulated by ascorbate at the post-transcriptional levels. Our extensive omics analyses provide a novel resource of altered gene expression patterns at the transcriptional and post-transcriptional levels under ascorbate deficiency.

Genome-wide study of gene-by-sex interactions identifies risks for cleft palate

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (GxS) interaction testing. There were 13 loci significant for GxS interactions, with the top finding in LTBP1 (RR=3.37 [2.04 - 5.56], p=1.93x10 -6 ). LTBP1 plays a role in regulating TGF-B bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant GxS interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p=0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins

BACKGROUND

The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remains underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers.

METHODS

Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N = 651) and the Netherlands Twin Register (NTR) (N = 665). Follow-up comprised 4 BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated in latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. In FinnTwin12, the sources of genetic and environmental variation underlying the protein abundances were quantified by twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) applying mixed-effects models and correlation networks.

RESULTS

We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 7 and 3 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers.

CONCLUSIONS

Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

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.

Female cardiovascular biology and resilience in the setting of physiological and pathological stress

For years, females were thought of as smaller men with complex hormonal cycles; as a result, females have been largely excluded from preclinical and clinical research. However, in the last ten years, with the increased focus on sex as a biological variable, it has become clear that this is not the case, and in fact, male and female cardiovascular biology and cardiac stress responses differ substantially. Premenopausal women are protected from cardiovascular diseases, such as myocardial infarction and resultant heart failure, having preserved cardiac function, reduced adverse remodeling, and increased survival. Many underlying biological processes that contribute to ventricular remodeling differ between the sexes, such as cellular metabolism; immune cell responses; cardiac fibrosis and extracellular matrix remodeling; cardiomyocyte dysfunction; and endothelial biology; however, it is unclear how these changes afford protection to the female heart. Although many of these changes are dependent on protection provided by female sex hormones, several of these changes occur independent of sex hormones, suggesting that the nature of these changes is more complex than initially thought. This may be why studies focused on the cardiovascular benefits of hormone replacement therapy in post-menopausal women have provided mixed results. Some of the complexity likely stems from the fact that the cellular composition of the heart is sexually dimorphic and that in the setting of MI, different subpopulations of these cell types are apparent. Despite the documented sex-differences in cardiovascular (patho)physiology, the underlying mechanisms that contribute are largely unknown due to inconsistent findings amongst investigators and, in some cases, lack of rigor in reporting and consideration of sex-dependent variables. Therefore, this review aims to describe current understanding of the sex-dependent differences in the myocardium in response to physiological and pathological stressors, with a focus on the sex-dependent differences that contribute to post-infarction remodeling and resultant functional decline.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins

Background

The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remain underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers.

Methods

Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N=651) and the Netherlands Twin Register (NTR) (N=665). Follow-up comprised four BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated using latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. The sources of genetic and environmental variation underlying the protein abundances were quantified using twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) using mixed-effect models and correlation networks.

Results

We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 6 and 4 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with many metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers.

Conclusions

Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Association of Circadian Clock Gene Expression with Pediatric/Adolescent Asthma and Its Comorbidities

The pathology of asthma is characterized by marked day-night variation, which is likely controlled by circadian clock activity. This study aimed to clarify the association of core circadian clock gene expression with clinical features of asthma. For this purpose, we accessed the National Center for Biotechnology Information database and analyzed transcriptomes of peripheral blood mononuclear cells and clinical characteristics of 134 pediatric/adolescent patients with asthma. Based on the expression patterns of seven core circadian clock genes (CLOCK, BMAL1, PER1-3, CRY1-2), we identified three circadian clusters (CCs) with distinct comorbidities and transcriptomic expressions. In the three CC subtypes, allergic rhinitis, and atopic dermatitis, both asthma comorbidities occurred in different proportions: CC1 had a high proportion of allergic rhinitis and atopic dermatitis; CC2 had a high proportion of atopic dermatitis but a low proportion of allergic rhinitis; and CC3 had a high proportion of allergic rhinitis but a low proportion of atopic dermatitis. This might be associated with the low activity of the FcεRI signaling pathway in CC2 and the cytokine-cytokine receptor interaction pathways in CC3. This is the first report to consider circadian clock gene expression in subcategories of patients with asthma and to explore their contribution to pathophysiology and comorbidity.

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 differences in neuroimmunoendocrine communication. Involvement on longevity

Endocrine, nervous, and immune systems work coordinately to maintain the global homeostasis of the organism. They show sex differences in their functions that, in turn, contribute to sex differences beyond reproductive function. Females display a better control of the energetic metabolism and improved neuroprotection and have more antioxidant defenses and a better inflammatory status than males, which is associated with a more robust immune response than that of males. These differences are present from the early stages of life, being more relevant in adulthood and influencing the aging trajectory in each sex and may contribute to the different life lifespan between sexes.

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

Blood transcriptome comparison between sexes and their function in 4-week Rhode Island red chickens

Sex is a major biological factor in the development and physiology of a sexual reproductive organism, and its role in the growing process is needed to be investigated in various species. We compare blood transcriptome between 5 males and 5 females in 4-week-old Rhode Island Red chickens and perform functional annotation of differentially expressed genes (DEGs). The results are as follows. 141 and 109 DEGs were located in autosomes and sex chromosomes, respectively. The gene ontology (GO) terms are significantly (p < 0.05) enriched, which were limb development, inner ear development, positive regulation of dendrite development, the KEGG pathway the TGF-beta signaling pathway, and melanogenesis (p < 0.05). These pathways are related to morphological maintenance and growth of the tissues. In addition, the SMAD2W and the BMP5 were involved in the TGF-beta signaling pathway, and both play an important role in maintaining tissue development. The major DEGs related to the development of neurons and synapses include the up-regulated NRN1, GDF10, SLC1A1, BMP5, NBEA, and NRXN1. Also, 7 DEGs were validated using RT-qPCR with high correlation (r ² = 0.74). In conclusion, the differential expression of blood tissue in the early growing chicken was enriched in TGF-beta signaling and related to the development of neurons and synapses including SMAD2W and BMP5. These results suggest that blood in the early growing stage is differentially affected in tissue development, nervous system, and pigmentation by sex. For future research, experimental characterization of DEGs and a holistic investigation of various tissues and growth stages will be required.

Alzheimer's disease-related transcriptional sex differences in myeloid cells

Sex differences have been identified in many diseases associated with dysregulated immune responses, including Alzheimer's disease (AD), for which approximately two-thirds of patients are women. An accumulating body of research indicates that microglia may play a causal role in the pathogenesis of this disease. We hypothesised that sex differences in the transcriptome of human myeloid cells may contribute to the sex difference observed in AD prevalence. To explore this, we assessed bulk and single-nuclear RNA sequencing data sets generated from four human derived myeloid cell populations: post-mortem microglial nuclei, peripheral monocytes, monocyte-derived macrophages (MDMs) and induced pluripotent stem cell derived microglial-like cells (MGLs). We found that expression of AD risk genes, gene signatures associated with the inflammatory response in AD, and genes related to proinflammatory immune responses were enriched in microglial nuclei isolated from aged female donors without ante-mortem neurological disease, relative to those from males. In addition, these inflammation-associated gene sets were found to be enriched in peripheral monocytes isolated from postmenopausal women and in MDMs obtained from premenopausal individuals relative to age-matched males. Expression of these gene sets did not differ in MDMs derived from women whose blood was sampled across the menstrual cycle or in MGLs cultured with 17β-oestradiol. This suggests that the observed gene set enrichments in myeloid cells from women were not being driven by acute hormonal influences. Together, these data support the hypothesis that the increased prevalence of AD in women may be partly explained by a myeloid cell phenotype biased towards expression of biological processes relevant to AD.

Impact of long-term exposure to PM2.5 on peripheral blood gene expression pathways involved in cell signaling and immune response

BACKGROUND

Exposure to ambient air pollution, even at low levels, is a major environmental health risk. The peripheral blood transcriptome provides a potential avenue for the elucidation of ambient air pollution related biological perturbations. We assessed the association between long-term estimates for seven priority air pollutants and perturbations in peripheral blood transcriptomics data collected in the Dutch National Twin Register (NTR) and Netherlands Study of Depression and Anxiety (NESDA) cohorts.

METHODS

In both the discovery (n = 2438) and replication (n = 1567) cohort, outdoor concentration of 7 air pollutants (NO₂, NO_x, particulate matter (PM_2.5, PM_2.5abs, PM₁₀, PM_coarse)_, and ultrafine particles) was predicted with land use regression models. Gene expression was assessed by Affymetrix U219 arrays. Multi-variable univariate mixed-effect models were applied to test for an association between the air pollutants and the transcriptome. Functional analysis was conducted in DAVID.

RESULTS

In the discovery cohort, we observed for 335 genes (374 probes with FDR < 5 %) a perturbation in peripheral blood gene expression that was associated with long-term average levels of PM_2.5. For 69 genes pooled effect estimates from the NTR and NESDA cohorts were significant. Identified genes play a role in biological pathways related to cell signaling and immune response. Sixty-two out of 69 genes had a similar direction of effect in an analysis in which we regressed the probes on differential PM_2.5 exposure within monozygotic twin pairs, indicating that the observed differences in gene expression were likely driven by differences in air pollution, rather than by confounding by genetic factors.

CONCLUSION

Our results indicate that PM_2.5 can elicit a response in cell signaling and the immune system, both hallmarks of environmental diseases. The differential effect that we observed between air pollutants may aid in the understanding of differential health effects that have been observed with these exposures.

17β-Oestradiol Protects from Hepatitis C Virus Infection through Induction of Type I Interferon

Background and Aims: Sex hormones are widely recognised to act as protective factors against several viral infections. Specifically, females infected by the hepatitis C virus display higher clearance rates and reduced disease progression than those found in males. Through modulation of particle release and spread, 17β-oestradiol controls HCV’s life cycle. We investigated the mechanism(s) behind oestrogen’s antiviral effect. Methods: We used cell culture-derived hepatitis C virus in in vitro assays to evaluate the effect of 17β-oestradiol on the innate immune response. Host immune responses were evaluated by enumerating gene transcripts via RT-qPCR in cells exposed to oestrogen in the presence or absence of viral infection. Antiviral effects were determined by focus-forming unit assay or HCV RNA quantification. Results: Stimulation of 17β-oestradiol triggers a pre-activated antiviral state in hepatocytes, which can be maintained for several hours after the hormone is removed. This induction results in the elevation of several innate immune genes, such as interferon alpha and beta, tumour necrosis factor, toll-like receptor 3 and interferon regulatory factor 5. We demonstrated that this pre-activation of immune response signalling is not affected by a viral presence, and the antiviral state can be ablated using an interferon-alpha/beta receptor alpha inhibitor. Finally, we proved that the oestrogen-induced stimulation is essential to generate an antiviral microenvironment mediated by activation of type I interferons. Conclusion: Resulting in viral control and suppression, 17β-oestradiol induces an interferon-mediated antiviral state in hepatocytes. Oestrogen-stimulated cells modulate the immune response through secretion of type I interferon, which can be countered by blocking interferon-alpha/beta receptor alpha signalling.

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.

The influence of age and sex in sarcoidosis

PURPOSE OF REVIEW

This review aims to describe how the clinical manifestations of sarcoidosis may be shaped by the effects of sex hormones and by age dependent changes in immune functions and physiology This review is intended to highlight the need to consider the effects of sex and sex in future studies of sarcoidosis.

RECENT FINDINGS

The clinical manifestations of sarcoidosis differ based on sex and gender There is emerging evidence that female and male hormones and X-linked genes are important determinants of immune responses to environmental antigens, which has important implications for granuloma formation in the context of sarcoidosis Furthermore, sex hormone levels predictably change throughout adolescence and adulthood, and this occurs in parallel with the onset immune senescence and changes in physiology with advanced age.

SUMMARY

Recent studies indicate that sex and age are important variables shaping the immune response of humans to environmental antigens We posit herein that sex and age are important determinants of sarcoidosis clinical phenotypes Many gaps in our understanding of the roles played by sex and gender in sarcoidosis, and these need to be considered in future studies.

Sex disparities in influenza: A multiscale network analysis

Sex differences in the pathogenesis of infectious diseases because of differential immune responses between females and males have been well-documented for multiple pathogens. However, the molecular mechanism underlying the observed sex differences in influenza virus infection remains poorly understood. In this study, we used a network-based approach to characterize the blood transcriptome collected over the course of infection with influenza A virus from female and male ferrets to dissect sex-biased gene expression. We identified significant differences in the temporal dynamics and regulation of immune responses between females and males. Our results elucidate sex-differentiated pathways involved in the unfolded protein response (UPR), lipid metabolism, and inflammatory responses, including a female-biased IRE1/XBP1 activation and male-biased crosstalk between metabolic reprogramming and IL-1 and AP-1 pathways. Overall, our study provides molecular insights into sex differences in transcriptional regulation of immune responses and contributes to a better understanding of sex biases in influenza pathogenesis.

•

Regulation of immune responses between females and males is significantly different

•

Rapid activation of UPR in females triggers potent immune and inflammatory responses

•

Male-specific regulatory pattern in the AP1 pathway indicate a bias in immune response

Serum starvation induces sexual dimorphisms in secreted proteins of human umbilical vein endothelial cells (HUVECs) from twin pairs

There is growing evidence for sex and gender differences in the clinical manifestation and outcomes of human diseases. Human primary endothelial cells represent a useful cardiovascular model to study sexual dimorphisms at the cellular level. Here, we analyzed sexual dimorphisms of the secretome after serum starvation using HUVECs from twin pairs of the opposite sex to minimize the impact of varying genetic background. HUVECs were starved for 5 and 16 h, respectively, and proteins of the cell culture supernatants were analyzed by tandem mass spectrometry. Altogether, 960 extracellular proteins were identified of which 683 were amendable to stringent quantification. Significant alterations were observed for 324 proteins between long-term and short-term starvation and the majority were similar in both sexes. Only 5 proteins showed significant sex-specific regulation between long- versus short-term starvation. Furthermore, 19 unique proteins with significant sexual dimorphisms at the same time points of serum starvation were observed. A larger number of proteins, e.g. tissue factor inhibitor 2 (TFPI2), displayed higher levels in the supernatants of females compared to male cells after long term serum starvation that might point to higher adaptation capacity of female cells. The overall results demonstrate that male and female cells differ in their secretome. This article is protected by copyright. All rights reserved.

Long non-coding RNA Xist contribution in systemic lupus erythematosus and rheumatoid arthritis

Growing evidence points towards the role of the long non-coding (lnc)-RNA Xist expressed in female cells as a predominant key actor for the sex bias observed in systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Indeed, in female cells, lnc-Xist controls transcription directly by spreading across the inactivated X chromosome (Xi) and indirectly by sequestring miRNAs as a sponge. The inactivation process at Xi is altered in lymphocytes from SLE women and associated with important variations in ribonucleoproteins (RNP) associated with lnc-Xist. In fibroblast-like synoviocytes (FLS) and osteoclasts from RA women, proinflammatory and proliferative pathways are upregulated due to the sequestration effect exerted by lnc-Xist overexpression on miRNAs. The key role played by lnc-Xist in SLE and RA is further supported by it's knock down that recapitulates the SLE B cell extrafollicular profile and controls RA associated FLS proinflammatory cytokine production and proliferation.

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.

Gender Variations in the Oral Microbiomes of Elderly Patients with Initial Periodontitis

Periodontitis is a globally prevalent disease that imposes a functional and aesthetic burden on patients. The oral microbiome influences human health. The aim of this study was at assessing gender variation in the subgingival bacterial microbiome of elderly patients with initial periodontitis and to determine the causes of this variation. Twelve males and twenty females (range 50-68 years old) with initial periodontitis provided subgingival plaque samples. 16S rRNA gene sequencing, QIIME-based data processing, and statistical analyses were carried out using several different analytical approaches to detect differences in the oral microbiome between the two groups. Males had higher Chao1 index, observed species, and phylogenetic diversity whole tree values than females. Analysis of β-diversity indicated that the samples were reasonably divided by the gender. The linear discriminant analysis effect size showed that the most representative biomarkers were the genus Haemophilus in males, whereas the dominant bacteria in females were Campylobacter. Kyoto Encyclopedia of Genes and Genomes analysis showed that predicting changes in the female oral microbiota may be related to the immune system and immune system diseases are the main factor in males. These data suggest that gender may be a differentiating factor in the microbial composition of subgingival plaques in elderly patients with initial periodontitis. These results could deepen our understanding of the role of gender in the oral microbiota present during initial periodontitis.

Inflammation and Mitochondrial Dysfunction in Autism Spectrum Disorder

Autism Spectrum Disorders (ASD) is a severe childhood psychiatric condition with an array of cognitive, language and social impairments that can significantly impact family life. ASD is classically characterized by reduced communication skills and social interactions, with limitations imposed by repetitive patterns of behavior, interests, and activities. The pathophysiology of ASD is thought to arise from complex interactions between environmental and genetic factors within the context of individual development. A growing body of research has raised the possibility of identifying the aetiological causes of the disorder. This review highlights the roles of immune-inflammatory pathways, nitro-oxidative stress and mitochondrial dysfunctions in ASD pathogenesis and symptom severity. The role of NK-cells, T helper, T regulatory and B-cells, coupled with increased inflammatory cytokines, lowered levels of immune-regulatory cytokines, and increased autoantibodies and microglial activation is elucidated. It is proposed that alterations in mitochondrial activity and nitrooxidative stress are intimately associated with activated immune-inflammatory pathways. Future research should determine as to whether the mitochondria, immune-inflammatory activity and nitrooxidative stress changes in ASD affect the development of amygdala-frontal cortex interactions. A number of treatment implications may arise, including prevention-orientated prenatal interventions, treatment of pregnant women with vitamin D, and sodium butyrate. Treatments of ASD children and adults with probiotics, sodium butyrate and butyrate-inducing diets, antipurinergic therapy with suramin, melatonin, oxytocin and taurine are also discussed.

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 of monocyte transcriptome in individuals with chronic low-grade inflammation

Sexual dimorphism in the immune system is evidenced by a higher prevalence of autoimmune diseases in women and higher susceptibility to infectious diseases in men. However, the molecular basis of these sex-based differences is not fully understood. We have characterized the transcriptome profiles of peripheral blood monocytes from males and postmenopausal females with chronic low-grade inflammation. We identified 41 sexually differentially expressed genes [adjusted p value (FDR) < 0.1], including genes involved in immune cell activation (e.g., CEACAM1, FCGR2B, and SLAMF7) and antigen presentation (e.g., AIM2, CD1E, and UBA1) with a higher expression in females than males. Moreover, signaling pathways of immune or inflammatory responses, including interferon (IFN) signaling [z-score = 2.45, -log(p) = 3.88], were found to be more upregulated in female versus male monocytes, based on a set of genes exhibiting sex-biased expression (p < 0.03). The contribution of IFN signaling to the sexual transcriptional differences was further confirmed by direct comparisons of the monocyte sex-biased genes with IFN signature genes (ISGs) that were previously curated in mouse macrophages. ISGs showed a greater overlap with female-biased genes than male-biased genes and a higher overall expression in female than male monocytes, particularly for the genes of antiviral and inflammatory responses to IFN. Given the role of IFN in immune defense and autoimmunity, our results suggest that sexual dimorphism in immune functions may be associated with more priming of innate immune pathways in female than male monocytes. These findings highlight the role of sex on the human immune transcriptome.

No evidence for differential gene expression in major depressive disorder PBMCs, but robust evidence of elevated biological ageing

The increasingly compelling data supporting the involvement of immunobiological mechanisms in Major Depressive Disorder (MDD) might provide some explanation forthe variance in this heterogeneous condition. Peripheral blood measures of cytokines and chemokines constitute the bulk of evidence, with consistent meta-analytic data implicating raised proinflammatory cytokines such as IL6, IL1β and TNF. Among the potential mechanisms linking immunobiological changes to affective neurobiology is the accelerated biological ageing seen in MDD, particularly via the senescence associated secretory phenotype (SASP). However, the cellular source of immunobiological markers remains unclear. Pre-clinical evidence suggests a role for peripheral blood mononuclear cells (PBMC), thus here we aimed to explore the transcriptomic profile using RNA sequencing in PBMCs in a clinical sample of people with various levels of depression and treatment response comparing it with that in healthy controls (HCs). There were three groups with major depressive disorder (MDD): treatment-resistant (n = 94), treatment-responsive (n = 47) and untreated (n = 46). Healthy controls numbered 44. Using PBMCs gene expression analysis was conducted using RNAseq to a depth of 54.5 million reads. Differential gene expression analysis was performed using DESeq2. The data showed no robust signal differentiating MDD and HCs. There was, however, significant evidence of elevated biological ageing in MDD vs HC. Biological ageing was evident in these data as a transcriptional signature of 888 age-associated genes (adjusted p < 0.05, absolute log2fold > 0.6) that also correlated strongly with chronological age (spearman correlation coefficient of 0.72). Future work should expand clinical sample sizes and reduce clinical heterogeneity. Exploration of RNA-seq signatures in other leukocyte populations and single cell RNA sequencing may help uncover more subtle differences. However, currently the subtlety of any PBMC signature mitigates against its convincing use as a diagnostic or predictive biomarker.

Y disruption, autosomal hypomethylation and poor male lung cancer survival

Lung cancer is the most frequent cause of cancer death worldwide. It affects more men than women, and men generally have worse survival outcomes. We compared gene co-expression networks in affected and unaffected lung tissue from 126 consecutive patients with Stage IA-IV lung cancer undergoing surgery with curative intent. We observed marked degradation of a sex-associated transcription network in tumour tissue. This disturbance, detected in 27.7% of male tumours in the discovery dataset and 27.3% of male tumours in a further 123-sample replication dataset, was coincident with partial losses of the Y chromosome and extensive autosomal DNA hypomethylation. Central to this network was the epigenetic modifier and regulator of sexually dimorphic gene expression, KDM5D. After accounting for prognostic and epidemiological covariates including stage and histology, male patients with tumour KDM5D deficiency showed a significantly increased risk of death (Hazard Ratio [HR] 3.80, 95% CI 1.40-10.3, P = 0.009). KDM5D deficiency was confirmed as a negative prognostic indicator in a further 1100 male lung tumours (HR 1.67, 95% CI 1.4-2.0, P = 1.2 × 10^-10). Our findings identify tumour deficiency of KDM5D as a prognostic marker and credible mechanism underlying sex disparity in lung cancer.

Sex-specific gene expression in the blood of four primates

Blood is an important non-reproductive tissue, but little is known about the sex-specific gene expressions in the blood. Therefore, we investigated sex-specific gene expression differences in the blood tissues of four primates, rhesus macaques (Macaca mulatta), Tibetan macaques (M. thibetana), yellow baboons (Papio cynocephalus), and humans. We identified seven sex-specific differentially expressed genes (SDEGs) in each non-human primate and 31 SDEGs in humans. The four primates had only one common SDEG, MAP7D2. In humans, immune-related SDEGs were identified as up-regulated, but also down-regulated in females. We also found that most of the X-Y gene pairs had similar expression levels between species, except pair EIF1AY/EIF1AX. The expression level of X-Y gene pairs of rhesus and Tibetan macaques showed no significant differential expression levels, while humans had six significant XY-biased and three XX-biased X-Y gene pairs. Our observed sex differences in blood should increase understanding of sex differences in primate blood tissue.

Mapping the transcriptomics landscape of post-traumatic stress disorder symptom dimensions in World Trade Center responders

Gene expression has provided promising insights into the pathophysiology of post-traumatic stress disorder (PTSD); however, specific regulatory transcriptomic mechanisms remain unknown. The present study addressed this limitation by performing transcriptome-wide RNA-Seq of whole-blood samples from 226 World Trade Center responders. The investigation focused on differential expression (DE) at the gene, isoform, and for the first time, alternative splicing (AS) levels associated with the symptoms of PTSD: total burden, re-experiencing, avoidance, numbing, and hyperarousal subdimensions. These symptoms were associated with 76, 1, 48, 15, and 49 DE genes, respectively (FDR < 0.05). Moreover, they were associated with 103, 11, 0, 43, and 32 AS events. Avoidance differed the most from other dimensions with respect to DE genes and AS events. Gene set enrichment analysis (GSEA) identified pathways involved in inflammatory and metabolic processes, which may have implications in the treatment of PTSD. Overall, the findings shed a novel light on the wide range of transcriptomic alterations associated with PTSD at the gene and AS levels. The results of DE analysis associated with PTSD subdimensions highlights the importance of studying PTSD symptom heterogeneity.

Transcriptomic analyses of humans and mice provide insights into depression

Accumulating studies have been conducted to identify risk genes and relevant biological mechanisms underlying major depressive disorder (MDD). In particular, transcriptomic analyses in brain regions engaged in cognitive and emotional processes, e.g., the dorsolateral prefrontal cortex (DLPFC), have provided essential insights. Based on three independent DLPFC RNA-seq datasets of 79 MDD patients and 75 healthy controls, we performed differential expression analyses using two alternative approaches for cross-validation. We also conducted transcriptomic analyses in mice undergoing chronic variable stress (CVS) and chronic social defeat stress (CSDS). We identified 12 differentially expressed genes (DEGs) through both analytical methods in MDD patients, the majority of which were also dysregulated in stressed mice. Notably, the mRNA level of the immediate early gene FOS ( Fos proto-oncogene) was significantly decreased in both MDD patients and CVS-exposed mice, and CSDS-susceptible mice exhibited a greater reduction in Fos expression compared to resilient mice. These findings suggest the potential key roles of this gene in the pathogenesis of MDD related to stress exposure. Altered transcriptomes in the DLPFC of MDD patients might be, at least partially, the result of stress exposure, supporting that stress is a primary risk factor for MDD.

Single-cell transcriptomics analysis of mild cognitive impairment in World Trade Center disaster responders

INTRODUCTION

Recent research has found that World Trade Center (WTC) responders in their mid-50s have an elevated prevalence of mild cognitive impairment (MCI) that is associated with neural degeneration and subcortical thinning. This article extends our understanding of the molecular complexity of MCI through gene expression profiling of blood.

METHODS

The transcriptomics of 40 male WTC responders were profiled across two cohorts (discovery: nine MCI and nine controls; replication: 11 MCI and 11 controls) using CITE-Seq at single-cell resolution in blood.

RESULTS

Comparing the transcriptomic signatures across seven major cell subpopulations, the largest differences were observed in monocytes in which 226 genes were differentially expressed. Pathway analysis on the genes unique to monocytes identified processes associated with cerebral immune response.

DISCUSSION

Our findings suggested monocytes may constitute a key cell type to target in blood-based biomarker studies for early detection of risk of MCI and development of new interventions.

An integrative study of five biological clocks in somatic and mental health

Biological clocks have been developed at different molecular levels and were found to be more advanced in the presence of somatic illness and mental disorders. However, it is unclear whether different biological clocks reflect similar aging processes and determinants. In ~3000 subjects, we examined whether five biological clocks (telomere length, epigenetic, transcriptomic, proteomic, and metabolomic clocks) were interrelated and associated to somatic and mental health determinants. Correlations between biological aging indicators were small (all r < 0.2), indicating little overlap. The most consistent associations of advanced biological aging were found for male sex, higher body mass index (BMI), metabolic syndrome, smoking, and depression. As compared to the individual clocks, a composite index of all five clocks showed most pronounced associations with health determinants. The large effect sizes of the composite index and the low correlation between biological aging indicators suggest that one's biological age is best reflected by combining aging measures from multiple cellular levels.

Integrating biological knowledge and gene expression data using pathway-guided random forests: a benchmarking study

MOTIVATION

High-throughput technologies allow comprehensive characterization of individuals on many molecular levels. However, training computational models to predict disease status based on omics data is challenging. A promising solution is the integration of external knowledge about structural and functional relationships into the modeling process. We compared four published random forest-based approaches using two simulation studies and nine experimental datasets.

RESULTS

The self-sufficient prediction error approach should be applied when large numbers of relevant pathways are expected. The competing methods hunting and learner of functional enrichment should be used when low numbers of relevant pathways are expected or the most strongly associated pathways are of interest. The hybrid approach synthetic features is not recommended because of its high false discovery rate.

AVAILABILITY AND IMPLEMENTATION

An R package providing functions for data analysis and simulation is available at GitHub (https://github.com/szymczak-lab/PathwayGuidedRF). An accompanying R data package (https://github.com/szymczak-lab/DataPathwayGuidedRF) stores the processed and quality controlled experimental datasets downloaded from Gene Expression Omnibus (GEO).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Sexual Dimorphism in Innate Immunity: The Role of Sex Hormones and Epigenetics

Sexual dimorphism refers to differences between biological sexes that extend beyond sexual characteristics. In humans, sexual dimorphism in the immune response has been well demonstrated, with females exhibiting lower infection rates than males for a variety of bacterial, viral, and parasitic pathogens. There is also a substantially increased incidence of autoimmune disease in females compared to males. Together, these trends indicate that females have a heightened immune reactogenicity to both self and non-self-molecular patterns. However, the molecular mechanisms driving the sexually dimorphic immune response are not fully understood. The female sex hormones estrogen and progesterone, as well as the male androgens, such as testosterone, elicit direct effects on the function and inflammatory capacity of immune cells. Several studies have identified a sex-specific transcriptome and methylome, independent of the well-described phenomenon of X-chromosome inactivation, suggesting that sexual dimorphism also occurs at the epigenetic level. Moreover, distinct alterations to the transcriptome and epigenetic landscape occur in synchrony with periods of hormonal change, such as puberty, pregnancy, menopause, and exogenous hormone therapy. These changes are also mirrored by changes in immune cell function. This review will outline the evidence for sex hormones and pregnancy-associated hormones as drivers of epigenetic change, and how this may contribute to the sexual dimorphism. Determining the effects of sex hormones on innate immune function is important for understanding sexually dimorphic autoimmune diseases, sex-specific responses to pathogens and vaccines, and how innate immunity is altered during periods of hormonal change (endogenous or exogenous).

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.

The Effects of Age, Cigarette Smoking, Sex, and Race on the Qualitative Characteristics of Lung Transcriptome

The within-sample relative expression orderings (REOs) of genes, which are stable qualitative transcriptional characteristics, can provide abundant information for a disease. Methods based on REO comparisons have been proposed for identifying differentially expressed genes (DEGs) at the individual level and for detecting disease-associated genes based on one-phenotype disease data by reusing data of normal samples from other sources. Here, we evaluated the effects of common potential confounding factors, including age, cigarette smoking, sex, and race, on the REOs of gene pairs within normal lung tissues transcriptome. Our results showed that age has little effect on REOs within lung tissues. We found that about 0.23% of the significantly stable REOs of gene pairs in nonsmokers' lung tissues are reversed in smokers' lung tissues, introduced by 344 DEGs between the two groups of samples (RankCompV2, FDR <0.05), which are enriched in metabolism of xenobiotics by cytochrome P450, glutathione metabolism, and other pathways (hypergeometric test, FDR <0.05). Comparison between the normal lung tissue samples of males and females revealed fewer reversal REOs introduced by 24 DEGs between the sex groups, among which 19 DEGs are located on sex chromosomes and 5 DEGs involving in spermatogenesis and regulation of oocyte are located on autosomes. Between the normal lung tissue samples of white and black people, we identified 22 DEGs (RankCompV2, FDR <0.05) which introduced a few reversal REOs between the two races. In summary, the REO-based study should take into account the confounding factors of cigarette smoking, sex, and race.

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.

Intrinsic transcriptomic sex differences in human endothelial cells at birth and in adults are associated with coronary artery disease targets

Sex differences in endothelial cell (EC) biology may reflect intrinsic differences driven by chromosomes or sex steroid exposure and gender differences accumulated over life. We analysed EC gene expression data from boy-girl twins at birth and in non-twin adults to detect sex differences at different stages of life, and show that 14-25% of the EC transcriptome is sex-biased. By combining data from both stages of life, we identified sex differences that are present at birth and maintained throughout life, and those that are acquired over life. Promisingly, we found that genes that present with an acquired sex difference in ECs are more likely to be targets of sex steroids. Annotating both gene sets with data from multiple genome-wide association studies (GWAS) revealed that genes with an intrinsic sex difference in ECs are enriched for coronary artery disease GWAS hits. This study underscores the need for treating sex as a biological variable.

Sex Chromosomes and Sex Phenotype Contribute to Biased DNA Methylation in Mouse Liver

Sex biases in the genome-wide distribution of DNA methylation and gene expression levels are some of the manifestations of sexual dimorphism in mammals. To advance our understanding of the mechanisms that contribute to sex biases in DNA methylation and gene expression, we conducted whole genome bisulfite sequencing (WGBS) as well as RNA-seq on liver samples from mice with different combinations of sex phenotype and sex-chromosome complement. We compared groups of animals with different sex phenotypes, but the same genetic sexes, and vice versa, same sex phenotypes, but different sex-chromosome complements. We also compared sex-biased DNA methylation in mouse and human livers. Our data show that sex phenotype, X-chromosome dosage, and the presence of Y chromosome shape the differences in DNA methylation between males and females. We also demonstrate that sex bias in autosomal methylation is associated with sex bias in gene expression, whereas X-chromosome dosage-dependent methylation differences are not, as expected for a dosage-compensation mechanism. Furthermore, we find partial conservation between the repertoires of mouse and human genes that are associated with sex-biased methylation, an indication that gene function is likely to be an important factor in this phenomenon.

Distinctions between sex and time in patterns of DNA methylation across puberty

BACKGROUND

There are significant sex differences in human physiology and disease; the genomic sources of these differences, however, are not well understood. During puberty, a drastic neuroendocrine shift signals physical changes resulting in robust sex differences in human physiology. Here, we explore how shifting patterns of DNA methylation may inform these pathways of biological plasticity during the pubertal transition. In this study we analyzed DNA methylation (DNAm) in saliva at two time points across the pubertal transition within the same individuals. Our purpose was to compare two domains of DNAm patterns that may inform processes of sexual differentiation 1) sex related sites, which demonstrated differences between males from females and 2) time related sites in which DNAm shifted significantly between timepoints. We further explored the correlated network structure sex and time related DNAm networks and linked these patterns to pubertal stage, assays of salivary testosterone, a reliable diagnostic of free, unbound hormone that is available to act on target tissues, and overlap with androgen response elements.

RESULTS

Sites that differed by biological sex were largely independent of sites that underwent change across puberty. Time-related DNAm sites, but not sex-related sites, formed correlated networks that were associated with pubertal stage. Both time and sex DNAm networks reflected salivary testosterone levels that were enriched for androgen response elements, with sex-related DNAm networks being informative of testosterone levels above and beyond biological sex later in the pubertal transition.

CONCLUSIONS

These results inform our understanding of the distinction between sex- and time-related differences in DNAm during the critical period of puberty and highlight a novel linkage between correlated patterns of sex-related DNAm and levels of salivary testosterone.

The impact of sex differences on genomic research

Sex and gender differences affect all dimensions of human health ranging from the biological basis of disease to therapeutic access, choice and response. Genomics research has long ignored the role of sex differences as potential modulators and the concept is gaining more attention only recently. In the present review we summarize the current knowledge of the impact of sex differences on genomic and epigenomic research, the potential interaction of genomics and gender and the role of these differences in disease etiopathogenesis. Sex differences can emerge from differences in the sex chromosomes themselves, from their interaction with the genome and from the influence of hormones on genomic processes. The impact of these processes on the incidence of autoimmune and oncologic disease is well documented. The growing field of systems biology, which aims at integrating information from different networks of the human body, could also greatly benefit from this approach. In the present review we summarize the current knowledge and provide recommendations for the future performance of sex-sensitive genomics research.

Male predominance in reported Visceral Leishmaniasis cases: Nature or nurture? A comparison of population-based with health facility-reported data

Background

Bangladesh, India, and Nepal aim for the elimination of Visceral Leishmaniasis (VL), a systemic parasitic infectious disease, as a public health problem by 2020. For decades, male patients have comprised the majority of reported VL cases in this region. By comparing this reported VL sex ratio to the one observed in population-based studies conducted in the Indian subcontinent, we tested the working hypothesis that mainly socio-cultural gender differences in healthcare-seeking behavior explain this gender imbalance.

Methodology/Principal findings

We compared the observed sex ratio of male versus female among all VL cases reported by the health system in Nepal and in the two most endemic states in India with that observed in population-based cohort studies in India and Nepal. Also, we assessed male sex as a potential risk factor for seroprevalence at baseline, seroconversion, and VL incidence in the same population-based data. The male/female ratio among VL cases reported by the health systems was 1.40 (95% CI 1.37–1.43). In the population cohort data, the age- and study site-adjusted male to female risk ratio was 1.27 (95% CI 1.08–1.51). Also, males had a 19% higher chance of being seropositive at baseline in the population surveys (RR 1.19; 95% CI 1.11–1.27), while we observed no significant difference in seroconversion rate between both sexes at the DAT cut-off titer defined as the primary endpoint.

Conclusions/Significance

Our population-based data show that male sex is a risk factor for VL, and not only as a socio-cultural determinant. Biological sex-related differences likely play an important role in the pathogenesis of this disease.

Visceral Leishmaniasis (VL) is a parasitic disease that is lethal if not treated timely and mainly affects impoverished populations. Bangladesh, India, and Nepal have targeted the elimination of this disease as a public health problem by 2020. The majority of VL patients attending the health services are male, and this is usually attributed to unequal access to health care for men and women in this sociocultural context. We analyzed two large datasets obtained in population surveillance projects, including regular door-to-door screening for VL, in India and Nepal. Thereby we minimized any potential differences in access to health care between both sexes, as every suspect VL case occurring in the community received a full diagnostic work-up, and the research project facilitated transport to treatment centers if needed. By comparing the observed sex ratio in the health services records with those of the population surveillance records, we aimed to reach meaningful conclusions about the pathway through which male gender exerts its leverage for increasing VL risk: socio-cultural determinants blocking women from accessing care (nurture) or biological factors (nature) making men more vulnerable to VL. Because in the population-based age adjusted-data, male VL cases were significantly more frequent than female, our findings strongly suggest that in the Indian subcontinent, biological differences between men and women play a more critical role in the pathogenesis of VL than previously assumed, and the observed male predominance in VL cases in health services cannot be explained by socio-cultural factors only. Moreover, data show that above the age of 14 years, males are seropositive more often than females and are at higher risk to develop VL disease.

Sex-Based Differences in Human Immunodeficiency Virus Type 1 Reservoir Activity and Residual Immune Activation

Plasma human immunodeficiency virus type 1 (HIV-1) RNA levels in women are lower early in untreated HIV-1 infection compared with those in men, but women have higher T-cell activation and faster disease progression when adjusted for viral load. It is not known whether these sex differences persist during effective antiretroviral therapy (ART), or whether they would be relevant for the evaluation and implementation of HIV-1 cure strategies. We prospectively enrolled a cohort of reproductive-aged women and matched men on suppressive ART and measured markers of HIV-1 persistence, residual virus activity, and immune activation. The frequency of CD4+ T cells harboring HIV-1 DNA was comparable between the sexes, but there was higher cell-associated HIV-1 RNA, higher plasma HIV-1 (single copy assay), and higher T-cell activation and PD-1 expression in men compared with women. These sex-related differences in immune phenotype and HIV-1 persistence on ART have significant implications for the design and measurement of curative interventions.

Sex Differences in the Blood Transcriptome Identify Robust Changes in Immune Cell Proportions with Aging and Influenza Infection

Sex differences in autoimmunity and infection suggest that a better understanding of molecular sex differences will improve the diagnosis and treatment of immune-related disease. We identified 144 differentially expressed genes, referred to as immune sex expression signature (iSEXS), between human males and females using an integrated multi-cohort analysis of blood transcriptome profiles from six discovery cohorts from five continents with 458 healthy individuals. We validated iSEXS in 11 additional cohorts of 524 peripheral blood samples. When we separated iSEXS into genes located on sex chromosomes (XY-iSEXS) or autosomes (autosomal-iSEXS), both modules distinguished males and females. iSEXS reflects sex differences in immune cell proportions, with female-associated genes showing higher expression by CD4+ T cells and male-associated genes showing higher expression by myeloid cells. Autosomal-iSEXS detected an increase in monocytes with age in females, reflected sex-differential immune cell dynamics during influenza infection, and predicted antibody response in males, but not females.

The Netherlands Twin Register: Longitudinal Research Based on Twin and Twin-Family Designs

The Netherlands Twin Register (NTR) is a national register in which twins, multiples and their parents, siblings, spouses and other family members participate. Here we describe the NTR resources that were created from more than 30 years of data collections; the development and maintenance of the newly developed database systems, and the possibilities these resources create for future research. Since the early 1980s, the NTR has enrolled around 120,000 twins and a roughly equal number of their relatives. The majority of twin families have participated in survey studies, and subsamples took part in biomaterial collection (e.g., DNA) and dedicated projects, for example, for neuropsychological, biomarker and behavioral traits. The recruitment into the NTR is all inclusive without any restrictions on enrollment. These resources - the longitudinal phenotyping, the extended pedigree structures and the multigeneration genotyping - allow for future twin-family research that will contribute to gene discovery, causality modeling, and studies of genetic and cultural inheritance.

A characterization of cis- and trans-heritability of RNA-Seq-based gene expression

Insights into individual differences in gene expression and its heritability (h²) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h²_total, composed of cis-heritability (h²_cis, the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h²_res, the residual variance explained by all other genome-wide variants). Mean h²_total was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h² = 0.14, p = 6.15 × 10^-258). Mean h²_cis was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, ρ = 0.76, p < 10^-308) and with estimates from earlier RNA-Seq-based studies. Mean h²_res was 0.20 and correlated with the beta of the corresponding trans-eQTL (ρ = 0.04, p < 1.89 × 10^-3) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 × 10^-15), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h² estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies.

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.

Sexual dimorphism is observed frequently in immune disorders, but the underlying insights are still unclear. Here the authors analyze transcriptome and epigenome changes induced by interferon in various mouse immune cell types, and find only a restricted set of sexual dimorphism genes in innate immunity and macrophages.

Gene dysregulation in peripheral blood of moyamoya disease and comparison with other vascular disorders

Objective

Moyamoya disease (MMD) is a chronic occlusive cerebrovascular disease with unknown etiology, sharing many similar clinical symptoms with other vascular disorders. This study aimed to investigate gene dysregulation in peripheral blood of MMD and compare it with other vascular disorders.

Methods

Transcriptomic profiles of 12 MMD patients and 8 healthy controls were obtained using RNA sequencing. Differentially expressed genes (DEGs) were identified and several were validated by quantitative real-time PCR in independent samples. Biological pathway enrichment analysis of DEGs and deconvolution of leukocyte subsets in peripheral blood were performed. Expression profiles for other vascular diseases were downloaded from public database and consistent DEGs were calculated. Gene set enrichment analysis (GSEA) was conducted to compare gene dysregulation pattern between MMD and other vascular diseases.

Results

A total of 533 DEGs were identified for MMD. Up-regulated genes were mainly involved in extracellular matrix (ECM) organization, whereas down-regulated genes were primarily associated with inflammatory and immune responses. As for cell populations, significantly increased naïve B cells and naïve CD4 cells as well as obviously decreased resting natural killer cells were observed in peripheral blood of MMD patients. GSEA analysis indicated that only up-regulated genes of ischemic stroke and down-regulated genes of coronary artery disease and myocardial infarction were enriched in up-regulated and down-regulated genes of MMD, respectively.

Conclusion

Dysregulated genes in peripheral blood of MMD mainly played key roles in ECM organization, inflammatory and immune responses. This gene dysregulation pattern was specific compared with other vascular diseases. Besides, naïve B cells, naïve CD4 cells and resting natural killer cells were aberrantly disrupted in peripheral blood of MMD patients. These results will help elucidate the complicated pathogenic mechanism of MMD.

SCN11A mRNA levels in female bipolar disorder PBMCs as tentative biomarker for distinct patient sub-phenotypes

Bipolar disorder (BD) is a complex neuropsychiatric disorder characterized by recurrent mania and depression episodes and requiring lifelong treatment with mood stabilizing drugs. Several lines of evidence, including with BD patient iPSC-derived neurons, suggest that neuronal hyperexcitability may underlie the key clinical symptoms of BD. Indeed, higher mRNA levels of SCN11A, coding for the voltage-gated sodium channel Na_V 1.9 implicated in nociception, were detected in iPSC-derived neurons from BD patients, and were normalized by in vitro lithium. Here we studied SCN11A expression in peripheral blood mononuclear cells (PBMCs) from well-phenotyped female BD patients and controls and evaluated their association with several clinical sub-phenotypes. We observed higher mRNA levels of SCN11A in PBMCs from female BD patients with no records of alcohol dependence (p = .0050), no records of psychosis (p = .0097), or no records of suicide attempts (p = .0409). A trend was observed for higher SCN11A expression (FD = 1.91; p = .052) in BD PBMCs compared with controls. Datamining of published postmortem gene expression datasets indicated higher SCN11A expression in dorsolateral prefrontal cortex and orbitofrontal cortex tissues from BD patients compared with controls. Higher phenotype-associated expression levels in PBMC from BD patients were also observed for ID2 (alcohol dependence, suicide attempts) and HDGFRP3 (seasonal BD pattern). Our findings suggest that higher PBMC SCN11A expression levels may be associated with certain behavioral BD sub-phenotypes, including lack of alcohol dependence and psychosis, among BD patients. The Na_V 1.9 voltage-gated sodium channel thus deserves consideration as a tentative phenotype modifier in BD.

Paradoxical role of lipocalin-2 in metabolic disorders and neurological complications

Lipocalin-2 (LCN2), also known as 24p3 and neutrophil gelatinase-associated lipocalin (NGAL), is a 25-kDa secreted protein implicated in various metabolic and inflammatory diseases. Early studies suggest the protective function of LCN2 in which it acts as a bacteriostatic agent that competes with bacteria for iron-bound siderophores. However, both detrimental and beneficial roles of LCN2 have recently been documented in metabolic and neuroinflammatory diseases. Metabolic inflammation, as observed in diabetes and obesity, has been closely associated with the upregulation of LCN2 in blood plasma and several tissues in both humans and rodents, suggesting its pro-diabetic and pro-obesogenic role. On the contrary, other studies imply an anti-diabetic and anti-obesogenic role of LCN2 whereby a deficiency in the Lcn2 gene results in the impairment of insulin sensitivity and enhances the high-fat-diet-induced expansion of fat. A similar dual role of LCN2 has also been reported in various animal models for neurological disorders. In the midst of these mixed findings, there is no experimental evidence to explain why LCN2 shows such a contrasting role in the various studies. This debate needs to be resolved (or reconciled) and an integrated view on the topic is desirable. Herein, we attempt to address this issue by reviewing the recent findings on LCN2 in metabolic disorders and assess the potential cellular or molecular mechanisms underlying the dual role of LCN2. We further discuss the possibilities and challenges of targeting LCN2 as a potential therapeutic strategy for metabolic disorders and neurological complications.