PLZF/ZBTB16, a glucocorticoid response gene in acute lymphoblastic leukemia, interferes with glucocorticoid-induced apoptosis

Enhanced ZBTB16 Levels by Progestin-Only Contraceptives Induces Decidualization and Inflammation

Progestin-only long-acting reversible-contraceptive (pLARC)-exposed endometria displays decidualized human endometrial stromal cells (HESCs) and hyperdilated thin-walled fragile microvessels. The combination of fragile microvessels and enhanced tissue factor levels in decidualized HESCs generates excess thrombin, which contributes to abnormal uterine bleeding (AUB) by inducing inflammation, aberrant angiogenesis, and proteolysis. The- zinc finger and BTB domain containing 16 (ZBTB16) has been reported as an essential regulator of decidualization. Microarray studies have demonstrated that ZBTB16 levels are induced by medroxyprogesterone acetate (MPA) and etonogestrel (ETO) in cultured HESCs. We hypothesized that pLARC-induced ZBTB16 expression contributes to HESC decidualization, whereas prolonged enhancement of ZBTB16 levels triggers an inflammatory milieu by inducing pro-inflammatory gene expression and tissue-factor-mediated thrombin generation in decidualized HESCs. Thus, ZBTB16 immunostaining was performed in paired endometria from pre- and post-depo-MPA (DMPA)-administrated women and oophorectomized guinea pigs exposed to the vehicle, estradiol (E2), MPA, or E2 + MPA. The effect of progestins including MPA, ETO, and levonorgestrel (LNG) and estradiol + MPA + cyclic-AMP (E2 + MPA + cAMP) on ZBTB16 levels were measured in HESC cultures by qPCR and immunoblotting. The regulation of ZBTB16 levels by MPA was evaluated in glucocorticoid-receptor-silenced HESC cultures. ZBTB16 was overexpressed in cultured HESCs for 72 h followed by a ± 1 IU/mL thrombin treatment for 6 h. DMPA administration in women and MPA treatment in guinea pigs enhanced ZBTB16 immunostaining in endometrial stromal and glandular epithelial cells. The in vitro findings indicated that: (1) ZBTB16 levels were significantly elevated by all progestin treatments; (2) MPA exerted the greatest effect on ZBTB16 levels; (3) MPA-induced ZBTB16 expression was inhibited in glucocorticoid-receptor-silenced HESCs. Moreover, ZBTB16 overexpression in HESCs significantly enhanced prolactin (PRL), insulin-like growth factor binding protein 1 (IGFBP1), and tissue factor (F3) levels. Thrombin-induced interleukin 8 (IL-8) and prostaglandin-endoperoxide synthase 2 (PTGS2) mRNA levels in control-vector-transfected HESCs were further increased by ZBTB16 overexpression. In conclusion, these results supported that ZBTB16 is enhanced during decidualization, and long-term induction of ZBTB16 expression by pLARCs contributes to thrombin generation through enhancing tissue factor expression and inflammation by enhancing IL-8 and PTGS2 levels in decidualized HESCs.

Super enhancers targeting ZBTB16 in osteogenesis protect against osteoporosis

As the major cell precursors in osteogenesis, mesenchymal stem cells (MSCs) are indispensable for bone homeostasis and development. However, the primary mechanisms regulating osteogenic differentiation are controversial. Composed of multiple constituent enhancers, super enhancers (SEs) are powerful cis-regulatory elements that identify genes that ensure sequential differentiation. The present study demonstrated that SEs were indispensable for MSC osteogenesis and involved in osteoporosis development. Through integrated analysis, we identified the most common SE-targeted and osteoporosis-related osteogenic gene, ZBTB16. ZBTB16, positively regulated by SEs, promoted MSC osteogenesis but was expressed at lower levels in osteoporosis. Mechanistically, SEs recruited bromodomain containing 4 (BRD4) at the site of ZBTB16, which then bound to RNA polymerase II-associated protein 2 (RPAP2) that transported RNA polymerase II (POL II) into the nucleus. The subsequent synergistic regulation of POL II carboxyterminal domain (CTD) phosphorylation by BRD4 and RPAP2 initiated ZBTB16 transcriptional elongation, which facilitated MSC osteogenesis via the key osteogenic transcription factor SP7. Bone-targeting ZBTB16 overexpression had a therapeutic effect on the decreased bone density and remodeling capacity of Brd4^fl/fl Prx1-cre mice and osteoporosis (OP) models. Therefore, our study shows that SEs orchestrate the osteogenesis of MSCs by targeting ZBTB16 expression, which provides an attractive focus and therapeutic target for osteoporosis. Without SEs located on osteogenic genes, BRD4 is not able to bind to osteogenic identity genes due to its closed structure before osteogenesis. During osteogenesis, histones on osteogenic identity genes are acetylated, and OB-gain SEs appear, enabling the binding of BRD4 to the osteogenic identity gene ZBTB16. RPAP2 transports RNA Pol II from the cytoplasm to the nucleus and guides Pol II to target ZBTB16 via recognition of the navigator BRD4 on SEs. After the binding of the RPAP2-Pol II complex to BRD4 on SEs, RPAP2 dephosphorylates Ser5 at the Pol II CTD to terminate the transcriptional pause, and BRD4 phosphorylates Ser2 at the Pol II CTD to initiate transcriptional elongation, which synergistically drives efficient transcription of ZBTB16, ensuring proper osteogenesis. Dysregulation of SE-mediated ZBTB16 expression leads to osteoporosis, and bone-targeting ZBTB16 overexpression is efficient in accelerating bone repair and treating osteoporosis.

Chronic stress-driven glucocorticoid receptor activation programs key cell phenotypes and functional epigenomic patterns in human fibroblasts

Chronic environmental stress can profoundly impact cell and body function. Although the underlying mechanisms are poorly understood, epigenetics has emerged as a key link between environment and health. The genomic effects of stress are thought to be mediated by the action of glucocorticoid stress hormones, primarily cortisol in humans, which act via the glucocorticoid receptor (GR). To dissect how chronic stress-driven GR activation influences epigenetic and cell states, human fibroblasts underwent prolonged exposure to physiological stress levels of cortisol and/or a selective GR antagonist. Cortisol was found to drive robust changes in cell proliferation, migration, and morphology, which were abrogated by concomitant GR blockade. The GR-driven cell phenotypes were accompanied by widespread, yet genomic context-dependent, changes in DNA methylation and mRNA expression, including gene loci with known roles in cell proliferation and migration. These findings provide insights into how chronic stress-driven functional epigenomic patterns become established to shape key cell phenotypes.

•

Physiological stress levels of cortisol drive robust changes in key cell phenotypes

•

Stress-driven changes in cell phenotypes are abrogated by concomitant GR blockade

•

GR activation induces functional and phenotypically relevant epigenomic changes

Short and long-term effect of dexamethasone on the transcriptome profile of primary human trabecular meshwork cells in vitro

In the quest of identifying newer molecular targets for the management of glucocorticoid-induced ocular hypertension (GC-OHT) and glaucoma (GCG), several microarray studies have attempted to investigate the genome-wide transcriptome profiling of primary human trabecular meshwork (TM) cells in response to dexamethasone (DEX). However, no studies are reported so far to demonstrate the temporal changes in the expression of genes in the cultured human TM cells in response to DEX treatment. Therefore, in the present study, the time-dependent changes in the genome-wide expression of genes in primary human TM cells after short (16 hours: 16 h) and long exposure (7 days: 7 d) of DEX was investigated using RNA sequencing. There were 199 (118 up-regulated; 81 down-regulated) and 525 (119 up-regulated; 406 down-regulated) DEGs in 16 h and 7 d treatment groups respectively. The unique genes identified in 16 h and 7 d treatment groups were 152 and 478 respectively. This study found a distinct gene signature and pathways between two treatment regimes. Longer exposure of DEX treatment showed a dys-regulation of Wnt and Rap1 signaling and so highlighted potential therapeutic targets for pharmacological management of GC-OHT/glaucoma.

The Flavonoid Baicalein Negatively Regulates Progesterone Target Genes in the Uterus in Vivo

Baicalein is a flavonoid extracted from the root of Scutellaria baicalensis (Chinese skullcap) and is consumed as part of this botanical dietary supplement to reduce oxidative stress, pain, and inflammation. We previously reported that baicalein can also modify receptor signaling through the progesterone receptor (PR) and glucocorticoid receptor (GR) in vitro, which is interesting due to the well-established roles of both PR and GR in reducing inflammation. To understand the effects of baicalein on PR and GR signaling in vivo in the uterus, ovariectomized CD-1 mice were treated with DMSO, progesterone (P4), baicalein, P4 with baicalein, and P4 with RU486, a PR antagonist, for a week. The uteri were collected for histology and RNA sequencing. Our results showed that baicalein attenuated the antiproliferative effect of P4 on luminal epithelium as well as on the PR target genes HAND2 and ZBTB16. Baicalein did not change levels of PR or GR RNA or protein in the uterus. RNA sequencing data indicated that many transcripts significantly altered by baicalein were regulated in the opposite direction by P4. Similarly, a large portion of GO/KEGG terms and GSEA gene sets were altered in the opposite direction by baicalein as compared to P4 treatment. Treatment of baicalein did not change body weight, organ weight, or blood glucose level. In summary, baicalein functioned as a PR antagonist in vivo and therefore may oppose P4 action under certain conditions such as uterine hyperplasia, fibroids, and uterine cancers.

Glucocorticoid signaling induces transcriptional memory and universally reversible chromatin changes

Glucocorticoids are stress hormones that elicit various cellular responses. These responses are typically reversible; however, in some instances, a previous exposure is “remembered” and influences the response to a subsequent hormone encounter.

Glucocorticoids are stress hormones that elicit cellular responses by binding to the glucocorticoid receptor, a ligand-activated transcription factor. The exposure of cells to this hormone induces wide-spread changes in the chromatin landscape and gene expression. Previous studies have suggested that some of these changes are reversible whereas others persist even when the hormone is no longer around. However, when we examined chromatin accessibility in human airway epithelial cells after hormone washout, we found that the hormone-induced changes were universally reversed after 1 d. Moreover, priming of cells by a previous exposure to hormone, in general, did not alter the transcriptional response to a subsequent encounter of the same cue except for one gene, ZBTB16, that displays transcriptional memory manifesting itself as a more robust transcriptional response upon repeated hormone stimulation. Single-cell analysis revealed that the more robust response is driven by a higher probability of primed cells to activate ZBTB16 and by a subset of cells that express the gene at levels that are higher than the induction levels observed for naïve cells.

Glucocorticoid signaling in pancreatic islets modulates gene regulatory programs and genetic risk of type 2 diabetes

Glucocorticoids are key regulators of glucose homeostasis and pancreatic islet function, but the gene regulatory programs driving responses to glucocorticoid signaling in islets and the contribution of these programs to diabetes risk are unknown. In this study we used ATAC-seq and RNA-seq to map chromatin accessibility and gene expression from eleven primary human islet samples cultured in vitro with the glucocorticoid dexamethasone at multiple doses and durations. We identified thousands of accessible chromatin sites and genes with significant changes in activity in response to glucocorticoids. Chromatin sites up-regulated in glucocorticoid signaling were prominently enriched for glucocorticoid receptor binding sites and up-regulated genes were enriched for ion transport and lipid metabolism, whereas down-regulated chromatin sites and genes were enriched for inflammatory, stress response and proliferative processes. Genetic variants associated with glucose levels and T2D risk were enriched in glucocorticoid-responsive chromatin sites, including fine-mapped variants at 51 known signals. Among fine-mapped variants in glucocorticoid-responsive chromatin, a likely casual variant at the 2p21 locus had glucocorticoid-dependent allelic effects on beta cell enhancer activity and affected SIX2 and SIX3 expression. Our results provide a comprehensive map of islet regulatory programs in response to glucocorticoids through which we uncover a role for islet glucocorticoid signaling in mediating genetic risk of T2D.

Differential DNA Methylation by Hispanic Ethnicity Among Firefighters in the United States

Firefighters are exposed to a variety of environmental hazards and are at increased risk for multiple cancers. There is evidence that risks differ by ethnicity, yet the biological or environmental differences underlying these differences are not known. DNA methylation is one type of epigenetic regulation that is altered in cancers. In this pilot study, we profiled DNA methylation with the Infinium MethylationEPIC in blood leukocytes from 31 Hispanic white and 163 non-Hispanic white firefighters. We compared DNA methylation (1) at 12 xenobiotic metabolizing genes and (2) at all loci on the array (>740 000), adjusting for confounders. Five of the xenobiotic metabolizing genes were differentially methylated at a raw P-value <.05 when comparing the 2 ethnic groups, yet were not statistically significant at a 5% false discovery rate (q-value <.05). In the epigenome-wide analysis, 76 loci exhibited DNA methylation differences at q < .05. Among these, 3 CpG sites in the promoter region of the biotransformation gene SULT1C2 had lower methylation in Hispanic compared to non-Hispanic firefighters. Other differentially methylated loci included genes that have been implicated in carcinogenesis in published studies (FOXK2, GYLTL1B, ZBTB16, ARHGEF10, and more). In this pilot study, we report differential DNA methylation between Hispanic and non-Hispanic firefighters in xenobiotic metabolism genes and other genes with functions related to cancer. Epigenetic susceptibility by ethnicity merits further study as this may alter risk for cancers linked to toxic exposures.

Identification of hub genes and small molecule therapeutic drugs related to breast cancer with comprehensive bioinformatics analysis

Breast cancer is one of the most common malignant tumors among women worldwide and has a high morbidity and mortality. This research aimed to identify hub genes and small molecule drugs for breast cancer by integrated bioinformatics analysis. After downloading multiple gene expression datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, 283 overlapping differentially expressed genes (DEGs) significantly enriched in different cancer-related functions and pathways were obtained using LIMMA, VennDiagram and ClusterProfiler packages of R. We then analyzed the topology of protein–protein interaction (PPI) network with overlapping DEGs and further obtained six hub genes (RRM2, CDC20, CCNB2, BUB1B, CDK1, and CCNA2) from the network via STRING and Cytoscape. Subsequently, we conducted genes expression verification, genetic alterations evaluation, immune infiltration prediction, clinicopathological parameters analysis, identification of transcriptional and post-transcriptional regulatory molecules, and survival analysis for these hub genes. Meanwhile, 29 possible drug candidates (e.g., Cladribine, Gallium nitrate, Alvocidib, 1β-hydroxyalantolactone, Berberine hydrochloride, Nitidine chloride) were identified from the DGIdb database and the GSE85871 dataset. In addition, some transcription factors and miRNAs (e.g., E2F1, PTTG1, TP53, ZBTB16, hsa-miR-130a-3p, hsa-miR-204-5p) targeting hub genes were identified as key regulators in the progression of breast cancer. In conclusion, our study identified six hub genes and 29 potential drug candidates for breast cancer. These findings may advance understanding regarding the diagnosis, prognosis and treatment of breast cancer.

BAL Cell Gene Expression in Severe Asthma Reveals Mechanisms of Severe Disease and Influences of Medications

Rationale: Gene expression of BAL cells, which samples the cellular milieu within the lower respiratory tract, has not been well studied in severe asthma.Objectives: To identify new biomolecular mechanisms underlying severe asthma by an unbiased, detailed interrogation of global gene expression.Methods: BAL cell expression was profiled in 154 asthma and control subjects. Of these participants, 100 had accompanying airway epithelial cell gene expression. BAL cell expression profiles were related to participant (age, sex, race, and medication) and sample traits (cell proportions), and then severity-related gene expression determined by correlating transcripts and coexpression networks to lung function, emergency department visits or hospitalizations in the last year, medication use, and quality-of-life scores.Measurements and Main Results: Age, sex, race, cell proportions, and medications strongly influenced BAL cell gene expression, but leading severity-related genes could be determined by carefully identifying and accounting for these influences. A BAL cell expression network enriched for cAMP signaling components most differentiated subjects with severe asthma from other subjects. Subsequently, an in vitro cellular model showed this phenomenon was likely caused by a robust upregulation in cAMP-related expression in nonsevere and β-agonist-naive subjects given a β-agonist before cell collection. Interestingly, ELISAs performed on BAL lysates showed protein levels may partly disagree with expression changes.Conclusions: Gene expression in BAL cells is influenced by factors seldomly considered. Notably, β-agonist exposure likely had a strong and immediate impact on cellular gene expression, which may not translate to important disease mechanisms or necessarily match protein levels. Leading severity-related genes were discovered in an unbiased, system-wide analysis, revealing new targets that map to asthma susceptibility loci.

Altered expression of dopaminergic cell fate regulating genes prior to manifestation of symptoms in a transgenic rat model of Huntington's disease

Hyperactivity of the dopaminergic pathway is thought to contribute to clinical symptoms in the early stages of Huntington's disease (HD). It is suggested to be result of a reduced dopaminergic inhibition by degeneration of medium spiny neurons in the striatum. Previously, we have shown that the number of dopaminergic cells is increased in the dorsal raphe nucleus (DRN) of HD patients and transgenic HD (tgHD) rats during the manifestation phase of the disease; as well as in the substantia nigra pars compacta (SNc) and ventral tegmental area (VTA) of tgHD rats. To address whether these changes are secondary to neurodegeneration or take place in the pre-manifest phase of the disease, we examined the expression of genes controlling neuronal cell fate and genes that define dopaminergic cell phenotype. In the SNc-VTA of tgHD rats, Msx1 was upregulated, which correlated with an altered expression of transcription factors Zbtb16 and Tcf12. Zbtb16 was upregulated in the DRN and it was the only gene that showed a correlated expression in the tgHD rats between SNc-VTA and DRN. Zbtb16 may be a candidate for regionally tuning its cell populations, resulting in the increase in dopaminergic cells observed in our previous studies. Here, we demonstrated an altered expression of genes related to dopaminergic cell fate regulation in the brainstem of 6 months-old tgHD rats. This suggests that changes in dopaminergic system in HD precede the manifestation of clinical symptoms, contradicting the theory that hyperdopaminergic status in HD is a consequence of neurodegeneration in the striatum.

Genetic and Epigenetic Perturbations by DNMT3A-R882 Mutants Impaired Apoptosis through Augmentation of PRDX2 in Myeloid Leukemia Cells

DNA methyltransferase 3A (DNMT3A) is mutated in various myeloid neoplasms including acute myeloid leukemia (AML), especially at the Arg882 and associated with inferior outcomes. Here, we report that the DNMT3A-Arg882His/Cys (R882H/C) mutations led to inactivation of apoptosis through DNA damage signaling following the impairment of differentiation of myeloid leukemia cells. Gene expression profiling analysis revealed aberrant expression of several cell-cycle and apoptosis-related genes, and the DNA methylation assay identified both hypo- and hypermethylation features in different regions throughout the whole genome of DNMT3A mutants-transduced myeloid cells. We found that DNMT3A-R882H/C mutations upregulated the expression of an antioxidant protein, pyroxiredoxin-2 (PRDX2), at the mRNA and protein levels with decreased accumulation of reactive oxygen species (ROS). Augmentation of ROS generation by ROS accumulating agent or by knockdown of PRDX2 from myeloid cells effectively increased drug sensitivity and apoptosis as a consequence of reduced cell proliferation. DNMT3A-R882C/H mutations decreased apoptosis induction in part by increasing the antioxidant capacity of the cell owing to upregulation of PRDX2. Molecularly, both DNMT3A-WT and R882H/C mutants interacted with PRDX2; and R882C/H mutation-induced hypomethylation increased PRDX2 expression which enhanced cell proliferation and growth with impairment of apoptosis, thereby contributing to leukemogenesis.

Cortisol limits selected actions of synthetic glucocorticoids in the airway epithelium

Cortisol, a physiologic glucocorticoid (GC), is essential for growth and differentiation of the airway epithelium. Epithelial function influences inflammation in chronic respiratory diseases. Synthetic GCs, including inhaled corticosteroids, exert anti-inflammatory effects in airway epithelium by transactivation of genes and by inhibition of proinflammatory cytokine release. We examined the effect of cortisol on the actions of synthetic GCs in the airway epithelium, demonstrating that cortisol acts like a partial agonist at the GC receptor (GR), limiting GC-induced GR-dependent transcription in the BEAS-2B human bronchial epithelial cell line. Cortisol also limited the inhibition of granulocyte macrophage colony-stimulating factor release by synthetic GCs in TNF-α-activated BEAS-2B cells. The relevance of these findings is supported by observations on tracheal epithelium obtained from mice treated for 5 d with systemic GC, showing limitations in selected GC effects, including inhibition of IL-6. Moreover, gene transactivation by synthetic GCs was compromised by standard air-liquid interface (ALI) growth medium cortisol concentration (1.4 μM) in the ALI-differentiated organotypic culture of primary human airway epithelial cells. These findings suggest that endogenous corticosteroids may limit certain actions of synthetic pharmacological GCs and contribute to GC insensitivity, particularly when corticosteroid levels are elevated by stress.-Prodanovic, D., Keenan, C. R., Langenbach, S., Li, M., Chen, Q., Lew, M. J., Stewart, A. G. Cortisol limits selected actions of synthetic glucocorticoids in the airway epithelium.

Gene expression changes in tumor free tongue tissue adjacent to tongue squamous cell carcinoma

Due to the high frequency of loco-regional recurrences, which could be explained by changes in the field surrounding the tumor, patients with squamous cell carcinoma of head and neck show poor survival. Here we identified a total of 554 genes as dysregulated in clinically tumor free tongue tissue in patients with tongue tumors when compared to healthy control tongue tissue. Among the top dysregulated genes when comparing control and tumor free tissue were those involved in apoptosis (CIDEC, MUC1, ZBTB16, PRNP, ECT2), immune response (IFI27) and differentiation (KRT36). Data suggest that these are important findings which can aid in earlier diagnosis of tumor development, a relapse or a novel squamous cell carcinoma of the tongue, in the absence of histological signs of a tumor.

Gene Expression Correlated with Severe Asthma Characteristics Reveals Heterogeneous Mechanisms of Severe Disease

RATIONALE

Severe asthma (SA) is a heterogeneous disease with multiple molecular mechanisms. Gene expression studies of bronchial epithelial cells in individuals with asthma have provided biological insight and underscored possible mechanistic differences between individuals.

OBJECTIVES

Identify networks of genes reflective of underlying biological processes that define SA.

METHODS

Airway epithelial cell gene expression from 155 subjects with asthma and healthy control subjects in the Severe Asthma Research Program was analyzed by weighted gene coexpression network analysis to identify gene networks and profiles associated with SA and its specific characteristics (i.e., pulmonary function tests, quality of life scores, urgent healthcare use, and steroid use), which potentially identified underlying biological processes. A linear model analysis confirmed these findings while adjusting for potential confounders.

MEASUREMENTS AND MAIN RESULTS

Weighted gene coexpression network analysis constructed 64 gene network modules, including modules corresponding to T1 and T2 inflammation, neuronal function, cilia, epithelial growth, and repair mechanisms. Although no network selectively identified SA, genes in modules linked to epithelial growth and repair and neuronal function were markedly decreased in SA. Several hub genes of the epithelial growth and repair module were found located at the 17q12-21 locus, near a well-known asthma susceptibility locus. T2 genes increased with severity in those treated with corticosteroids but were also elevated in untreated, mild-to-moderate disease compared with healthy control subjects. T1 inflammation, especially when associated with increased T2 gene expression, was elevated in a subgroup of younger patients with SA.

CONCLUSIONS

In this hypothesis-generating analysis, gene expression networks in relation to asthma severity provided potentially new insight into biological mechanisms associated with the development of SA and its phenotypes.

Gene Expression Correlated with Severe Asthma Characteristics Reveals Heterogeneous Mechanisms of Severe Disease

RATIONALE

Severe asthma (SA) is a heterogeneous disease with multiple molecular mechanisms. Gene expression studies of bronchial epithelial cells in individuals with asthma have provided biological insight and underscored possible mechanistic differences between individuals.

OBJECTIVES

Identify networks of genes reflective of underlying biological processes that define SA.

METHODS

Airway epithelial cell gene expression from 155 subjects with asthma and healthy control subjects in the Severe Asthma Research Program was analyzed by weighted gene coexpression network analysis to identify gene networks and profiles associated with SA and its specific characteristics (i.e., pulmonary function tests, quality of life scores, urgent healthcare use, and steroid use), which potentially identified underlying biological processes. A linear model analysis confirmed these findings while adjusting for potential confounders.

MEASUREMENTS AND MAIN RESULTS

Weighted gene coexpression network analysis constructed 64 gene network modules, including modules corresponding to T1 and T2 inflammation, neuronal function, cilia, epithelial growth, and repair mechanisms. Although no network selectively identified SA, genes in modules linked to epithelial growth and repair and neuronal function were markedly decreased in SA. Several hub genes of the epithelial growth and repair module were found located at the 17q12-21 locus, near a well-known asthma susceptibility locus. T2 genes increased with severity in those treated with corticosteroids but were also elevated in untreated, mild-to-moderate disease compared with healthy control subjects. T1 inflammation, especially when associated with increased T2 gene expression, was elevated in a subgroup of younger patients with SA.

CONCLUSIONS

In this hypothesis-generating analysis, gene expression networks in relation to asthma severity provided potentially new insight into biological mechanisms associated with the development of SA and its phenotypes.

Blood transcriptome based biomarkers for human circadian phase

Diagnosis and treatment of circadian rhythm sleep-wake disorders both require assessment of circadian phase of the brain’s circadian pacemaker. The gold-standard univariate method is based on collection of a 24-hr time series of plasma melatonin, a suprachiasmatic nucleus-driven pineal hormone. We developed and validated a multivariate whole-blood mRNA-based predictor of melatonin phase which requires few samples. Transcriptome data were collected under normal, sleep-deprivation and abnormal sleep-timing conditions to assess robustness of the predictor. Partial least square regression (PLSR), applied to the transcriptome, identified a set of 100 biomarkers primarily related to glucocorticoid signaling and immune function. Validation showed that PLSR-based predictors outperform published blood-derived circadian phase predictors. When given one sample as input, the R² of predicted vs observed phase was 0.74, whereas for two samples taken 12 hr apart, R² was 0.90. This blood transcriptome-based model enables assessment of circadian phase from a few samples.

DOI:http://dx.doi.org/10.7554/eLife.20214.001

Genome-Wide Association between Transcription Factor Expression and Chromatin Accessibility Reveals Regulators of Chromatin Accessibility

To better understand genome regulation, it is important to uncover the role of transcription factors in the process of chromatin structure establishment and maintenance. Here we present a data-driven approach to systematically characterise transcription factors that are relevant for this process. Our method uses a linear mixed modelling approach to combine datasets of transcription factor binding motif enrichments in open chromatin and gene expression across the same set of cell lines. Applying this approach to the ENCODE dataset, we confirm already known and imply numerous novel transcription factors that play a role in the establishment or maintenance of open chromatin. In particular, our approach rediscovers many factors that have been annotated as pioneer factors.

Transcription factor binding occurs mainly in regions of open chromatin. For many transcription factors, it is unclear whether binding is the cause or the consequence of open chromatin. Here, we used datasets on open chromatin and gene expression provided by the ENCODE project to predict which transcription factors drive transitions between open and closed states. A signature of such a factor is that its expression values are correlated to chromatin accessibility at its motif across the same set of cell lines. Our method assesses this correlation while accounting for the fact that some tested cell lines are more related than others. We find many transcription factors showing evidence of driving transitions and a high proportion of these transcription factors are known pioneer factors, i.e., they play a role in opening up closed chromatin.

Regulation of Bim in Health and Disease

The BH3-only Bim protein is a major determinant for initiating the intrinsic apoptotic pathway under both physiological and pathophysiological conditions. Tight regulation of its expression and activity at the transcriptional, translational and post-translational levels together with the induction of alternatively spliced isoforms with different pro-apoptotic potential, ensure timely activation of Bim. Under physiological conditions, Bim is essential for shaping immune responses where its absence promotes autoimmunity, while too early Bim induction eliminates cytotoxic T cells prematurely, resulting in chronic inflammation and tumor progression. Enhanced Bim induction in neurons causes neurodegenerative disorders including Alzheimer's, Parkinson's and Huntington's diseases. Moreover, type I diabetes is promoted by genetically predisposed elevation of Bim in β-cells. On the contrary, cancer cells have developed mechanisms that suppress Bim expression necessary for tumor progression and metastasis. This review focuses on the intricate network regulating Bim activity and its involvement in physiological and pathophysiological processes.

ZBTB16 induces osteogenic differentiation marker genes in dental follicle cells independent from RUNX2

BACKGROUND

Dental follicle cells (DFCs) are neural crest cell-derived cells and the genuine precursor cells of cementoblast and alveolar osteoblasts. After osteogenic differentiation, expression levels of the transcription factor zinc factor and BTB domain containing 16 (ZBTB16) were significantly increased. ZBTB16 is associated with the process of osteogenic differentiation in bone marrow-derived mesenchymal stem cells and crucial for the expression of the osteogenic transcription factor runt-related transcription factor 2 (RUNX2). It is proposed that ZBTB16 plays also a crucial role for the differentiation of DFCs into osteoblasts.

METHODS

In this study, the differentiation of DFCs by alkaline phosphatase (ALP) activity measurement, alizarin red staining, and electron-dispersive x-ray spectrometry (EDX) analysis is investigated. The expression of ZBTB16 during osteogenic differentiation and the expression of osteogenic differentiation markers were assessed by real-time reverse transcription polymerase chain reaction. Glucocorticoid stimulation was inhibited using RU486 (11β-[p-(Dimethylamino)phenyl]-17β-hydroxy-17-(1-propynyl)estra-4,9-dien-3-one), and ZBTB16 was overexpressed via transient transfection of an expression vector.

RESULTS

After the initiation of osteogenic differentiation, ZBTB16 levels were increased highly in DFCs, whereas RUNX2 was expressed constitutively only. An EDX analysis verified the differentiation of DFCs into osteoblast-like cells because clusters of mineralization consisted of hydroxyapatite. ZBTB16 induced the expression of nuclear receptor subfamily 4, group A, member 3; osteocalcin; and stanniocalcin 1 (STC1) but not of RUNX2 and ALP in DFCs. STC1 was upregulated in DFCs downstream of ZBTB16 and after the osteogenic differentiation. The overexpression of STC1 in DFCs increased the expression of ZBTB16 and specific markers for biomineralization.

CONCLUSIONS

The present study shows that ZBTB16 induced the expression of osteogenic differentiation markers independently of RUNX2. Moreover, STC1 is a new candidate for the evaluation of late mechanisms of osteogenic differentiation downstream of ZBTB16.

Problems encountered in bicistronic IRES-GFP expression vectors employed in functional analyses of GC-induced genes

Our laboratory has developed a series of Gateway(®) compatible lentiviral expression systems for constitutive and conditional gene knock-down and over-expression. For tetracycline-regulated transgenic expression, we constructed a lentiviral "DEST" plasmid (pHR-TetCMV-Dest-IRES-GFP5) containing a tetracycline-responsive minimal CMV promoter, followed by an attP site-flanked DEST cassette (for efficient cloning of cDNAs by "Gateway(®)" recombination cloning) and green fluorescent protein (GFP) driven by an internal ribosomal entry site (IRES).This lentiviral bicistronic plasmid allows immediate FACS identification and characterization of successfully transfected cell lines. Although this system worked well with several cDNAs, we experienced serious problems with SLA, Bam and BMF. Particularly, we cloned the cDNA for human SLA (Src-like adapter), a candidate gene in GC-induced apoptosis, into this plasmid. The resulting construct (pHR-TetCMV-SLA-IRES-GFP5) was transfected into HEK 293-T packaging cells to produce viral particles for transduction of CEM-C7H2-2C8 cells. Although the construct produced many green fluorescent colonies at the HEK 293-T and the CEM-C7H2-2C8 level, we could not detect any SLA protein with α-SLA antibody from corresponding cell lysates. In contrast, the antibody readily detected SLA in whole cell lysate of HEK 293-T cells transfected with a GST-flagged SLA construct lacking IRES-GFP. To directly address the potential role of the IRES-GFP sequence, we cloned the SLA coding region into pHR-TetCMV-Dest, a vector that differs from pHR-TetCMV-Dest-IRES-GFP5 just by the absence of the IRES-GFP cassette. The resulting pHR-TetCMV-SLA construct was used for transfection of HEK 293-T cells. Corresponding lysates were assayed with α-SLA antibody and found positive. These data, in concert with previous findings, suggest that the IRES-GFP cassette may interfere with translation of certain smaller size cDNAs (like SLA) or generate fusion proteins and entail defective virus production in an unpredictable manner.

Hepatic gene expression analysis of mice exposed to raw water from Meiliang Bay, Lake Taihu, China

Lake water is a micro-polluted water system, and characterization of its toxicity remains difficult. Microarray-based determination of altered gene expression might be an alterative approach. We chose raw water from Meiliang Bay, Lake Taihu, China as the target water. Male mice were exposed to the lake water for 90 days. Total hepatic RNA was applied to interrogate the Affymetrix Mouse Genome 430A 2.0 array. Gene ontology analysis, pathway analysis and gene network analysis were used to identify biological effects of differently expressed genes. The results showed that the expressions of 170 genes were altered. Nine biological processes and nine biological pathways were significantly perturbed (P ≤ 0.01), mainly linked to the regulation of cell processes, DNA repair, chromatin modification, oxidative reduction and carbohydrate metabolism. Important genes, such as Prkca, Pik3r1, Fgfr1 and Zbtb16, were identified by gene network analysis. This study provided excellent insights into early toxicological effects related to raw Lake Taihu water, and illustrated that the toxicogenomic approach might be a useful tool to evaluate potential environmental health effects of raw lake water.

Evaluation of the NOD/SCID xenograft model for glucocorticoid-regulated gene expression in childhood B-cell precursor acute lymphoblastic leukemia

Background

Glucocorticoids such as prednisolone and dexamethasone are critical drugs used in multi-agent chemotherapy protocols used to treat acute lymphoblastic leukemia (ALL), and response to glucocorticoids is highly predictive of outcome. The NOD/SCID xenograft mouse model of ALL is a clinically relevant model in which the mice develop a systemic leukemia which retains the fundamental biological characteristics of the original disease. Here we report a study evaluating the NOD/SCID xenograft mouse model to investigate glucocorticoid-induced gene expression. Cells from a glucocorticoid-sensitive xenograft derived from a child with B-cell precursor ALL were inoculated into NOD/SCID mice. When highly engrafted the mice were randomized into groups of 4 to receive dexamethasone 15 mg/kg by intraperitoneal injection or vehicle control. Leukemia cells were harvested from mice spleens at 0, 8, 24 or 48 hours thereafter, and gene expression analyzed on Illumina WG-6_V3 chips, comparing all groups to time 0 hours.

Results

The 8 hour dexamethasone-treated timepoint had the highest number of significantly differentially expressed genes, with fewer observed at the 24 and 48 hour timepoints, and with minimal changes seen across the time-matched controls. When compared to publicly available datasets of glucocorticoid-induced gene expression from an in vitro cell line study and from an in vivo study of patients with ALL, at the level of pathways, expression changes in the 8 hour xenograft samples showed a similar response to patients treated with glucocorticoids. Replicate analysis revealed that at the 8 hour timepoint, a dataset with high signal and differential expression, using data from 3 replicates instead of 4 resulted in excellent recovery scores of > 0.9. However at other timepoints with less signal very poor recovery scores were obtained with 3 replicates.

Conclusions

The NOD/SCID xenograft mouse model provides a reproducible experimental system in which to investigate clinically-relevant mechanisms of drug-induced gene regulation in ALL; the 8 hour timepoint provides the highest number of significantly differentially expressed genes; time-matched controls are redundant and excellent recovery scores can be obtained with 3 replicates.

Research resource: transcriptional response to glucocorticoids in childhood acute lymphoblastic leukemia

Glucocorticoids (GC) induce apoptosis in lymphoblasts and are thus essential in the treatment of acute lymphoblastic leukemia (ALL). Their effects result from gene regulations via the GC receptor (NR3C1/GR), but it is unknown how these changes evolve, what the primary GR targets are, and to what extent responses differ between ALL subtypes and nonlymphoid malignancies. We delineated the transcriptional response to GC on the exon level in a time-resolved manner in a precursor B- and a T childhood ALL model employing Exon microarrays and combined this with genome-wide NR3C1-binding site detection using chromatin immunoprecipitation-on-chip technology. This integrative approach showed that the response was strongly influenced by kinetics and extent of GR autoinduction in both models. Although remarkable differences between the ALL systems were apparent, we defined a set of common response genes enriched in apoptosis-related processes. Globally, GR binding was higher for GC-induced vs. -repressed genes, suggesting that GR mediates gene repression by interaction with distant enhancers or by cross talk with other transcription factors. Exon level analysis defined several new GC-regulated transcript variants of genes, including ATP4B, GPR98, TBCD, and ZBTB16. Our study provides unprecedented insight into the transcriptional response to GC in ALL cells, essential to understand this biologically and clinically important phenomenon. We found evidence of cell type-specific as well as common responses, possibly related to apoptosis induction, and detected induction of novel transcript variants by GC in the investigated systems. Finally, we implemented a bioinformatic framework that might be useful for high-density microarray analyses to identify alternative transcript variant expression.

Promyelocytic leukemia zinc finger turns on the effector T cell program without requirement for agonist TCR signaling

Thymocytes expressing the NKT cell semi-invariant αβ TCR are thought to undergo agonist interactions with CD1d ligands prior to expressing promyelocytic leukemia zinc finger (PLZF), a broad complex, tramtrack, bric-a-brac, poxvirus, and zinc finger transcription factor that directs acquisition of the effector program of these innate-like T cells. Whether PLZF can mediate this effector conversion independently of agonist signaling has not been investigated. We demonstrated that transgenic (Tg) expression of PLZF under the CD4 promoter induced the innate effector program in two different MHC class II-restricted TCR-Tg Rag1(-/-) models examined. In CD4 thymocytes expressing a fixed Tg TCR β-chain, the associated TCRα sequences in wild-type and PLZF-Tg mice overlapped extensively, further demonstrating that PLZF could induce the effector program in most CD4 T cells that would normally be selected as naive cells. In contrast, PLZF altered the negative selection of thymocytes expressing TCR β-chains reactive against several retroviral superantigens. Thus, PLZF is remarkable in that it is a transcription factor capable of inducing an effector program in the absence of T cell agonist interactions or cell division. Its expression may also enhance the survival of agonist-signaled thymocytes.

Protective effect of dexamethasone against hypoxia-induced disruption of barrier function in human corneal epithelial cells

The corneal epithelium functions as a barrier to protect the cornea from external agents such as infectious organisms and toxins and thereby contributes to corneal homeostasis. The barrier function of epithelia is dependent on the formation of tight and adherens junctions between adjacent epithelial cells. We have previously shown that hypoxia disrupts the barrier function of cultured human corneal epithelial (HCE) cells by affecting tight junctions. We have now examined the effect of dexamethasone on this barrier disruption induced by hypoxia in HCE cells. Measurement of transepithelial electrical resistance revealed that the hypoxia-induced decrease in the barrier function of HCE cells was inhibited by dexamethasone in a concentration-dependent manner. The hypoxia-induced loss of the tight junction protein ZO-1 from the borders of adjacent HCE cells (as revealed by immunofluorescence analysis) as well as the hypoxia-induced down-regulation of ZO-1 expression (as revealed by immunoblot analysis) were also inhibited by dexamethasone, whereas this drug had no effect on the expression or distribution of the tight junction protein occludin or of the adherens junction proteins E-cadherin and β-catenin. Moreover, dexamethasone attenuated the reorganization of the actin cytoskeleton, the formation of focal adhesions, and the up-regulation of myosin light chain kinase expression induced by hypoxia in HCE cells. Our results thus suggest that dexamethasone protects corneal epithelial cells from the hypoxia-induced disruption of barrier function by maintaining the distribution and expression of ZO-1 as well as the organization of the actin cytoskeleton.

A corticosteroid-responsive transcription factor, promyelocytic leukemia zinc finger protein, mediates protection of the cochlea from acoustic trauma

Animals can be induced to resist cochlear damage associated with acoustic trauma by exposure to a variety of "conditioning" stimuli, including restraint stress, moderate level sound, heat stress, hypoxia, and corticosteroids. Here we identify in mice a corticosteroid-responsive transcription factor, PLZF (promyelocytic leukemia zinc finger protein), which mediates conditioned protection of the cochlea from acoustic trauma. PLZF mRNA levels in the cochlea are increased following conditioning stimuli, including restraint stress, dexamethasone administration, and moderate-to-high level acoustic stimulation. Heterozygous mutant (luxoid.Zbtb16(LU)/J) mice deficient in PLZF have hearing and responses to acoustic trauma similar to their wild type littermates but are unable to generate conditioning-induced protection from acoustic trauma. PLZF immunoreactivity is present in the spiral ganglion, lateral wall of the cochlea, and organ of Corti, all targets for acoustic trauma. PLZF is also present in the brain and PLZF mRNA in brain is elevated following conditioning stimuli. The identification of a transcription factor that mediates conditioned protection from trauma provides a tool for understanding the protective action of corticosteroids, which are widely used in treating acute hearing loss, and has relevance to understanding the role of corticosteroids in trauma protection.