Dynamic transcription programs during ES cell differentiation towards mesoderm in serum versus serum-freeBMP4 culture

Expression profiling of embryonic stem (ES) cell differentiation in the presence of serum has been performed previously. It remains unclear if transcriptional activation is dependent on complext growth factor mixtures in serum or whether this process is intrinsic to ES cells once the stem cell program has been inactivated. The aims of this study were to determine the transcriptional programs associated with the stem cell state and to characterize mesoderm differentiation between serum and serum-free culture

A global role for KLF1 in erythropoiesis revealed by ChIP-seq in primary erythroid cells

KLF1 regulates a diverse suite of genes to direct erythroid cell differentiation from bipotent progenitors. To determine the local cis-regulatory contexts and transcription factor networks in which KLF1 operates, we performed KLF1 ChIP-seq in the mouse. We found at least 945 sites in the genome of E14.5 fetal liver erythroid cells which are occupied by endogenous KLF1. Many of these recovered sites reside in erythroid gene promoters such as Hbb-bl, but the majority are distant to any known gene. Our data suggests KLF1 directly regulates most aspects of terminal erythroid differentiation including production of alpha- and beta-globin protein chains, heme biosynthesis, coordination of proliferation and anti-apoptotic pathways, and construction of the red cell membrane and cytoskeleton by functioning primarily as a transcriptional activator. Additionally, we suggest new mechanisms for KLF1 cooperation with other transcription factors, in particular the erythroid transcription factor GATA1, to maintain homeostasis in the erythroid compartment

Ampullary cancers harbor ELF3 tumor suppressor gene mutations and exhibit frequent WNT dysregulation

The ampulla of Vater is a complex cellular environment from which adenocarcinomas arise to form a group of histopathologically heterogenous tumors. To evaluate the molecular features of these tumors, 98 ampullary adenocarcinomas were evaluated and compared to 44 distal bile duct and 18 duodenal adenocarcinomas. Genomic analyses revealed mutations in the WNT signaling pathway among half of the patients and in all three adenocarcinomas irrespective of their origin and histological morphology. These tumors were characterized by a high frequency of inactivating mutations of ELF3, a high rate of microsatellite instability, and common focal deletions and amplifications, suggesting common attributes in the molecular pathogenesis are at play in these tumors. The high frequency of WNT pathway activating mutation, coupled with small-molecule inhibitors of β-catenin in clinical trials, suggests future treatment decisions for these patients may be guided by genomic analysis

Eddy covariance measurements highlight sources of nitrogen oxide emissions missing from inventories for central London

During March–June 2017 emissions of nitrogen oxides were measured via eddy covariance at the British Telecom Tower in central London, UK. Through the use of a footprint model the expected emissions were simulated from the spatially resolved National Atmospheric Emissions Inventory for 2017 and compared with the measured emissions. These simulated emissions were shown to underestimate measured emissions during the daytime by a factor of 1.48, but they agreed well overnight. Furthermore, underestimations were spatially mapped, and the areas around the measurement site responsible for differences in measured and simulated emissions were inferred. It was observed that areas of higher traffic, such as major roads near national rail stations, showed the greatest underestimation by the simulated emissions. These discrepancies are partially attributed to a combination of the inventory not fully capturing traffic conditions in central London and both the spatial and temporal resolution of the inventory not fully describing the high heterogeneity of the urban centre. Understanding of this underestimation may be further improved with longer measurement time series to better understand temporal variation and improved temporal scaling factors to better simulate sub-annual emissions

A Continuum of Cell States Spans Pluripotency and Lineage Commitment in Human Embryonic Stem Cells

Background: Commitment in embryonic stem cells is often depicted as a binary choice between alternate cell states, pluripotency and specification to a particular germ layer or extraembryonic lineage. However, close examination of human ES cell cultures has revealed significant heterogeneity in the stem cell compartment. Methodology/Principal Findings: We isolated subpopulations of embryonic stem cells using surface markers, then examined their expression of pluripotency genes and lineage specific transcription factors at the single cell level, and tested their ability to regenerate colonies of stem cells. Transcript analysis of single embryonic stem cells showed that there is a gradient and a hierarchy of expression of pluripotency genes in the population. Even cells at the top of the hierarchy generally express only a subset of the stem cell genes studied. Many cells co-express pluripotency and lineage specific genes. Cells along the continuum show a progressively decreasing likelihood of self renewal as their expression of stem cell surface markers and pluripotency genes wanes. Most cells that are positive for stem cell surface markers express Oct-4, but only those towards the top of the hierarchy express the nodal receptor TDGF-1 and the growth factor GDF3. Significance: These findings on gene expression in single embryonic stem cells are in concert with recent studies of early mammalian development, which reveal molecular heterogeneity and a stochasticity of gene expression in blastomeres. Our work indicates that only a small fraction of the population resides at the top of the hierarchy, that lineage priming (co-expression of stem cell and lineage specific genes) characterizes pluripotent stem cell populations, and that extrinsic signaling pathways are upstream of transcription factor networks that control pluripotency

Use of DNA–Damaging Agents and RNA Pooling to Assess Expression Profiles Associated with BRCA1 and BRCA2 Mutation Status in Familial Breast Cancer Patients

A large number of rare sequence variants of unknown clinical significance have been identified in the breast cancer susceptibility genes, BRCA1 and BRCA2. Laboratory-based methods that can distinguish between carriers of pathogenic mutations and non-carriers are likely to have utility for the classification of these sequence variants. To identify predictors of pathogenic mutation status in familial breast cancer patients, we explored the use of gene expression arrays to assess the effect of two DNA–damaging agents (irradiation and mitomycin C) on cellular response in relation to BRCA1 and BRCA2 mutation status. A range of regimes was used to treat 27 lymphoblastoid cell-lines (LCLs) derived from affected women in high-risk breast cancer families (nine BRCA1, nine BRCA2, and nine non-BRCA1/2 or BRCAX individuals) and nine LCLs from healthy individuals. Using an RNA–pooling strategy, we found that treating LCLs with 1.2 µM mitomycin C and measuring the gene expression profiles 1 hour post-treatment had the greatest potential to discriminate BRCA1, BRCA2, and BRCAX mutation status. A classifier was built using the expression profile of nine QRT–PCR validated genes that were associated with BRCA1, BRCA2, and BRCAX status in RNA pools. These nine genes could distinguish BRCA1 from BRCA2 carriers with 83% accuracy in individual samples, but three-way analysis for BRCA1, BRCA2, and BRCAX had a maximum of 59% prediction accuracy. Our results suggest that, compared to BRCA1 and BRCA2 mutation carriers, non-BRCA1/2 (BRCAX) individuals are genetically heterogeneous. This study also demonstrates the effectiveness of RNA pools to compare the expression profiles of cell-lines from BRCA1, BRCA2, and BRCAX cases after treatment with irradiation and mitomycin C as a method to prioritize treatment regimes for detailed downstream expression analysis

HNF4A and GATA6 loss reveals therapeutically actionable subtypes in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) can be divided into transcriptomic subtypes with two broad lineages referred to as classical (pancreatic) and squamous. We find that these two subtypes are driven by distinct metabolic phenotypes. Loss of genes that drive endodermal lineage specification, HNF4A and GATA6, switch metabolic profiles from classical (pancreatic) to predominantly squamous, with glycogen synthase kinase 3 beta (GSK3β) a key regulator of glycolysis. Pharmacological inhibition of GSK3β results in selective sensitivity in the squamous subtype; however, a subset of these squamous patient-derived cell lines (PDCLs) acquires rapid drug tolerance. Using chromatin accessibility maps, we demonstrate that the squamous subtype can be further classified using chromatin accessibility to predict responsiveness and tolerance to GSK3β inhibitors. Our findings demonstrate that distinct patterns of chromatin accessibility can be used to identify patient subgroups that are indistinguishable by gene expression profiles, highlighting the utility of chromatin-based biomarkers for patient selection in the treatment of PDAC

Transcriptome-Wide Prediction of miRNA Targets in Human and Mouse Using FASTH

Transcriptional regulation by microRNAs (miRNAs) involves complementary base-pairing at target sites on mRNAs, yielding complex secondary structures. Here we introduce an efficient computational approach and software (FASTH) for genome-scale prediction of miRNA target sites based on minimizing the free energy of duplex structure. We apply our approach to identify miRNA target sites in the human and mouse transcriptomes. Our results show that short sequence motifs in the 5′ end of miRNAs frequently match mRNAs perfectly, not only at validated target sites but additionally at many other, energetically favourable sites. High-quality matching regions are abundant and occur at similar frequencies in all mRNA regions, not only the 3′UTR. About one-third of potential miRNA target sites are reassigned to different mRNA regions, or gained or lost altogether, among different transcript isoforms from the same gene. Many potential miRNA target sites predicted in human are not found in mouse, and vice-versa, but among those that do occur in orthologous human and mouse mRNAs most are situated in corresponding mRNA regions, i.e. these sites are themselves orthologous. Using a luciferase assay in HEK293 cells, we validate four of six predicted miRNA-mRNA interactions, with the mRNA level reduced by an average of 73%. We demonstrate that a thermodynamically based computational approach to prediction of miRNA binding sites on mRNAs can be scaled to analyse complete mammalian transcriptome datasets. These results confirm and extend the scope of miRNA-mediated species- and transcript-specific regulation in different cell types, tissues and developmental conditions