The Krüppel-like factor 9 cistrome in mouse hippocampal neurons reveals predominant transcriptional repression via proximal promoter binding

Comparative analysis of Krüppel-like factors expression in the retinas of zebrafish and mice during development and after injury

The Krüppel-like factors (KLFs) have emerged as important transcriptional regulators of various cellular processes, including neural development. Some of them have been described as intrinsic factors involved in axon regeneration in the central nervous system (CNS) of vertebrates. Zebrafish are known for their ability to regenerate several tissues in adulthood, including the CNS, a capability lost during vertebrate evolution and absent in adult mammals. The role that KLFs could play in this differential ability remains unknown. Therefore, in this study, we analyzed the endogenous response of certain KLFs implicated in axon regeneration (KLFs 6, 7, 9, and 13) during retina development and after axon injury. The results showed that the expression of Klfs 6, 7, and 13 decreases in the developing retina of mice but not in zebrafish, while the mRNA levels of Klf9 strongly increase in both species. The response to injury was further analyzed using optic nerve crush (ONC) as a model of lesion. Our analysis during the acute phase (hours) demonstrated an induction of Klfs 6 and 7 expression exclusively in the zebrafish retina, while Klfs 9 and 13 mRNA levels increased in both species. Further analysis of the chronic response (days) showed that mRNA levels of Klf6 transiently increase in the retinas of both zebrafish and mice, whereas those of Klf7 decrease later after optic nerve injury. In addition, the analysis revealed that the expression of Klf9 decreases, while that of Klf13 increases in the retinas of zebrafish in response to optic nerve injury but remains unaltered in mice. Altogether, these findings support the hypothesis that KLFs may play a role in the differential axon regeneration abilities exhibited by fish and mice.

Krüppel-like Factor-9 and Krüppel-like Factor-13: Highly Related, Multi-Functional, Transcriptional Repressors and Activators of Oncogenesis

Specificity Proteins/Krüppel-like Factors (SP/KLF family) are a conserved family of transcriptional regulators. These proteins share three highly conserved, contiguous zinc fingers in their carboxy-terminus, requisite for binding to cis elements in DNA. Each SP/KLF protein has unique primary sequence within its amino-terminal and carboxy-terminal regions, and it is these regions which interact with co-activators, co-repressors, and chromatin-modifying proteins to support the transcriptional activation and repression of target genes. Krüppel-like Factor 9 (KLF9) and Krüppel-like Factor 13 (KLF13) are two of the smallest members of the SP/KLF family, are paralogous, emerged early in metazoan evolution, and are highly conserved. Paradoxically, while most similar in primary sequence, KLF9 and KLF13 display many distinct roles in target cells. In this article, we summarize the work that has identified the roles of KLF9 (and to a lesser degree KLF13) in tumor suppression or promotion via unique effects on differentiation, pro- and anti-inflammatory pathways, oxidative stress, and tumor immune cell infiltration. We also highlight the great diversity of miRNAs, lncRNAs, and circular RNAs which provide mechanisms for the ubiquitous tumor-specific suppression of KLF9 mRNA and protein. Elucidation of KLF9 and KLF13 in cancer biology is likely to provide new inroads to the understanding of oncogenesis and its prevention and treatments.

Klf9 plays a critical role in GR -dependent metabolic adaptations in cardiomyocytes

Glucocorticoids through activation of the Glucocorticoid receptor (GR) play an essential role in cellular homeostasis during physiological variations and in response to stress. Our genomic GR binding and transcriptome data from Dexamethasone (Dex) treated cardiomyocytes showed an early differential regulation of mostly transcription factors, followed by sequential change in genes involved in downstream functional pathways. We examined the role of Krüppel-like factor 9 (Klf9), an early direct target of GR in cardiomyocytes. Klf9-ChIPseq identified 2150 genes that showed an increase in Klf9 binding in response to Dex. Transcriptome analysis of Dex treated cardiomyocytes with or without knockdown of Klf9 revealed differential regulation of 1777 genes, of which a reversal in expression is seen in 1640 genes with knockdown of Klf9 compared to Dex. Conversely, only 137 (∼8%) genes show further dysregulation in expression with siKLf9, as seen with Dex treated cardiomyocytes. Functional annotation identified genes of metabolic pathways on the top of differentially expressed genes, including those involved in glycolysis and oxidative phosphorylation. Knockdown of Klf9 in cardiomyocytes inhibited Dex induced increase in glycolytic function and mitochondrial spare respiratory capacity, as measured by glycolysis and mito stress tests, respectively. Thus, we conclude that cyclic, diurnal GR activation, through Klf9 -dependent feedforward signaling plays a central role in maintaining cellular homeostasis through metabolic adaptations in cardiomyocytes.

Loss of Krüppel-like factor 9 deregulates both physiological gene expression and development

Krüppel-like factor 9 (Klf9) is a ubiquitously expressed transcription factor that is a feedforward regulator of multiple stress-responsive and endocrine signaling pathways. We previously described how loss of Klf9 function affects the transcriptome of zebrafish larvae sampled at a single time point 5 days post-fertilization (dpf). However, klf9 expression oscillates diurnally, and the sampled time point corresponded to its expression nadir. To determine if the transcriptomic effects of the klf9-/- mutation vary with time of day, we performed bulk RNA-seq on 5 dpf zebrafish embryos sampled at three timepoints encompassing the predawn peak and midmorning nadir of klf9 expression. We found that while the major effects of the klf9-/- mutation that we reported previously are robust to time of day, the mutation has additional effects that manifest only at the predawn time point. We used a published single-cell atlas of zebrafish development to associate the effects of the klf9-/- mutation with different cell types and found that the mutation increased mRNA associated with digestive organs (liver, pancreas, and intestine) and decreased mRNA associated with differentiating neurons and blood. Measurements from confocally-imaged larvae suggest that overrepresentation of liver mRNA in klf9-/- mutants is due to development of enlarged livers.

KLF13 Regulates the Activity of the GH-Induced JAK/STAT Signaling by Targeting Genes Involved in the Pathway

The Krüppel-like factor 13 (KLF13) has emerged as an important transcription factor involved in essential processes of the central nervous system (CNS). It predominantly functions as a transcriptional repressor, impacting the activity of several signaling pathways with essential roles in the CNS, including the JAK/STAT pathway, which is the canonical mediator of growth hormone (GH) signaling. It is now recognized that GH has important actions as a neurotrophic factor. Therefore, we analyzed the effects of KLF13 on the activity of the JAK/STAT signaling pathway in the hippocampus-derived cell line HT22. Results showed that KLF13 directly regulates the expression of several genes involved in the JAK-STAT pathway, including Jak1, Jak2, Jak3, and Socs1, by associating with their proximal gene promoters. In addition, it was found that in KLF13-deficient HT22 neurons, the expression of Jak1, Stat3, Socs1, Socs3, and Igf1 was dysregulated, exhibiting mRNA levels that went up to 7-fold higher than the control cell line. KLF13 displayed a differential effect on the GH-induced JAK/STAT pathway activity, decreasing the STAT3 branch while enhancing the STAT5 branch. In KLF13-deficient HT22 cells, the activity of the STAT3 branch was enhanced, mediating the GH-dependent augmented expression of the JAK/STAT output genes Socs1, Socs3, Igf1, and Bdnf. Furthermore, GH treatment increased both the nuclear content of KLF13 and Klf13 mRNA levels, suggesting that KLF13 could be part of the mechanisms that maintain the homeostatic state of this pathway. These findings support the notion that KLF13 is a regulator of JAK/STAT activity.

The Role of Krüppel-like Factors in Pancreatic Physiology and Pathophysiology

Krüppel-like factors (KLFs) belong to the family of transcription factors with three highly conserved zinc finger domains in the C-terminus. They regulate homeostasis, development, and disease progression in many tissues. It has been shown that KLFs play an essential role in the endocrine and exocrine compartments of the pancreas. They are necessary to maintain glucose homeostasis and have been implicated in the development of diabetes. Furthermore, they can be a vital tool in enabling pancreas regeneration and disease modeling. Finally, the KLF family contains proteins that act as tumor suppressors and oncogenes. A subset of members has a biphasic function, being upregulated in the early stages of oncogenesis and stimulating its progression and downregulated in the late stages to allow for tumor dissemination. Here, we describe KLFs' function in pancreatic physiology and pathophysiology.

Epigenetic and transcriptomic characterization reveals progression markers and essential pathways in clear cell renal cell carcinoma

Identifying tumor-cell-specific markers and elucidating their epigenetic regulation and spatial heterogeneity provides mechanistic insights into cancer etiology. Here, we perform snRNA-seq and snATAC-seq in 34 and 28 human clear cell renal cell carcinoma (ccRCC) specimens, respectively, with matched bulk proteogenomics data. By identifying 20 tumor-specific markers through a multi-omics tiered approach, we reveal an association between higher ceruloplasmin (CP) expression and reduced survival. CP knockdown, combined with spatial transcriptomics, suggests a role for CP in regulating hyalinized stroma and tumor-stroma interactions in ccRCC. Intratumoral heterogeneity analysis portrays tumor cell-intrinsic inflammation and epithelial-mesenchymal transition (EMT) as two distinguishing features of tumor subpopulations. Finally, BAP1 mutations are associated with widespread reduction of chromatin accessibility, while PBRM1 mutations generally increase accessibility, with the former affecting five times more accessible peaks than the latter. These integrated analyses reveal the cellular architecture of ccRCC, providing insights into key markers and pathways in ccRCC tumorigenesis.

Krüppel-like factor 9 (KLF9) links hormone dysregulation and circadian disruption to breast cancer pathogenesis

BACKGROUND

Circadian disruption is an emerging driver of breast cancer (BCa), with epidemiological studies linking shift work and chronic jet lag to increased BCa risk. Indeed, several clock genes participate in the gating of mitotic entry, regulation of DNA damage response, and epithelial-to-mesenchymal transition, thus impacting BCa etiology. Dysregulated estrogen (17β-estradiol, E2) and glucocorticoid (GC) signaling prevalent in BCa may further contribute to clock desynchrony by directly regulating the expression and cycling dynamics of genes comprising the local breast oscillator. In this study, we investigated the tumor suppressor gene, Krüppel-like factor 9 (KLF9), as an important point of crosstalk between hormone signaling and the circadian molecular network, and further examine its functional role in BCa.

METHODS

Through meta-analysis of publicly available RNA- and ChIP-sequencing datasets from BCa tumor samples and cell lines, and gene expression analysis by RT-qPCR and enhancer- reporter assays, we elucidated the molecular mechanism behind the clock and hormone regulation of KLF9. Lentiviral knockdown and overexpression of KLF9 in three distinct breast epithelial cell lines (MCF10A, MCF7 and MDA-MB-231) was generated to demonstrate the role of KLF9 in orthogonal assays on breast epithelial survival, proliferation, apoptosis, and migration.

RESULTS

We determined that KLF9 is a direct GC receptor target in mammary epithelial cells, and that induction is likely mediated through coordinate transcriptional activation from multiple GC-responsive enhancers in the KLF9 locus. More interestingly, rhythmic expression of KLF9 in MCF10A cells was abolished in the highly aggressive MDA-MB-231 line. In turn, forced expression of KLF9 altered the baseline and GC/E2-responsive expression of several clock genes, indicating that KLF9 may function as a regulator of the core clock machinery. Characterization of the role of KLF9 using complementary cancer hallmark assays in the context of the hormone-circadian axis revealed that KLF9 plays a tumor-suppressive role in BCa regardless of molecular subtype. KLF9 potentiated the anti-tumorigenic effects of GC in E2 receptor + luminal MCF7 cells, while it restrained GC-enhanced oncogenicity in triple-negative MCF10A and MDA-MB-231 cells.

CONCLUSIONS

Taken together, our findings support that dysregulation of KLF9 expression and oscillation in BCa impinges on circadian network dynamics, thus ultimately affecting the BCa oncogenic landscape.

Krüppel-like factors in tumors: Key regulators and therapeutic avenues

Krüppel-like factors (KLFs) are a group of DNA-binding transcriptional regulators with multiple essential functions in various cellular processes, including proliferation, migration, inflammation, and angiogenesis. The aberrant expression of KLFs is often found in tumor tissues and is essential for tumor development. At the molecular level, KLFs regulate multiple signaling pathways and mediate crosstalk among them. Some KLFs may also be molecular switches for specific biological signals, driving their transition from tumor suppressors to promoters. At the histological level, the abnormal expression of KLFs is closely associated with tumor cell stemness, proliferation, apoptosis, and alterations in the tumor microenvironment. Notably, the role of each KLF in tumors varies according to tumor type and different stages of tumor development rather than being invariant. In this review, we focus on the advances in the molecular biology of KLFs, particularly the regulations of several classical signaling pathways by these factors, and the critical role of KLFs in tumor development. We also highlight their strong potential as molecular targets in tumor therapy and suggest potential directions for clinical translational research.

Iodide intake during pregnancy and lactation stimulates KLF9, BDNF expression in offspring brain with elevated DHA, EPA metabolites

To investigate the effect of different iodide intake during pregnancy and lactation on thyroid function, docosahexaenoic acid (DHA), Eicosapentaenoic acid (EPA) metabolites, the expression of Krüppel-like factor KLF9 (KLF9), brain-derived neurotrophic factor (BDNF) in brain in offspring rats. In both male and female offspring rats, serum FT3, FT4 levels and the expression of KLF9, thyroid hormone receptors (TR)α, TRβ and BDNF in the hippocampal region and cerebellum were significantly increased in 5 times higher-than-normal pregnant iodide intake (5 HI) and 10 times higher-than-normal pregnant iodide intake (10 HI) group. The median levels of DHA metabolite (17-HDoHE) and EPA metabolites (15-HEPE, 17,18-EEQ, 9-HEPE and 14,15-DiHETE) were significantly increased in 5 HI and 10 HI group of offspring rats. Serum DHA, EPA metabolites and KLF9 as well as BDNF expression in brain might be potential iodine status biomarkers to reflect brain development in offspring.

Mutated KLF4(K409Q) in meningioma binds STRs and activates FGF3 gene expression

Krüppel-like factor 4 (KLF4) is a transcription factor that has been proven necessary for both induction and maintenance of pluripotency and self-renewal. Whole-genome sequencing defined a unique mutation in KLF4 (KLF4K409Q) in human meningiomas. However, the molecular mechanism of this tumor-specific KLF4 mutation is unknown. Using genome-wide high-throughput and focused quantitative transcriptional approaches in human cell lines, primary meningeal cells, and meningioma tumor tissue, we found that a change in the evolutionarily conserved DNA-binding domain of KLF4 alters its DNA recognition preference, resulting in a shift in downstream transcriptional activity. In the KLF4K409Q-specific targets, the normally silent fibroblast growth factor 3 (FGF3) is activated. We demonstrated a neomorphic function of KLF4K409Q in stimulating FGF3 transcription through binding to its promoter and in using short tandem repeats (STRs) located within the locus as enhancers.

SNX5 suppresses clear cell renal cell carcinoma progression by inducing CD44 internalization and epithelial-to-mesenchymal transition

Aberrant expression of SNX5 can contribute to tumorigenesis, invasion, and metastasis of several human cancers. However, the clinicopathological and biological significance of SNX5 in clear cell renal cell carcinoma (ccRCC) remain unclear. In this study, we found that SNX5 expression was downregulated and negatively correlated with tumor size, American Joint Committee on Cancer stage, tumor thrombus of inferior vena cava, and poor prognosis in human ccRCC. Ectopic expression of SNX5 inhibited ccRCC cell proliferation and metastasis, whereas knockdown of SNX5 increased these activities both in vitro and in vivo. Mechanistically, overexpression of SNX5 blocked internalization and intracellular trafficking of CD44 in ccRCC cells. Knockdown of SNX5 was associated with epithelial-to-mesenchymal transition (EMT) in ccRCC cells. Overexpression of SNX5 inhibited TGF-β-induced migration, invasion, and EMT in ccRCC cells. KLF9 directly bound to the SNX5 promoter and increased SNX5 transcription. Moreover, we found that the combination of SNX5 and CD44 or E-cadherin or KLF9 was a more powerful predictor of poor prognosis than either parameter alone. Collectively, our data reveal a mechanism that KLF9-mediated SNX5 expression was associated with poor prognosis via trafficking of CD44 and promoting EMT in ccRCC. SNX5 may be a potential prognostic biomarker and therapeutic target for patients with ccRCC.

Transcriptional regulation of neural stem cell expansion in the adult hippocampus

Experience governs neurogenesis from radial-glial neural stem cells (RGLs) in the adult hippocampus to support memory. Transcription factors (TFs) in RGLs integrate physiological signals to dictate self-renewal division mode. Whereas asymmetric RGL divisions drive neurogenesis during favorable conditions, symmetric divisions prevent premature neurogenesis while amplifying RGLs to anticipate future neurogenic demands. The identities of TFs regulating RGL symmetric self-renewal, unlike those that regulate RGL asymmetric self-renewal, are not known. Here, we show in mice that the TF Kruppel-like factor 9 (Klf9) is elevated in quiescent RGLs and inducible, deletion of Klf9 promotes RGL activation state. Clonal analysis and longitudinal intravital two-photon imaging directly demonstrate that Klf9 functions as a brake on RGL symmetric self-renewal. In vivo translational profiling of RGLs lacking Klf9 generated a molecular blueprint for RGL symmetric self-renewal that was characterized by upregulation of genetic programs underlying Notch and mitogen signaling, cell cycle, fatty acid oxidation, and lipogenesis. Together, these observations identify Klf9 as a transcriptional regulator of neural stem cell expansion in the adult hippocampus.

In humans and other mammals, a region of the brain known as the hippocampus plays important roles in memory. New experiences guide cells in the hippocampus known as radial-glial neural stem cells (RGLs) to divide to make new neurons and other types of cells involved in forming memories.

Each time an RGL divides, it can choose to divide asymmetrically to maintain a copy of itself and make a new cell of another type, or divide symmetrically (a process known as symmetric self-renewal) to produce two RGLs. Symmetric self-renewal helps to restore and replenish the pool of stem cells in the hippocampus that are lost due to injury or age, allowing us to continue making new neurons.

Proteins known as transcription factors are believed to control how RGLs divide. Previous studies have identified several transcription factors that regulate the RGLs splitting asymmetrically to make neurons and other cells. But the identities of the transcription factors that regulate symmetric self-renewal in the adult hippocampus have remained elusive.

Here, Guo et al. searched for transcription factors that regulate symmetric self-renewal of RGLs in mice. The experiments found that RGLs that are resting and not dividing (referred to as ‘quiescent’) have higher levels of a transcription factor called Klf9 than RGLs that are actively dividing. Loss of the gene encoding Klf9 triggered quiescent RGLs to start dividing, and further experiments showed that Klf9 directly inhibited symmetric self-renewal. Guo et al. then used an approach called in vivo translational profiling to generate a blueprint that revealed new insights into the molecular processes involved in this symmetric division.

These findings pave the way for researchers to develop strategies that may expand the numbers of stem cells in the hippocampus. This could eventually be used to help replenish brain circuits with neurons and improve the memory of individuals with Alzheimer’s disease or other conditions that cause memory loss.

Glucocorticoid-Mediated Developmental Programming of Vertebrate Stress Responsivity

Glucocorticoids, vertebrate steroid hormones produced by cells of the adrenal cortex or interrenal tissue, function dynamically to maintain homeostasis under constantly changing and occasionally stressful environmental conditions. They do so by binding and thereby activating nuclear receptor transcription factors, the Glucocorticoid and Mineralocorticoid Receptors (MR and GR, respectively). The GR, by virtue of its lower affinity for endogenous glucocorticoids (cortisol or corticosterone), is primarily responsible for transducing the dynamic signals conveyed by circadian and ultradian glucocorticoid oscillations as well as transient pulses produced in response to acute stress. These dynamics are important determinants of stress responsivity, and at the systemic level are produced by feedforward and feedback signaling along the hypothalamus-pituitary-adrenal/interrenal axis. Within receiving cells, GR signaling dynamics are controlled by the GR target gene and negative feedback regulator fkpb5. Chronic stress can alter signaling dynamics via imperfect physiological adaptation that changes systemic and/or cellular set points, resulting in chronically elevated cortisol levels and increased allostatic load, which undermines health and promotes development of disease. When this occurs during early development it can "program" the responsivity of the stress system, with persistent effects on allostatic load and disease susceptibility. An important question concerns the glucocorticoid-responsive gene regulatory network that contributes to such programming. Recent studies show that klf9, a ubiquitously expressed GR target gene that encodes a Krüppel-like transcription factor important for metabolic plasticity and neuronal differentiation, is a feedforward regulator of GR signaling impacting cellular glucocorticoid responsivity, suggesting that it may be a critical node in that regulatory network.

Synergistic gene regulation by thyroid hormone and glucocorticoid in the hippocampus

The hippocampus is considered the center for learning and memory in the brain, and its development and function is greatly affected by the thyroid and stress axes. Thyroid hormone (TH) and glucocorticoids (GC) are known to have a synergistic effect on developmental programs across several vertebrate species, and their effects on hippocampal structure and function are well-documented. However, there are few studies that focus on the processes and genes that are cooperatively regulated by the two hormone axes. Cross-regulation of the thyroid and stress axes in the hippocampus occurs on multiple levels such that TH can regulate the expression of the GC receptor (GR) while GC can modulate tissue sensitivity to TH by controlling the expression of TH receptor (TR) and enzymes involved in TH biosynthesis. Thyroid hormone and GC are also known to synergistically regulate the transcription of genes associated with neuronal function and development. Synergistic gene regulation by TH and GC may occur through the direct, cooperative action of TR and GR on common target genes, or by indirect mechanisms involving gene regulatory cascades activated by TR and GR. In this chapter, we describe the known physiological effects and underlying molecular mechanisms of TH and GC synergistic gene regulation in the hippocampus.

Glucocorticoid-Responsive Transcription Factor Krüppel-Like Factor 9 Regulates fkbp5 and Metabolism

Krüppel-like factor 9 (Klf9) is a feedforward regulator of glucocorticoid receptor (GR) signaling. Here we show that in zebrafish klf9 is expressed with GR-dependent oscillatory dynamics in synchrony with fkbp5, a GR target that encodes a negative feedback regulator of GR signaling. We found that fkbp5 transcript levels are elevated in klf9 -/- mutants and that Klf9 associates with chromatin at the fkbp5 promoter, which becomes hyperacetylated in klf9 -/ - mutants, suggesting that the GR regulates fkbp5 via an incoherent feedforward loop with klf9. As both the GR and Fkbp5 are known to regulate metabolism, we asked how loss of Klf9 affects metabolic rate and gene expression. We found that klf9 -/- mutants have a decreased oxygen consumption rate (OCR) and upregulate glycolytic genes, the promoter regions of which are enriched for potential Klf9 binding motifs. Our results suggest that Klf9 functions downstream of the GR to regulate cellular glucocorticoid responsivity and metabolic homeostasis.

An Intact Krüppel-like factor 9 Gene Is Required for Acute Liver Period 1 mRNA Response to Restraint Stress

The clock protein period 1 (PER1) is a central component of the core transcription-translation feedback loop governing cell-autonomous circadian rhythms in animals. Transcription of Per1 is directly regulated by the glucocorticoid (GC) receptor (GR), and Per1 mRNA is induced by stressors or injection of GC. Circulating GCs may synchronize peripheral clocks with the central pacemaker located in the suprachiasmatic nucleus of the brain. Krüppel-like factor 9 (KLF9) is a zinc finger transcription factor that, like Per1, is directly regulated by liganded GR, and it associates in chromatin at clock and clock-output genes, including at Per1. We hypothesized that KLF9 modulates stressor-dependent Per1 transcription. We exposed wild-type (WT) and Klf9 null mice (Klf9-/-) of both sexes to 1 hour restraint stress, which caused similar 2- to 2.5-fold increases in plasma corticosterone (B) in each genotype and sex. Although WT mice of both sexes showed a 2-fold increase in liver Per1 mRNA level after restraint stress, this response was absent in Klf9-/- mice. However, injection of B in WT and Klf9-/- mice induced similar increases in Per1 mRNA. Our findings support that an intact Klf9 gene is required for liver Per1 mRNA responses to an acute stressor, but a possible role for GCs in this response requires further investigation.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L

CONTEXT

Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited.

OBJECTIVE

To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts.

METHOD

We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed.

RESULTS

We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH.

CONCLUSIONS

This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Genomic determinants implicated in the glucocorticoid-mediated induction of KLF9 in pulmonary epithelial cells

Ligand-activated glucocorticoid receptor (GR) elicits variable glucocorticoid-modulated transcriptomes in different cell types. However, some genes, including Krüppel-like factor 9 (KLF9), a putative transcriptional repressor, demonstrate conserved responses. We show that glucocorticoids induce KLF9 expression in the human airways in vivo and in differentiated human bronchial epithelial (HBE) cells grown at air-liquid interface (ALI). In A549 and BEAS-2B pulmonary epithelial cells, glucocorticoids induce KLF9 expression with similar kinetics to primary HBE cells in submersion culture. A549 and BEAS-2B ChIP-seq data reveal four common glucocorticoid-induced GR binding sites (GBSs). Two GBSs mapped to the 5'-proximal region relative to KLF9 transcription start site (TSS) and two occurred at distal sites. These were all confirmed in primary HBE cells. Global run-on (GRO) sequencing indicated robust enhancer RNA (eRNA) production from three of these GBSs in BEAS-2B cells. This was confirmed in A549 cells, plus submersion, and ALI culture of HBE cells. Cloning each GBS into luciferase reporters revealed glucocorticoid-induced activity requiring a glucocorticoid response element (GRE) within each distal GBS. While the proximal GBSs drove modest reporter induction by glucocorticoids, this region exhibited basal eRNA production, RNA polymerase II enrichment, and looping to the TSS, plausibly underlying constitutive KLF9 expression. Post glucocorticoid treatment, interactions between distal and proximal GBSs and the TSS correlated with KLF9 induction. CBP/P300 silencing reduced proximal GBS activity, but negligibly affected KLF9 expression. Overall, a model for glucocorticoid-mediated regulation of KLF9 involving multiple GBSs is depicted. This work unequivocally demonstrates that mechanistic insights gained from cell lines can translate to physiologically relevant systems.

Krüppel-Like Factors 9 and 13 Block Axon Growth by Transcriptional Repression of Key Components of the cAMP Signaling Pathway

Krüppel-like factors (KLFs) are zinc finger transcription factors implicated in diverse biological processes, including differentiation of neural cells. The ability of mammalian neurons to elongate axons decreases during postnatal development in parallel with a decrease in cAMP, and increase in expression of several Klf genes. The paralogous KLFs 9 and 13 inhibit neurite outgrowth, and we hypothesized that their actions are mediated through repression of cAMP signaling. To test this we used the adult mouse hippocampus-derived cell line HT22 engineered to control expression of Klf9 or Klf13 with doxycycline, or made deficient for these Klfs by CRISPR/Cas9 genome editing. We also used primary hippocampal cells isolated from wild type, Klf9^–/– and Klf13^–/– mice. Forced expression of Klf9 or Klf13 in HT22 changed the mRNA levels of several genes involved with cAMP signaling; the predominant action was gene repression, and KLF13 influenced ∼4 times more genes than KLF9. KLF9 and KLF13 repressed promoter activity of the protein kinase a catalytic subunit alpha gene in transfection-reporter assays; KLF13, but not KLF9 repressed the calmodulin 3 promoter. Forskolin activation of a cAMP-dependent promoter was reduced after forced expression of Klf9 or Klf13, but was enhanced in Klf gene knockout cells. Forced expression of Klf9 or Klf13 blocked cAMP-dependent neurite outgrowth in HT22 cells, and axon growth in primary hippocampal neurons, while Klf gene knockout enhanced the effect of elevated cAMP. Taken together, our findings show that KLF9 and KLF13 inhibit neurite/axon growth in hippocampal neurons, in part, by inhibiting the cAMP signaling pathway.

Klf9 is a key feedforward regulator of the transcriptomic response to glucocorticoid receptor activity

The zebrafish has recently emerged as a model system for investigating the developmental roles of glucocorticoid signaling and the mechanisms underlying glucocorticoid-induced developmental programming. To assess the role of the Glucocorticoid Receptor (GR) in such programming, we used CRISPR-Cas9 to produce a new frameshift mutation, GR^369-, which eliminates all potential in-frame initiation codons upstream of the DNA binding domain. Using RNA-seq to ask how this mutation affects the larval transcriptome under both normal conditions and with chronic cortisol treatment, we find that GR mediates most of the effects of the treatment, and paradoxically, that the transcriptome of cortisol-treated larvae is more like that of larvae lacking a GR than that of larvae with a GR, suggesting that the cortisol-treated larvae develop GR resistance. The one transcriptional regulator that was both underexpressed in GR^369- larvae and consistently overexpressed in cortisol-treated larvae was klf9. We therefore used CRISPR-Cas9-mediated mutation of klf9 and RNA-seq to assess Klf9-dependent gene expression in both normal and cortisol-treated larvae. Our results indicate that Klf9 contributes significantly to the transcriptomic response to chronic cortisol exposure, mediating the upregulation of proinflammatory genes that we reported previously.

Molecular Mechanisms for Krüppel-Like Factor 13 Actions in Hippocampal Neurons

Krüppel-like factors (KLFs) play key roles in nervous system development and function. Several KLFs are known to promote, and then maintain neural cell differentiation. Our previous work focused on the actions of KLF9 in mouse hippocampal neurons. Here we investigated genomic targets and functions of KLF9's paralog KLF13, with the goal of understanding how these two closely related transcription factors influence hippocampal cell function, proliferation, survival, and regeneration. We engineered the adult mouse hippocampus-derived cell line HT22 to control Klf13 expression with doxycycline. We also generated HT22 Klf13 knock out cells, and we analyzed primary hippocampal cells from wild type and Klf13^-/- mice. RNA sequencing showed that KLF13, like KLF9, acts predominantly as a transcriptional repressor in hippocampal neurons and can regulate other Klf genes. Pathway analysis revealed that genes regulated by KLF13 are involved in cell cycle, cell survival, cytoarchitecture regulation, among others. Chromatin-streptavidin sequencing conducted on chromatin isolated from HT22 cells expressing biotinylated KLF13 identified 9506 genomic targets; 79% were located within 1-kb upstream of transcription start sites. Transfection-reporter assays confirmed that KLF13 can directly regulate transcriptional activity of its target genes. Comparison of the target genes of KLF9 and KLF13 found that they share some functions that were likely present in their common ancestor, but they have also acquired distinct functions during evolution. Flow cytometry showed that KLF13 promotes cell cycle progression, and it protects cells from glutamate-induced excitotoxic damage. Taken together, our findings establish novel roles and molecular mechanisms for KLF13 actions in mammalian hippocampal neurons.

The Paralogous Krüppel-like Factors 9 and 13 Regulate the Mammalian Cellular Circadian Clock Output Gene Dbp

An intricate transcription-translation feedback loop (TTFL) governs cellular circadian rhythms in mammals. Here, we report that the zinc finger transcription factor Krüppel-like factor 9 (KLF9) is regulated by this TTFL, it associates in chromatin at the core circadian clock and clock-output genes, and it acts to modulate transcription of the clock-output gene Dbp. Our earlier genome-wide analysis of the mouse hippocampus-derived cell line HT22 showed that KLF9 associates in chromatin with Per1, Per3, Dbp, Tef, Bhlhe40, Bhlhe41, Nr1d1, and Nr1d2. Of the 3514 KLF9 peaks identified in HT22 cells, 1028 contain E-box sequences to which the transcriptional activators CLOCK and BMAL1 may bind, a frequency significantly greater than expected by chance. Klf9 mRNA showed circadian oscillation in synchronized HT22 cells, mouse hippocampus, and liver. At the clock-output gene Dbp, KLF9 exhibited circadian rhythmicity in its association in chromatin in HT22 cells and hippocampus. Forced expression of KLF9 in HT22 cells repressed basal Dbp transcription and strongly inhibited CLOCK+BMAL1-dependent transcriptional activation of a transfected Dbp reporter. Mutational analysis showed that this action of KLF9 depended on 2 intact KLF9-binding motifs within the Dbp locus that are in close proximity to E-boxes. Knockout of Klf9 or the paralogous gene Klf13 using CRISPR/Cas9 genome editing in HT22 cells had no effect on Dbp expression, but combined knockout of both genes strongly impaired circadian Dbp mRNA oscillation. Like KLF9, KLF13 also showed association in chromatin with clock- and clock-output genes, and forced expression of KLF13 inhibited the actions of CLOCK+BMAL1 on Dbp transcription. Our results suggest novel and partly overlapping roles for KLF9 and KLF13 in modulating cellular circadian clock output by a mechanism involving direct interaction with the core TTFL.

Thyroid Hormone Protects Primary Cortical Neurons Exposed to Hypoxia by Reducing DNA Methylation and Apoptosis

Traumatic brain injury (TBI) is associated with disruption of cerebral blood flow leading to localized brain hypoxia. Thyroid hormone (TH) treatment, administered shortly after injury, has been shown to promote neural protection in rodent TBI models. The mechanism of TH protection, however, is not established. We used mouse primary cortical neurons to investigate the effectiveness and possible pathways of T3-promoted cell survival after exposure to hypoxic injury. Cultured primary cortical neurons were exposed to hypoxia (0.2% oxygen) for 7 hours with or without T3 (5 nM). T3 treatment enhanced DNA 5-hydroxymethylcytosine levels and attenuated the hypoxia-induced increase in DNA 5-methylcytosine (5-mc). In the presence of T3, mRNA expression of Tet family genes was increased and DNA methyltransferase (Dnmt) 3a and Dnmt3b were downregulated, compared with conditions in the absence of T3. These T3-induced changes decreased hypoxia-induced DNA de novo methylation, which reduced hypoxia-induced neuronal damage and apoptosis. We used RNA sequencing to characterize T3-regulated genes in cortical neurons under hypoxic conditions and identified 22 genes that were upregulated and 15 genes that were downregulated. Krüppel-like factor 9 (KLF9), a multifunctional transcription factor that plays a key role in central nervous system development, was highly upregulated by T3 treatment in hypoxic conditions. Knockdown of the KLF9 gene resulted in early apoptosis and abolished the beneficial role of T3 in neuronal survival. KLF9 mediates, in part, the neuronal protective role of T3. T3 treatment reduces hypoxic damage, although pathways that reduce DNA methylation and apoptosis remain to be elucidated.

Coordinated transcriptional regulation by thyroid hormone and glucocorticoid interaction in adult mouse hippocampus-derived neuronal cells

The hippocampus is a well-known target of thyroid hormone (TH; e.g., 3,5,3’-triiodothyronine—T₃) and glucocorticoid (GC; e.g., corticosterone—CORT) action. Despite evidence that TH and GC play critical roles in neural development and function, few studies have identified genes and patterns of gene regulation influenced by the interaction of these hormones at a genome-wide scale. In this study we investigated gene regulation by T₃, CORT, and T₃ + CORT in the mouse hippocampus-derived cell line HT-22. We treated cells with T₃, CORT, or T₃ + CORT for 4 hr before cell harvest and RNA isolation for microarray analysis. We identified 9 genes regulated by T₃, 432 genes by CORT, and 412 genes by T₃ + CORT. Among the 432 CORT-regulated genes, there were 203 genes that exhibited an altered CORT response in the presence of T₃, suggesting that T₃ plays a significant role in modulating CORT-regulated genes. We also found 80 genes synergistically induced, and 73 genes synergistically repressed by T₃ + CORT treatment. We performed in silico analysis using publicly available mouse neuronal chromatin immunoprecipitation-sequencing datasets and identified a considerable number of synergistically regulated genes with TH receptor and GC receptor peaks mapping within 1 kb of chromatin marks indicative of hormone-responsive enhancer regions. Functional annotation clustering of synergistically regulated genes reveal the relevance of proteasomal-dependent degradation, neuroprotective effect of growth hormones, and neuroinflammatory responses as key pathways to how TH and GC may coordinately influence learning and memory. Taken together, our transcriptome data represents a promising exploratory dataset for further study of common molecular mechanisms behind synergistic TH and GC gene regulation, and identify specific genes and their role in processes mediated by cross-talk between the thyroid and stress axes in a mammalian hippocampal model system.

Glucocorticoid-driven transcriptomes in human airway epithelial cells: commonalities, differences and functional insight from cell lines and primary cells

Background

Glucocorticoids act on the glucocorticoid receptor (GR; NR3C1) to resolve inflammation and, as inhaled corticosteroids (ICS), are the cornerstone of treatment for asthma. However, reduced efficacy in severe disease or exacerbations indicates a need to improve ICS actions.

Methods

Glucocorticoid-driven transcriptomes were compared using PrimeView microarrays between primary human bronchial epithelial (HBE) cells and the model cell lines, pulmonary type II A549 and bronchial epithelial BEAS-2B cells.

Results

In BEAS-2B cells, budesonide induced (≥2-fold, P ≤ 0.05) or, in a more delayed fashion, repressed (≤0.5-fold, P ≤ 0.05) the expression of 63, 133, 240, and 257 or 15, 56, 236, and 344 mRNAs at 1, 2, 6, and 18 h, respectively. Within the early-induced mRNAs were multiple transcriptional activators and repressors, thereby providing mechanisms for the subsequent modulation of gene expression. Using the above criteria, 17 (BCL6, BIRC3, CEBPD, ERRFI1, FBXL16, FKBP5, GADD45B, IRS2, KLF9, PDK4, PER1, RGCC, RGS2, SEC14L2, SLC16A12, TFCP2L1, TSC22D3) induced and 8 (ARL4C, FLRT2, IER3, IL11, PLAUR, SEMA3A, SLC4A7, SOX9) repressed mRNAs were common between A549, BEAS-2B and HBE cells at 6 h. As absolute gene expression change showed greater commonality, lowering the cut-off (≥1.25 or ≤ 0.8-fold) within these groups produced 93 induced and 82 repressed genes in common. Since large changes in few mRNAs and/or small changes in many mRNAs may drive function, gene ontology (GO)/pathway analyses were performed using both stringency criteria. Budesonide-induced genes showed GO term enrichment for positive and negative regulation of transcription, signaling, proliferation, apoptosis, and movement, as well as FOXO and PI3K-Akt signaling pathways. Repressed genes were enriched for inflammatory signaling pathways (TNF, NF-κB) and GO terms for cytokine activity, chemotaxis and cell signaling. Reduced growth factor expression and effects on proliferation and apoptosis were highlighted.

Conclusions

While glucocorticoids repress mRNAs associated with inflammation, prior induction of transcriptional activators and repressors may explain longer-term responses to these agents. Furthermore, positive and negative effects on signaling, proliferation, migration and apoptosis were revealed. Since many such gene expression changes occurred in human airways post-ICS inhalation, the effects observed in cell lines and primary HBE cells in vitro may be relevant to ICS in vivo.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0467-2) contains supplementary material, which is available to authorized users.

Serine Threonine Kinase Receptor-Associated Protein Deficiency Impairs Mouse Embryonic Stem Cells Lineage Commitment Through CYP26A1-Mediated Retinoic Acid Homeostasis

Retinoic acid (RA) signaling is essential for the differentiation of embryonic stem cells (ESCs) and vertebrate development. RA biosynthesis and metabolism are controlled by a series of enzymes, but the molecular regulators of these enzymes remain largely obscure. In this study, we investigated the functional role of the WD-domain protein STRAP (serine threonine kinase receptor-associated protein) in the pluripotency and lineage commitment of murine ESCs. We generated Strap knockout (KO) mouse ESCs and subjected them to spontaneous differentiation. We observed that, despite the unchanged characteristics of ESCs, Strap KO ESCs exhibited defects for lineage differentiation. Signature gene expression analyses revealed that Strap deletion attenuated intracellular RA signaling in embryoid bodies (EBs), and exogenous RA significantly rescued this deficiency. Moreover, loss of Strap selectively induced Cyp26A1 expression in mouse EBs, suggesting a potential role of STRAP in RA signaling. Mechanistically, we identified putative Krüppel-like factor 9 (KLF9) binding motifs to be critical in the enhancement of non-canonical RA-induced transactivation of Cyp26A1. Increased KLF9 expression in the absence of STRAP is partially responsible for Cyp26A1 induction. Interestingly, STRAP knockdown in Xenopus embryos influenced anterior-posterior neural patterning and impaired the body axis and eye development during early Xenopus embryogenesis. Taken together, our study reveals an intrinsic role for STRAP in the regulation of RA signaling and provides new molecular insights for ESC fate determination. Stem Cells 2018;36:1368-1379.

Glucocorticoids, genes and brain function

The identification of key genes in transcriptomic data constitutes a huge challenge. Our review of microarray reports revealed 88 genes whose transcription is consistently regulated by glucocorticoids (GCs), such as cortisol, corticosterone and dexamethasone, in the brain. Replicable transcriptomic data were combined with biochemical and physiological data to create an integrated view of the effects induced by GCs. The most frequently reported genes were Errfi1 and Ddit4. Their up-regulation was associated with the altered transcription of genes regulating growth factor and mTORC1 signaling (Gab1, Tsc22d3, Dusp1, Ndrg2, Ppp5c and Sesn1) and progression of the cell cycle (Ccnd1, Cdkn1a and Cables1). The GC-induced reprogramming of cell function involves changes in the mRNA level of genes responsible for the regulation of transcription (Klf9, Bcl6, Klf15, Tle3, Cxxc5, Litaf, Tle4, Jun, Sox4, Sox2, Sox9, Irf1, Sall2, Nfkbia and Id1) and the selective degradation of mRNA (Tob2). Other genes are involved in the regulation of metabolism (Gpd1, Aldoc and Pdk4), actin cytoskeleton (Myh2, Nedd9, Mical2, Rhou, Arl4d, Osbpl3, Arhgef3, Sdc4, Rdx, Wipf3, Chst1 and Hepacam), autophagy (Eva1a and Plekhf1), vesicular transport (Rhob, Ehd3, Vps37b and Scamp2), gap junctions (Gjb6), immune response (Tiparp, Mertk, Lyve1 and Il6r), signaling mediated by thyroid hormones (Thra and Sult1a1), calcium (Calm2), adrenaline/noradrenaline (Adcy9 and Adra1d), neuropeptide Y (Npy1r) and histamine (Hdc). GCs also affected genes involved in the synthesis of polyamines (Azin1) and taurine (Cdo1). The actions of GCs are restrained by feedback mechanisms depending on the transcription of Sgk1, Fkbp5 and Nr3c1. A side effect induced by GCs is increased production of reactive oxygen species. Available data show that the brain's response to GCs is part of an emergency mode characterized by inactivation of non-core activities, restrained inflammation, restriction of investments (growth), improved efficiency of energy production and the removal of unnecessary or malfunctioning cellular components to conserve energy and maintain nutrient supply during the stress response.