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.

Cancer cells adopt various modes of migration during metastasis. How the ubiquitination machinery contributes to cancer cell motility remains underexplored. Here, we report that tripartite motif (TRIM) 59 is frequently up-regulated in metastatic breast cancer, which is correlated with advanced clinical stages and reduced survival among breast cancer patients. TRIM59 knockdown (KD) promoted apoptosis and inhibited tumor growth, while TRIM59 overexpression led to the opposite effects. Importantly, we uncovered TRIM59 as a key regulator of cell contractility and adhesion to control the plasticity of metastatic tumor cells. At the molecular level, we identified programmed cell death protein 10 (PDCD10) as a target of TRIM59. TRIM59 stabilized PDCD10 by suppressing RING finger and transmembrane domain-containing protein 1 (RNFT1)-induced lysine 63 (K63) ubiquitination and subsequent phosphotyrosine-independent ligand for the Lck SH2 domain of 62 kDa (p62)-selective autophagic degradation. TRIM59 promoted PDCD10-mediated suppression of Ras homolog family member A (RhoA)-Rho-associated coiled-coil kinase (ROCK) 1 signaling to control the transition between amoeboid and mesenchymal invasiveness. PDCD10 overexpression or administration of a ROCK inhibitor reversed TRIM59 loss-induced contractile phenotypes, thereby accelerating cell migration, invasion, and tumor formation. These findings establish the rationale for targeting deregulated TRIM59/PDCD10 to treat breast cancer.

The protein translational system, including transfer RNAs (tRNAs) and several categories of enzymes, plays a key role in regulating cell proliferation. Translation dysregulation also contributes to cancer development, though relatively little is known about the changes that occur to the translational system in cancer. Here, we present global analyses of tRNAs and three categories of enzymes involved in translational regulation in ~10,000 cancer patients across 31 cancer types from The Cancer Genome Atlas. By analyzing the expression levels of tRNAs at the gene, codon, and amino acid levels, we identified unequal alterations in tRNA expression, likely due to the uneven distribution of tRNAs decoding different codons. We find that overexpression of tRNAs recognizing codons with a low observed-over-expected ratio may overcome the translational bottleneck in tumorigenesis. We further observed overall overexpression and amplification of tRNA modification enzymes, aminoacyl-tRNA synthetases, and translation factors, which may play synergistic roles with overexpression of tRNAs to activate the translational systems across multiple cancer types.

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.

Here, we show that tumor ADORA1 deletion suppresses cell growth in human melanoma cell lines in vitro and tumor development in vivo in immune-deficient xenografts. However, this deletion induces the upregulation of PD-L1 levels, which inactivates cocultured T cells in vitro, compromises anti-tumor immunity in vivo, and reduces anti-tumor efficacy in an immune-competent mouse model. Functionally, PD-1 mAb treatment enhances the efficacy of ADORA1-deficient or ADORA1 antagonist-treated melanoma and NSCLC immune-competent mouse models. Mechanistically, we identify ATF3 as the factor transcriptionally upregulating PD-L1 expression. Tumor ATF3 deletion improves the effect of ADORA1 antagonist treatment of melanoma and NSCLC xenografts. We observe higher ADORA1, lower ATF3, and lower PD-L1 expression levels in tumor tissues from nonresponders among PD-1 mAb-treated NSCLC patients.

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.

RNA-binding proteins (RBPs) can recognize thousands of RNAs that help to maintain cell homeostasis, and RBP dysfunction is frequently observed in various cancers. However, whether specific RBPs are involved in tumor immune evasion by regulating programmed death ligand-1 (PD-L1) is unclear. Here, we perform targeted RBP CRISPR/Cas9 screening and identify density regulated re-initiation and release factor (DENR) as a PD-L1 regulator. DENR-depleted cancer cells exhibit reduced PD-L1 expression in vitro and in vivo. DENR depletion significantly suppresses tumor growth and enhances the tumor-killing activity of CD8+ T cells. Mechanistically, DENR antagonizes the translational repression of three consecutive upstream open reading frames (uORFs) upstream of Janus kinase 2 (Jak2); thus, DENR deficiency impairs JAK2 translation and the IFNγ-JAK-STAT signaling pathway, resulting in reduced PD-L1 expression in tumors. Overall, we discover an RBP DENR that could regulate PD-L1 expression for tumor immune evasion, and highlight the potential of DENR as a therapeutic target for immunotherapy.

One of the major obstacles to treating pancreatic ductal adenocarcinoma (PDAC) is its immunoresistant microenvironment. The functional importance and molecular mechanisms of Schwann cells in PDAC remains largely elusive. We characterized the gene signature of tumor-associated nonmyelinating Schwann cells (TASc) in PDAC and indicated that the abundance of TASc was correlated with immune suppressive tumor microenvironment and the unfavorable outcome of patients with PDAC. Depletion of pancreatic-specific TASc promoted the tumorigenesis of PDAC tumors. TASc-expressed long noncoding RNA (lncRNA) plasmacytoma variant translocation 1 (PVT1) was triggered by the tumor cell-produced interleukin-6. Mechanistically, PVT1 modulated RAF proto-oncogene serine/threonine protein kinase-mediated phosphorylation of tryptophan 2,3-dioxygenase in TASc, facilitating its enzymatic activities in catalysis of tryptophan to kynurenine. Depletion of TASc-expressed PVT1 suppressed PDAC tumor growth. Furthermore, depletion of TASc using a small-molecule inhibitor effectively sensitized PDAC to immunotherapy, signifying the important roles of TASc in PDAC immune resistance.

Cardiac differentiation from human pluripotent stem cells provides a unique opportunity to study human heart development in vitro and offers a potential cell source for cardiac regeneration. Compared to the large body of studies investigating cardiac maturation and cardiomyocyte subtype-specific induction, molecular events underlying cardiac lineage commitment from pluripotent stem cells at early stage remain poorly characterized.

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.