Hypersensitive Sites Research Articles

Principles of nucleosome organization revealed by single-cell micrococcal nuclease sequencing.

Nucleosome positioning is critical to chromatin accessibility and is associated with gene expression programs in cells1-3. Previous nucleosome mapping methods assemble profiles from cell populations and reveal a cell-averaged pattern: nucleosomes are positioned and form a phased array that surrounds the transcription start sites of active genes3-6 and DNase I hypersensitive sites7. However, even in a homogenous population of cells, cells exhibit heterogeneity in expression in response to active signalling8,9 that may be related to heterogeneity in chromatin accessibility10-12. Here we report a technique, termed single-cell micrococcal nuclease sequencing (scMNase-seq), that can be used to simultaneously measure genome-wide nucleosome positioning and chromatin accessibility in single cells. Application of scMNase-seq to NIH3T3 cells, mouse primary naive CD4 T cells and mouse embryonic stem cells reveals two principles of nucleosome organization: first, nucleosomes in heterochromatin regions, or that surround the transcription start sites of silent genes, show large variation in positioning across different cells but are highly uniformly spaced along the nucleosome array; and second, nucleosomes that surround the transcription start sites of active genes and DNase I hypersensitive sites show little variation in positioning across different cells but are relatively heterogeneously spaced along the nucleosome array. We found a bimodal distribution of nucleosome spacing at DNase I hypersensitive sites, which corresponds to inaccessible and accessible states and is associated with nucleosome variation and variation in accessibility across cells. Nucleosome variation is smaller within single cells than across cells, and smaller within the same cell type than across cell types. A large fraction of naive CD4 T cells and mouse embryonic stem cells shows depleted nucleosome occupancy at the de novo enhancers detected in their respective differentiated lineages, revealing the existence of cells primed for differentiation to specific lineages in undifferentiated cell populations.

Genome-wide analysis of DNA methylation in buccal cells: a study of monozygotic twins and mQTLs

BackgroundDNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA.ResultsGood-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1–10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10−16) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10−16). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3.ConclusionsWe conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.

Prospective Isolation and Characterization of Genetically and Functionally Distinct AML Subclones

Intra-tumor heterogeneity caused by clonal evolution is a major problem in cancer treatment. To address this problem, we performed label-free quantitative proteomics on primary acute myeloid leukemia (AML) samples. We identified 50 leukemia-enriched plasma membrane proteins enabling the prospective isolation of genetically distinct subclones from individual AML patients. Subclones differed in their regulatory phenotype, drug sensitivity, growth, and engraftment behavior, as determined by RNA sequencing, DNase I hypersensitive site mapping, transcription factor occupancy analysis, invitro culture, and xenograft transplantation. Finally, we show that these markers can be used to identify and longitudinally track distinct leukemic clones in patients in routine diagnostics. Our study describes a strategy for a major improvement in stratifying cancer diagnosis and treatment.

CoTRaCTE predicts co-occurring transcription factors within cell-type specific enhancers

Cell-type specific gene expression is regulated by the combinatorial action of transcription factors (TFs). In this study, we predict transcription factor (TF) combinations that cooperatively bind in a cell-type specific manner. We first divide DNase hypersensitive sites into cell-type specifically open vs. ubiquitously open sites in 64 cell types to describe possible cell-type specific enhancers. Based on the pattern contrast between these two groups of sequences we develop “co-occurring TF predictor on Cell-Type specific Enhancers” (coTRaCTE) - a novel statistical method to determine regulatory TF co-occurrences. Contrasting the co-binding of TF pairs between cell-type specific and ubiquitously open chromatin guarantees the high cell-type specificity of the predictions. coTRaCTE predicts more than 2000 co-occurring TF pairs in 64 cell types. The large majority (70%) of these TF pairs is highly cell-type specific and overlaps in TF pair co-occurrence are highly consistent among related cell types. Furthermore, independently validated co-occurring and directly interacting TFs are significantly enriched in our predictions. Focusing on the regulatory network derived from the predicted co-occurring TF pairs in embryonic stem cells (ESCs) we find that it consists of three subnetworks with distinct functions: maintenance of pluripotency governed by OCT4, SOX2 and NANOG, regulation of early development governed by KLF4, STAT3, ZIC3 and ZNF148 and general functions governed by MYC, TCF3 and YY1. In summary, coTRaCTE predicts highly cell-type specific co-occurring TFs which reveal new insights into transcriptional regulatory mechanisms.

C3D: a tool to predict 3D genomic interactions between cis-regulatory elements.

The 3D genome architecture influences the regulation of genes by facilitating chromatin interactions between distal cis-regulatory elements and gene promoters. We implement Cross Cell-type Correlation based on DNA accessibility (C3D), a customizable computational tool that predicts chromatin interactions using an unsupervised algorithm that utilizes correlations in chromatin measurements, such as DNaseI hypersensitivity signals. C3D accurately predicts 32.7%, 18.3% and 24.1% of interactions, validated by ChIA-PET assays, between promoters and distal regions that overlie DNaseI hypersensitive sites in K562, MCF-7 and GM12878 cells, respectively. Source code is open-source and freely available on GitHub (https://github.com/LupienLabOrganization/C3D) under the GNU GPLv3 license. C3D is implemented in Bash and R; it runs on any platform with Bash (≥4.0), R (≥3.1.1) and BEDTools (≥2.19.0). It requires the following R packages: GenomicRanges, Sushi, data.table, preprocessCore and dynamicTreeCut. Supplementary data are available at Bioinformatics online.

DNase-seq to Study Chromatin Accessibility in Early Xenopus tropicalis Embryos.

Transcriptional regulatory elements are typically found in relatively nucleosome-free genomic regions, often referred to as "open chromatin." Deoxyribonuclease I (DNase I) can digest nucleosome-depleted DNA (presumably bound by transcription factors), but DNA in nucleosomes or higher-order chromatin fibers is less accessible to the nuclease. The DNase-seq method uses high-throughput sequencing to permit the interrogation of DNase hypersensitive sites (DHSs) across the entire genome and does not require prior knowledge of histone modifications, transcription factor binding sites, or high quality antibodies to identify potentially active regions of chromatin. Here, discontinuous iodixanol gradients are used as a gentle preparation of the nuclei from Xenopus embryos. Short DNase I digestion times are followed by size selection of digested genomic DNA, yielding DHS fragments. These DNA fragments are subjected to real-time quantitative polymerase chain reaction (qPCR) and sequencing library construction. A library generation method and pipeline for analyzing DNase-seq data are also described.

The bithorax complex iab-7 Polycomb response element has a novel role in the functioning of the Fab-7 chromatin boundary.

Expression of the three bithorax complex homeotic genes is orchestrated by nine parasegment-specific regulatory domains. Autonomy of each domain is conferred by boundary elements (insulators). Here, we have used an in situ replacement strategy to reanalyze the sequences required for the functioning of one of the best-characterized fly boundaries, Fab-7. It was initially identified by a deletion, Fab-71, that transformed parasegment (PS) 11 into a duplicate copy of PS12. Fab-71 deleted four nuclease hypersensitive sites, HS*, HS1, HS2, and HS3, located between the iab-6 and iab-7 regulatory domains. Transgenic and P-element excision experiments mapped the boundary to HS*+HS1+HS2, while HS3 was shown to be the iab-7 Polycomb response element (PRE). Recent replacement experiments showed that HS1 is both necessary and sufficient for boundary activity when HS3 is also present in the replacement construct. Surprisingly, while HS1+HS3 combination has full boundary activity, we discovered that HS1 alone has only minimal function. Moreover, when combined with HS3, only the distal half of HS1, dHS1, is needed. A ~1,000 kD multiprotein complex containing the GAF protein, called the LBC, binds to the dHS1 sequence and we show that mutations in dHS1, that disrupt LBC binding in nuclear extracts, eliminate boundary activity and GAF binding in vivo. HS3 has binding sites for GAF and Pho proteins that are required for PRE silencing. In contrast, HS3 boundary activity only requires the GAF binding sites. LBC binding with HS3 in nuclear extracts, and GAF association in vivo, depend upon the HS3 GAF sites, but not the Pho sites. Consistent with a role for the LBC in HS3 boundary activity, the boundary function of the dHS1+HS3mPho combination is lost when the flies are heterozygous for a mutation in the GAF gene. Taken together, these results reveal a novel function for the iab-7 PREs in chromosome architecture.

An Optimized Lentiviral Vector Efficiently Corrects the Human Sickle Cell Disease Phenotype

Autologous transplantation of hematopoietic stem cells transduced with a lentiviral vector (LV) expressing an anti-sickling HBB variant is a potential treatment for sickle cell disease (SCD). With a clinical trial as our ultimate goal, we generated LV constructs containing an anti-sickling HBB transgene (HBBAS3), a minimal HBB promoter, and different combinations of DNase I hypersensitive sites (HSs) from the locus control region (LCR). Hematopoietic stem progenitor cells (HSPCs) from SCD patients were transduced with LVs containing either HS2 and HS3 (β-AS3) or HS2, HS3, and HS4 (β-AS3 HS4). The inclusion of the HS4 element drastically reduced vector titer and infectivity in HSPCs, with negligible improvement of transgene expression. Conversely, the LV containing only HS2 and HS3 was able to efficiently transduce SCD bone marrow and Plerixafor-mobilized HSPCs, with anti-sickling HBB representing up to ∼60% of the total HBB-like chains. The expression of the anti-sickling HBB and the reduced incorporation of the βS-chain in hemoglobin tetramers allowed up to 50% reduction in the frequency of RBC sickling under hypoxic conditions. Together, these results demonstrate the ability of a high-titer LV to express elevated levels of a potent anti-sickling HBB transgene ameliorating the SCD cell phenotype.

Characterization of the accessible genome in the human malaria parasite Plasmodium falciparum.

Human malaria is a devastating disease and a major cause of poverty in resource-limited countries. To develop and adapt within hosts Plasmodium falciparum undergoes drastic switches in gene expression. To identify regulatory regions in the parasite genome, we performed genome-wide profiling of chromatin accessibility in two culture-adapted isogenic subclones at four developmental stages during the intraerythrocytic cycle by using the Assay for Transposase-Accessible Chromatin by sequencing (ATAC-seq). Tn5 transposase hypersensitivity sites (THSSs) localize preferentially at transcriptional start sites (TSSs). Chromatin accessibility by ATAC-seq is predictive of active transcription and of the levels of histone marks H3K9ac and H3K4me3. Our assay allows the identification of novel regulatory regions including TSS and enhancer-like elements. We show that the dynamics in the accessible chromatin profile matches temporal transcription during development. Motif analysis of stage-specific ATAC-seq sites predicts the in vivo binding sites and function of multiple ApiAP2 transcription factors. At last, the alternative expression states of some clonally variant genes (CVGs), including eba, phist, var and clag genes, associate with a differential ATAC-seq signal at their promoters. Altogether, this study identifies genome-wide regulatory regions likely to play an essential function in the developmental transitions and in CVG expression in P. falciparum.

Novel long-range regulatory mechanisms controlling PKD2 gene expression

BackgroundCis-regulatory elements control gene expression over large distances through the formation of chromatin loops, which allow contact between enhancers and gene promoters. Alterations in cis-acting regulatory systems could be linked to human genetic diseases. Here, we analyse the spatial organization of a large region spanning the polycystic kidney disease 2 (PKD2) gene, one of the genes responsible of autosomal dominant polycystic kidney disease (ADPKD).ResultsBy using chromosome conformation capture carbon copy (5C) technology in primary human renal cyst epithelial cells, we identify novel contacts of the PKD2 promoter with chromatin regions, which display characteristics of regulatory elements. In parallel, by using functional analysis with a reporter assay, we demonstrate that three DNAse I hypersensitive sites regions are involved in the regulation of PKD2 gene expression.ConclusionsFinally, through alignment of CCCTC-binding factor (CTCF) sites, we suggest that these novel enhancer elements are brought to the PKD2 promoter by chromatin looping via the recruitment of CTCF.

Prenatal immune activation alters the adult neural epigenome but can be partly stabilised by a n-3 polyunsaturated fatty acid diet

An unstable epigenome is implicated in the pathophysiology of neurodevelopmental disorders such as schizophrenia and autism. This is important because the epigenome is potentially modifiable. We have previously reported that adult offspring exposed to maternal immune activation (MIA) prenatally have significant global DNA hypomethylation in the hypothalamus. However, what genes had altered methylation state, their functional effects on gene expression and whether these changes can be moderated, have not been addressed. In this study, we used next-generation sequencing (NGS) for methylome profiling in a MIA rodent model of neurodevelopmental disorders. We assessed whether differentially methylated regions (DMRs) affected the chromatin state by mapping known DNase I hypersensitivity sites (DHSs), and selected overlapping genes to confirm a functional effect of MIA on gene expression using qPCR. Finally, we tested whether methylation differences elicited by MIA could be limited by post-natal dietary (omega) n-3 polyunsaturated fatty acid (PUFA) supplementation. These experiments were conducted using hypothalamic brain tissue from 12-week-old offspring of mice injected with viral analogue PolyI:C on gestation day 9 of pregnancy or saline on gestation day 9. Half of the animals from each group were fed a diet enriched with n-3 PUFA from weaning (MIA group, n = 12 units, n = 39 mice; Control group, n = 12 units, n = 38 mice). The results confirmed our previous finding that adult offspring exposed to MIA prenatally had significant global DNA hypomethylation. Furthermore, genes linked to synaptic plasticity were over-represented among differentially methylated genes following MIA. More than 80% of MIA-induced hypomethylated sites, including those affecting chromatin state and MECP2 binding, were stabilised by the n-3 PUFA intervention. MIA resulted in increased expression of two of the ‘top five’ genes identified from an integrated analysis of DMRs, DHSs and MECP2 binding sites, namely Abat (t = 2.46, p < 0.02) and Gnas9 (t = 2.96, p < 0.01), although these changes were not stabilised by dietary intervention. Thus, prenatal MIA exposure impacts upon the epigenomic regulation of gene pathways linked to neurodevelopmental conditions; and many of the changes can be attenuated by a low-cost dietary intervention.

Abstract 3173: Lymphocyte-specific chromatin accessibility predetermines glucocorticoid resistance in acute lymphoblastic leukemia

Abstract Glucocorticoids are well known for their immunosuppressive activity, and play a critical role in the treatment of lymphoid malignancies. However, the development of resistance remains a significant barrier to cure and the mechanisms are poorly defined. Moreover, glucocorticoids are rarely efficacious in treating myeloid and other non-lymphoid malignancies, and the mechanisms of lymphocyte-specific efficacy are unclear. To address these issues, we first carried out a global analysis of DNase I hypersensitive sites in 18 lymphoid and 64 non-lymphoid cell types to map lymphocyte-specific open chromatin domains (LSOs). We then integrated these domains with genome-wide glucocorticoid-induced gene transcription and epigenetic modulation in an in vivo patient-derived xenograft (PDX) model of acute lymphoblastic leukemia (ALL). Performing chromatin immunoprecipitation and assays for transposase-accessible chromatin (ATAC), we determined a strong correlation between glucocorticoid receptor (GR) binding and chromatin accessibility, acetylated histone marks and binding of a DNA structural protein (CTCF), and identified 1,536 GR bound LSOs. We next analyzed RNA-seq data in glucocorticoid sensitive and resistant ALL PDXs for expression changes in genes located within 100 kb of the LSOs after glucocorticoid treatment in vivo, and identified four groups comprising 389 genes that were significantly differentially expressed. Of the 198 up-regulated genes, 143 showed increased H3K27Ac enrichments at 177 LSOs. Forty-two LSOs showed the increase only in glucocorticoid sensitive but not resistant ALLs. Applying this to an extended panel of ALL PDXs, basal DNA methylation at the 42 LSOs was significantly higher in resistant ALLs despite no difference within the gene bodies, whereas the basal chromatin accessibility indicated by ATAC abundance was diminished. Pathway analysis indicated that the LSO regulated genes were involved in repressing B- and T- cell receptor signaling pathways and activating the apoptotic pathway. One such LSO was at the pro-apoptotic BIM gene locus, where CTCF binding was found only in lymphocytes but not in other cell types. The GR cooperated with CTCF to mediate interactions between the promoter and the BIM LSO to direct DNA looping, thus triggering BIM transcription. Importantly, this LSO was heavily methylated in resistant ALLs as well as non-lymphoid cells. Azacitidine, a DNA demethylating drug that is routinely used in the clinic, could partially reverse these changes and restore glucocorticoid sensitivity. Taken together, this study demonstrated for the first time that lymphocyte-specific epigenetic modifications pre-determine glucocorticoid resistance in ALL and may account for the lack of glucocorticoid sensitivity in other cell type. Reversal of these epigenetic changes may lead to improvements in the use of glucocorticoids in the clinic. Citation Format: Duohui Jing, Yizhou Huang, Xiaoyun Liu, Keith Sia, Rebecca C. Poulos, Miriam Span, Chao Zhang, Jianqing Mi, Jason WH Wong, Dominik Beck, John E. Pimanda, Richard B. Lock. Lymphocyte-specific chromatin accessibility predetermines glucocorticoid resistance in acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3173.

Finding Genetic Regulators, Forecasting Function

Finding Genetic Regulators, Forecasting Function

Abstract 224: Evidence for a shared genetic basis between breast cancer and mammographic density phenotypes

Abstract Background: Radiologically dense area (DA) on a mammogram represents the stromal and epithelial tissue in the breast, while non-dense area (NDA) represents adipose tissue. Quantitative measurements of DA, NDA as well as the proportion of dense area in the breast (percent density, PD) adjusted for age and body mass index (BMI) have all been associated with breast cancer (BC) risk. Previous genome-wide association studies (GWAS) have identified genetic loci that show association with both mammographic density (MD) phenotypes and BC, but little is known about the genome-wide and local genetic correlation between MD phenotypes and BC. Methods: We estimated the genome-wide and local genetic correlation between MD phenotypes adjusted for age and BMI and BC, leveraging GWAS summary statistics for MD phenotypes based on 12,866 women (MODE consortium) and for BC based on 122,977 cases and 105,974 controls (BCAC consortium). We also conducted partitioned heritability analysis to identify specific functional annotations that contribute to MD heritability. Results: We observed a negative genetic correlation between DA and NDA (-0.62, p=0.002). Both DA and PD showed positive genetic correlations with overall BC (DA: rg=0.24, p=0.004; PD: rg=0.27, p=0.0001), estrogen receptor-positive (ER+) BC (DA: rg=0.23, p=0.005; PD: rg=0.23, p=0.001), and ER- BC (DA: rg=0.26, p=0.01; PD: rg=0.32, p=0.0006). In contrast, NDA showed a nonsignificant negative genetic correlation with BC (rg=-0.11, p=0.21). Local genetic correlation analyses between MD phenotypes and BC reflected known findings including positive correlation at 10q21 region between both PD and DA and BC (both p&amp;lt;5x10-7), negative correlation at 22q13.1 between DA and BC (p=4.6x10-6), and positive correlation at 8p11.23 between NDA and BC (p=2.8x10-6). Two known BC regions that have not previously been associated with MD showed suggestive genetic correlations between BC and DA: 1q21.2 (negative genetic correlation, p=1.3x10-4) and 12q21.31 (positive genetic correlation, p=6.7x10-4). In exploratory analyses estimating the genome-wide genetic correlation between MD and hormone-related phenotypes in UK Biobank (n=337,199), we observed nominally significant positive genetic correlations between PD and benign mammary dysplasia (rg=0.45, p=0.01), NDA and endometriosis (rg=0.45, p=0.02), and negative genetic correlation between NDA and menstrual cycle length (rg=-0.30, p=0.05). We identified multiple functional annotations that contribute to MD heritability including histone mark H3K27ac that marks active enhancers (~2,5-fold enrichment in heritability for both DA and NDA) and DNAsel Hypersensitive Sites (DHS, 9.1-fold enrichment in heritability for DA). Conclusion: Our results provide evidence of shared genetic basis between mammographic density and breast cancer. Citation Format: Hongjie Chen, Julie Douglas, William Gordon, Deborah Thompson, Fergus Couch, Peter Fasching, Per Hall, John Hopper, Gertraud Maskarenic, Marina Pollan, Jennifer Stone, Rulla Tamimi, Sara Lindstroem. Evidence for a shared genetic basis between breast cancer and mammographic density phenotypes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 224.

PO-348 The other 98%: making sense of non-coding variation

IntroductionDespite the ever-increasing availability of whole genome sequencing datasets, our ability to understand functional impact of genetic variation remains largely restricted to the coding portion of the genome. Although a selection of methods exist to estimate variant damage potential, these methods are heavily biassed towards coding variation and have limited ability to detect potentially functional variants in non-coding regions. Furthermore, these methods either ignore tissue context of the variant or combine epigenetic annotations from multiple tissues, often dominated by blood cell lines, and are therefore poorly suited for predicting functional impact of variants in regulatory elements that act in a tissue-specific manner.Material and methodsIn order to address these limitations, we have developed a machine learning classifier aimed at prioritising non-coding variants while taking into account relevant tissue context. We defined ’functional impact’ as the propensity of a variant to cause allele-specific chromatin accessibility. DNase Hypersensitive Site (DHS) variants found to display preferential accessibility, as well as a matched set of negative variants, were annotated with a range of features including a set of core epigenetic marks from the ROADMAP Epigenomics project matched to the tissue of origin of the positive variant set. We then combined several machine learning algorithms to train an ensemble classifier. On a balanced test set (equal number of positive and negative variants) our model achieved area under curve (AUC) of ~75%. On a realistic dataset, where the negative variants outnumbered the positive ones 100:1, we achieved AUC of 90%.Results and discussionsTo demonstrate the importance of tissue context when estimating functional impact of non-coding variants, we compared the classifier’s performance on the same set of variants annotated with epigenetic information from either a closely matched or an unrelated tissue. We observed a notable drop in performance when using epigenetic context from a mismatched tissue. We show that our method drastically outperforms existing damage estimation tools in its ability to predict allele-specific chromatin accessibility. Finally, we demonstrate the utility of our method by successfully prioritising experimentally validated regulatory variants among the large number of variants within the same linkage disequilibrium block.ConclusionThe above method will aid in interpretation of whole genome sequencing datasets generated in cancer and rare disease studies.

Identification of open chromosomal regions and key genes in prostate cancer via integrated analysis of DNase‑seq and RNA‑seq data.

Prostate cancer is a type of adenocarcinoma arising from the peripheral zone of the prostate gland, and metastasized prostate cancer is incurable with the current available therapies. The present study aimed to identify open chromosomal regions and differentially expressed genes (DEGs) associated with prostate cancer development. The DNase sequencing data (GSE33216) and RNA sequencing data (GSE22260) were downloaded from the Gene Expression Omnibus database. DNaseI hypersensitive sites were detected and analyzed. Subsequently, DEGs were identified and their potential functions were enriched. Finally, upstream regulatory elements of DEGs were predicted. In LNCaP cells, following androgen receptor activation, 244 upregulated and 486 downregulated open chromosomal regions were identified. However, only 1% of the open chromosomal regions were dynamically altered. The 41genes with upregulated open chromosomal signals within their promoter regions were primarily enriched in biological processes. Additionally, 211upregulated and 150downregulated DEGs were identified in prostate cancer, including eight transcription factors (TFs). Finally, nine regulatory elements associated with prostate cancer were predicted. In particular, inhibitor of DNA binding1 HLH protein (ID1) was the only significantly upregulated TF which exhibited motif enrichment in the promoter regions of upregulated genes. CCCTC‑binding factor (CTCF) and ELK1 ETS transcription factor (ELK1), enriched in the open promoter regions of downregulated genes, were potential upstream regulatory elements. Furthermore, reverse transcription‑quantitative polymerase chain reaction analysis confirmed that ID1 expression was significantly upregulated in LNCaP cells and 5α‑dihydrotestosterone (DHT)‑treated LNCaP cells compared with that in BPH1 cells, while CTCF and ELK1 expression was significantly downregulated in LNCaP cells and DHT‑treated LNCaP cells. In conclusion, ID1, CTCF and ELK1 may be associated with prostate cancer, and may be potential therapeutic targets for the treatment of this disease.

Histone variants H3.3 and H2A.Z are incorporated into the β-globin locus during transcription activation via different mechanisms

Histone variants H3.3 and H2A.Z are often enriched in enhancers and transcriptionally active genes. However, the incorporation dynamics of these variants and the mechanisms of their incorporation are unclear. Here, we examined the distribution of H3.3 and H2A.Z in the human β-globin locus and analyzed their incorporation dynamics during transcription activation. Locus control region hypersensitive sites (LCR HSs), acting as enhancers, and active globin genes were enriched by H3.3 and H2A.Z in erythroid K562 cells, but inactive globin genes were not. Both variants were dynamically incorporated into the β-globin locus after transcription induction in MEL/ch11 cells, and prior to gene transcription the LCR HSs became enriched with the variants. In the activated β-globin gene, H3.3 was highly incorporated during transcription, whereas H2A.Z incorporation appeared to precede it. To further explore the relationship between gene transcription and variant incorporation, we deleted the LCR HS3 enhancer or the β-globin proximal promoter from the β-globin locus using the CRISPR-Cas9 genome editing system. H2A.Z was incorporated into the β-globin gene in the locus lacking promoter, even though the β-globin gene transcription was abolished by these deletions. However, H3.3 incorporation was reduced in the untranscribed β-globin gene. These results suggest that H3.3 and H2A.Z are systematically incorporated into the LCR enhancer and β-globin gene as part of transcription activation, but that their incorporation is carried out via different mechanisms.

Chromosomal architecture and placental expression of the human growth hormone gene family are targeted by pre-pregnancy maternal obesity.

The human (h) placental lactogenic hormone chorionic somatomammotropin (CS) is highly produced during pregnancy and acts as a metabolic adaptor in response to maternal insulin resistance. Maternal obesity can exacerbate this "resistance", and a >75% decrease in CS RNA levels was observed in term placentas from obese vs. lean women. The genes coding for hCS ( hCS-A and hCS-B) and placental growth hormone ( hGH-V) as well as the hCS-L pseudogene and pituitary growth hormone (GH) gene ( hGH-N) are located at a single locus on chromosome 17. Three remote hypersensitive sites (HS III-V) located >28 kb upstream of hGH-N as well as local hCS gene promoter and enhancer regions are implicated in hCS gene expression. A placenta-specific chromosomal architecture, including interaction between HS III-V and hCS but not hGH gene promoters, was detected in placentas from lean women (BMI <25 kg/m2) by using the chromosome conformation capture assay. This architecture was disrupted by pre-pregnancy maternal obesity (BMI >35 kg/m2), resulting in a predominant interaction between HS III and the hGH-N promoter, which was also observed in nonplacental tissues. This was accompanied by a decrease in hCS levels, which was consistent with reduced RNA polymerase II and CCAAT/enhancer-binding protein-β association with individual hCS promoter and enhancer sequences, respectively. Thus, pre-pregnancy maternal obesity disrupts the placental hGH/CS gene locus chromosomal architecture. However, based on data from obese women who develop GDM, insulin treatment partially recapitulates the chromosomal architecture seen in lean women and positively affects hCS production, presumably facilitating prolactin receptor-related signaling by hCS.

Prediction and identification of transcriptional regulatory elements at the lung cancer-specific DKK1 locus.

The glycoprotein dickkopf 1 (DKK1) is highly expressed in lung cancer cell lines and tissues. Our previous study demonstrated that DKK1 promoter activity is low in lung cancer cell lines. This may be because it lacks the necessary transcriptional regulatory elements (TREs) required for higher activity levels. However, it is difficult to computationally predict functionally significant TREs, as TREs from different locations can affect large segments of distant DNA. The Encyclopedia of DNA Elements project features multiple integrated technologies and approaches for the discovery and definition of functional elements, including enhancer elements and enhancer-blocking insulators. In the present study, DNase I hypersensitive sites and histone modifications of DKK1 were investigated in the A549 lung cancer cell line using the UCSC Genome Browser. A set of cis-acting enhancer elements were identified by a dual-luciferase reporter gene assay system to increase activity of the DKK1 promoter with lung cancer specificity. To the best of our knowledge, these data provide the first insight into the role of the DKK1 locus in lung cancer, and confirm the contribution of intronic cis-acting elements to the regulation of DKK1 expression, providing a new insight into gene regulation in lung cancer, which could inform the development of targeted therapy.

Reprogramming of Chromatin Accessibility in Somatic Cell Nuclear Transfer Is DNA Replication Independent

SUMMARYMammalian oocytes have the ability to reset the transcriptional program of differentiated somatic cells into that of totipotent embryos through somatic cell nuclear transfer (SCNT). However, the mechanisms underlying SCNT-mediated reprogramming are largely unknown. To understand the mechanisms governing chromatin reprogramming during SCNT, we profiled DNase I hypersensitive sites (DHSs) in donor cumulus cells and one-cell stage SCNT embryos. To our surprise, the chromatin accessibility landscape of the donor cells is drastically changed to recapitulate that of the in vitro fertilization (IVF)-derived zygotes within 12 hr. Interestingly, this DHS reprogramming takes place even in the presence of a DNA replication inhibitor, suggesting that SCNT-mediated DHS reprogramming is independent of DNA replication. Thus, this study not only reveals the rapid and drastic nature of the changes in chromatin accessibility through SCNT but also establishes a DNA replication-independent model for studying cellular reprogramming.