Search

Somatic cell reprogramming and oncogenic transformation share surprisingly similar features, yet transformed cells are highly resistant to reprogramming. There must be barriers that block transformed cells from reprogramming, but the nature of those barriers is unclear. In this study, we generated a systematic panel of transformed mouse embryonic fibroblasts (MEFs) using a variety of oncogenic transgenes, and discovered transformed cell lines that remain compatible with reprogramming when transfected withOct4/Sox2/Klf4/Myc. By comparing the reprogramming-capable and incapable transformed lines we identified multiple stages of failure in the reprogramming process. Some transformed lines failed very early, whilst other lines seemed to progress through a normal-looking reprogramming process. Finally, we show that MEK inhibition overcomes one critical reprogramming barrier by indirectly suppressing a hyperactive epigenetic state in some of the transformed cells. This study reveals that the barriers underlying resistance to reprogramming vary between the different transformation methods.Key findingsSomatic cell reprogramming of transformed cells is context-specificInhibition of MEK converts some cell lines to reprogramming-capableTransformed cell lines are characterized by a hyperactive chromatin stateMEK inhibition indirectly affects chromatin to enable reprogramming

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

Capybara (Hydrochoerus hydrochaeri) is the largest species among the extant rodents. The draft genome of capybara was sequenced with the estimated genome size of 2.6 Gb. Although capybara is about 60 times larger than guinea pig, comparative analyses revealed that the neutral evolutionary rates of the two species were not substantially different. However, analyses of 39 mammalian genomes revealed very heterogeneous evolutionary rates. The highest evolutionary rate, 8.5 times higher than the human rate, was found in the Cricetidae–Muridae common ancestor after the divergence of Spalacidae. Muridae, the family with the highest number of species among mammals, emerged after the rate acceleration. Factors responsible for the evolutionary rate heterogeneity were investigated through correlations between the evolutionary rate and longevity, gestation length, litter frequency, litter size, body weight, generation interval, age at maturity, and taxonomic order. The regression analysis of these factors showed that the model with three factors (taxonomic order, generation interval, and litter size) had the highest predictive power (R2 = 0.74). These three factors determine the number of meiosis per unit time. We also conducted transcriptome analysis and found that the evolutionary rate dynamics affects the evolution of gene expression patterns.

Species of the genus Pagurus have diversified into a wide variety of marine habitats across the world. Despite their worldwide abundance, the genus diversity and biogeographical relationship are relatively less understood at species-level. We evaluated the phylogenetic relationship and genetic diversity among the Pagurus species based on publicly available mitochondrial and nuclear markers. While independent analyses of different markers allowed for larger coverage of taxa and produced largely consistent results, the concatenation of 16S and COI partial sequences led to higher confidence in the phylogenetic relationships. Our analyses established several monophyletic species clusters, substantially corresponding to the previously established morphology-based species groups. The comprehensive species inclusion in the molecular phylogeny resolved the taxonomic position of a number of recently described species that had not been assigned to any morpho-group. In mitochondrial markers-based phylogenies, the “Provenzanoi” group was identified as the basal lineage of Pagurus. The divergence time estimation of the major groups of Pagurus revealed that the Pacific species originated and diversified from the Atlantic lineages around 25–51 MYA. The molecular results suggested a higher inter-regional species diversity and complex phylogenetic relationships within the diverse and heterogeneous members of the genus Pagurus. The study presents a comprehensive snapshot of the diversity of pagurid hermit crabs across multiple geographic regions.

BackgroundInvestigating the functions and activities of genes requires proper annotation of the transcribed units. However, transcript assembly efforts have produced a surprisingly large variation in the number of transcripts, and especially so for noncoding transcripts. The heterogeneity of the assembled transcript sets might be partially explained by sequencing depth.ResultsHere, we used real and simulated short-read sequencing data as well as long-read data to systematically investigate the impact of sequencing depths on the accuracy of assembled transcripts. We assembled and analyzed transcripts from 671 human short-read data sets and four long-read data sets. At the first level, there is a positive correlation between the number of reads and the number of recovered transcripts. However, the effect of the sequencing depth varied based on cell or tissue type, the type of read considered and the nature and expression levels of the transcripts. The detection of coding transcripts saturated rapidly for both short-read and long-reads, however, there was no sign of saturation for noncoding transcripts at any sequencing depth. Increasing long-read sequencing depth specifically benefited transcripts containing transposable elements. Finally, we show how single-cell RNA-seq can be guided by transcripts assembled from bulk long-read samples, and demonstrate that noncoding transcripts are expressed at similar levels to coding transcripts but are expressed in fewer cells.ConclusionsThis study shows the impact of sequencing depth on transcript assembly. Sequencing read depth has a relatively minor impact on coding transcript assembly, but a major effect on the assembly of noncoding transcripts. This study highlights important factors to consider when deciding the sequencing read depths to be used for transcript assembly.

Around 60% of in vitro fertilized (IVF) human embryos irreversibly arrest before compaction between the 3- to 8-cell stage, posing a significant clinical problem. The mechanisms behind this arrest are unclear. Here, we show that the arrested embryos enter a senescent-like state, marked by cell cycle arrest, the down-regulation of ribosomes and histones and down-regulation of MYC and p53 activity. The arrested embryos can be divided into 3 types. Type I embryos fail to complete the maternal-zygotic transition, and Type II/III embryos have low levels of glycolysis and either high (Type II) or low (Type III) levels of oxidative phosphorylation. Treatment with the SIRT agonist resveratrol or nicotinamide riboside (NR) can partially rescue the arrested phenotype, which is accompanied by changes in metabolic activity. Overall, our data suggests metabolic and epigenetic dysfunctions underlie the arrest of human embryos.

Around 60% of in vitro fertilized (IVF) human embryos irreversibly arrest before compaction between the 3-8-cell stage, posing a significant clinical problem. The mechanisms behind this arrest are unclear. Here, we show that the arrested embryos enter a quiescent-like state, marked by cell cycle arrest, the downregulation of ribosomes and histones and downregulation of MYC and p53 activity. Mechanistically, the arrested embryos can be divided into three types. Type I embryos fail to complete the maternal-zygotic transition, and type II/III embryos have erroneously low levels of glycolysis and variable levels of oxidative phosphorylation. Treatment with resveratrol or nicotinamide riboside (NR) can partially rescue the arrested phenotype. The mechanism of reactivation involves the upregulation of SIRT1, and activation of glycolysis and fatty acid oxidation which forces the embryos out of a quiescent state. Overall, our data reveal how human embryo arrest can be overcome by modulating metabolic pathways.

To satisfy both the safety and rapidity of glycemic control, muscles' insulin response must be bistable, as theoretically predicted. Here, we test the bistability hypothesis by combining cellular experiments (to measure the threshold values

CTCF plays key roles in gene regulation, chromatin insulation, imprinting, X chromosome inactivation and organizing the higher-order chromatin architecture of mammalian genomes. Previous studies have mainly focused on the roles of the canonical CTCF isoform. Here, we explore the functions of an alternatively spliced human CTCF isoform in which exons 3 and 4 are skipped, producing a shorter isoform (CTCF-s). Functionally, we find that CTCF-s competes with the genome binding of canonical CTCF and binds a similar DNA sequence. CTCF-s binding disrupts CTCF/cohesin binding, alters CTCF-mediated chromatin looping and promotes the activation of IFI6 that leads to apoptosis. This effect is caused by an abnormal long-range interaction at the IFI6 enhancer and promoter. Taken together, this study reveals a non-canonical function for CTCF-s that antagonizes the genomic binding of canonical CTCF and cohesin, and that modulates chromatin looping and causes apoptosis by stimulating IFI6 expression., CTCF plays key roles in gene regulation, chromatin insulation and organizing the higher-order chromatin architecture of mammalian genomes. Here the authors investigate the function an alternatively spliced shorter CTCF isoform, finding that this isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis.

High-throughput sequencing technologies are increasingly used in molecular cell biology to assess genome-wide chromatin dynamics of proteins bound to DNA, through techniques such as chromatin immunoprecipitation sequencing (ChIP-seq). These techniques often rely on an analysis strategy based on identifying genomic regions with increased sequencing signal to infer the binding location or chemical modifications of proteins bound to DNA. Peak calling within individual samples has been well described, however relatively little attention has been devoted to the merging of replicate samples, and the cross-comparison of many samples. Here, we present a generalized strategy to enable the unification of ChIP-seq datasets, enabling enhanced cross-comparison of binding patterns. The strategy works by merging peak data between different (even unrelated) samples, and then using a local background to recalculate enrichment. This strategy redefines the peaks within each experiment, allowing for more accurate cross-comparison of datasets.

High-throughput sequencing technologies are increasingly used in molecular cell biology to assess genome-wide chromatin dynamics of proteins bound to DNA, through techniques such as chromatin immunoprecipitation sequencing (ChIP-seq). These techniques often rely on an analysis strategy based on identifying genomic regions with increased sequencing signal to infer the binding location or chemical modifications of proteins bound to DNA. Peak calling within individual samples has been well described, however relatively little attention has been devoted to the merging of replicate samples, and the cross-comparison of many samples. Here, we present a generalized strategy to enable the unification of ChIP-seq datasets, enabling enhanced cross-comparison of binding patterns. The strategy works by merging peak data between different (even unrelated) samples, and then using a local background to recalculate enrichment. This strategy redefines the peaks within each experiment, allowing for more accurate cross-comparison of datasets.

Transposable elements (TEs) occupy nearly 50% of mammalian genomes and are both potential dangers to genome stability and functional genetic elements. TEs can be expressed and exonised as part of a transcript, however, their full contribution to the transcript splicing remains unresolved. Here, guided by long and short read sequencing of RNAs, we show that 26% of coding and 65% of non-coding transcripts of human pluripotent stem cells (hPSCs) contain TEs. Different TE families have unique integration patterns with diverse consequences on RNA expression and function. We identify hPSC-specific splicing of endogenous retroviruses (ERVs) as well as LINE L1 elements into protein coding genes that generate TE-derived peptides. Finally, single cell RNA-seq reveals that proliferating hPSCs are dominated by ERV-containing transcripts, and subpopulations express SINE or LINE-containing transcripts. Overall, we demonstrate that TE splicing modulates the pluripotency transcriptome by enhancing and impairing transcript expression and generating novel transcripts and peptides.

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity and fate in unclear ways. Single cell-sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE activity has not been addressed. Here, we developed a single-cell TE processing pipeline, scTE, and report the activity of TEs in single cells in a range of biological contexts. Specific TE types were expressed in subpopulations of embryonic stem cells and were dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs were expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we applied scTE to single cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results reveal the dynamic patterns of TEs in single cells and their contributions to cell fate and heterogeneity.

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity., How transposable elements (TE) contribute to cell fate changes is unclear. Here, the authors generate a pipeline to quantify TE expression from single cell data. They show the dynamic expression of TEs from gastrulation to somatic cell reprogramming and human disease

β-Catenin recruits BRD4 and other coregulators to protect pluripotency gene transcription against network perturbation., Mouse embryonic stem cells cultured with MEK (mitogen-activated protein kinase kinase) and GSK3 (glycogen synthase kinase 3) inhibitors (2i) more closely resemble the inner cell mass of preimplantation blastocysts than those cultured with SL [serum/leukemia inhibitory factor (LIF)]. The transcriptional mechanisms governing this pluripotent ground state are unresolved. Release of promoter-proximal paused RNA polymerase II (Pol2) is a multistep process necessary for pluripotency and cell cycle gene transcription in SL. We show that β-catenin, stabilized by GSK3 inhibition in medium with 2i, supplies transcriptional coregulators at pluripotency loci. This selectively strengthens pluripotency loci and renders them addicted to transcription initiation for productive gene body elongation in detriment to Pol2 pause release. By contrast, cell cycle genes are not bound by β-catenin, and proliferation/self-renewal remains tightly controlled by Pol2 pause release under 2i conditions. Our findings explain how pluripotency is reinforced in the ground state and also provide a general model for transcriptional resilience/adaptation upon network perturbation in other contexts.

Summary Cells are generally resistant to cell type conversions, but can be converted by the application of growth factors, chemical inhibitors and ectopic expression of genes. However, it remains difficult to accurately identify the destination cell type or differentiation bias when these techniques are used to alter cell type. Consequently, there is demand for computational techniques that can help researchers understand both the cell type and differentiation bias. While advanced tools identifying cell types exist for single cell data and the deconvolution of mixed cell populations, the problem of exploring partially differentiated cells of indeterminate transcriptional identity has not been addressed. To fill this gap, we developed driver-predictor, which relies on scoring per gene transcriptional similarity between RNA-Seq datasets to reveal directional bias of differentiation. By comparing against large cell type transcriptome libraries or a desired target expression profile, the tool enables the user to visualize both the changes in transcriptional identity as well as the genes accounting for the cell type changes. This software will be a powerful tool for researchers to explore in vitro experiments that involve cell type conversions. Availability and implementation Source code is open source under the MIT license and is freely available on https://github.com/LoaloaF/DPre. Supplementary information Supplementary data are available at Bioinformatics online.

The interplay between the Yamanaka factors (OCT4, SOX2, KLF4 and c-MYC) and transcriptional/epigenetic co-regulators in somatic cell reprogramming is incompletely understood. Here, we demonstrate that the histone H3 lysine 27 trimethylation (H3K27me3) demethylase JMJD3 plays conflicting roles in mouse reprogramming. On one side, JMJD3 induces the pro-senescence factor Ink4a and degrades the pluripotency regulator PHF20 in a reprogramming factor-independent manner. On the other side, JMJD3 is specifically recruited by KLF4 to reduce H3K27me3 at both enhancers and promoters of epithelial and pluripotency genes. JMJD3 also promotes enhancer-promoter looping through the cohesin loading factor NIPBL and ultimately transcriptional elongation. This competition of forces can be shifted towards improved reprogramming by using early passage fibroblasts or boosting JMJD3’s catalytic activity with vitamin C. Our work, thus, establishes a multifaceted role for JMJD3, placing it as a key partner of KLF4 and a scaffold that assists chromatin interactions and activates gene transcription., Previous work suggested that histone demethylase JMJD3 is detrimental to somatic cell reprogramming. Here, the authors show that while JMJD3 has a context-independent detrimental effect on early stages of reprogramming, during late stages it activates epithelial and pluripotency genes together with Klf4.

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs., Transposable elements (TEs) fulfill essential but poorly understood roles in genome organization and gene expression control. Here the authors show that the regulation of TEs occurs through overlapping epigenetic mechanisms that control the expression and chromatin signatures at TEs.

Family Hominidae, which includes humans and great apes, is recognized for unique complex social behavior and intellectual abilities. Despite the increasing genome data, however, the genomic origin of its phenotypic uniqueness has remained elusive. Clade-specific genes and highly conserved noncoding sequences (HCNSs) are among the high-potential evolutionary candidates involved in driving clade-specific characters and phenotypes. On this premise, we analyzed whole genome sequences along with gene orthology data retrieved from major DNA databases to find Hominidae-specific (HS) genes and HCNSs. We discovered that Down syndrome critical region 4 (DSCR4) is the only experimentally verified gene uniquely present in Hominidae. DSCR4 has no structural homology to any known protein and was inferred to have emerged in several steps through LTR/ERV1, LTR/ERVL retrotransposition, and transversion. Using the genomic distance as neutral evolution threshold, we identified 1,658 HS HCNSs. Polymorphism coverage and derived allele frequency analysis of HS HCNSs showed that these HCNSs are under purifying selection, indicating that they may harbor important functions. They are overrepresented in promoters/untranslated regions, in close proximity of genes involved in sensory perception of sound and developmental process, and also showed a significantly lower nucleosome occupancy probability. Interestingly, many ancestral sequences of the HS HCNSs showed very high evolutionary rates. This suggests that new functions emerged through some kind of positive selection, and then purifying selection started to operate to keep these functions.

Experimental studies have found the involvement of certain conserved noncoding sequences (CNSs) in the regulation of the proximal protein-coding genes in mammals. However, reported cases of long range enhancer activities and inter-chromosomal regulation suggest that proximity of CNSs to protein-coding genes might not be important for regulation. To test the importance of the CNS genomic location, we extracted the CNSs conserved between chicken and four mammalian species (human, mouse, dog, and cattle). These CNSs were confirmed to be under purifying selection. The intergenic CNSs are often found in clusters in gene deserts, where protein-coding genes are in paucity. The distribution pattern, ChIP-Seq, and RNA-Seq data suggested that the CNSs are more likely to be regulatory elements and not corresponding to long intergenic noncoding RNAs. Physical distances between CNS and their nearest protein coding genes were well conserved between human and mouse genomes, and CNS-flanking genes were often found in evolutionarily conserved genomic neighborhoods. ChIP-Seq signal and gene expression patterns also suggested that CNSs regulate nearby genes. Interestingly, genes with more CNSs have more evolutionarily conserved expression than those with fewer CNSs. These computationally obtained results suggest that the genomic locations of CNSs are important for their regulatory functions. In fact, various kinds of evolutionary constraints may be acting to maintain the genomic locations of CNSs and protein-coding genes in mammals to ensure proper regulation.

Cynomolgus macaques (Macaca fascicularis) were introduced to the island of Mauritius by humans around the 16th century. The unique demographic history of the Mauritian cynomolgus macaques provides the opportunity to not only examine the genetic background of well-established nonhuman primates for biomedical research but also understand the effect of an extreme population bottleneck on the pattern of polymorphisms in genomes. We sequenced the whole genomes of six Mauritian cynomolgus macaques and obtained an average of 20-fold coverage of the genome sequences for each individual. The overall level of nucleotide diversity was 23% smaller than that of the Malaysian cynomolgus macaques, and a reduction of low-frequency polymorphisms was observed. In addition, we also confirmed that the Mauritian cynomolgus macaques were genetically closer to a representative of the Malaysian population than to a representative of the Indochinese population. Excess of nonsynonymous polymorphisms in low frequency, which has been observed in many other species, was not very strong in the Mauritian samples, and the proportion of heterozygous nonsynonymous polymorphisms relative to synonymous polymorphisms is higher within individuals in Mauritian than Malaysian cynomolgus macaques. Those patterns indicate that the extreme population bottleneck made purifying selection overwhelmed by the power of genetic drift in the population. Finally, we estimated the number of founding individuals by using the genome-wide site frequency spectrum of the six samples. Assuming a simple demographic scenario with a single bottleneck followed by exponential growth, the estimated number of founders (∼20 individuals) is largely consistent with previous estimates.

Conserved noncoding sequences (CNSs) of vertebrates are considered to be closely linked with protein-coding gene regulatory functions. We examined the abundance and genomic distribution of CNSs in four mammalian orders: primates, rodents, carnivores, and cetartiodactyls. We defined the two thresholds for CNS using conservation level of coding genes; using all the three coding positions and using only first and second codon positions. The abundance of CNSs varied among lineages, with primates and rodents having highest and lowest number of CNSs, respectively, whereas carnivores and cetartiodactyls had intermediate values. These CNSs cover 1.3–5.5% of the mammalian genomes and have signatures of selective constraints that are stronger in more ancestral than the recent ones. Evolution of new CNSs as well as retention of ancestral CNSs contribute to the differences in abundance. The genomic distribution of CNSs is dynamic with higher proportions of rodent and primate CNSs located in the introns compared with carnivores and cetartiodactyls. In fact, 19% of orthologous single-copy CNSs between human and dog are located in different genomic regions. If CNSs can be considered as candidates of gene expression regulatory sequences, heterogeneity of CNSs among the four mammalian orders may have played an important role in creating the order-specific phenotypes. Fewer CNSs in rodents suggest that rodent diversity is related to lower regulatory conservation. With CNSs shown to cluster around genes involved in nervous systems and the higher number of primate CNSs, our result suggests that CNSs may be involved in the higher complexity of the primate nervous system.