Adaptive variation regulates the expression of the human SGK1 gene in response to stress

Genotype × environment interactions in gene regulation and complex traits

Genotype × environment interactions (GxE) have long been recognized as a key mechanism underlying human phenotypic variation. Technological developments over the past 15 years have dramatically expanded our appreciation of the role of GxE in both gene regulation and complex traits. The richness and complexity of these datasets also required parallel efforts to develop robust and sensitive statistical and computational approaches. Although our understanding of the genetic architecture of molecular and complex traits has been maturing, a large proportion of complex trait heritability remains unexplained. Furthermore, there are increasing efforts to characterize the effect of environmental exposure on human health. We therefore review GxE in human gene regulation and complex traits, advocating for a comprehensive approach that jointly considers genetic and environmental factors in human health and disease.

Dexmedetomidine improves the acute stress reactivity of male rat through interventions of serum- and glucocorticoid-inducible kinase 1 and nNOS in the bed nucleus of the stria terminalis

Moderate acute stress responses are beneficial for adaptation and maintenance of homeostasis. Exposure of male rat to stress induces effects in the bed nucleus of the stria terminalis (BNST), for it can be activated by the same stimuli that induce activation of the hypothalamic-pituitary-adrenal axis. However, the underlying mechanism of the BNST on male stress reactivity remains unclear. In this study, we explored whether systematic administration of dexmedetomidine (DEXM) altered the acute stress reactivity through its effect on the BNST. Male Sprague-Dawley rats in the stress (STRE) group, DEXM group, and the DEXM + GSK-650394 (GSK, an antagonist of serum- and glucocorticoid-inducible kinase 1 (SGK1)) group, except those in the vehicle (VEH) group, underwent 1-h restraint plus water-immersion (RPWI) exposure. All the rats proceeded the open field test (OFT) 24 h before RPWI and 1 h after RPWI. After the second OFT, the rats received VEH, DEXM (75 μg/kg i.p.), or were pretreated with GSK (2 μM i.p.) 0.5 h ahead of DEXM respectively. The third OFT was conducted 6 h after drug administration and then the rats were sacrificed. The rats that experienced RPWI showed dramatically elevated serum corticosterone (CORT), multiplied neuronal nitric oxide synthase (nNOS) and SGK1 in the BNST, and terrible OFT behavior. We discovered when the nNOS and SGK1 were decreased in the rat BNST through DEXM treatment, the serum CORT was reduced and the OFT manifestation was ameliorated, whereas these were restrained by GSK application. Our results reveal that modest interventions to SGK1 and nNOS in the BNST improve the male rat reactivity to acute stress, and DEXM was one modulator of these effects.

Emerging Role of Serum Glucocorticoid-Regulated Kinase 1 in Pathological Pain

Currently, the management of acute and chronic pain in clinical practice remains unsatisfactory due to the existence of limited effective treatments, and novel therapeutic strategies for pathological pain are urgently needed. In the past few decades, the role of serum and glucocorticoid-inducible kinase 1 (SGK1) in the development of pain and diurnal rhythms has been implicated in numerous studies. The expression levels of SGK1 mRNA and protein were found to be elevated in the spinal cord and brain in various pathological pain models. Blocking SGK1 significantly attenuated pain-like responses and the development of pathological pain. These studies provide strong evidence that SGK1 plays a role in the development of various types of pathological pain and that targeting SGK1 may be a novel therapeutic strategy for pain management. In this review article, we provide evidence from animal models for the potential role of SGK1 in the regulation of pathological pain caused by inflammation, nerve injury, psychiatric disorders, and chronic opioid exposure.

The Emerging Role of SGK1 (Serum- and Glucocorticoid-Regulated Kinase 1) in Major Depressive Disorder: Hypothesis and Mechanisms

Major depressive disorder (MDD) is a heterogeneous psychiatric disease characterized by persistent low mood, diminished interests, and impaired cognitive and social functions. The multifactorial etiology of MDD is still largely unknown because of the complex genetic and environmental interactions involved. Therefore, no established mechanism can explain all the aspects of the disease. In this light, an extensive research about the pathophysiology of MDD has been carried out. Several pathogenic hypotheses, such as monoamines deficiency and neurobiological alterations in the stress-responsive system, including the hypothalamic-pituitary-adrenal (HPA) axis and the immune system, have been proposed for MDD. Over time, remarkable studies, mainly on preclinical rodent models, linked the serum- and glucocorticoid-regulated kinase 1 (SGK1) to the main features of MDD. SGK1 is a serine/threonine kinase belonging to the AGK Kinase family. SGK1 is ubiquitously expressed, which plays a pivotal role in the hormonal regulation of several ion channels, carriers, pumps, and transcription factors or regulators. SGK1 expression is modulated by cell stress and hormones, including gluco- and mineralocorticoids. Compelling evidence suggests that increased SGK1 expression or function is related to the pathogenic stress hypothesis of major depression. Therefore, the first part of the present review highlights the putative role of SGK1 as a critical mediator in the dysregulation of the HPA axis, observed under chronic stress conditions, and its controversial role in the neuroinflammation as well. The second part depicts the negative regulation exerted by SGK1 in the expression of both the brain-derived neurotrophic factor (BDNF) and the vascular endothelial growth factor (VEGF), resulting in an anti-neurogenic activity. Finally, the review focuses on the antidepressant-like effects of anti-oxidative nutraceuticals in several preclinical model of depression, resulting from the restoration of the physiological expression and/or activity of SGK1, which leads to an increase in neurogenesis. In summary, the purpose of this review is a systematic analysis of literature depicting SGK1 as molecular junction of the complex mechanisms underlying the MDD in an effort to suggest the kinase as a potential biomarker and strategic target in modern molecular antidepressant therapy.

Leonurine promotes neurite outgrowth and neurotrophic activity by modulating the GR/SGK1 signaling pathway in cultured PC12 cells

Depression is a common psychiatric disorder that affects almost 10% of children and adolescents worldwide. Numerous synthetic chemical antidepressants used to treat depression have adverse side effects. Therefore, new therapeutic approaches for depression treatment are urgently needed. Leonurus cardiaca has recently been shown to be effective for the treatment of nervous system diseases such as depression, but its mechanism is not clear. In this study, we aimed to reveal the mechanism underlying leonurine’s antidepressant activity. Leonurine was used to treat corticosterone-induced PC12 cells to examine its effect on neurite outgrowth and neurotrophic factors after treatment with the inhibitor of glucocorticoid receptor (GR) and serum-inducible and glucocorticoid-inducible kinase 1 (SGK1). Methyl thiazolyl tetrazolium assays were used to evaluate the viability of cells. High content analysis was used to detect cell area, total neurite length, maximum neurite length, and expression of GR, SGK1, brain-derived neurotrophic factor (BDNF), neurotrophic factor-3 (NT-3), and B-cell lymphoma-2 (BCL-2). The results showed that leonurine increased cell viability in a concentration-dependent manner, with the maximal prosurvival effect at 60 μM. Leonurine increased cell area, total neurite length, and maximum neurite length of corticosterone-induced PC12 cells, increased the expression of GR, BDNF, NT-3, and BCL-2, and decreased the expression of SGK1. After treatment with GR inhibitor RU486, the expressions of GR, BDNF, NT-3, and BCL-2 were significantly decreased and SGK1 was increased. In contrast, treatment with GSK650394 had the opposite effect of RU486. Our data indicate that leonurine promotes neurite outgrowth and neurotrophic activity in cultured PC12 cells, and its potential mechanism may involve the GR/SGK1 signaling pathway.

Implication of Spiritual Network Support System in Epigenomic Modulation and Health Trajectory

With challenges in understanding the multifactorial etiologies of disease and individual treatment effect heterogeneities over the past four decades, much has been acquired on how physical, chemical and social environments affect human health, predisposing certain subpopulations to adverse health outcomes, especially the socio-environmentally disadvantaged (SED). Current translational data on gene and adverse environment interaction have revealed how adverse gene-environment interaction, termed aberrant epigenomic modulation, translates into impaired gene expression via messenger ribonucleic acid (mRNA) dysregulation, reflecting abnormal protein synthesis and hence dysfunctional cellular differentiation and maturation. The environmental influence on gene expression observed in most literature includes physical, chemical, physicochemical and recently social environment. However, data are limited on spiritual or religious environment network support systems, which reflect human psychosocial conditions and gene interaction. With this limited information, we aimed to examine the available data on spiritual activities characterized by prayers and meditation for a possible explanation of the nexus between the spiritual network support system (SNSS) as a component of psychosocial conditions, implicated in social signal transduction, and the gene expression correlate. With the intent to incorporate SNSS in human psychosocial conditions, we assessed the available data on bereavement, loss of spouse, loneliness, social isolation, low socio-economic status (SES), chronic stress, low social status, social adversity (SA) and early life stress (ELS), as surrogates for spiritual support network connectome. Adverse human psychosocial conditions have the tendency for impaired gene expression through an up-regulated conserved transcriptional response to adversity (CTRA) gene expression via social signal transduction, involving the sympathetic nervous system (SNS), beta-adrenergic receptors, the hypothalamus-pituitary-adrenal (HPA) axis and the glucocorticoid response. This review specifically explored CTRA gene expression and the nuclear receptor subfamily 3 group C member 1 (NR3C1) gene, a glucocorticoid receptor gene, in response to stress and the impaired negative feedback, given allostatic overload as a result of prolonged and sustained stress and social isolation as well as the implied social interaction associated with religiosity. While more remains to be investigated on psychosocial and immune cell response and gene expression, current data on human models do implicate appropriate gene expression via the CTRA and NR3C1 gene in the SNSS as observed in meditation, yoga and thai-chi, implicated in malignant neoplasm remission. However, prospective epigenomic studies in this context are required in the disease causal pathway, prognosis and survival, as well as cautious optimism in the application of these findings in clinical and public health settings, due to unmeasured and potential confoundings implicated in these correlations.

Association of SGK1 Polymorphisms With Susceptibility to Coronary Heart Disease in Chinese Han Patients With Comorbid Depression

There is a strong link between heart disease and depression, both of which are closely related to lifetime stress exposure. Serum/glucocorticoid-regulated kinase 1 (SGK1) is a stress-responsive gene with a pivotal role in both the heart and brain. To determine the role of SGK1 polymorphisms (rs2758151, rs1743963, rs9493857, rs1763509, rs9376026, and rs9389154) in susceptibility to comorbid coronary heart disease (CHD) and depression, we conducted a hospital-based case–control study involving 257 CHD cases (including 69 cases with depression and 188 cases without depression) and 107 controls in a Chinese Han population. Six single-nucleotide polymorphisms (SNPs) in the SGK1 gene were successfully genotyped by polymerase chain reaction–ligase detection reaction (PCR-LDR) assay. Our results showed no significant differences in SGK1 genetic polymorphisms between CHD patients and controls, whereas significant associations were observed between SGK1 SNPs (rs1743963 and rs1763509) and the development of depression in CHD patients (P = 0.018 by genotype, P = 0.032 by allele; P = 0.017 by genotype, P = 0.003 by allele, respectively). However, none of these associations remained significant after Bonferroni correction (P = 0.054 for rs1743963; P = 0.051 for rs1763509). Interestingly, both the GG genotype of SGK1 rs1743963 and AA genotype of SGK1 rs1763509 were associated with a higher risk of depression in CHD patients; for rs1763509, the Patient Health Questionnaire-9 (PHQ-9) scores in the carriers of the risk genotype for comorbid depression, AA, were significantly higher than in GG and AG carriers (P = 0.008). Notably, haplotype analysis indicated that haplotype GGA significantly increased the risk of depression in CHD patients (P = 0.011, odds ratio (OR) = 1.717, 95% confidence interval (CI) = 1.132–2.605), whereas haplotype AAG may be a protective factor for CHD patients with comorbid depression (P = 0.038, OR = 0.546, 95% CI = 0.307–0.972). It should be noted that only the significance of haplotype GGA survived after Bonferroni adjustment (P = 0.044) and that no significant differences were found for other SGK1 SNPs (rs2758151, rs9493857, rs9376026, and rs9389154) between CHD patients with and without depression. These findings, for the first time, elucidate the important role of SGK1 variants in the comorbidity of CHD and depression.

Region- and time-dependent gene regulation in the amygdala and anterior cingulate cortex of a PTSD-like mouse model

Posttraumatic stress disorder is developed by exposure to a threatening and/or a horrifying event and characterized by the presence of anxiety, hyperarousal, avoidance, and sleep abnormality for a prolonged period of time. To elucidate the potential molecular mechanisms, we constructed a mouse model by electric foot shock followed by situational reminders and performed transcriptome analysis in brain tissues. The stressed mice acquired anxiety-like behavior after 2 weeks and exaggerated startle response after 4 weeks. Avoidance latency and freezing behavior were sustained up to 5 weeks post stress and abnormal static behavior was observed during the sleep period. RNA sequencing was performed in two of the emotional regulatory regions, anterior cingulate cortex and amygdala, at 2 and 5 weeks post stress. More than 1000 differentially expressed genes were identified at 2 weeks in both regions. The number of the regulated genes remained constant in amygdala at 5 weeks post stress, whereas those in anterior cingulate cortex were plummeted. Although synaptic remodeling and endocrine system were the most enriched signaling pathways in both anterior cingulate cortex and amygdala, the individual gene expression profile was regulated in a region- and time-dependent manner. In addition, several genes associated with PTSD involved in Hypothalamic-Pituitary-Adrenal axis were differentially regulated. These findings suggested that global gene expression profile was dynamically regulated in accordance with the disease development stage, and therefore targeting the distinct signaling molecules in different region and development stage might be critical for effective treatment to PTSD.

Social status alters chromatin accessibility and the gene regulatory response to glucocorticoid stimulation in rhesus macaques

Low social status is an important predictor of disease susceptibility and mortality risk in humans and other social mammals. These effects are thought to stem in part from dysregulation of the glucocorticoid (GC)-mediated stress response. However, the molecular mechanisms that connect low social status and GC dysregulation to downstream health outcomes remain elusive. Here, we used an in vitro GC challenge to investigate the consequences of experimentally manipulated social status (i.e., dominance rank) for immune cell gene regulation in female rhesus macaques, using paired control and GC-treated peripheral blood mononuclear cell samples. We show that social status not only influences immune cell gene expression but also chromatin accessibility at hundreds of regions in the genome. Social status effects on gene expression were less pronounced following GC treatment than under control conditions. In contrast, social status effects on chromatin accessibility were stable across conditions, resulting in an attenuated relationship between social status, chromatin accessibility, and gene expression after GC exposure. Regions that were more accessible in high-status animals and regions that become more accessible following GC treatment were enriched for a highly concordant set of transcription factor binding motifs, including motifs for the GC receptor cofactor AP-1. Together, our findings support the hypothesis that social status alters the dynamics of GC-mediated gene regulation and identify chromatin accessibility as a mechanism involved in social stress-driven GC resistance. More broadly, they emphasize the context-dependent nature of social status effects on gene regulation and implicate epigenetic remodeling of chromatin accessibility as a contributing factor.

Maternal mRNA input of growth and stress-response-related genes in cichlids in relation to egg size and trophic specialization

BACKGROUND

Egg size represents an important form of maternal effect determined by a complex interplay of long-term adaptation and short-term plasticity balancing egg size with brood size. Haplochromine cichlids are maternal mouthbrooders showing differential parental investment in different species, manifested in great variation in egg size, brood size and duration of maternal care. Little is known about maternally determined molecular characters of eggs in fishes and their relation to egg size and trophic specialization. Here we investigate maternal mRNA inputs of selected growth- and stress-related genes in eggs of mouthbrooding cichlid fishes adapted to different trophic niches from Lake Tanganyika, Lake Malawi, Lake Victoria and compare them to their riverine allies.

RESULTS

We first identified two reference genes, atf7ip and mid1ip1, to be suitable for cross-species quantification of mRNA abundance via qRT-PCR in the cichlid eggs. Using these reference genes, we found substantial variation in maternal mRNA input for a set of candidate genes related to growth and stress response across species and lakes. We observed negative correlation of mRNA abundance between two of growth hormone receptor paralogs (ghr1 and ghr2) across all haplochromine cichlid species which also differentiate the species in the two younger lakes, Malawi and Lake Victoria, from those in Lake Tanganyika and ancestral riverine species. Furthermore, we found correlations between egg size and maternal mRNA abundance of two growth-related genes igf2 and ghr2 across the haplochromine cichlids as well as distinct clustering of the species based on their trophic specialization using maternal mRNA abundance of five genes (ghr1, ghr2, igf2, gr and sgk1).

CONCLUSIONS

These findings indicate that variations in egg size in closely related cichlid species can be linked to differences in maternal RNA deposition of key growth-related genes. In addition, the cichlid species with contrasting trophic specialization deposit different levels of maternal mRNAs in their eggs for particular growth-related genes; however, it is unclear whether such differences contribute to differential morphogenesis at later stages of development. Our results provide first insights into this aspect of gene activation, as a basis for future studies targeting their role during ecomorphological specialization and adaptive radiation.

Vitamin D Regulation of the Uridine Phosphorylase 1 Gene and Uridine-Induced DNA Damage in Colon in African Americans and European Americans

BACKGROUND & AIMS

African Americans have the greatest colorectal cancer (CRC) burden in the United States; interethnic differences in protective effects of vitamin D might contribute to disparities. 1α,25(OH)₂D₃ vitamin D (the active form of vitamin D) induces transcription of the uridine phosphorylase gene (UPP1) in colon tissues of European Americans but to a lesser extent in colon tissues of African Americans. UPP1-knockout mice have increased intestinal concentrations of uridine and Deoxyuridine triphosphate (dUTP), have increased uridine-induced DNA damage, and develop colon tumors. We studied 1α,25(OH)₂D₃ regulation of UPP1 and uridine-induced DNA damage in the colon and differences in these processes between African and European Americans.

METHODS

We quantified expression and activity of UPP1 in response to 1α,25(OH)₂D₃ in young adult mouse colonic cells, human CRC cells (LS174T), and organoids (derived from rectosigmoid biopsy samples of healthy individuals undergoing colonoscopies) using quantitative polymerase chain reaction, immunoblot, and immunocytochemistry assays. Binding of the vitamin D receptor to UPP1 was tested by chromatin immunoprecipitation. Uridine-induced DNA damage was measured by fragment-length analysis in repair enzyme assays. Allele-specific 1α,25(OH)₂D₃ responses were tested using luciferase assays.

RESULTS

Vitamin D increased levels of UPP1 mRNA, protein, and enzymatic activity and increased vitamin D receptor binding to the UPP1 promoter in young adult mouse colonic cells, LS174T cells, and organoids. 1α,25(OH)₂D₃ significantly reduced levels of uridine and uridine-induced DNA damage in these cells, which required UPP1 expression. Organoids derived from colon tissues of African Americans expressed lower levels of UPP1 after exposure to 1α,25(OH)₂D₃ and had increased uridine-induced DNA damage compared with organoids derived from tissues of European Americans. Luciferase assays with the T allele of single nucleotide polymorphism rs28605337 near UPP1, which is found more frequently in African Americans than European Americans, expressed lower levels of UPP1 after exposure to 1α,25(OH)₂D₃ than assays without this variant.

CONCLUSIONS

We found vitamin D to increase expression of UPP1, leading to reduce uridine-induced DNA damage, in colon cells and organoids. A polymorphism in UPP1 found more frequently in African Americans than European Americans reduced UPP1 expression upon cell exposure to 1α,25(OH)₂D₃. Differences in expression of UPP1 in response to vitamin D could contribute to the increased risk of CRC in African Americans.

Dietary factors in rheumatic autoimmune diseases: a recipe for therapy?

Today, we are facing a new era of digitization in the health care system, and with increased access to health care information has come a growing demand for safe, cost-effective and easy to administer therapies. Dietary habits have a crucial influence on human health, affecting an individual's risk for hypertension, heart disease and stroke, as well as influencing the risk of developing of cancer. Moreover, an individual's lifestyle choices can greatly influence the progression and manifestation of chronic autoimmune rheumatic diseases. In light of these effects, it makes sense that the search for additional therapies to attenuate such diseases would include investigations into lifestyle modifications. When considering the complex web of factors that influence autoimmunity, it is not surprising to find that several dietary elements are involved in disease progression or prevention. In this Review, several common nutritional components of the human diet are presented, and the evidence for their effects on rheumatic diseases is discussed.

SGK-1 protects kidney cells against apoptosis induced by ceramide and TNF-α

Ceramide regulates several different cellular responses including mechanisms leading to apoptosis. Serum- and glucocorticoid-inducible protein kinase (SGK)-1 is a serine threonine kinase, which activates survival pathways in response to stress stimuli. Recently, we demonstrated an anti-apoptotic role of SGK-1 in human umbilical endothelial cells treated with high glucose. In the present study, since ceramide induces apoptosis by multiple mechanisms in diabetes and its complication such as nephropathy, we aimed to investigate whether SGK-1 may protect even against apoptosis induced by ceramide in kidney cells. Human embryonic kidney (HEK)-293 cells stable transfected with SGK-1 wild type (SGK-1wt) and its dominant negative gene (SGK-1dn) have been used in this study. Apoptotic stimuli were induced by C₂-ceramide and TNF-α to increase endogenous synthesis of ceramide. Upon activation with these stimuli, SGK-1wt transfected cells have a statistically significant reduction of apoptosis compared with SGK-1dn cells (P<0.001). This protection was dependent on activation of caspase-3 and Poly-ADP-ribose-polymerase-1 (PARP-1) cleavage. SGK-1 and AKT-1 two highly homologous kinases differently reacted to ceramide treatment, since SGK-1 increases in response to apoptotic stimulus while AKT-1 decreases. This enhancement of SGK-1 was dependent on p38-mitogen-activated-protein kinases (p38MAPK), cyclic-adenosine-monophosphate/protein kinase A (cAMP/PKA) and phosphoinositide-3-kinase (PI3K) pathways. Especially, by using selective LY294002 inhibitor, we demonstrated that the most involved pathway in the SGK-1 mediated process of protection was PI3K. Treatment with inhibitor of SGK-1 (GSK650394) significantly enhanced TNF-α-dependent apoptosis in HEK-293 cells overexpressing SGK-1wt. Caspase-3, -8 and -9 selective inhibitors confirmed that SGK-1 reduced the activation of caspase-dependent apoptosis, probably by both intrinsic and extrinsic pathways. In conclusion, we demonstrated that in kidney cells, overexpression of SGK-1 is protective against ceramide-induced apoptosis and the role of SGK-1 can be potentially explored as a therapeutic target in conditions like diabetes, where ceramide levels are increased.

Human social genomics

A growing literature in human social genomics has begun to analyze how everyday life circumstances influence human gene expression. Social-environmental conditions such as urbanity, low socioeconomic status, social isolation, social threat, and low or unstable social status have been found to associate with differential expression of hundreds of gene transcripts in leukocytes and diseased tissues such as metastatic cancers. In leukocytes, diverse types of social adversity evoke a common conserved transcriptional response to adversity (CTRA) characterized by increased expression of proinflammatory genes and decreased expression of genes involved in innate antiviral responses and antibody synthesis. Mechanistic analyses have mapped the neural "social signal transduction" pathways that stimulate CTRA gene expression in response to social threat and may contribute to social gradients in health. Research has also begun to analyze the functional genomics of optimal health and thriving. Two emerging opportunities now stand to revolutionize our understanding of the everyday life of the human genome: network genomics analyses examining how systems-level capabilities emerge from groups of individual socially sensitive genomes and near-real-time transcriptional biofeedback to empirically optimize individual well-being in the context of the unique genetic, geographic, historical, developmental, and social contexts that jointly shape the transcriptional realization of our innate human genomic potential for thriving.

Genetic analysis of reproductive traits and antibody response in a PRRS outbreak herd

Porcine reproductive and respiratory syndrome (PRRS) is the most economically significant disease impacting pig production in North America, Europe, and Asia, causing reproductive losses such as increased rates of stillbirth and mummified piglets. The objective of this study was to explore the genetic basis of host response to the PRRS virus (PRRSV) in a commercial multiplier sow herd before and after a PRRS outbreak, using antibody response and reproductive traits. Reproductive data comprising number born alive (NBA), number alive at 24 h (NA24), number stillborn (NSB), number born mummified (NBM), proportion born dead (PBD), number born dead (NBD), number weaned (NW), and number of mortalities through weaning (MW) of 5,227 litters from 1,967 purebred Landrace sows were used along with a pedigree comprising 2,995 pigs. The PRRS outbreak date was estimated from rolling averages of farrowing traits and was used to split the data into a pre-PRRS phase and a PRRS phase. All 641 sows in the herd during the outbreak were blood sampled 46 d after the estimated outbreak date and were tested for anti-PRRSV IgG using ELISA (sample-to-positive [S/P] ratio). Genetic parameters of traits were estimated separately for the pre-PRRS and PRRS phase data sets. Sows were genotyped using the PorcineSNP60 BeadChip, and genome-wide association studies (GWAS) were performed using method Bayes B. Heritability estimates for reproductive traits ranged from 0.01 (NBM) to 0.12 (NSB) and from 0.01 (MW) to 0.12 (NBD) for the pre-PRRS and PRRS phases, respectively. S/P ratio had heritability (0.45) and strong genetic correlations with most traits, ranging from -0.72 (NBM) to 0.73 (NBA). In the pre-PRRS phase, regions associated with NSB and PBD explained 1.6% and 3% of the genetic variance, respectively. In the PRRS phase, regions associated with NBD, NSB, and S/P ratio explained 0.8%, 11%, and 50.6% of the genetic variance, respectively. For S/P ratio, 2 regions on SSC 7 (SSC7) separated by 100 Mb explained 40% of the genetic variation, including a region encompassing the major histocompatibility complex, which explained 25% of the genetic variance. These results indicate a significant genomic component associated with PRRSV antibody response and NSB in this data set. Also, the high heritability and genetic correlation estimates for S/P ratio during the PRRS phase suggest that S/P ratio could be used as an indicator of the impact of PRRS on reproductive traits.

Human expression QTLs are enriched in signals of environmental adaptation

Expression quantitative trait loci (eQTLs) have been found to be enriched in trait-associated single-nucleotide polymorphisms (SNPs). However, whether eQTLs are adaptive to different environmental factors and its relative evolutionary significance compared with nonsynonymous SNPs (NS SNPs) are still elusive. Compiling environmental correlation data from three studies for more than 500,000 SNPs and 42 environmental factors, including climate, subsistence, pathogens, and dietary patterns, we performed a systematic examination of the adaptive patterns of eQTLs to local environment. Compared with intergenic SNPs, eQTLs are significantly enriched in the lower tail of a transformed rank statistic in the environmental correlation analysis, indicating possible adaptation of eQTLs to the majority of 42 environmental factors. The mean enrichment of eQTLs across 42 environmental factors is as great as, if not greater than, that of NS SNPs. The enrichment of eQTLs, although significant across all levels of recombination rate, is inversely correlated with recombination rate, suggesting the presence of selective sweep or background selection. Further pathway enrichment analysis identified a number of pathways with possible environmental adaption from eQTLs. These pathways are mostly related with immune function and metabolism. Our results indicate that eQTLs might have played an important role in recent and ongoing human adaptation and are of special importance for some environmental factors and biological pathways.

Quantitative trait loci for magnitude of the plasma cortisol response to confinement in rainbow trout

Better understanding of the mechanisms underlying interindividual variation in stress responses and their links with production traits is a key issue for sustainable animal breeding. In this study, we searched for quantitative trait loci (QTL) controlling the magnitude of the plasma cortisol stress response and compared them to body size traits in five F2 full-sib families issued from two rainbow trout lines divergently selected for high or low post-confinement plasma cortisol level. Approximately 1000 F2 individuals were individually tagged and exposed to two successive acute confinement challenges (1 month interval). Post-stress plasma cortisol concentrations were determined for each fish. A medium density genome scan was carried out (268 markers, overall marker spacing less than 10 cM). QTL detection was performed using qtlmap software, based on an interval mapping method (http://www.inra.fr/qtlmap). Overall, QTL of medium individual effects on cortisol responsiveness (<10% of phenotypic variance) were detected on 18 chromosomes, strongly supporting the hypothesis that control of the trait is polygenic. Although a core array of QTL controlled cortisol concentrations at both challenges, several QTL seemed challenge specific, suggesting that responses to the first and to a subsequent exposure to the confinement stressor are distinct traits sharing only part of their genetic control. Chromosomal location of the steroidogenic acute regulatory protein (STAR) makes it a good potential candidate gene for one of the QTL. Finally, comparison of body size traits QTL (weight, length and body conformation) with cortisol-associated QTL did not support evidence for negative genetic relationships between the two types of traits.

Role for the kinase SGK1 in stress, depression, and glucocorticoid effects on hippocampal neurogenesis

Stress and glucocorticoid hormones regulate hippocampal neurogenesis, but the molecular mechanisms mediating these effects are poorly understood. Here we identify the glucocorticoid receptor (GR) target gene, serum- and glucocorticoid-inducible kinase 1 (SGK1), as one such mechanism. Using a human hippocampal progenitor cell line, we found that a small molecule inhibitor for SGK1, GSK650394, counteracted the cortisol-induced reduction in neurogenesis. Moreover, gene expression and pathway analysis showed that inhibition of the neurogenic Hedgehog pathway by cortisol was SGK1-dependent. SGK1 also potentiated and maintained GR activation in the presence of cortisol, and even after cortisol withdrawal, by increasing GR phosphorylation and GR nuclear translocation. Experiments combining the inhibitor for SGK1, GSK650394, with the GR antagonist, RU486, demonstrated that SGK1 was involved in the cortisol-induced reduction in progenitor proliferation both downstream of GR, by regulating relevant target genes, and upstream of GR, by increasing GR function. Corroborating the relevance of these findings in clinical and rodent settings, we also observed a significant increase of SGK1 mRNA in peripheral blood of drug-free depressed patients, as well as in the hippocampus of rats subjected to either unpredictable chronic mild stress or prenatal stress. Our findings identify SGK1 as a mediator for the effects of cortisol on neurogenesis and GR function, with particular relevance to stress and depression.

Predicting PTSD: pre-existing vulnerabilities in glucocorticoid-signaling and implications for preventive interventions

Posttraumatic stress disorder (PTSD) is an anxiety disorder that may develop in response to a traumatic event. Approximately 10% of trauma-exposed individuals subsequently develop PTSD. It is hypothesized that the development of PTSD is associated with biological vulnerability factors, which are already present prior to the onset of symptoms. In this review we present an overview of currently identified vulnerability factors in the glucocorticoid (GC) signaling pathway for the development of PTSD. In addition, the implications of the identified vulnerability factors for potential preventive intervention strategies, including glucocorticoid receptor (GR) agonists and oxytocin, are discussed. Summarized, the findings of these studies indicate that individuals vulnerable for development of PTSD have dysregulations on various levels of the GC-signaling cascade: i.e. low levels of circulating levels of cortisol shortly after trauma, high GR number in peripheral blood mononuclear cells (PBMCs), high GILZ mRNA expression and low FKBP5 expression in PBMCs prior to trauma, and high sensitivity of T-cells for regulation by GCs prior to trauma. Furthermore, single nucleotide polymorphisms in the GR and FKBP5 genes have been found to be associated with increased risk for PTSD. Collectively, the identified vulnerability factors tentatively suggest that the development of PTSD may be preceded by a high sensitivity of various cells for regulation by GCs. The identification of these vulnerability factors may ultimately aid selective targeting of preventive interventions towards individuals at risk for PTSD. In addition, the identification of these vulnerability factors may eventually result in new preventive pharmacological strategies for PTSD.

Genetic, functional and molecular features of glucocorticoid receptor binding

Glucocorticoids (GCs) are key mediators of stress response and are widely used as pharmacological agents to treat immune diseases, such as asthma and inflammatory bowel disease, and certain types of cancer. GCs act mainly by activating the GC receptor (GR), which interacts with other transcription factors to regulate gene expression. Here, we combined different functional genomics approaches to gain molecular insights into the mechanisms of action of GC. By profiling the transcriptional response to GC over time in 4 Yoruba (YRI) and 4 Tuscans (TSI) lymphoblastoid cell lines (LCLs), we suggest that the transcriptional response to GC is variable not only in time, but also in direction (positive or negative) depending on the presence of specific interacting transcription factors. Accordingly, when we performed ChIP-seq for GR and NF-κB in two YRI LCLs treated with GC or with vehicle control, we observed that features of GR binding sites differ for up- and down-regulated genes. Finally, we show that eQTLs that affect expression patterns only in the presence of GC are 1.9-fold more likely to occur in GR binding sites, compared to eQTLs that affect expression only in its absence. Our results indicate that genetic variation at GR and interacting transcription factors binding sites influences variability in gene expression, and attest to the power of combining different functional genomic approaches.

Therapeutic potential of serum and glucocorticoid inducible kinase inhibition

INTRODUCTION

Expression of serum-and-glucocorticoid-inducible kinase-1 (SGK1) is low in most cells, but dramatically increases under certain pathophysiological conditions, such as glucocorticoid or mineralocorticoid excess, inflammation with TGFβ release, hyperglycemia, cell shrinkage and ischemia. SGK1 is activated by insulin and growth factors via phosphatidylinositide-3-kinase, 3-phosphoinositide-dependent kinase and mammalian target of rapamycin. SGK1 sensitive functions include activation of ion channels (including epithelial Na(+) channel ENaC, voltage gated Na(+) channel SCN5A transient receptor potential channels TRPV4 - 6, Ca(2+) release activated Ca(2+) channel Orai1/STIM1, renal outer medullary K(+) channel ROMK, voltage gated K(+) channels KCNE1/KCNQ1, kainate receptor GluR6, cystic fibrosis transmembrane regulator CFTR), carriers (including Na(+),Cl(-) symport NCC, Na(+),K(+),2Cl(-) symport NKCC, Na(+)/H(+) exchangers NHE1 and NHE3, Na(+), glucose symport SGLT1, several amino acid transporters), and Na(+)/K(+)-ATPase. SGK1 regulates several enzymes (e.g., glycogen synthase kinase-3, ubiquitin-ligase Nedd4-2) and transcription factors (e.g., forkhead transcription factor 3a, β-catenin, nuclear factor kappa B).

AREAS COVERED

The phenotype of SGK1 knockout mice is mild and SGK1 is apparently dispensible for basic functions. Excessive SGK1 expression and activity, however, contributes to the pathophysiology of several disorders, including hypertension, obesity, diabetes, thrombosis, stroke, fibrosing disease, infertility and tumor growth. A SGK1 gene variant (prevalence ∼ 3 - 5% in Caucasians and ∼ 10% in Africans) is associated with hypertension, stroke, obesity and type 2 diabetes. SGK1 inhibitors have been developed and shown to reduce blood pressure of hyperinsulinemic mice and to counteract tumor cell survival.

EXPERT OPINION

Targeting SGK1 may be a therapeutic option in several clinical conditions, including metabolic syndrome and tumor growth.

Allelic imbalance assays to quantify allele-specific gene expression and transcription factor binding

A growing number of noncoding variants are found to influence the susceptibility to common diseases and interindividual variation in drug response. However, the mechanisms by which noncoding variation affects cellular and clinical phenotypes remain to be elucidated. Allele-specific assays allow testing directly the differential properties of the alleles at a regulatory variant, which are detected as an allelic imbalance. Two widely used allelic imbalance assays target cDNA and DNA from chromatin immunoprecipitation (ChIP) experiments, and therefore revealing allele-specific gene expression and transcription factor binding, respectively. The throughput of allelic imbalance assays ranges from single variant to the genome scale, which are made possible by the recent advances in genotyping and sequencing technologies (e.g., genome-wide quantitative cDNA genotyping, ChIP-seq).

Mapping gene-environment interactions at regulatory polymorphisms: insights into mechanisms of phenotypic variation

Genetic effects on gene regulation make a substantial contribution to phenotypic diversity, yet their mechanisms remain elusive. Here, we discuss the potential insights to be gained from mapping gene-environment interactions at regulatory polymorphisms (i.e., genetic variation that affects gene expression under specific environmental conditions). We highlight a novel statistical method to identify specific patterns of gene-environment interaction at these regulatory polymorphisms. Reviewing its application to a study that mapped gene expression in the presence and absence of glucocorticoids, we discuss the mechanistic insights that this approach provides.

Glucocorticoid receptor pathway components predict posttraumatic stress disorder symptom development: a prospective study

BACKGROUND

Biological correlates of posttraumatic stress disorder (PTSD) have mostly been studied using cross-sectional or posttrauma prospective designs. Therefore, it remains largely unknown whether previously observed biological correlates of PTSD precede trauma exposure. We investigated whether glucocorticoid receptor (GR) pathway components assessed in leukocytes before military deployment represent preexisting vulnerability factors for development of PTSD symptoms.

METHODS

Four hundred forty-eight male soldiers were assessed before and 6 months after deployment to a combat zone. Participants were assigned to the PTSD or comparison group based on Self-Rating Inventory for PTSD scores after deployment. Logistic regression analysis was applied to predict development of a high level of PTSD symptoms based on predeployment GR number, messenger (m)RNA expression of GR target genes FKBP5, GILZ, and SGK1, plasma cortisol, and childhood trauma. We also investigated whether predeployment GR number and FKBP5 mRNA expression were associated with single nucleotide polymorphisms in the GR and FKBP5 genes, either alone or in interaction with childhood trauma.

RESULTS

Several GR pathway components predicted subsequent development of a high level of PTSD symptoms: predeployment high GR number, low FKBP5 mRNA expression, and high GILZ mRNA expression were independently associated with increased risk for a high level of PTSD symptoms. Childhood trauma also independently predicted development of a high level of PTSD symptoms. Additionally, we observed a significant interaction effect of GR haplotype BclI and childhood trauma on GR number.

CONCLUSIONS

Collectively, our results indicate that predeployment GR pathway components are vulnerability factors for subsequent development of a high level of PTSD symptoms.

Inter-ethnic differences in lymphocyte sensitivity to glucocorticoids reflect variation in transcriptional response

Glucocorticoids (GCs) are steroid hormones widely used as pharmaceutical interventions, which act mainly by regulating gene expression levels. A large fraction of patients (~30%), especially those of African descent, show a weak response to treatment. To interrogate the contribution of variable transcriptional response to inter-ethnic differences, we measured in vitro lymphocyte GC sensitivity (LGS) and transcriptome-wide response to GCs in peripheral blood mononuclear cells (PBMCs) from African-American and European-American healthy donors. We found that transcriptional response after 8hrs treatment was significantly correlated with variation in LGS within and between populations. We found that NFKB1, a gene previously found to predict LGS within populations, was more strongly down-regulated in European-Americans on average. NFKB1 could not completely explain population differences, however, and we found an additional 177 genes with population differences in the average log₂ fold change (FDR<0.05), most of which also showed a weaker transcriptional response in African-Americans. These results suggest that inter-ethnic differences in GC sensitivity reflect variation in transcriptional response at many genes, including regulators with large effects (e.g. NFKB1) and numerous other genes with smaller effects.

Activation of the glucocorticoid receptor is associated with poor prognosis in estrogen receptor-negative breast cancer

Estrogen receptor-negative (ER(-)) breast cancers have limited treatment options and are associated with earlier relapses. Because glucocorticoid receptor (GR) signaling initiates antiapoptotic pathways in ER(-) breast cancer cells, we hypothesized that activation of these pathways might be associated with poor prognosis in ER(-) disease. Here we report findings from a genome-wide study of GR transcriptional targets in a premalignant ER(-) cell line model of early breast cancer (MCF10A-Myc) and in primary early-stage ER(-) human tumors. Chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) coupled to time-course expression profiling led us to identify epithelial-to-mesenchymal transition (EMT) pathways as an important aspect associated with GR activation. We validated these findings by carrying out a meta-analysis of primary breast tumor gene expression from 1,378 early-stage breast cancer patients with long-term clinical follow-up, confirming that high levels of GR expression significantly correlated with shorter relapse-free survival in ER(-) patients who were treated or untreated with adjuvant chemotherapy. Notably, in ER(+) breast cancer patients, high levels of GR expression in tumors were significantly associated with better outcome relative to low levels of GR expression. Gene expression analysis revealed that ER(-) tumors expressing high GR levels exhibited differential activation of EMT, cell adhesion, and inflammation pathways. Our findings suggest a direct transcriptional role for GR in determining the outcome of poor-prognosis ER(-) breast cancers.

Interactions between glucocorticoid treatment and cis-regulatory polymorphisms contribute to cellular response phenotypes

Glucocorticoids (GCs) mediate physiological responses to environmental stress and are commonly used as pharmaceuticals. GCs act primarily through the GC receptor (GR, a transcription factor). Despite their clear biomedical importance, little is known about the genetic architecture of variation in GC response. Here we provide an initial assessment of variability in the cellular response to GC treatment by profiling gene expression and protein secretion in 114 EBV-transformed B lymphocytes of African and European ancestry. We found that genetic variation affects the response of nearby genes and exhibits distinctive patterns of genotype-treatment interactions, with genotypic effects evident in either only GC-treated or only control-treated conditions. Using a novel statistical framework, we identified interactions that influence the expression of 26 genes known to play central roles in GC-related pathways (e.g. NQO1, AIRE, and SGK1) and that influence the secretion of IL6.

Effect of chronic mild stress on hippocampal transcriptome in mice selected for high and low stress-induced analgesia and displaying different emotional behaviors

There is increasing evidence that mood disorders may derive from the impact of environmental pressure on genetically susceptible individuals. Stress-induced hippocampal plasticity has been implicated in depression. We studied hippocampal transcriptomes in strains of mice that display high (HA) and low (LA) swim stress-induced analgesia and that differ in emotional behaviors and responses to different classes of antidepressants. Chronic mild stress (CMS) affected expression of a number of genes common for both strains. CMS also produced strain specific changes in expression suggesting that hippocampal responses to stress depend on genotype. Considerably larger number of genes, biological processes, molecular functions, biochemical pathways, and gene networks were affected by CMS in LA than in HA mice. The results suggest that potential drug targets against detrimental effects of stress include glutamate transporters, and cholinergic, cholecystokinin (CCK), glucocorticoids, and thyroid hormones receptors. Furthermore, some biological processes evoked by stress and different between the strains, such as apoptosis, neurogenesis and chromatin modifications, may be responsible for the long-term, irreversible effects of stress and suggest a role for epigenetic regulation of mood related stress responses.

Allele-specific down-regulation of RPTOR expression induced by retinoids contributes to climate adaptations

The mechanistic target of rapamycin (MTOR) pathway regulates cell growth, energy homeostasis, apoptosis, and immune response. The regulatory associated protein of MTOR encoded by the RPTOR gene is a key component of this pathway. A previous survey of candidate genes found that RPTOR contains multiple SNPs with strong correlations between allele frequencies and climate variables, consistent with the action of selective pressures that vary across environments. Using data from a recent genome scan for selection signals, we honed in on a SNP (rs11868112) 26 kb upstream to the transcription start site of RPTOR that exhibits the strongest association with temperature variables. Transcription factor motif scanning and mining of recently mapped transcription factor binding sites identified a binding site for POU class 2 homeobox 1 (POU2F1) spanning the SNP and an adjacent retinoid acid receptor (RAR) binding site. Using expression quantification, chromatin immunoprecipitation (ChIP), and reporter gene assays, we demonstrate that POU2F1 and RARA do bind upstream of the RPTOR gene to regulate its expression in response to retinoids; this regulation is affected by the allele status at rs11868112 with the derived allele resulting in lower expression levels. We propose a model in which the derived allele influences thermogenesis or immune response by altering MTOR pathway activity and thereby increasing fitness in colder climates. Our results show that signatures of genetic adaptations can identify variants with functional effects, consistent with the idea that selection signals may be used for SNP annotation.

Climate has exerted strong selective pressures in human populations during their dispersal, and signatures of these adaptations are still detectable in the geographic distribution of polymorphisms. RPTOR is a key component of the mechanistic target of rapamycin pathway, which regulates cell growth, metabolism, and immune response; and its deregulation is associated with human diseases, including cancer and diabetes. Previous studies showed that variation in RPTOR carry strong signatures of adaptations to different climates. Here, we used evolutionary genetics approaches coupled with transcription factor motif data mining to refine the location of the selection target. We then used functional assays to show that the selected polymorphism resides in a sequence element that regulates gene expression levels in response to retinoids. The derived allele at this SNP, which results in lower expression levels, increases in frequency with decreasing temperatures, consistent with the notion that it confers a selective advantage in colder climates possibly through its effects on energy metabolism or immune response. These results suggest a novel regulatory role for retinoids in MTOR signaling. Moreover, they support the proposal that evolutionary approaches can be informative for SNP functional annotation.

The dissection of transcriptional modules regulated by various drugs of abuse in the mouse striatum

BACKGROUND

Various drugs of abuse activate intracellular pathways in the brain reward system. These pathways regulate the expression of genes that are essential to the development of addiction. To reveal genes common and distinct for different classes of drugs of abuse, we compared the effects of nicotine, ethanol, cocaine, morphine, heroin and methamphetamine on gene expression profiles in the mouse striatum.

RESULTS

We applied whole-genome microarray profiling to evaluate detailed time-courses (1, 2, 4 and 8 hours) of transcriptome alterations following acute drug administration in mice. We identified 42 drug-responsive genes that were segregated into two main transcriptional modules. The first module consisted of activity-dependent transcripts (including Fos and Npas4), which are induced by psychostimulants and opioids. The second group of genes (including Fkbp5 and S3-12), which are controlled, in part, by the release of steroid hormones, was strongly activated by ethanol and opioids. Using pharmacological tools, we were able to inhibit the induction of particular modules of drug-related genomic profiles. We selected a subset of genes for validation by in situ hybridization and quantitative PCR. We also showed that knockdown of the drug-responsive genes Sgk1 and Tsc22d3 resulted in alterations to dendritic spines in mice, possibly reflecting an altered potential for plastic changes.

CONCLUSIONS

Our study identified modules of drug-induced genes that share functional relationships. These genes may play a critical role in the early stages of addiction.

Multiple Functional Variants in cis Modulate PDYN Expression

Understanding genetic variation and its functional consequences within cis-regulatory regions remains an important challenge in human genetics and evolution. Here, we present a fine-scale functional analysis of segregating variation within the cis-regulatory region of prodynorphin, a gene that encodes an endogenous opioid precursor with roles in cognition and disease. In order to characterize the functional consequences of segregating variation in cis in a region under balancing selection in different human populations, we examined associations between specific polymorphisms and gene expression in vivo and in vitro. We identified five polymorphisms within the 5' flanking region that affect transcript abundance: a 68-bp repeat recognized in prior studies, as well as two microsatellites and two single nucleotide polymorphisms not previously implicated as functional variants. The impact of these variants on transcription differs by brain region, sex, and cell type, implying interactions between cis genotype and the differentiated state of cells. The effects of individual variants on expression level are not additive in some combinations, implying epistatic interactions between nearby variants. These data reveal an unexpectedly complex relationship between segregating genetic variation and its expression-trait consequences and highlights the importance of close functional scrutiny of natural genetic variation within even relatively well-studied cis-regulatory regions.

Spatial patterns of variation due to natural selection in humans

Empowered by technology and sampling efforts designed to facilitate genome-wide association mapping, human geneticists are now studying the geography of genetic variation in unprecedented detail. With high genomic coverage and geographic resolution, these studies are identifying loci with spatial signatures of selection, such as extreme levels of differentiation and correlations with environmental variables. Collectively, patterns at these loci are beginning to provide new insights into the process of human adaptation. Here, we review the challenges of these studies and emerging results, including how human population structure has influenced the response to novel selective pressures.