The killifish germline regulates longevity and somatic repair in a sex-specific manner.

Classical evolutionary theories propose tradeoffs among reproduction, damage repair and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. In this study, we used the turquoise killifish (Nothobranchius furzeri) to genetically arrest germline development at discrete stages and examine how different modes of infertility impact life history. We first constructed a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. We show here that germline depletion-but not arresting germline differentiation-enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted Caenorhabditis elegans. Our results, therefore, demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

The killifish germline regulates longevity and somatic repair in a sex-specific manner.

Classical evolutionary theories propose tradeoffs between reproduction, damage repair, and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. Here, we use the turquoise killifish ( N. furzeri ) to genetically arrest germline development at discrete stages, and examine how different modes of infertility impact life-history. We first construct a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. Next, we show that germline depletion - but not arresting germline differentiation - enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted C. elegans . Our results therefore demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

Repetitive DNA symmetry elements negatively regulate gene expression in embryonic stem cells.

Transcription factor (TF) binding to genomic DNA elements constitutes one of the key mechanisms that regulates gene expression program in cells. Both consensus and nonconsensus DNA sequence elements influence the recognition specificity of TFs. Based on the analysis of experimentally determined c-Myc binding preferences to genomic DNA, here we statistically predict that certain repetitive, nonconsensus DNA symmetry elements can relatively reduce TF-DNA binding preferences. This is in contrast to a different set of repetitive, nonconsensus symmetry elements that can increase the strength of TF-DNA binding. Using c-Myc enhancer reporter system containing consensus motif flanked by nonconsensus sequences in embryonic stem cells, we directly demonstrate that the enrichment in such negatively regulating repetitive symmetry elements is sufficient to reduce the gene expression level compared with native genomic sequences. Negatively regulating repetitive symmetry elements around consensus c-Myc motif and DNA sequences containing consensus c-Myc motif flanked by entirely randomized sequences show similar expression baseline. A possible explanation for this observation is that rather than complete repression, negatively regulating repetitive symmetry elements play a regulatory role in fine-tuning the reduction of gene expression, most probably by binding TFs other than c-Myc.

Esrrb is a cell-cycle-dependent associated factor balancing pluripotency and XEN differentiation.

Cell cycle and differentiation decisions are linked; however, the underlying principles that drive these decisions are unclear. Here, we combined cell-cycle reporter system and single-cell RNA sequencing (scRNA-seq) profiling to study the transcriptomes of embryonic stem cells (ESCs) in the context of cell-cycle states and differentiation. By applying retinoic acid, to G1 and G2/M ESCs, we show that, while both populations can differentiate toward epiblast stem cells (EpiSCs), only G2/M ESCs could differentiate into extraembryonic endoderm cells. We identified Esrrb, a pluripotency factor that is upregulated during G2/M, as a driver of extraembryonic endoderm stem cell (XEN) differentiation. Furthermore, enhancer chromatin states based on wild-type (WT) and ESRRB knockout (KO) ESCs show association of ESRRB with XEN poised enhancers. G1 cells overexpressing Esrrb allow ESCs to produce XENs, while ESRRB-KO ESCs lost their potential to differentiate into XEN. Overall, this study reveals a vital link between Esrrb and cell-cycle states during the exit from pluripotency.

CloneSeq - Single-cell clonal 3D culture and analysis protocol.

This CloneSeq protocol combines clonal expansion inside 3D hydrogel spheres and droplet-based RNA sequencing to resolve the limited sensitivity of single-cell approaches. CloneSeq can reveal rare subpopulations and support cellular stemness. CloneSeq can be adapted to different biological systems to discover rare subpopulations by leveraging clonal enhanced sensitivity. Important considerations include the hydrogel composition, adaptation of 3D cultured clones to the inDrops system, and inherent adhesive properties of the cells. CloneSeq is only validated for cell lines so far. For complete details on the use and execution of this protocol, please refer to (Bavli et al., 2021).

Bisection of the X chromosome disrupts the initiation of chromosome silencing during meiosis in Caenorhabditis elegans.

During meiosis, gene expression is silenced in aberrantly unsynapsed chromatin and in heterogametic sex chromosomes. Initiation of sex chromosome silencing is disrupted in meiocytes with sex chromosome-autosome translocations. To determine whether this is due to aberrant synapsis or loss of continuity of sex chromosomes, we engineered Caenorhabditis elegans nematodes with non-translocated, bisected X chromosomes. In early meiocytes of mutant males and hermaphrodites, X segments are enriched with euchromatin assembly markers and active RNA polymerase II staining, indicating active transcription. Analysis of RNA-seq data showed that genes from the X chromosome are upregulated in gonads of mutant worms. Contrary to previous models, which predicted that any unsynapsed chromatin is silenced during meiosis, our data indicate that unsynapsed X segments are transcribed. Therefore, our results suggest that sex chromosome chromatin has a unique character that facilitates its meiotic expression when its continuity is lost, regardless of whether or not it is synapsed.

H3K27Ac modification and gene expression in psoriasis.

BACKGROUND: Numerous alterations in gene expression have been described in psoriatic lesions compared to uninvolved or healthy skin. However, the mechanisms which induce this altered expression remain unclear. Epigenetic modifications play a key role in regulating genes' expression. Only three studies compared the whole-genome DNA methylation of psoriasis versus healthy skin. The present is the first study of genome-wide comparison of histone modifications between psoriatic to healthy skins. OBJECTIVE: Our objective was to explore the pattern of H3K27Ac modifications in psoriatic lesions compared to uninvolved psoriatic and healthy skin, in order to identify new genes involved in the pathogenesis of psoriasis. METHOD: Using ChIP-seq with anti H3K27Ac we compared the acetylation of lysine 27 on histone 3 (H3K27Ac) modification between psoriatic to healthy skins, combined with mRNA array. RESULTS: We found a differential H3K27Ac pattern between psoriatic compared to uninvolved or healthy skins. We found that many of the overexpressed and H3K27Ac enriched genes in psoriasis, harbor a putative GRHL transcription factor-binding site. CONCLUSIONS: In the most overexpressed genes in psoriasis, there is an enrichment of H3K27Ac. However, the loss of H3K27 acetylation modification does not correlate with decreased gene expression. GRHL appears to play an important role in the pathogenesis of psoriasis and therefore, might be a new target for psoriasis therapeutics.

Delineating the heterogeneity of matrix-directed differentiation toward soft and stiff tissue lineages via single-cell profiling.

Mesenchymal stromal/stem cells (MSCs) form a heterogeneous population of multipotent progenitors that contribute to tissue regeneration and homeostasis. MSCs assess extracellular elasticity by probing resistance to applied forces via adhesion, cytoskeletal, and nuclear mechanotransducers that direct differentiation toward soft or stiff tissue lineages. Even under controlled culture conditions, MSC differentiation exhibits substantial cell-to-cell variation that remains poorly characterized. By single-cell transcriptional profiling of nonconditioned, matrix-conditioned, and early differentiating cells, we identified distinct MSC subpopulations with distinct mechanosensitivities, differentiation capacities, and cell cycling. We show that soft matrices support adipogenesis of multipotent cells and early endochondral ossification of nonadipogenic cells, whereas intramembranous ossification and preosteoblast proliferation are directed by stiff matrices. Using diffusion pseudotime mapping, we outline hierarchical matrix-directed differentiation and perform whole-genome screening of mechanoresponsive genes. Specifically, top-ranked tropomyosin-1 is highly sensitive to stiffness cues both at RNA and protein levels, and changes in TPM1 expression determine the differentiation toward soft versus stiff tissue lineage. Consistent with actin stress fiber stabilization, tropomyosin-1 overexpression maintains YAP1 nuclear localization, activates YAP1 target genes, and directs osteogenic differentiation. Knockdown of tropomyosin-1 reversed YAP1 nuclear localization consistent with relaxation of cellular contractility, suppressed osteogenesis, activated early endochondral ossification genes after 3 d of culture in induction medium, and facilitated adipogenic differentiation after 1 wk. Our results delineate cell-to-cell variation of matrix-directed MSC differentiation and highlight tropomyosin-mediated matrix sensing.

CloneSeq: A highly sensitive analysis platform for the characterization of 3D-cultured single-cell-derived clones.

Single-cell assays have revealed the importance of heterogeneity in many biological systems. However, limited sensitivity is a major hurdle for uncovering cellular variation. To overcome it, we developed CloneSeq, combining clonal expansion inside 3D hydrogel spheres and droplet-based RNA sequencing (RNA-seq). We show that clonal cells maintain similar transcriptional profiles and cell states. CloneSeq of lung cancer cells revealed cancer-specific subpopulations, including cancer stem-like cells, that were not revealed by scRNA-seq. Clonal expansion within 3D soft microenvironments supported cellular stemness of embryonic stem cells (ESCs) even without pluripotent media, and it improved epigenetic reprogramming efficiency of mouse embryonic fibroblasts. CloneSeq of ESCs revealed that the differentiation decision is made early during Oct4 downregulation and is maintained during early clonal expansion. Together, we show CloneSeq can be adapted to different biological systems to discover rare subpopulations by leveraging the enhanced sensitivity within clones.

DNA methylation patterns expose variations in enhancer-chromatin modifications during embryonic stem cell differentiation.

In mammals, cellular identity is defined through strict regulation of chromatin modifications and DNA methylation that control gene expression. Methylation of cytosines at CpG sites in the genome is mainly associated with suppression; however, the reason for enhancer-specific methylation is not fully understood. We used sequential ChIP-bisulfite-sequencing for H3K4me1 and H3K27ac histone marks. By collecting data from the same genomic region, we identified enhancers differentially methylated between these two marks. We observed a global gain of CpG methylation primarily in H3K4me1-marked nucleosomes during mouse embryonic stem cell differentiation. This gain occurred largely in enhancer regions that regulate genes critical for differentiation. The higher levels of DNA methylation in H3K4me1- versus H3K27ac-marked enhancers, despite it being the same genomic region, indicates cellular heterogeneity of enhancer states. Analysis of single-cell RNA-seq profiles demonstrated that this heterogeneity correlates with gene expression during differentiation. Furthermore, heterogeneity of enhancer methylation correlates with transcription start site methylation. Our results provide insights into enhancer-based functional variation in complex biological systems.

TrypOx, a Novel Eukaryotic Homolog of the Redox-Regulated Chaperone Hsp33 in Trypanosoma brucei.

ATP-independent chaperones are widespread across all domains of life and serve as the first line of defense during protein unfolding stresses. One of the known crucial chaperones for bacterial survival in a hostile environment (e.g., heat and oxidative stress) is the highly conserved, redox-regulated ATP-independent bacterial chaperone Hsp33. Using a bioinformatic analysis, we describe novel eukaryotic homologs of Hsp33 identified in eukaryotic pathogens belonging to the kinetoplastids, a family responsible for lethal human diseases such as Chagas disease as caused by Trypanosoma cruzi, African sleeping sickness caused by Trypanosoma brucei spp., and leishmaniasis pathologies delivered by various Leishmania species. During their pathogenic life cycle, kinetoplastids need to cope with elevated temperatures and oxidative stress, the same conditions which convert Hsp33 into a powerful chaperone in bacteria, thus preventing aggregation of a wide range of misfolded proteins. Here, we focused on a functional characterization of the Hsp33 homolog in one of the members of the kinetoplastid family, T. brucei, (Tb927.6.2630), which we have named TrypOx. RNAi silencing of TrypOx led to a significant decrease in the survival of T. brucei under mild oxidative stress conditions, implying a protective role of TrypOx during the Trypanosomes growth. We then adopted a proteomics-driven approach to investigate the role of TrypOx in defining the oxidative stress response. Depletion of TrypOx significantly altered the abundance of proteins mediating redox homeostasis, linking TrypOx with the antioxidant system. Using biochemical approaches, we identified the redox-switch domain of TrypOx, showing its modularity and oxidation-dependent structural plasticity. Kinetoplastid parasites such as T. brucei need to cope with high levels of oxidants produced by the innate immune system, such that parasite-specific antioxidant proteins like TrypOx - which are depleted in mammals - are highly promising candidates for drug targeting.

Mass spectrometry reveals the chemistry of formaldehyde cross-linking in structured proteins.

Whole-cell cross-linking coupled to mass spectrometry is one of the few tools that can probe protein-protein interactions in intact cells. A very attractive reagent for this purpose is formaldehyde, a small molecule which is known to rapidly penetrate into all cellular compartments and to preserve the protein structure. In light of these benefits, it is surprising that identification of formaldehyde cross-links by mass spectrometry has so far been unsuccessful. Here we report mass spectrometry data that reveal formaldehyde cross-links to be the dimerization product of two formaldehyde-induced amino acid modifications. By integrating the revised mechanism into a customized search algorithm, we identify hundreds of cross-links from in situ formaldehyde fixation of human cells. Interestingly, many of the cross-links could not be mapped onto known atomic structures, and thus provide new structural insights. These findings enhance the use of formaldehyde cross-linking and mass spectrometry for structural studies.

Transcription Factor Binding in Embryonic Stem Cells Is Constrained by DNA Sequence Repeat Symmetry.

Transcription factor (TF) recognition is dictated by the underlying DNA motif sequence specific for each TF. Here, we reveal that DNA sequence repeat symmetry plays a central role in defining TF-DNA-binding preferences. In particular, we find that different TFs bind similar symmetry patterns in the context of different developmental layers. Most TFs possess dominant preferences for similar DNA repeat symmetry types. However, in some cases, preferences of specific TFs are changed during differentiation, suggesting the importance of information encoded outside of known motif regions. Histone modifications also exhibit strong preferences for similar DNA repeat symmetry patterns unique to each type of modification. Next, using an in vivo reporter assay, we show that gene expression in embryonic stem cells can be positively modulated by the presence of genomic and computationally designed DNA oligonucleotides containing identified nonconsensus-repetitive sequence elements. This supports the hypothesis that certain nonconsensus-repetitive patterns possess a functional ability to regulate gene expression. We also performed a solution NMR experiment to probe the stability of double-stranded DNA via imino proton resonances for several double-stranded DNA sequences characterized by different repetitive patterns. We suggest that such local stability might play a key role in determining TF-DNA binding preferences. Overall, our findings show that despite the enormous sequence complexity of the TF-DNA binding landscape in differentiating embryonic stem cells, this landscape can be quantitatively characterized in simple terms using the notion of DNA sequence repeat symmetry.

Direct Induction of the Three Pre-implantation Blastocyst Cell Types from Fibroblasts.

Following fertilization, totipotent cells undergo asymmetric cell divisions, resulting in three distinct cell types in the late pre-implantation blastocyst: epiblast (Epi), primitive endoderm (PrE), and trophectoderm (TE). Here, we aim to understand whether these three cell types can be induced from fibroblasts by one combination of transcription factors. By utilizing a sophisticated fluorescent knockin reporter system, we identified a combination of five transcription factors, Gata3, Eomes, Tfap2c, Myc, and Esrrb, that can reprogram fibroblasts into induced pluripotent stem cells (iPSCs), induced trophoblast stem cells (iTSCs), and induced extraembryonic endoderm stem cells (iXENs), concomitantly. In-depth transcriptomic, chromatin, and epigenetic analyses provide insights into the molecular mechanisms that underlie the reprogramming process toward the three cell types. Mechanistically, we show that the interplay between Esrrb and Eomes during the reprogramming process determines cell fate, where high levels of Esrrb induce a XEN-like state that drives pluripotency and high levels of Eomes drive trophectodermal fate.

The G-rich Repeats in FMR1 and C9orf72 Loci Are Hotspots for Local Unpairing of DNA.

Pathological mutations involving noncoding microsatellite repeats are typically located near promoters in CpG islands and are coupled with extensive repeat instability when sufficiently long. What causes these regions to be prone to repeat instability is not fully understood. There is a general consensus that instability results from the induction of unusual structures in the DNA by the repeats as a consequence of mispairing between complementary strands. In addition, there is some evidence that repeat instability is mediated by RNA transcription through the formation of three-stranded nucleic structures composed of persistent DNA:RNA hybrids, concomitant with single-strand DNA displacements (R-loops). Using human embryonic stem cells with wild-type and repeat expanded alleles in the FMR1 (CGGs) and C9orf72 (GGGGCCs) genes, we show that these loci constitute preferential sites (hotspots) for DNA unpairing. When R-loops are formed, DNA unpairing is more extensive, and is coupled with the interruptions of double-strand structures by the nontranscribing (G-rich) DNA strand. These interruptions are likely to reflect unusual structures in the DNA that drive repeat instability when the G-rich repeats considerably expand. Further, we demonstrate that when the CGGs in FMR1 are hyper-methylated and transcriptionally inactive, local DNA unpairing is abolished. Our study thus takes one more step toward the identification of dynamic, unconventional DNA structures across the G-rich repeats at FMR1 and C9orf72 disease-associated loci.

Fine-Resolution Mapping of TF Binding and Chromatin Interactions.

Transcription factor (TF) binding to DNA is crucial for transcriptional regulation. There are multiple methods for mapping such binding. These methods balance between input requirements, spatial resolution, and compatibility with high-throughput automation. Here, we describe SLIM-ChIP (short-fragment-enriched, low-input, indexed MNase ChIP), which combines enzymatic fragmentation of chromatin and on-bead indexing to address these desiderata. SLIM-ChIP reproduces a high-resolution binding map of yeast Reb1 comparable with existing methods, yet with less input material and full compatibility with high-throughput procedures. We demonstrate the robustness and flexibility of SLIM-ChIP by probing additional factors in yeast and mouse. Finally, we show that SLIM-ChIP provides information on the chromatin landscape surrounding the bound transcription factor. We identify a class of Reb1 sites where the proximal -1 nucleosome tightly interacts with Reb1 and maintains unidirectional transcription. SLIM-ChIP is an attractive solution for mapping DNA binding proteins and charting the surrounding chromatin occupancy landscape at a single-cell level.

Single-molecule decoding of combinatorially modified nucleosomes.

Different combinations of histone modifications have been proposed to signal distinct gene regulatory functions, but this area is poorly addressed by existing technologies. We applied high-throughput single-molecule imaging to decode combinatorial modifications on millions of individual nucleosomes from pluripotent stem cells and lineage-committed cells. We identified definitively bivalent nucleosomes with concomitant repressive and activating marks, as well as other combinatorial modification states whose prevalence varies with developmental potency. We showed that genetic and chemical perturbations of chromatin enzymes preferentially affect nucleosomes harboring specific modification states. Last, we combined this proteomic platform with single-molecule DNA sequencing technology to simultaneously determine the modification states and genomic positions of individual nucleosomes. This single-molecule technology has the potential to address fundamental questions in chromatin biology and epigenetic regulation.

A Multiplexed System for Quantitative Comparisons of Chromatin Landscapes.

Genome-wide profiling of histone modifications can provide systematic insight into the regulatory elements and programs engaged in a given cell type. However, conventional chromatin immunoprecipitation and sequencing (ChIP-seq) does not capture quantitative information on histone modification levels, requires large amounts of starting material, and involves tedious processing of each individual sample. Here, we address these limitations with a technology that leverages DNA barcoding to profile chromatin quantitatively and in multiplexed format. We concurrently map relative levels of multiple histone modifications across multiple samples, each comprising as few as a thousand cells. We demonstrate the technology by monitoring dynamic changes following inhibition of p300, EZH2, or KDM5, by linking altered epigenetic landscapes to chromatin regulator mutations, and by mapping active and repressive marks in purified human hematopoietic stem cells. Hence, this technology enables quantitative studies of chromatin state dynamics across rare cell types, genotypes, environmental conditions, and drug treatments.

Single-cell ChIP-seq reveals cell subpopulations defined by chromatin state.

Chromatin profiling provides a versatile means to investigate functional genomic elements and their regulation. However, current methods yield ensemble profiles that are insensitive to cell-to-cell variation. Here we combine microfluidics, DNA barcoding and sequencing to collect chromatin data at single-cell resolution. We demonstrate the utility of the technology by assaying thousands of individual cells and using the data to deconvolute a mixture of ES cells, fibroblasts and hematopoietic progenitors into high-quality chromatin state maps for each cell type. The data from each single cell are sparse, comprising on the order of 1,000 unique reads. However, by assaying thousands of ES cells, we identify a spectrum of subpopulations defined by differences in chromatin signatures of pluripotency and differentiation priming. We corroborate these findings by comparison to orthogonal single-cell gene expression data. Our method for single-cell analysis reveals aspects of epigenetic heterogeneity not captured by transcriptional analysis alone.

High-Throughput Single-Cell Labeling (Hi-SCL) for RNA-Seq Using Drop-Based Microfluidics.

The importance of single-cell level data is increasingly appreciated, and significant advances in this direction have been made in recent years. Common to these technologies is the need to physically segregate individual cells into containers, such as wells or chambers of a micro-fluidics chip. High-throughput Single-Cell Labeling (Hi-SCL) in drops is a novel method that uses drop-based libraries of oligonucleotide barcodes to index individual cells in a population. The use of drops as containers, and a microfluidics platform to manipulate them en-masse, yields a highly scalable methodological framework. Once tagged, labeled molecules from different cells may be mixed without losing the cell-of-origin information. Here we demonstrate an application of the method for generating RNA-sequencing data for multiple individual cells within a population. Barcoded oligonucleotides are used to prime cDNA synthesis within drops. Barcoded cDNAs are then combined and subjected to second generation sequencing. The data are deconvoluted based on the barcodes, yielding single-cell mRNA expression data. In a proof-of-concept set of experiments we show that this method yields data comparable to other existing methods, but with unique potential for assaying very large numbers of cells.