Epigenomic analysis of multilineage differentiation of human embryonic stem cells.

Epigenetic mechanisms have been proposed to play crucial roles in mammalian development, but their precise functions are only partially understood. To investigate epigenetic regulation of embryonic development, we differentiated human embryonic stem cells into mesendoderm, neural progenitor cells, trophoblast-like cells, and mesenchymal stem cells and systematically characterized DNA methylation, chromatin modifications, and the transcriptome in each lineage. We found that promoters that are active in early developmental stages tend to be CG rich and mainly engage H3K27me3 upon silencing in nonexpressing lineages. By contrast, promoters for genes expressed preferentially at later stages are often CG poor and primarily employ DNA methylation upon repression. Interestingly, the early developmental regulatory genes are often located in large genomic domains that are generally devoid of DNA methylation in most lineages, which we termed DNA methylation valleys (DMVs). Our results suggest that distinct epigenetic mechanisms regulate early and late stages of ES cell differentiation.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography.

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

MicroRNA-145 regulates OCT4, SOX2, and KLF4 and represses pluripotency in human embryonic stem cells.

MicroRNAs (miRNAs) are posttranscriptional modulators of gene expression and play an important role in many developmental processes. We report here that expression of microRNA-145 (miR-145) is low in self-renewing human embryonic stem cells (hESCs) but highly upregulated during differentiation. We identify the pluripotency factors OCT4, SOX2, and KLF4 as direct targets of miR-145 and show that endogenous miR-145 represses the 3' untranslated regions of OCT4, SOX2, and KLF4. Increased miR-145 expression inhibits hESC self-renewal, represses expression of pluripotency genes, and induces lineage-restricted differentiation. Loss of miR-145 impairs differentiation and elevates OCT4, SOX2, and KLF4. Furthermore, we find that the miR-145 promoter is bound and repressed by OCT4 in hESCs. This work reveals a direct link between the core reprogramming factors and miR-145 and uncovers a double-negative feedback loop involving OCT4, SOX2, KLF4, and miR-145.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization.

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Serial KinderMiner (SKiM) discovers and annotates biomedical knowledge using co-occurrence and transformer models.

BACKGROUND: The PubMed archive contains more than 34 million articles; consequently, it is becoming increasingly difficult for a biomedical researcher to keep up-to-date with different knowledge domains. Computationally efficient and interpretable tools are needed to help researchers find and understand associations between biomedical concepts. The goal of literature-based discovery (LBD) is to connect concepts in isolated literature domains that would normally go undiscovered. This usually takes the form of an A-B-C relationship, where A and C terms are linked through a B term intermediate. Here we describe Serial KinderMiner (SKiM), an LBD algorithm for finding statistically significant links between an A term and one or more C terms through some B term intermediate(s). The development of SKiM is motivated by the observation that there are only a few LBD tools that provide a functional web interface, and that the available tools are limited in one or more of the following ways: (1) they identify a relationship but not the type of relationship, (2) they do not allow the user to provide their own lists of B or C terms, hindering flexibility, (3) they do not allow for querying thousands of C terms (which is crucial if, for instance, the user wants to query connections between a disease and the thousands of available drugs), or (4) they are specific for a particular biomedical domain (such as cancer). We provide an open-source tool and web interface that improves on all of these issues. RESULTS: We demonstrate SKiM's ability to discover useful A-B-C linkages in three control experiments: classic LBD discoveries, drug repurposing, and finding associations related to cancer. Furthermore, we supplement SKiM with a knowledge graph built with transformer machine-learning models to aid in interpreting the relationships between terms found by SKiM. Finally, we provide a simple and intuitive open-source web interface ( https://skim.morgridge.org ) with comprehensive lists of drugs, diseases, phenotypes, and symptoms so that anyone can easily perform SKiM searches. CONCLUSIONS: SKiM is a simple algorithm that can perform LBD searches to discover relationships between arbitrary user-defined concepts. SKiM is generalized for any domain, can perform searches with many thousands of C term concepts, and moves beyond the simple identification of an existence of a relationship; many relationships are given relationship type labels from our knowledge graph.

Pollen carryover, pollinator movement, and spatial context impact the delivery of pollination services in apple orchards.

Assessing the relative contributions of different pollinator taxa to pollination services is a central task in both basic eco-evolutionary research and applied conservation and agriculture. To that end, many studies have quantified single-visit pollen deposition and visitation frequency, which together determine a pollinator species' rate of conspecific pollen delivery. However, for plant species that require or benefit from outcrossing, pollination service quality further depends upon the ratio of outcross to self-pollen deposited, which is determined by two additional pollinator traits: pollen carryover and movement patterns among genetically compatible plant individuals. Here, we compare the pollination capacities of managed honey bees, native bumble bees, and native mining bees in apple-a varietally self-incompatible commercial crop-when pollen carryover and pollinator movement patterns are considered. We constructed simulation models of outcross pollen deposition parameterized using empirically measured single-visit pollen deposition, visitation frequency, and probabilities of intertree movement exhibited by each pollinator type, as well as pollen carryover patterns simulated based on parameters reported in the literature. In these models, we also explicitly specified the spatial relationships among cross-compatible trees based on field-realistic orchard layout schemes. We found that estimated pollination service delivery was considerably reduced for all pollinator types when pollen carryover and pollinator movement patterns were considered, as compared to when only single-visit pollen deposition and visitation frequency were considered. We also found that the performance of different pollinator types varied greatly across simulated orchard layout schemes and pollen carryover scenarios, including one instance where bumble and mining bees reversed their relative rankings. In all simulations, native bumble and mining bees outperformed managed honey bees in terms of both outcross pollen delivery per unit time and per flower visited, with disparities being greatest under scenarios of low pollen carryover. We demonstrate the degree to which pollination studies may reach inaccurate conclusions regarding pollination service delivery when pollen carryover and pollinator movement patterns are ignored. Our finding of the strong context dependence of pollination efficiency, even within a single plant-pollinator taxon pair, cautions that future studies in both basic and applied pollination biology should explicitly consider the ecological context in which pollination interactions take place.

General Approach to Enantiopure 1-Aminopyrrolizidines: Application to the Asymmetric Synthesis of the Loline Alkaloids.

The synthesis of a range of loline alkaloids is reported. The C(7) and C(7a) stereogenic centers for the targets were formed by the established conjugate addition of lithium (S)-N-benzyl-N-(α-methylbenzyl)amide to tert-butyl 5-benzyloxypent-2-enoate, ensuing enolate oxidation to give an α-hydroxy-ß-amino ester, and then formal exchange of the resultant amino and hydroxyl functionalities (via the intermediacy of the corresponding aziridinium ion) to give an α-amino-ß-hydroxy ester. Subsequent transformation gave a 3-hydroxyprolinal derivative which was converted to the corresponding N-tert-butylsulfinylimine. Mannich-type reaction with the enolate derived from O-Boc protected methyl glycolate then formed the remaining C(1) and C(2) stereogenic centers for the targets. The 2,7-ether bridge was formed by a displacement reaction, completing construction of the loline alkaloid core. Facile manipulations then gave a range of loline alkaloids, including loline itself.

Ronin is essential for embryogenesis and the pluripotency of mouse embryonic stem cells.

Pluripotency is a unique biological state that allows cells to differentiate into any tissue type. Here we describe a candidate pluripotency factor, Ronin, that possesses a THAP domain, which is associated with sequence-specific DNA binding and epigenetic silencing of gene expression. Ronin is expressed primarily during the earliest stages of murine embryonic development, and its deficiency in mice produces periimplantational lethality and defects in the inner cell mass. Conditional knockout of Ronin prevents the growth of ES cells while forced expression of Ronin allows ES cells to proliferate without differentiation under conditions that normally do not promote self-renewal. Ectopic expression also partly compensates for the effects of Oct4 knockdown. We demonstrate that Ronin binds directly to HCF-1, a key transcriptional regulator. Our findings identify Ronin as an essential factor underlying embryogenesis and ES cell pluripotency. Its association with HCF-1 suggests an epigenetic mechanism of gene repression in pluripotent cells.

A haplotype-resolved genome assembly of the Nile rat facilitates exploration of the genetic basis of diabetes.

BACKGROUND: The Nile rat (Avicanthis niloticus) is an important animal model because of its robust diurnal rhythm, a cone-rich retina, and a propensity to develop diet-induced diabetes without chemical or genetic modifications. A closer similarity to humans in these aspects, compared to the widely used Mus musculus and Rattus norvegicus models, holds the promise of better translation of research findings to the clinic. RESULTS: We report a 2.5 Gb, chromosome-level reference genome assembly with fully resolved parental haplotypes, generated with the Vertebrate Genomes Project (VGP). The assembly is highly contiguous, with contig N50 of 11.1 Mb, scaffold N50 of 83 Mb, and 95.2% of the sequence assigned to chromosomes. We used a novel workflow to identify 3613 segmental duplications and quantify duplicated genes. Comparative analyses revealed unique genomic features of the Nile rat, including some that affect genes associated with type 2 diabetes and metabolic dysfunctions. We discuss 14 genes that are heterozygous in the Nile rat or highly diverged from the house mouse. CONCLUSIONS: Our findings reflect the exceptional level of genomic resolution present in this assembly, which will greatly expand the potential of the Nile rat as a model organism.

Rapid changes in chromatin structure during dedifferentiation of primary hepatocytes in vitro.

Primary hepatocytes are widely used in the pharmaceutical industry to screen drug candidates for hepatotoxicity, but hepatocytes quickly dedifferentiate and lose their mature metabolic function in culture. Attempts have been made to better recapitulate the in vivo liver environment in culture, but the full spectrum of signals required to maintain hepatocyte function ex vivo remains elusive. To elucidate molecular changes that accompany, and may contribute to dedifferentiation of hepatocytes ex vivo, we performed lineage tracing and comprehensive profiling of alterations in their gene expression profiles and chromatin landscape during culture. First, using genetically tagged hepatocytes we demonstrate that expression of the fetal gene alpha-fetoprotein in cultured hepatocytes comes from cells that previously expressed the mature gene albumin, and not from a population of albumin-negative precursor cells, proving mature hepatocytes undergo true dedifferentiation in culture. Next we studied the dedifferentiation process in detail through bulk RNA-sequencing of hepatocytes cultured over an extended period. We identified three distinct phases of dedifferentiation: an early phase, where mature hepatocyte genes are rapidly downregulated in a matter of hours; a middle phase, where fetal genes are activated; and a late phase, where initially rare contaminating non-parenchymal cells proliferate, taking over the culture. Lastly, to better understand the signaling events that result in the rapid downregulation of mature genes in hepatocytes, we examined changes in chromatin accessibility in these cells during the first 24 h of culture using Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). We find that drastic and rapid changes in chromatin accessibility occur immediately upon the start of culture. Using binding motif analysis of the areas of open chromatin sharing similar temporal profiles, we identify several candidate transcription factors potentially involved in the dedifferentiation of primary hepatocytes in culture.

Simulated driving in the epilepsy monitoring unit: Effects of seizure type, consciousness, and motor impairment.

People with epilepsy face serious driving restrictions, determined using retrospective studies. To relate seizure characteristics to driving impairment, we aimed to study driving behavior during seizures with a simulator. Patients in the Yale New Haven Hospital undergoing video-electroencephalographic monitoring used a laptop-based driving simulator during ictal events. Driving function was evaluated by video review and analyzed in relation to seizure type, impairment of consciousness/responsiveness, or motor impairment during seizures. Fifty-one seizures in 30 patients were studied. In terms of seizure type, we found that focal to bilateral tonic-clonic or myoclonic seizures (5/5) and focal seizures with impaired consciousness/responsiveness (11/11) always led to driving impairment; focal seizures with spared consciousness/responsiveness (0/10) and generalized nonmotor (generalized spike-wave bursts; 1/19) usually did not lead to driving impairment. Regardless of seizure type, we found that seizures with impaired consciousness (15/15) or with motor involvement (13/13) always led to impaired driving, but those with spared consciousness (0/20) or spared motor function (5/38) usually did not. These results suggest that seizure types with impaired consciousness/responsiveness and abnormal motor function contribute to impaired driving. Expanding this work in a larger cohort could further determine how results with a driving simulator may translate into real world driving safety.

Time-dependent Pax3-mediated chromatin remodeling and cooperation with Six4 and Tead2 specify the skeletal myogenic lineage in developing mesoderm.

The transcriptional mechanisms driving lineage specification during development are still largely unknown, as the interplay of multiple transcription factors makes it difficult to dissect these molecular events. Using a cell-based differentiation platform to probe transcription function, we investigated the role of the key paraxial mesoderm and skeletal myogenic commitment factors-mesogenin 1 (Msgn1), T-box 6 (Tbx6), forkhead box C1 (Foxc1), paired box 3 (Pax3), Paraxis, mesenchyme homeobox 1 (Meox1), sine oculis-related homeobox 1 (Six1), and myogenic factor 5 (Myf5)-in paraxial mesoderm and skeletal myogenesis. From this study, we define a genetic hierarchy, with Pax3 emerging as the gatekeeper between the presomitic mesoderm and the myogenic lineage. By assaying chromatin accessibility, genomic binding and transcription profiling in mesodermal cells from mouse and human Pax3-induced embryonic stem cells and Pax3-null embryonic day (E)9.5 mouse embryos, we identified conserved Pax3 functions in the activation of the skeletal myogenic lineage through modulation of Hedgehog, Notch, and bone morphogenetic protein (BMP) signaling pathways. In addition, we demonstrate that Pax3 molecular function involves chromatin remodeling of its bound elements through an increase in chromatin accessibility and cooperation with sine oculis-related homeobox 4 (Six4) and TEA domain family member 2 (Tead2) factors. To our knowledge, these data provide the first integrated analysis of Pax3 function, demonstrating its ability to remodel chromatin in mesodermal cells from developing embryos and proving a mechanistic footing for the transcriptional hierarchy driving myogenesis.

Interspecies chimeric conditions affect the developmental rate of human pluripotent stem cells.

Human pluripotent stem cells hold significant promise for regenerative medicine. However, long differentiation protocols and immature characteristics of stem cell-derived cell types remain challenges to the development of many therapeutic applications. In contrast to the slow differentiation of human stem cells in vitro that mirrors a nine-month gestation period, mouse stem cells develop according to a much faster three-week gestation timeline. Here, we tested if co-differentiation with mouse pluripotent stem cells could accelerate the differentiation speed of human embryonic stem cells. Following a six-week RNA-sequencing time course of neural differentiation, we identified 929 human genes that were upregulated earlier and 535 genes that exhibited earlier peaked expression profiles in chimeric cell cultures than in human cell cultures alone. Genes with accelerated upregulation were significantly enriched in Gene Ontology terms associated with neurogenesis, neuron differentiation and maturation, and synapse signaling. Moreover, chimeric mixed samples correlated with in utero human embryonic samples earlier than human cells alone, and acceleration was dose-dependent on human-mouse co-culture ratios. The altered gene expression patterns and developmental rates described in this report have implications for accelerating human stem cell differentiation and the use of interspecies chimeric embryos in developing human organs for transplantation.

Synthesis and Configuration of O-Acetyl Microgrewiapine A: Phantomization of O-Acetyl 6-epi-Microgrewiapine A.

The formation of O-acetyl microgrewiapine A is investigated. NMR data for the authentic sample derived from the natural product are corrected. Wholly synthetic samples, produced from reductive N-methylation of synthetic microcosamine A (to give synthetic microgrewiapine A) followed by O-acetylation, exhibit NMR data that are identical to those of the authentic sample. The previous report that this two-step transformation proceeds with epimerization at C-6 is thus shown to be in error: the purported sample of O-acetyl 6-epi-microgrewiapine A is structurally misassigned and is, in fact, O-acetyl microgrewiapine A. A plausible rationale for the structural misassignment is advanced.

Microgrewiapine C: Asymmetric Synthesis, Spectroscopic Data, and Configuration Assignment.

The first asymmetric synthesis of microgrewiapine C, a piperidine alkaloid isolated from Microcos paniculata, is reported. This synthesis prompted correction of the 1H and 13C NMR data for the natural sample of the alkaloid, which was achieved by reanalysis of the original spectra. The corrected data for the natural product were found to be identical to those of the synthetic sample prepared herein, thus confirming the structural and relative configurational assignment of microgrewiapine C. Although comparison of specific rotation values indicates that the (1R,2S,3S,6S) absolute configuration should be assigned to the alkaloid, consideration of potential common biosynthetic origins of microgrewiapine C and congeners suggests that further phytochemical investigations are warranted.

Anthropogenic disturbance in tropical forests can double biodiversity loss from deforestation.

Concerted political attention has focused on reducing deforestation, and this remains the cornerstone of most biodiversity conservation strategies. However, maintaining forest cover may not reduce anthropogenic forest disturbances, which are rarely considered in conservation programmes. These disturbances occur both within forests, including selective logging and wildfires, and at the landscape level, through edge, area and isolation effects. Until now, the combined effect of anthropogenic disturbance on the conservation value of remnant primary forests has remained unknown, making it impossible to assess the relative importance of forest disturbance and forest loss. Here we address these knowledge gaps using a large data set of plants, birds and dung beetles (1,538, 460 and 156 species, respectively) sampled in 36 catchments in the Brazilian state of Pará. Catchments retaining more than 69­80% forest cover lost more conservation value from disturbance than from forest loss. For example, a 20% loss of primary forest, the maximum level of deforestation allowed on Amazonian properties under Brazil's Forest Code, resulted in a 39­54% loss of conservation value: 96­171% more than expected without considering disturbance effects. We extrapolated the disturbance-mediated loss of conservation value throughout Pará, which covers 25% of the Brazilian Amazon. Although disturbed forests retained considerable conservation value compared with deforested areas, the toll of disturbance outside Pará's strictly protected areas is equivalent to the loss of 92,000­139,000 km2 of primary forest. Even this lowest estimate is greater than the area deforested across the entire Brazilian Amazon between 2006 and 2015 (ref. 10). Species distribution models showed that both landscape and within-forest disturbances contributed to biodiversity loss, with the greatest negative effects on species of high conservation and functional value. These results demonstrate an urgent need for policy interventions that go beyond the maintenance of forest cover to safeguard the hyper-diversity of tropical forest ecosystems.

TimeMeter assesses temporal gene expression similarity and identifies differentially progressing genes.

Comparative time series transcriptome analysis is a powerful tool to study development, evolution, aging, disease progression and cancer prognosis. We develop TimeMeter, a statistical method and tool to assess temporal gene expression similarity, and identify differentially progressing genes where one pattern is more temporally advanced than the other. We apply TimeMeter to several datasets, and show that TimeMeter is capable of characterizing complicated temporal gene expression associations. Interestingly, we find: (i) the measurement of differential progression provides a novel feature in addition to pattern similarity that can characterize early developmental divergence between two species; (ii) genes exhibiting similar temporal patterns between human and mouse during neural differentiation are under strong negative (purifying) selection during evolution; (iii) analysis of genes with similar temporal patterns in mouse digit regeneration and axolotl blastema differentiation reveals common gene groups for appendage regeneration with potential implications in regenerative medicine.

Pluripotent stem cell-derived myogenic progenitors remodel their molecular signature upon in vivo engraftment.

Optimal cell-based therapies for the treatment of muscle degenerative disorders should not only regenerate fibers but provide a quiescent satellite cell pool ensuring long-term maintenance and regeneration. Conditional expression of Pax3/Pax7 in differentiating pluripotent stem cells (PSCs) allows the generation of myogenic progenitors endowed with enhanced regenerative capacity. To identify the molecular determinants underlying their regenerative potential, we performed transcriptome analyses of these cells along with primary myogenic cells from several developmental stages. Here we show that in vitro-generated PSC-derived myogenic progenitors possess a molecular signature similar to embryonic/fetal myoblasts. However, compared with fetal myoblasts, following transplantation they show superior myofiber engraftment and ability to seed the satellite cell niche, respond to multiple reinjuries, and contribute to long-term regeneration. Upon engraftment, the transcriptome of reisolated Pax3/Pax7-induced PSC-derived myogenic progenitors changes toward a postnatal molecular signature, particularly in genes involved in extracellular matrix remodeling. These findings demonstrate that Pax3/Pax7-induced myogenic progenitors remodel their molecular signature and functionally mature upon in vivo exposure to the adult muscle environment.

Mutual kinetic resolution: probing enantiorecognition phenomena and screening for kinetic resolution with racemic reagents.

Enantiorecognition between a racemic reagent and a racemic substrate can be a valuable process in organic synthesis. This review highlights representative examples of this phenomenon and the use of mutual kinetic resolution as a method for screening of kinetic and/or parallel kinetic resolutions.