Age-related and species-specific methylation changes in the protein-coding marmoset sperm epigenome.

The sperm epigenome is thought to affect the developmental programming of the resulting embryo, influencing health and disease in later life. Age-related methylation changes in the sperm of old fathers may mediate the increased risks for reproductive and offspring medical problems. The impact of paternal age on sperm methylation has been extensively studied in humans and, to a lesser extent, in rodents and cattle. Here, we performed a comparative analysis of paternal age effects on protein-coding genes in the human and marmoset sperm methylomes. The marmoset has gained growing importance as a non-human primate model of aging and age-related diseases. Using reduced representation bisulfite sequencing, we identified age-related differentially methylated transcription start site (ageTSS) regions in 204 marmoset and 27 human genes. The direction of methylation changes was the opposite, increasing with age in marmosets and decreasing in humans. None of the identified ageTSS was differentially methylated in both species. Although the average methylation levels of all TSS regions were highly correlated between marmosets and humans, with the majority of TSS being hypomethylated in sperm, more than 300 protein-coding genes were endowed with species-specifically (hypo)methylated TSS. Several genes of the glycosphingolipid (GSL) biosynthesis pathway, which plays a role in embryonic stem cell differentiation and regulation of development, were hypomethylated (<5%) in human and fully methylated (>95%) in marmoset sperm. The expression levels and patterns of defined sets of GSL genes differed considerably between human and marmoset pre-implantation embryo stages and blastocyst tissues, respectively.

Rat post-implantation epiblast-derived pluripotent stem cells produce functional germ cells.

In mammals, pluripotent cells transit through a continuum of distinct molecular and functional states en route to initiating lineage specification. Capturing pluripotent stem cells (PSCs) mirroring in vivo pluripotent states provides accessible in vitro models to study the pluripotency program and mechanisms underlying lineage restriction. Here, we develop optimal culture conditions to derive and propagate post-implantation epiblast-derived PSCs (EpiSCs) in rats, a valuable model for biomedical research. We show that rat EpiSCs (rEpiSCs) can be reset toward the naive pluripotent state with exogenous Klf4, albeit not with the other five candidate genes (Nanog, Klf2, Esrrb, Tfcp2l1, and Tbx3) effective in mice. Finally, we demonstrate that rat EpiSCs retain competency to produce authentic primordial germ cell-like cells that undergo functional gametogenesis leading to the birth of viable offspring. Our findings in the rat model uncover principles underpinning pluripotency and germline competency across species.

Origin and segregation of the human germline.

Human germline-soma segregation occurs during weeks 2-3 in gastrulating embryos. Although direct studies are hindered, here, we investigate the dynamics of human primordial germ cell (PGCs) specification using in vitro models with temporally resolved single-cell transcriptomics and in-depth characterisation using in vivo datasets from human and nonhuman primates, including a 3D marmoset reference atlas. We elucidate the molecular signature for the transient gain of competence for germ cell fate during peri-implantation epiblast development. Furthermore, we show that both the PGCs and amnion arise from transcriptionally similar TFAP2A-positive progenitors at the posterior end of the embryo. Notably, genetic loss of function experiments shows that TFAP2A is crucial for initiating the PGC fate without detectably affecting the amnion and is subsequently replaced by TFAP2C as an essential component of the genetic network for PGC fate. Accordingly, amniotic cells continue to emerge from the progenitors in the posterior epiblast, but importantly, this is also a source of nascent PGCs.

Epigenetic resetting in the human germ line entails histone modification remodeling.

Epigenetic resetting in the mammalian germ line entails acute DNA demethylation, which lays the foundation for gametogenesis, totipotency, and embryonic development. We characterize the epigenome of hypomethylated human primordial germ cells (hPGCs) to reveal mechanisms preventing the widespread derepression of genes and transposable elements (TEs). Along with the loss of DNA methylation, we show that hPGCs exhibit a profound reduction of repressive histone modifications resulting in diminished heterochromatic signatures at most genes and TEs and the acquisition of a neutral or paused epigenetic state without transcriptional activation. Efficient maintenance of a heterochromatic state is limited to a subset of genomic loci, such as evolutionarily young TEs and some developmental genes, which require H3K9me3 and H3K27me3, respectively, for efficient transcriptional repression. Accordingly, transcriptional repression in hPGCs presents an exemplary balanced system relying on local maintenance of heterochromatic features and a lack of inductive cues.

Microgel culture and spatial identity mapping elucidate the signalling requirements for primate epiblast and amnion formation.

The early specification and rapid growth of extraembryonic membranes are distinctive hallmarks of primate embryogenesis. These complex tasks are resolved through an intricate combination of signals controlling the induction of extraembryonic lineages and, at the same time, safeguarding the pluripotent epiblast. Here, we delineate the signals orchestrating primate epiblast and amnion identity. We encapsulated marmoset pluripotent stem cells into agarose microgels and identified culture conditions for the development of epiblast- and amnion-spheroids. Spatial identity mapping authenticated spheroids generated in vitro by comparison with marmoset embryos in vivo. We leveraged the microgel system to functionally interrogate the signalling environment of the post-implantation primate embryo. Single-cell profiling of the resulting spheroids demonstrated that activin/nodal signalling is required for embryonic lineage identity. BMP4 promoted amnion formation and maturation, which was counteracted by FGF signalling. Our combination of microgel culture, single-cell profiling and spatial identity mapping provides a powerful approach to decipher the essential cues for embryonic and extraembryonic lineage formation in primate embryogenesis.

Spatial profiling of early primate gastrulation in utero.

Gastrulation controls the emergence of cellular diversity and axis patterning in the early embryo. In mammals, this transformation is orchestrated by dynamic signalling centres at the interface of embryonic and extraembryonic tissues1-3. Elucidating the molecular framework of axis formation in vivo is fundamental for our understanding of human development4-6 and to advance stem-cell-based regenerative approaches7. Here we illuminate early gastrulation of marmoset embryos in utero using spatial transcriptomics and stem-cell-based embryo models. Gaussian process regression-based 3D transcriptomes delineate the emergence of the anterior visceral endoderm, which is hallmarked by conserved (HHEX, LEFTY2, LHX1) and primate-specific (POSTN, SDC4, FZD5) factors. WNT signalling spatially coordinates the formation of the primitive streak in the embryonic disc and is counteracted by SFRP1 and SFRP2 to sustain pluripotency in the anterior domain. Amnion specification occurs at the boundaries of the embryonic disc through ID1, ID2 and ID3 in response to BMP signalling, providing a developmental rationale for amnion differentiation of primate pluripotent stem cells (PSCs). Spatial identity mapping demonstrates that primed marmoset PSCs exhibit the highest similarity to the anterior embryonic disc, whereas naive PSCs resemble the preimplantation epiblast. Our 3D transcriptome models reveal the molecular code of lineage specification in the primate embryo and provide an in vivo reference to decipher human development.

Nuclear transfer system for the direct induction of embryonic transcripts from intra- and cross-species nuclei using mouse 4-cell embryos.

Reprogramming of somatic nuclei toward the embryonic state has been studied using nuclear transfer (NT) to an oocyte at the metaphase II (MII) stage. In this NT, a somatic nucleus transplanted into an MII oocyte of the same species undergoes DNA replication and cell division before activating embryonic genes. Here, we describe a direct NT protocol using 4-cell stage mouse embryos that enables reprogramming of intra- and cross-species nuclei to express embryonic genes without requiring DNA replication and cell division. For complete details on the use and execution of this protocol, please refer to Tomikawa et al. (2021).

Sequential enhancer state remodelling defines human germline competence and specification.

Germline-soma segregation is a fundamental event during mammalian embryonic development. Here we establish the epigenetic principles of human primordial germ cell (hPGC) development using in vivo hPGCs and stem cell models recapitulating gastrulation. We show that morphogen-induced remodelling of mesendoderm enhancers transiently confers the competence for hPGC fate, but further activation favours mesoderm and endoderm fates. Consistently, reducing the expression of the mesendodermal transcription factor OTX2 promotes the PGC fate. In hPGCs, SOX17 and TFAP2C initiate activation of enhancers to establish a core germline programme, including the transcriptional repressor PRDM1 and pluripotency factors POU5F1 and NANOG. We demonstrate that SOX17 enhancers are the critical components in the regulatory circuitry of germline competence. Furthermore, activation of upstream cis-regulatory elements by an optimized CRISPR activation system is sufficient for hPGC specification. We reveal an enhancer-linked germline transcription factor network that provides the basis for the evolutionary divergence of mammalian germlines.

Cell division- and DNA replication-free reprogramming of somatic nuclei for embryonic transcription.

Nuclear transfer systems represent the efficient means to reprogram a cell and in theory provide a basis for investigating the development of endangered species. However, conventional nuclear transfer using oocytes of laboratory animals does not allow reprogramming of cross-species nuclei owing to defects in cell divisions and activation of embryonic genes. Here, we show that somatic nuclei transferred into mouse four-cell embryos arrested at the G2/M phase undergo reprogramming toward the embryonic state. Remarkably, genome-wide transcriptional reprogramming is induced within a day, and ZFP281 is important for this replication-free reprogramming. This system further enables transcriptional reprogramming of cells from Oryx dammah, now extinct in the wild. Thus, our findings indicate that arrested mouse embryos are competent to induce intra- and cross-species reprogramming. The direct induction of embryonic transcripts from diverse genomes paves a unique approach for identifying mechanisms of transcriptional reprogramming and genome activation from a diverse range of species.

Building a stem cell-based primate uterus.

The uterus is the organ for embryo implantation and fetal development. Most current models of the uterus are centred around capturing its function during later stages of pregnancy to increase the survival in pre-term births. However, in vitro models focusing on the uterine tissue itself would allow modelling of pathologies including endometriosis and uterine cancers, and open new avenues to investigate embryo implantation and human development. Motivated by these key questions, we discuss how stem cell-based uteri may be engineered from constituent cell parts, either as advanced self-organising cultures, or by controlled assembly through microfluidic and print-based technologies.

Time-series transcriptomics reveals a BBX32-directed control of acclimation to high light in mature Arabidopsis leaves.

The photosynthetic capacity of mature leaves increases after several days' exposure to constant or intermittent episodes of high light (HL) and is manifested primarily as changes in chloroplast physiology. How this chloroplast-level acclimation to HL is initiated and controlled is unknown. From expanded Arabidopsis leaves, we determined HL-dependent changes in transcript abundance of 3844 genes in a 0-6 h time-series transcriptomics experiment. It was hypothesized that among such genes were those that contribute to the initiation of HL acclimation. By focusing on differentially expressed transcription (co-)factor genes and applying dynamic statistical modelling to the temporal transcriptomics data, a regulatory network of 47 predominantly photoreceptor-regulated transcription (co-)factor genes was inferred. The most connected gene in this network was B-BOX DOMAIN CONTAINING PROTEIN32 (BBX32). Plants overexpressing BBX32 were strongly impaired in acclimation to HL and displayed perturbed expression of photosynthesis-associated genes under LL and after exposure to HL. These observations led to demonstrating that as well as regulation of chloroplast-level acclimation by BBX32, CRYPTOCHROME1, LONG HYPOCOTYL5, CONSTITUTIVELY PHOTOMORPHOGENIC1 and SUPPRESSOR OF PHYA-105 are important. In addition, the BBX32-centric gene regulatory network provides a view of the transcriptional control of acclimation in mature leaves distinct from other photoreceptor-regulated processes, such as seedling photomorphogenesis.

Inferring Gene Regulatory Networks from Multiple Datasets.

Gaussian process dynamical systems (GPDS) represent Bayesian nonparametric approaches to inference of nonlinear dynamical systems, and provide a principled framework for the learning of biological networks from multiple perturbed time series measurements of gene or protein expression. Such approaches are able to capture the full richness of complex ODE models, and can be scaled for inference in moderately large systems containing hundreds of genes. Related hierarchical approaches allow for inference from multiple datasets in which the underlying generative networks are assumed to have been rewired, either by context-dependent changes in network structure, evolutionary processes, or synthetic manipulation. These approaches can also be used to leverage experimentally determined network structures from one species into another where the network structure is unknown. Collectively, these methods provide a comprehensive and flexible platform for inference from a diverse range of data, with applications in systems and synthetic biology, as well as spatiotemporal modelling of embryo development. In this chapter we provide an overview of GPDS approaches and highlight their applications in the biological sciences, with accompanying tutorials available as a Jupyter notebook from https://github.com/cap76/GPDS .

Branch-recombinant Gaussian processes for analysis of perturbations in biological time series.

Motivation: A common class of behaviour encountered in the biological sciences involves branching and recombination. During branching, a statistical process bifurcates resulting in two or more potentially correlated processes that may undergo further branching; the contrary is true during recombination, where two or more statistical processes converge. A key objective is to identify the time of this bifurcation (branch or recombination time) from time series measurements, e.g. by comparing a control time series with perturbed time series. Gaussian processes (GPs) represent an ideal framework for such analysis, allowing for nonlinear regression that includes a rigorous treatment of uncertainty. Currently, however, GP models only exist for two-branch systems. Here, we highlight how arbitrarily complex branching processes can be built using the correct composition of covariance functions within a GP framework, thus outlining a general framework for the treatment of branching and recombination in the form of branch-recombinant Gaussian processes (B-RGPs). Results: We first benchmark the performance of B-RGPs compared to a variety of existing regression approaches, and demonstrate robustness to model misspecification. B-RGPs are then used to investigate the branching patterns of Arabidopsis thaliana gene expression following inoculation with the hemibotrophic bacteria, Pseudomonas syringae DC3000, and a disarmed mutant strain, hrpA. By grouping genes according to the number of branches, we could naturally separate out genes involved in basal immune response from those subverted by the virulent strain, and show enrichment for targets of pathogen protein effectors. Finally, we identify two early branching genes WRKY11 and WRKY17, and show that genes that branched at similar times to WRKY11/17 were enriched for W-box binding motifs, and overrepresented for genes differentially expressed in WRKY11/17 knockouts, suggesting that branch time could be used for identifying direct and indirect binding targets of key transcription factors. Availability and implementation: https://github.com/cap76/BranchingGPs. Supplementary information: Supplementary data are available at Bioinformatics online.

Comparative transcriptome analysis in Arabidopsis ein2/ore3 and ahk3/ore12 mutants during dark-induced leaf senescence.

Leaf senescence involves degenerative but active biological processes that require balanced regulation of pro- and anti-senescing activities. Ethylene and cytokinin are major antagonistic regulatory hormones that control the timing and progression rate of leaf senescence. To identify the roles of these hormones in the regulation of leaf senescence in Arabidopsis, global gene expression profiles in detached leaves of the wild type, an ethylene-insensitive mutant (ein2/ore3), and a constitutive cytokinin response mutant (ahk3/ore12) were investigated during dark-induced leaf senescence. Comparative transcriptome analyses revealed that genes involved in oxidative or salt stress response were preferentially altered in the ein2/ore3 mutant, whereas genes involved in ribosome biogenesis were affected in the ahk3/ore12 mutant during dark-induced leaf senescence. Similar results were also obtained for developmental senescence. Through extensive molecular and physiological analyses in ein2/ore3 and ahk3/ore12 during dark-induced leaf senescence, together with responses when treated with cytokinin and ethylene inhibitor, we conclude that ethylene acts as a senescence-promoting factor via the transcriptional regulation of stress-related responses, whereas cytokinin acts as an anti-senescing agent by maintaining cellular activities and preserving the translational machinery. These findings provide new insights into how plants utilize two antagonistic hormones, ethylene and cytokinin, to regulate the molecular programming of leaf senescence.

Stella modulates transcriptional and endogenous retrovirus programs during maternal-to-zygotic transition.

The maternal-to-zygotic transition (MZT) marks the period when the embryonic genome is activated and acquires control of development. Maternally inherited factors play a key role in this critical developmental process, which occurs at the 2-cell stage in mice. We investigated the function of the maternally inherited factor Stella (encoded by Dppa3) using single-cell/embryo approaches. We show that loss of maternal Stella results in widespread transcriptional mis-regulation and a partial failure of MZT. Strikingly, activation of endogenous retroviruses (ERVs) is significantly impaired in Stella maternal/zygotic knockout embryos, which in turn leads to a failure to upregulate chimeric transcripts. Amongst ERVs, MuERV-L activation is particularly affected by the absence of Stella, and direct in vivo knockdown of MuERV-L impacts the developmental potential of the embryo. We propose that Stella is involved in ensuring activation of ERVs, which themselves play a potentially key role during early development, either directly or through influencing embryonic gene expression.

Inferring the perturbation time from biological time course data.

MOTIVATION: Time course data are often used to study the changes to a biological process after perturbation. Statistical methods have been developed to determine whether such a perturbation induces changes over time, e.g. comparing a perturbed and unperturbed time course dataset to uncover differences. However, existing methods do not provide a principled statistical approach to identify the specific time when the two time course datasets first begin to diverge after a perturbation; we call this the perturbation time. Estimation of the perturbation time for different variables in a biological process allows us to identify the sequence of events following a perturbation and therefore provides valuable insights into likely causal relationships. RESULTS: We propose a Bayesian method to infer the perturbation time given time course data from a wild-type and perturbed system. We use a non-parametric approach based on Gaussian Process regression. We derive a probabilistic model of noise-corrupted and replicated time course data coming from the same profile before the perturbation time and diverging after the perturbation time. The likelihood function can be worked out exactly for this model and the posterior distribution of the perturbation time is obtained by a simple histogram approach, without recourse to complex approximate inference algorithms. We validate the method on simulated data and apply it to study the transcriptional change occurring in Arabidopsis following inoculation with Pseudomonas syringae pv. tomato DC3000 versus the disarmed strain DC3000hrpA AVAILABILITY AND IMPLEMENTATION: : An R package, DEtime, implementing the method is available at https://github.com/ManchesterBioinference/DEtime along with the data and code required to reproduce all the results. CONTACT: Jing.Yang@manchester.ac.uk or Magnus.Rattray@manchester.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Time-Series Transcriptomics Reveals That AGAMOUS-LIKE22 Affects Primary Metabolism and Developmental Processes in Drought-Stressed Arabidopsis.

In Arabidopsis thaliana, changes in metabolism and gene expression drive increased drought tolerance and initiate diverse drought avoidance and escape responses. To address regulatory processes that link these responses, we set out to identify genes that govern early responses to drought. To do this, a high-resolution time series transcriptomics data set was produced, coupled with detailed physiological and metabolic analyses of plants subjected to a slow transition from well-watered to drought conditions. A total of 1815 drought-responsive differentially expressed genes were identified. The early changes in gene expression coincided with a drop in carbon assimilation, and only in the late stages with an increase in foliar abscisic acid content. To identify gene regulatory networks (GRNs) mediating the transition between the early and late stages of drought, we used Bayesian network modeling of differentially expressed transcription factor (TF) genes. This approach identified AGAMOUS-LIKE22 (AGL22), as key hub gene in a TF GRN. It has previously been shown that AGL22 is involved in the transition from vegetative state to flowering but here we show that AGL22 expression influences steady state photosynthetic rates and lifetime water use. This suggests that AGL22 uniquely regulates a transcriptional network during drought stress, linking changes in primary metabolism and the initiation of stress responses.

Transcriptional Dynamics Driving MAMP-Triggered Immunity and Pathogen Effector-Mediated Immunosuppression in Arabidopsis Leaves Following Infection with Pseudomonas syringae pv tomato DC3000.

Transcriptional reprogramming is integral to effective plant defense. Pathogen effectors act transcriptionally and posttranscriptionally to suppress defense responses. A major challenge to understanding disease and defense responses is discriminating between transcriptional reprogramming associated with microbial-associated molecular pattern (MAMP)-triggered immunity (MTI) and that orchestrated by effectors. A high-resolution time course of genome-wide expression changes following challenge with Pseudomonas syringae pv tomato DC3000 and the nonpathogenic mutant strain DC3000hrpA- allowed us to establish causal links between the activities of pathogen effectors and suppression of MTI and infer with high confidence a range of processes specifically targeted by effectors. Analysis of this information-rich data set with a range of computational tools provided insights into the earliest transcriptional events triggered by effector delivery, regulatory mechanisms recruited, and biological processes targeted. We show that the majority of genes contributing to disease or defense are induced within 6 h postinfection, significantly before pathogen multiplication. Suppression of chloroplast-associated genes is a rapid MAMP-triggered defense response, and suppression of genes involved in chromatin assembly and induction of ubiquitin-related genes coincide with pathogen-induced abscisic acid accumulation. Specific combinations of promoter motifs are engaged in fine-tuning the MTI response and active transcriptional suppression at specific promoter configurations by P. syringae.

Inferring orthologous gene regulatory networks using interspecies data fusion.

MOTIVATION: The ability to jointly learn gene regulatory networks (GRNs) in, or leverage GRNs between related species would allow the vast amount of legacy data obtained in model organisms to inform the GRNs of more complex, or economically or medically relevant counterparts. Examples include transferring information from Arabidopsis thaliana into related crop species for food security purposes, or from mice into humans for medical applications. Here we develop two related Bayesian approaches to network inference that allow GRNs to be jointly inferred in, or leveraged between, several related species: in one framework, network information is directly propagated between species; in the second hierarchical approach, network information is propagated via an unobserved 'hypernetwork'. In both frameworks, information about network similarity is captured via graph kernels, with the networks additionally informed by species-specific time series gene expression data, when available, using Gaussian processes to model the dynamics of gene expression. RESULTS: Results on in silico benchmarks demonstrate that joint inference, and leveraging of known networks between species, offers better accuracy than standalone inference. The direct propagation of network information via the non-hierarchical framework is more appropriate when there are relatively few species, while the hierarchical approach is better suited when there are many species. Both methods are robust to small amounts of mislabelling of orthologues. Finally, the use of Saccharomyces cerevisiae data and networks to inform inference of networks in the budding yeast Schizosaccharomyces pombe predicts a novel role in cell cycle regulation for Gas1 (SPAC19B12.02c), a 1,3-beta-glucanosyltransferase. AVAILABILITY AND IMPLEMENTATION: MATLAB code is available from http://go.warwick.ac.uk/systemsbiology/software/.

CSI: a nonparametric Bayesian approach to network inference from multiple perturbed time series gene expression data.

Here we introduce the causal structure identification (CSI) package, a Gaussian process based approach to inferring gene regulatory networks (GRNs) from multiple time series data. The standard CSI approach infers a single GRN via joint learning from multiple time series datasets; the hierarchical approach (HCSI) infers a separate GRN for each dataset, albeit with the networks constrained to favor similar structures, allowing for the identification of context specific networks. The software is implemented in MATLAB and includes a graphical user interface (GUI) for user friendly inference. Finally the GUI can be connected to high performance computer clusters to facilitate analysis of large genomic datasets.