Polycomb Repressive Complex 2 Regulates Lineage Fidelity during Embryonic Stem Cell Differentiation

Polycomb Repressive Complex 2 (PRC2) catalyzes histone H3 lysine 27 tri-methylation (H3K27me3), an epigenetic modification associated with gene repression. H3K27me3 is enriched at the promoters of a large cohort of developmental genes in embryonic stem cells (ESCs). Loss of H3K27me3 leads to a failure of ESCs to properly differentiate, making it difficult to determine the precise roles of PRC2 during lineage commitment. Moreover, while studies suggest that PRC2 prevents DNA methylation, how these two epigenetic regulators coordinate to regulate lineage programs is poorly understood. Using several PRC2 mutant ESC lines that maintain varying levels of H3K27me3, we found that partial maintenance of H3K27me3 allowed for proper temporal activation of lineage genes during directed differentiation of ESCs to spinal motor neurons (SMNs). In contrast, genes that function to specify other lineages failed to be repressed in these cells, suggesting that PRC2 is also necessary for lineage fidelity. We also found that loss of H3K27me3 leads to a modest gain in DNA methylation at PRC2 target regions in both ESCs and in SMNs. Our study demonstrates a critical role for PRC2 in safeguarding lineage decisions and in protecting genes against inappropriate DNA methylation.National Cancer Institute (U.S.) (Cancer Center Support (Core) Grant P30-CA14051)National Institutes of Health (U.S.) (Training Grant T 32 GM007287)Smith Family Foundation (Contract LTR DATED 11/6/09

Pulsatile Blood Flow in Anatomically Accurate Vessels with Multiple Aneurysms: A Medical Intervention Planning Application of Computational Haemodynamics

In the present study, we are investigating computationally the pulsatile flow in an anatomically accurate cerebral arterial segment exhibiting two saccular aneurysms. Our focus is on the haemodynamic patterns observed within the two aneurysms, in terms of inflow-outflow regions, emergence and disappearance of coherent structures and mixing throughout the cardiac cycle. The results obtained carry interesting features, important for thrombosis, pharmacokinetics and particularly for interventional planning for aneurysm treatment. For the latter, being the center-point of this study, we show that the two aneurysms behave in a dissimilar manner, since the blood inflow region oscillates significantly for one of them and practically does not oscillate at all for the second. This information can guide the medical interventionist in designing the optimal approach, particularly in cases where total obliteration of the aneurysm neck opening is impossibl

The genetic basis for individual differences in mRNA splicing and APOBEC1 editing activity in murine macrophages

Alternative splicing and mRNA editing are known to contribute to transcriptome diversity. Although alternative splicing is pervasive and known to contribute to a variety of pathologies, including cancer, the genetic context for individual differences in isoform usage is still evolving. Similarly, although mRNA editing is ubiquitous and associated with important biological processes such as intracellular viral replication and cancer development, individual variations in and the genetic transmissibility of mRNA editing are equivocal. Here, we have used linkage analysis to show that both mRNA editing and alternative splicing are regulated by the macrophage genetic background and environmental cues. We show that distinct loci, potentially harboring variable splice factors, regulate the splicing of multiple transcripts. Additionally, we show that individual genetic variability at the Apobec1 locus results in differential rates of C-to-U editing in murine macrophages; with mouse strains expressing mostly a truncated isoform of Apobec1 exhibiting lower rates of editing. As a proof of concept, we have used linkage analysis to identify 36 high confidence novel edited sites. These results provide a novel and complementary method that can be used to identify C-to-U editing sites in individuals segregating at specific loci and show that, beyond individual DNA sequence and structural changes, differential isoform usage and mRNA editing can contribute to intra-species genomic and phenotypic diversity

Methods for time series analysis of RNA-seq data with application to human Th17 cell differentiation

Motivation: Gene expression profiling using RNA-seq is a powerful technique for screening RNA species’ landscapes and their dynamics in an unbiased way. While several advanced methods exist for differential expression analysis of RNA-seq data, proper tools to anal.yze RNA-seq time-course have not been proposed. Results: In this study, we use RNA-seq to measure gene expression during the early human T helper 17 (Th17) cell differentiation and T-cell activation (Th0). To quantify Th17-specific gene expression dynamics, we present a novel statistical methodology, DyNB, for analyzing time-course RNA-seq data. We use non-parametric Gaussian processes to model temporal correlation in gene expression and combine that with negative binomial likelihood for the count data. To account for experiment-specific biases in gene expression dynamics, such as differences in cell differentiation efficiencies, we propose a method to rescale the dynamics between replicated measurements. We develop an MCMC sampling method to make inference of differential expression dynamics between conditions. DyNB identifies several known and novel genes involved in Th17 differentiation. Analysis of differentiation efficiencies revealed consistent patterns in gene expression dynamics between different cultures. We use qRT-PCR to validate differential expression and differentiation efficiencies for selected genes. Comparison of the results with those obtained via traditional timepoint-wise analysis shows that time-course analysis together with time rescaling between cultures identifies differentially expressed genes which would not otherwise be detected. Availability: An implementation of the proposed computational methods will be available at http://research.ics.aalto.fi/csb/software/Academy of Finland (Centre of Excellence in Moleculary Systems Immunology and Physiology Research (2012-2017) Grant 135320)Seventh Framework Programme (European Commission) (Grant EC-FP7-SYBILLA-201106)EU ERASysBio ERA-NETSigrid Juslius FoundationFICS Graduate Schoo

Host proteostasis modulates influenza evolution

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host-pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance.National Institutes of Health (U.S.) (Award 1DP2GM119162)National Institutes of Health (U.S.) (Grant P30-ES002109

Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction

Background Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. Results We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. Conclusions These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them. Keywords: RNAseqr; RNA depletion; Illumina; NGS; ABRF; TranscriptomicsNational Cancer Institute (U.S.) (Grant P30-CA14051)National Institute of Environmental Health Sciences (Grant P30-ES002109

Genomic mapping of phosphorothioates reveals partial modification of short consensus sequences

Bacterial phosphorothioate (PT) DNA modifications are incorporated by Dnd proteins A-E and often function with DndF-H as a restriction-modification (R-M) system, as in Escherichia coli B7A. However, bacteria such as Vibrio cyclitrophicus FF75 lack dndF-H, which points to other PT functions. Here we report two novel, orthogonal technologies to map PTs across the genomes of B7A and FF75 with >90% agreement: single molecule, real-time sequencing and deep sequencing of iodine-induced cleavage at PT (ICDS). In B7A, we detect PT on both strands of G[subscript ps]AAC/G[subscript ps]TTC motifs, but with only 12% of 40,701 possible sites modified. In contrast, PT in FF75 occurs as a single-strand modification at C[subscript ps]CA, again with only 14% of 160,541 sites modified. Single-molecule analysis indicates that modification could be partial at any particular genomic site even with active restriction by DndF-H, with direct interaction of modification proteins with GAAC/GTTC sites demonstrated with oligonucleotides. These results point to highly unusual target selection by PT-modification proteins and rule out known R-M mechanisms.National Natural Science Foundation (China)Ministry of Science and Technology of the People's Republic of China (973 and 863 Programs)Shanghai Municipal Council of Science and Technology. Shanghai Pujiang ProgramNational Science Foundation (U.S.) (CHE-1019990)National Institute of Environmental Health Sciences (ES002109)Singapore. National Research Foundation (Singapore-MIT Alliance for Research and Technology

Braveheart, a Long Noncoding RNA Required for Cardiovascular Lineage Commitment

Long noncoding RNAs (lncRNAs) are often expressed in a development-specific manner, yet little is known about their roles in lineage commitment. Here, we identified Braveheart (Bvht), a heart-associated lncRNA in mouse. Using multiple embryonic stem cell (ESC) differentiation strategies, we show that Bvht is required for progression of nascent mesoderm toward a cardiac fate. We find that Bvht is necessary for activation of a core cardiovascular gene network and functions upstream of mesoderm posterior 1 (MesP1), a master regulator of a common multipotent cardiovascular progenitor. We also show that Bvht interacts with SUZ12, a component of polycomb-repressive complex 2 (PRC2), during cardiomyocyte differentiation, suggesting that Bvht mediates epigenetic regulation of cardiac commitment. Finally, we demonstrate a role for Bvht in maintaining cardiac fate in neonatal cardiomyocytes. Together, our work provides evidence for a long noncoding RNA with critical roles in the establishment of the cardiovascular lineage during mammalian development.Damon Runyon Cancer Research Foundation (DRG 2032-09)Damon Runyon Cancer Research Foundation (DFS 04-12)European Molecular Biology Organization (Long-term Fellowship)National Heart, Lung, and Blood Institute. Bench to Bassinet Program (U01HL098179)National Heart, Lung, and Blood Institute. Bench to Bassinet Program (U01HL098188)Smith Family FoundationPew Charitable Trusts. Program in the Biomedical Science

Emergent mechanical control of vascular morphogenesis

Vascularization is driven by morphogen signals and mechanical cues that coordinately regulate cellular force generation, migration, and shape change to sculpt the developing vascular network. However, it remains unclear whether developing vasculature actively regulates its own mechanical properties to achieve effective vascularization. We engineered tissue constructs containing endothelial cells and fibroblasts to investigate the mechanics of vascularization. Tissue stiffness increases during vascular morphogenesis resulting from emergent interactions between endothelial cells, fibroblasts, and ECM and correlates with enhanced vascular function. Contractile cellular forces are key to emergent tissue stiffening and synergize with ECM mechanical properties to modulate the mechanics of vascularization. Emergent tissue stiffening and vascular function rely on mechanotransduction signaling within fibroblasts, mediated by YAP1. Mouse embryos lacking YAP1 in fibroblasts exhibit both reduced tissue stiffness and develop lethal vascular defects. Translating our findings through biology-inspired vascular tissue engineering approaches will have substantial implications in regenerative medicine

Transcriptional Analysis of Murine Macrophages Infected with Different Toxoplasma Strains Identifies Novel Regulation of Host Signaling Pathways

Most isolates of Toxoplasma from Europe and North America fall into one of three genetically distinct clonal lineages, the type I, II and III lineages. However, in South America these strains are rarely isolated and instead a great variety of other strains are found. T. gondii strains differ widely in a number of phenotypes in mice, such as virulence, persistence, oral infectivity, migratory capacity, induction of cytokine expression and modulation of host gene expression. The outcome of toxoplasmosis in patients is also variable and we hypothesize that, besides host and environmental factors, the genotype of the parasite strain plays a major role. The molecular basis for these differences in pathogenesis, especially in strains other than the clonal lineages, remains largely unexplored. Macrophages play an essential role in the early immune response against T. gondii and are also the cell type preferentially infected in vivo. To determine if non-canonical Toxoplasma strains have unique interactions with the host cell, we infected murine macrophages with 29 different Toxoplasma strains, representing global diversity, and used RNA-sequencing to determine host and parasite transcriptomes. We identified large differences between strains in the expression level of known parasite effectors and large chromosomal structural variation in some strains. We also identified novel strain-specifically regulated host pathways, including the regulation of the type I interferon response by some atypical strains. IFNβ production by infected cells was associated with parasite killing, independent of interferon gamma activation, and dependent on endosomal Toll-like receptors in macrophages and the cytoplasmic receptor retinoic acid-inducible gene 1 (RIG-I) in fibroblasts.National Institutes of Health (U.S.) (R01-AI080621)New England Regional Center of Excellence for Biodefense and Emerging Infectious Diseases (Developmental Grant AIO57159)Pew Charitable Trusts (Biomedical Scholars Program)Robert A. Swanson Career Development awardThe Knights Templar Eye Foundation, Inc.Pre-Doctoral Grant in the Biological Sciences (5-T32-GM007287-33)Cleo and Paul Schimmel Foundatio