During early embryogenesis, mechanical constraints and localized biochemical signals co-occur around anteroposterior axis determination and symmetry breaking. Their relative roles, however, are hard to tease apart in vivo. Using brachyury (Bra), a primitive streak and mesendoderm marker in mouse embryoid bodies (EBs), we studied how contact, biochemical cues and neighboring cell cues affect the positioning of a primitive streak-like locus and thus determine the anteroposterior axis. We show that a Bra-competent layer must be formed in the EB before Bra expression initiates, and that Bra onset locus position is biased by contact points of the EB with its surrounding, probably through modulation of chemical cues rather than by mechanical signaling. We can push or pull Bra onset away from contact points by introducing a separate localized Wnt signal source, or maneuver Bra onset to a few loci or to an isotropic peripheral pattern. Furthermore, we show that Foxa2-positive cells are predictive of the future location of Bra onset, demonstrating an earlier symmetrybreaking event. Our analysis of factors affecting symmetry breaking and spatial fate choice during this developmental process could prove valuable for in vitro differentiation and organoid formation.

Embryonic stem cells can spontaneously differentiate into cell types of all germ layers within embryoid bodies (EBs) in a highly variable manner. Whether there exists an intrinsic differentiation program common to all EBs is unknown. Here, we present a novel combination of high-throughput live two-photon imaging and gene expression profiling to study early differentiation dynamics spontaneously occurring within developing EBs. Onset timing of Brachyury-GFP was highly variable across EBs, while the spatial patterns as well as the dynamics of mesendodermal progression following onset were remarkably similar. We therefore defined a ‘developmental clock’ using the Brachyury-GFP signal onset timing. Mapping snapshot gene expression measurements to this clock revealed their temporal trends, indicating that loss of pluripotency, formation of primitive streak and mesodermal lineage progression are synchronized in EBs. Exogenous activation of Wnt or BMP signaling accelerated the intrinsic clock. CHIR down-regulated Wnt3, allowing insights into dependency mechanisms between canonical Wnt signaling and multiple genes. Our findings reveal a developmental clock characteristic of an early differentiation program common to all EBs, further establishing them as an in vitro developmental model.

During early embryogenesis, mechanical signals, localized biochemical signals and neighboring cell layers interaction coordinate around anteroposterior axis determination and symmetry breaking. Deciphering their relative roles, which are hard to tease apart in vivo, will enhance our understanding of how these processes are driven. In recent years, in vitro 3D models of early mammalian development, such as embryoid bodies (EBs) and gastruloids, were successful in mimicking various aspects of the early embryo, providing high throughput accessible systems for studying the basic rules shaping cell fate and morphology during embryogenesis. Using Brachyury (Bry), a primitive streak and mesendoderm marker in EBs, we study how contact, biochemical and neighboring cell cues affect the positioning of a primitive streak-like locus, determining the AP axis. We show that a Bry-competent layer must be formed in the EB before Bry expression initiates, and that Bry onset locus selection depends on contact points of the EB with its surrounding. We can maneuver Bry onset to occur at a specific locus, a few loci, or in an isotropic peripheral pattern. By spatially separating contact and biochemical signal sources, we show these two modalities can be integrated by the EB to generate a single Bry locus. Finally, we show Foxa2+ cells are predictive of the future location of Bry onset, demonstrating an earlier symmetry-breaking event. By delineating the temporal signaling pathway dependencies of Bry and Foxa2, we were able to selectively abolish either, or spatially decouple the two cell types during EB differentiation. These findings demonstrate multiple inputs integration during an early developmental process, and may prove valuable in directing in vitro differentiation.

Chronological and biological age correlate with DNA methylation levels at specific sites in the genome. Linear combinations of multiple methylation sites, termed epigenetic clocks, can inform us of the chronological age and predict multiple health-related outcomes. However, why some sites correlate with lifespan, healthspan, or specific medical conditions remains poorly understood. Kidney fibrosis is the common pathway for Chronic Kidney Disease (CKD) which affects 10% of Europe and USA population. Currently estimated glomerular filtration rate (eGFR) is the common diagnostic measure. Here we identify epigenetic clocks and methylation sites that correlate with kidney function. Moreover, we identify methylation sites that have a unique methylation signature in the kidney. Methylation levels in the majority of these sites correlates with kidney state and function. When kidney function deteriorates, as measured by interstitial fibrosis (IF) on kidney biopsy and by eGFR, all of these sites regress towards the common methylation pattern observed in other tissues. Interestingly, while the majority of sites are less methylated in the kidney and become more methylated with loss of function, a fraction of the sites are highly methylated in the kidney and become less methylated when kidney function declines. These methylation sites are enriched for specific transcription factor binding sites. In a large subset of sites, changes in methylation pattern are accompanied by changes in gene expression in kidneys of chronic kidney disease patients. These results support the information theory of aging, and the hypothesis that the unique tissue identity, as captured by methylation patterns, is lost as tissue function declines. However, this information loss is not random, but guided towards a baseline that is dependent on the genomic loci.

Aging is a major risk factor for a plethora of diseases. The information theory of aging posits that epigenetic information loss, especially alterations in methylation patterns, serves as a principal driver of the aging process. Despite this, the connection between epigenetic information loss and disease has not been thoroughly investigated. In this study, we mapped tissue-unique methylation patterns in both healthy and pathologically diagnosed organs. Our findings revealed that in multiple diseases and tissues, including kidney in Chronic Kidney Disease (CKD), liver in liver diseases, and adipose in Type 2 Diabetes (T2D), methylation patterns degrade in a specific manner, regressing towards the mean form observed across the body. We interpret this as epigenetic information loss, where tissue-unique patterns erode. By contrast, in pancreas of T2D patients, methylation patterns diverge away from the mean. Information loss is not limited to diseases. Sun exposure, for instance, was associated with information loss in the epidermis, but not in the dermis. Age-related erosion of unique methylation patterns was also observed in brain and breast tissues, while the colon showed divergence. Our findings demonstrate that analyzing methylation patterns in tissue-unique sites can effectively distinguish between patients and healthy controls across a range of diseases. It also underscores the role of epigenetic information loss as a common feature in various pathological conditions.