Aflatoxin is a natural potent carcinogen and a major cause of liver cancer. However, the molecular mechanisms of hepatocellular carcinogenesis remain largely unexplored. In this study, we profiled global gene expression in liver tissues of rats that developed hepatocellular carcinoma (HCC) from aflatoxin B1 (AFB1) administration and those that were AFB1-resistant, as well as rats without AFB1 exposure as a control. AFB1 exposure resulted in extensive perturbation in gene expression with different functions in HCC and AFB1 resistance (AR) samples. The differentially expressed genes (DEGs) in HCC sample were enriched for cell proliferation, cell adhesion and vasculature development that largely contribute to carcinogenesis. Anti-apoptosis genes were up-regulated in HCC sample whereas apoptosis-induction genes were up-regulated in AR sample. AFB1 exposure also caused extensive alteration in expression level of lncRNAs. Among all the 4511 annotated lncRNAs, half of them were highly expressed only in HCC sample and up-regulated a group of protein-coding genes with cancer-related functions: apoptosis regulation, DNA repair, and cell cycle. Intriguingly, these genes were down-regulated by lncRNAs highly expressed in AR sample. Collectively, apoptosis is the critical biological process for carcinogenesis in response to AFB1 exposure through changes in expression level of both protein-coding and lncRNA genes.

The core pluripotency factor Oct4 plays key roles in somatic cell reprogramming through transcriptional control. Here, we profile Oct4 occupancy, epigenetic changes, and gene expression in reprogramming. We find that Oct4 binds in a hierarchical manner to target sites with primed epigenetic modifications. Oct4 binding is temporally continuous and seldom switches between bound and unbound. Oct4 occupancy in most of promoters is maintained throughout the entire reprogramming process. In contrast, somatic cell-specific enhancers are silenced in the early and intermediate stages, whereas stem cell-specific enhancers are activated in the late stage in parallel with cell fate transition. Both epigenetic remodeling and Oct4 binding contribute to the hyperdynamic enhancer signature transitions. The hierarchical Oct4 bindings are associated with distinct functional themes at different stages. Collectively, our results provide a comprehensive molecular roadmap of Oct4 binding in concert with epigenetic rearrangements and rich resources for future reprogramming studies.

Nucleosome organization determines the chromatin state, which in turn controls gene expression or silencing. Nucleosome remodeling occurs during somatic cell reprogramming, but it is still unclear to what degree the re-established nucleosome organization of induced pluripotent stem cells (iPSCs) resembles embryonic stem cells (ESCs), and whether the iPSCs inherit some residual gene expression from the parental fibroblast cells.

Enhancers are distal cis-regulatory elements that modulate gene expression. They are depleted of nucleosomes and enriched in specific histone modifications; thus, calling DNase-seq and histone mark ChIP-seq peaks can predict enhancers. We evaluated nine peak-calling algorithms for predicting enhancers validated by transgenic mouse assays. DNase and H3K27ac peaks were consistently more predictive than H3K4me1/2/3 and H3K9ac peaks. DFilter and Hotspot2 were the best DNase peak callers, while HOMER, MUSIC, MACS2, DFilter and F-seq were the best H3K27ac peak callers. We observed that the differential DNase or H3K27ac signals between two distant tissues increased the area under the precision-recall curve (PR-AUC) of DNase peaks by 17.5-166.7% and that of H3K27ac peaks by 7.1-22.2%. We further improved this differential signal method using multiple contrast tissues. Evaluated using a blind test, the differential H3K27ac signal method substantially improved PR-AUC from 0.48 to 0.75 for predicting heart enhancers. We further validated our approach using postnatal retina and cerebral cortex enhancers identified by massively parallel reporter assays, and observed improvements for both tissues. In summary, we compared nine peak callers and devised a superior method for predicting tissue-specific mouse developmental enhancers by reranking the called peaks.

Human adipose tissue-derived stem cells (ADSCs) are an attractive multipotent stem cell source with therapeutic applicability across diverse fields for the repair and regeneration of acute and chronically damaged tissues. In recent years, there has been increasing interest in ADSC for tissue engineering applications. However, the mechanisms underlying the regulation of ADSC proliferation are not fully understood. Here we show that 47 transcripts are up-regulated while 23 are down-regulated in ADSC compared to terminally differentiated cells based on global mRNA profiling and microRNA profiling. Among the up-regulated genes, the expression of vascular endothelial growth factor (VEGF) is fine-tuned by miR-199a-5p. Further investigation indicates that VEGF accelerates ADSC proliferation whereas the multipotency of ADSC remains stable in terms of adipogenic, chondrogenic and osteogenic potentials after VEGF treatment, suggesting that VEGF may serve as an excellent supplement for accelerating ADSC proliferation during in vitro expansion.

Bisphenol A (BPA) is an environmental endocrine disruptor which has been detected in human bodies. Many studies have implied that BPA exposure is harmful to human health. Previous studies mainly focused on BPA effects on estrogen receptor (ER)-positive cells. Genome-wide impacts of BPA on gene expression in ER-negative cells is unclear. In this study, we performed RNA-seq to characterize BPA-induced cellular and molecular impacts on ER-negative HEK293 cells. The microscopic observation showed that low-dose BPA exposure did not affect cell viability and morphology. Gene expression profiling analysis identified a list of differentially expressed genes in response to BPA exposure in HEK293 cells. These genes were involved in variable important biological processes including ion transport, cysteine metabolic process, apoptosis, DNA damage repair, etc. Notably, BPA up-regulated the expression of ERCC5 encoding a DNA endonuclease for nucleotide-excision repair. Further electrochemical experiment showed that BPA induced significant DNA damage in ER-positive MCF-7 cells but not in ER-negative HEK293 cells. Collectively, our study revealed that ER-negative HEK293 cells employed mechanisms in response to BPA exposure different from ER-positive cells.

Both environmental cues and intracellular bioenergetic states profoundly affect intracellular pH (pHi). How a cell responds to pHi changes to maintain bioenergetic homeostasis remains elusive. Here we show that Smad5, a well-characterized downstream component of bone morphogenetic protein (BMP) signaling responds to pHi changes. Cold, basic or hypertonic conditions increase pHi, which in turn dissociates protons from the charged amino acid clusters within the MH1 domain of Smad5, prompting its relocation from the nucleus to the cytoplasm. On the other hand, heat, acidic or hypotonic conditions decrease pHi, blocking the nuclear export of Smad5, and thus causing its nuclear accumulation. Active nucleocytoplasmic shuttling of Smad5 induced by environmental changes and pHi fluctuation is independent of BMP signaling, carboxyl terminus phosphorylation and Smad4. In addition, ablation of Smad5 causes chronic and irreversible dysregulation of cellular bioenergetic homeostasis and disrupted normal neural developmental processes as identified in a differentiation model of human pluripotent stem cells. Importantly, these metabolic and developmental deficits in Smad5-deficient cells could be rescued only by cytoplasmic Smad5. Cytoplasmic Smad5 physically interacts with hexokinase 1 and accelerates glycolysis. Together, our findings indicate that Smad5 acts as a pHi messenger and maintains the bioenergetic homeostasis of cells by regulating cytoplasmic metabolic machinery.

Epigenetic mechanisms play key roles in initiation and progression of prostate cancer by changing gene expression. The Prostate Epigenetic Database (PEpiD: http://wukong.tongji.edu.cn/pepid) archives the three extensively characterized epigenetic mechanisms DNA methylation, histone modification, and microRNA implicated in prostate cancer of human, mouse, and rat. PEpiD uses a distinct color scheme to present the three types of epigenetic data and provides a user-friendly interface for flexible query. The retrieved information includes Refseq ID, gene symbol, gene alias, genomic loci of epigenetic changes, tissue source, experimental method, and supportive references. The change of histone modification (hyper or hypo) and the corresponding gene expression change (up or down) are also indicated. A graphic view of DNA methylation with exon-intron structure and predicted CpG islands is provided as well. Moreover, the prostate-related ENCODE tracks (DNA methylation, histone modifications, chromatin remodelers), and other key transcription factors with reported roles in prostate are displayed in the browser as well. The reversibility of epigenetic aberrations has made them potential markers for diagnosis and prognosis, and targets for treatment of cancers. This curated information will improve our understanding of epigenetic mechanisms of gene regulation in prostate cancer, and serve as an important resource for epigenetic research in prostate cancer.

Dictyostelium discoideum is an amoebozoa that exists in both a free-living unicellular and a multicellular form. It is situated in a deep branch in the evolutionary tree and is particularly noteworthy in having a very A/T-rich genome. Dictyostelium provides an ideal system to examine the extreme to which nucleotide bias may be employed in organizing promoters, genes, and nucleosomes across a genome. We find that Dictyostelium genes are demarcated precisely at their 5' ends by poly-T tracts and precisely at their 3' ends by poly-A tracts. These tracts are also associated with nucleosome-free regions and are embedded with precisely positioned TATA boxes. Homo- and heteropolymeric tracts of A and T demarcate nucleosome border regions. Together, these findings reveal the presence of a variety of functionally distinct polymeric A/T elements. Strikingly, Dictyostelium chromatin may be organized in di-nucleosome units but is otherwise organized as in animals. This includes a +1 nucleosome in a position that predicts the presence of a paused RNA polymerase II. Indeed, we find a strong phylogenetic relationship between the presence of the NELF pausing factor and positioning of the +1 nucleosome. Pausing and +1 nucleosome positioning may have coevolved in animals.

The local environment of single nucleotide polymorphisms (SNPs) contains abundant genetic information for the study of mechanisms of mutation, genome evolution, and causes of diseases. Recent studies revealed that neighboring-nucleotide biases on SNPs were strong and the genome-wide bias patterns could be represented by a small subset of the total SNPs. It remains unsolved for the estimation of the effective SNP size, the number of SNPs that are sufficient to represent the bias patterns observed from the whole SNP data.

Induced pluripotent stem cells can be derived from somatic cells through ectopic expression of transcription factors or chemical cocktails. Chemical iPSCs (C-iPSCs) and OSKM-iPSCs (4F-iPSCs) have been suggested to have similar characteristics to mouse embryonic stem cells (mESCs). However, their epigenetic equivalence remains incompletely understood throughout the genome. In this study, we have generated mouse C-iPSCs and 4F-iPSCs, and further compared the genome-wide DNA methylomes of C-iPSCs, 4F-iPSCs, and mESCs that were maintained in 2i and LIF. Three pluripotent stem cells tend to be low methylated overall, however, DNA methylations in some specific regions (such as retrotransposons) are cell type-specific. Importantly, C-iPSCs are more hypomethylated than 4F-iPSCs. Bisulfite sequencing indicated that DNA methylation status in several known imprinted clusters, such as: Dlk1-Dio3 and Peg12-Ube3a, in C-iPSCs are closer to those of mESCs than 4F-iPSCs. Overall, our data demonstrate the reprogramming methods-dependent epigenetic differences of C-iPSCs and 4F-iPSCs and reveal that C-iPSCs are more hypomethylated than OSKM-integrated iPSCs.

Drosophila neural development undergoes extensive chromatin remodeling and precise epigenetic regulation. However, the roles of chromatin remodeling in establishment and maintenance of cell identity during cell fate transition remain enigmatic. Here, we compared the changes in gene expression, as well as the dynamics of nucleosome positioning and key histone modifications between the four major neural cell types during Drosophila neural development. We find that the neural progenitors can be separated from the terminally differentiated cells based on their gene expression profiles, whereas nucleosome distribution in the flanking regions of transcription start sites fails to identify the relationships between the progenitors and the differentiated cells. H3K27me3 signal in promoters and enhancers can not only distinguish the progenitors from the differentiated cells but also identify the differentiation path of the neural stem cells (NSCs) to the intermediate progenitor cells to the glial cells. In contrast, H3K9ac signal fails to identify the differentiation path, although it activates distinct sets of genes with neuron-specific and glia-related functions during the differentiation of the NSCs into neurons and glia, respectively. Together, our study provides novel insights into the crucial roles of chromatin remodeling in determining cell type during Drosophila neural development.

Identification of RNA targets of RNA-binding proteins (RBPs) is essential for complete understanding of their biological functions. However, it is still a challenge to identify the biologically relevant targets of RBPs through in vitro strategies of RIP-seq, HITS-CLIP, or GoldCLIP due to the potentially high background and complicated manipulation. In malaria parasites, RIP-seq and gene disruption are the few tools available currently for identification of RBP targets. Here, we have adopted the TRIBE (Targets of RNA binding proteins identified by editing) system to in vivo identify the RNA targets of PfDis3, a key exoribonuclease subunit of RNA exosome in Plasmodium falciparum. We generated a transgenic parasite line of PfDis3-ADARcd, which catalyzes an adenosine (A)-to-inosine (I) conversion at the potential interacting sites of PfDis3-targeting RNAs. Most of PfDis3 target genes contain one edit site. The majority of the edit sites detected by PfDis3-TRIBE locate in exons and spread across the entire coding regions. The nucleotides adjacent to the edit sites contain ∼75% of A + T. PfDis3-TRIBE target genes are biases toward higher RIP enrichment, suggesting that PfDis3-TRIBE preferentially detects stronger PfDis3 RIP targets. Collectively, PfDis3-TRIBE is a favorable tool to identify in vivo target genes of RBP with high efficiency and reproducibility. Additionally, the PfDis3-targeting genes are involved in stage-related biological processes during the blood-stage development. Thus PfDis3 appears to shape the dynamic transcriptional transcriptome of malaria parasites through post-transcriptional degradation of a variety of unwanted transcripts from both strands in the asexual blood stage.

Differentiated somatic cells can be reprogrammed to totipotent embryos through somatic cell nuclear transfer (SCNT) with low efficiency. The histone deacetylase inhibitor trichostatin A (TSA) has been found to improve SCNT efficiency, but the underlying mechanism remains undetermined. Here, we examined genome-wide H3K9ac during SCNT embryo development and found that aberrant H3K9ac regions resulted in reduced 2-cell genome activation. TSA treatment largely corrects aberrant acetylation in SCNT embryos with an efficiency that is dictated by the native epigenetic environment. We further identified that the overexpression of Dux greatly improves SCNT efficiency by correcting the aberrant H3K9ac signal at its target sites, ensuring appropriate 2-cell genome activation. Intriguingly, the improvement in development mediated by TSA and Kdm4b is impeded by Dux knockout in SCNT embryos. Together, our study reveals that reprogramming of H3K9ac is important for optimal SCNT efficiency and identifies Dux as a crucial transcription factor in this process.

Assisted reproductive technology has been widely applied in the treatment of human infertility. However, accumulating evidence indicates that in vitro fertilization (IVF) is associated with a low pregnancy rate, placental defects, and metabolic diseases in offspring. Here, we find that IVF manipulation notably disrupts extraembryonic tissue-specific gene expression, and 334 epiblast (Epi)-specific genes and 24 Epi-specific transcription factors are abnormally expressed in extraembryonic ectoderm (ExE) of IVF embryos at embryonic day 7.5. Combined histone modification analysis reveals that aberrant H3K4me3 modification at the Epi active promoters results in increased expression of these genes in ExE. Importantly, we demonstrate that knockdown of the H3K4me3-recruited regulator Kmt2e, which is highly expressed in IVF embryos, greatly improves the development of IVF embryos and reduces abnormal gene expression in ExE. Our study therefore identifies that abnormal H3K4me3 modification in extraembryonic tissue is a major cause of implantation failure and abnormal placental development of IVF embryos.

Haploid embryonic stem cells (haESCs) are derived from the inner cell mass of the haploid blastocyst, containing only one set of chromosomes. Extensive and accurate chromatin remodelling occurs during haESC derivation, but the intrinsic transcriptome profiles and chromatin structure of haESCs have not been fully explored. We profiled the transcriptomes, nucleosome positioning, and key histone modifications of four mouse haESC lines, and compared these profiles with those of other closely-related stem cell lines, MII oocytes, round spermatids, sperm, and mouse embryonic fibroblasts. haESCs had transcriptome profiles closer to those of naïve pluripotent stem cells. Consistent with the one X chromosome in haESCs, Xist was repressed, indicating no X chromosome inactivation. haESCs and ESCs shared a similar global chromatin structure. However, a nucleosome depletion region was identified in 2056 promoters in ESCs, which was absent in haESCs. Furthermore, three characteristic spatial relationships were formed between transcription factor motifs and nucleosomes in both haESCs and ESCs, specifically in the linker region, on the nucleosome central surface, and nucleosome borders. Furthermore, the chromatin state of 4259 enhancers was off in haESCs but active in ESCs. Functional annotation of these enhancers revealed enrichment in regulation of the cell cycle, a predominantly reported mechanism of haESC self-diploidization. Notably, the transcriptome profiles and chromatin structure of haESCs were highly preserved during passaging but different from those of differentiated cell types.

Neural crest defects lead to congenital heart disease involving outflow tract malformation. Integrin-linked-kinase (ILK) plays important roles in multiple cellular processes and embryogenesis. ILK is expressed in the neural crest, but its role in neural crest and outflow tract morphogenesis remains unknown.

Neuroblastoma is the most common solid tumor during early childhood. One of the key features of neuroblastoma is extensive tumor-driven angiogenesis due to hypoxia. However, the mechanism through which neuroblastoma cells drive angiogenesis is poorly understood. Here we show that the long noncoding RNA MALAT1 was upregulated in human neuroblastoma cell lines under hypoxic conditions. Conditioned media from neuroblastoma cells transfected with small interfering RNAs (siRNA) targeting MALAT1, compared with conditioned media from neuroblastoma cells transfected with control siRNAs, induced significantly less endothelial cell migration, invasion and vasculature formation. Microarray-based differential gene expression analysis showed that one of the genes most significantly down-regulated following MALAT1 suppression in human neuroblastoma cells under hypoxic conditions was fibroblast growth factor 2 (FGF2). RT-PCR and immunoblot analyses confirmed that MALAT1 suppression reduced FGF2 expression, and Enzyme-Linked Immunosorbent Assays revealed that transfection with MALAT1 siRNAs reduced FGF2 protein secretion from neuroblastoma cells. Importantly, addition of recombinant FGF2 protein to the cell culture media reversed the effects of MALAT1 siRNA on vasculature formation. Taken together, our data suggest that up-regulation of MALAT1 expression in human neuroblastoma cells under hypoxic conditions increases FGF2 expression and promotes vasculature formation, and therefore plays an important role in tumor-driven angiogenesis.

Neuroectoderm is an important neural precursor. However, chromatin remodeling and its epigenetic regulatory roles during the differentiation of human neuroectodermal cells (hNECs) from human embryonic stem cells (hESCs) remain largely unexplored. Here, we obtained hNECs through directed differentiation from hESCs, and determined chromatin states in the two cell types. Upon differentiation, H2A.Z-mediated nucleosome depletion leads to an open chromatin structure in promoters and upregulates expression of neuroectodermal genes. Increase in H3K9ac signals and decrease in H3K27me3 signals in promoters result in an active chromatin state and activate neuroectodermal genes. Conversely, decrease in H3K9ac signals and increase in H3K27me3 signals in promoters repress pluripotency genes. Moreover, H3K9ac signals facilitate the pluripotency factor Sox2 binding to target sites unique to hNECs. Knockdown of the acetyltransferase Kat2b erases H3K9ac signals, disrupts Sox2 binding, and fails the differentiation. Our results demonstrate a hierarchy of epigenetic regulation of gene expression during the differentiation of hNECs from hESCs through chromatin remodeling.

Nucleosome occupancy results in complex sequence variation rate heterogeneity by either increasing mutation rate or inhibiting DNA repair in yeast, fish, and human. H2A.Z nucleosome is extensively involved in gene transcription activation and regulation. To test whether H2A.Z nucleosome has the similar impact on sequence variability in the Drosophila genome, we profiled the H2A.Z nucleosome occupancy and sequence variation rate at gene ends and splicing sites. Consistent with previous studies, H2A.Z nucleosome positioning helps to demarcate the borders of exons. Nucleosome occupancy is anticorrelated with sequence divergence rate in the regions flanking transcription start sites and splicing sites. However, there is no rate heterogeneity between the linker DNA and H2A.Z nucleosomal DNA regardless of nucleosome occupancy, fuzziness, positioning in promoter, coding, and intergenic regions, young or old genes. But the rate at intergenic nucleosomes and the flanking linker regions is higher than that at the genic counterparts. Further analyses found that the high sequence divergence rate in the promoter regions that are usually nucleosome depleted regions may be likely resulted from the high mutation rate in the enriched tandem repeats. Interestingly, within nucleosomes spanning splicing sites, sequence variability of nucleosomal DNA significantly increases from the end within exons to the other end protruding into introns. The relaxed functional constraint in introns contributes to the high rate of nucleosomal DNA residing in introns while the strict functional constraint in exons maintains the low rate of nucleosomal DNA residing in exons. Taken together, H2A.Z nucleosome occupancy has no effect on sequence variability of Drosophila genome, which is likely determined by local sequence composition and the concomitant selection pressure.