The regulation of Dkk1 expression during embryonic development

First Transcriptome Analysis of Hepatoblastoma in Brazil: Unraveling the Pivotal Role of Noncoding RNAs and Metabolic Pathways

Hepatoblastoma stands as the most prevalent liver cancer in the pediatric population. Characterized by a low mutational burden, chromosomal and epigenetic alterations are key drivers of its tumorigenesis. Transcriptome analysis is a powerful tool for unraveling the molecular intricacies of hepatoblastoma, shedding light on the effects of genetic and epigenetic changes on gene expression. In this study conducted in Brazilian patients, an in-depth whole transcriptome analysis was performed on 14 primary hepatoblastomas, compared to control liver tissues. The analysis unveiled 1,492 differentially expressed genes (1,031 upregulated and 461 downregulated), including 920 protein-coding genes (62%). Upregulated biological processes were linked to cell differentiation, signaling, morphogenesis, and development, involving known hepatoblastoma-associated genes (DLK1, MEG3, HDAC2, TET1, HMGA2, DKK1, DKK4), alongside with novel findings (GYNG4, CDH3, and TNFRSF19). Downregulated processes predominantly centered around oxidation and metabolism, affecting amines, nicotinamides, and lipids, featuring novel discoveries like the repression of SYT7, TTC36, THRSP, CCND1, GCK and CAMK2B. Two genes, which displayed a concordant pattern of DNA methylation alteration in their promoter regions and dysregulation in the transcriptome, were further validated by RT-qPCR: the upregulated TNFRSF19, a key gene in the embryonic development, and the repressed THRSP, connected to lipid metabolism. Furthermore, based on protein-protein interaction analysis, we identified genes holding central positions in the network, such as HDAC2, CCND1, GCK, and CAMK2B, among others, that emerged as prime candidates warranting functional validation in future studies. Notably, a significant dysregulation of non-coding RNAs (ncRNAs), predominantly upregulated transcripts, was observed, with 42% of the top 50 highly expressed genes being ncRNAs. An integrative miRNA-mRNA analysis revealed crucial biological processes associated with metabolism, oxidation reactions of lipids and carbohydrates, and methylation-dependent chromatin silencing. In particular, four upregulated miRNAs (miR-186, miR-214, miR-377, and miR-494) played a pivotal role in the network, potentially targeting multiple protein-coding transcripts, including CCND1 and CAMK2B. In summary, our transcriptome analysis highlighted disrupted embryonic development as well as metabolic pathways, particularly those involving lipids, emphasizing the emerging role of ncRNAs as epigenetic regulators in hepatoblastomas. These findings provide insights into the complexity of the hepatoblastoma transcriptome and identify potential targets for future therapeutic interventions.

Transient suppression of Wnt signaling in poor-quality buffalo oocytes improves their developmental competence

Introduction

One of the most evolutionary conserved communication systems, the Wnt signaling pathway is a major gene regulatory pathway that affects the developmental competence of oocytes and regulates most embryonic developmental processes. The present study was undertaken to modulate the canonical Wnt (Wingless/integration) signaling pathway in the poor-quality (colorless cytoplasm after Brilliant Cresyl Blue staining, BCB-) buffalo cumulus-oocyte complexes (COCs) to improve their in vitro maturation (IVM) and embryo production (IVEP) rates.

Methods

The expression of key Wnt pathway genes was initially assessed in the good (blue cytoplasm after Brilliant Cresyl Blue staining, BCB+) and poor quality (BCB-) buffalo COCs to establish a differential activity of the Wnt pathway. The BCB- COCs were supplemented with the Wnt pathway inhibitor, Dickkopf-related protein 1 (DKK1) and later subjected to IVM and IVEP along with the BCB+ and BCB- controls. The cumulus expansion index (CEI), rate of nuclear maturation (mean percentage of oocytes in the MII stage) and embryo production, and the expression of developmentally important genes were evaluated to assess the effect of Wnt pathway inhibition on the development competence of these poor-quality oocytes.

Results

The Wnt pathway genes exhibited a significantly higher expression (p < 0.05) in the poor-quality BCB- oocytes compared to the good-quality BCB+ oocytes during the early maturation stages. The supplementation of BCB- COCs with 100 ng/mL DKK1 effectively inhibited the expression of the key mediators of the Wnt pathway (β-catenin and dishevelled homolog 1, DVL1). DKK1 supplemented BCB- COCs exhibited significantly improved cytoplasmic and nuclear maturation indices, development rates and significantly elevated expression (p < 0.05) of genes implicated in germinal vesicle breakdown (GVBD) and embryonic genome activation (EGA) vis-à-vis BCB- control COCs.

Conclusion

These data indicate that inhibition of the Wnt pathway during the initial course of oocyte maturation can improve the development competence of poor-quality buffalo oocytes.

Sex, hormones and cerebrovascular function: from development to disorder

Proper cerebrovascular development and neurogliovascular unit assembly are essential for brain growth and function throughout life, ensuring the continuous supply of nutrients and oxygen. This involves crucial events during pre- and postnatal stages through key pathways, including vascular endothelial growth factor (VEGF) and Wnt signaling. These pathways are pivotal for brain vascular growth, expansion, and blood-brain barrier (BBB) maturation. Interestingly, during fetal and neonatal life, cerebrovascular formation coincides with the early peak activity of the hypothalamic-pituitary-gonadal axis, supporting the idea of sex hormonal influence on cerebrovascular development and barriergenesis.Sex hormonal dysregulation in early development has been implicated in neurodevelopmental disorders with highly sexually dimorphic features, such as autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD). Both disorders show higher prevalence in men, with varying symptoms between sexes, with boys exhibiting more externalizing behaviors, such as aggressivity or hyperactivity, and girls displaying higher internalizing behaviors, including anxiety, depression, or attention disorders. Indeed, ASD and ADHD are linked to high prenatal testosterone exposure and reduced aromatase expression, potentially explaining sex differences in prevalence and symptomatology. In line with this, high estrogen levels seem to attenuate ADHD symptoms. At the cerebrovascular level, sex- and region-specific variations of cerebral blood flow perfusion have been reported in both conditions, indicating an impact of gonadal hormones on the brain vascular system, disrupting its ability to respond to neuronal demands.This review aims to provide an overview of the existing knowledge concerning the impact of sex hormones on cerebrovascular formation and maturation, as well as the onset of neurodevelopmental disorders. Here, we explore the concept of gonadal hormone interactions with brain vascular and BBB development to function, with a particular focus on the modulation of VEGF and Wnt signaling. We outline how these pathways may be involved in the underpinnings of ASD and ADHD. Outstanding questions and potential avenues for future research are highlighted, as uncovering sex-specific physiological and pathological aspects of brain vascular development might lead to innovative therapeutic approaches in the context of ASD, ADHD and beyond.

Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects

Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.

Diagnostic implications of pitfalls in causal variant identification based on 4577 molecularly characterized families

Despite large sequencing and data sharing efforts, previously characterized pathogenic variants only account for a fraction of Mendelian disease patients, which highlights the need for accurate identification and interpretation of novel variants. In a large Mendelian cohort of 4577 molecularly characterized families, numerous scenarios in which variant identification and interpretation can be challenging are encountered. We describe categories of challenges that cover the phenotype (e.g. novel allelic disorders), pedigree structure (e.g. imprinting disorders masquerading as autosomal recessive phenotypes), positional mapping (e.g. double recombination events abrogating candidate autozygous intervals), gene (e.g. novel gene-disease assertion) and variant (e.g. complex compound inheritance). Overall, we estimate a probability of 34.3% for encountering at least one of these challenges. Importantly, our data show that by only addressing non-sequencing-based challenges, around 71% increase in the diagnostic yield can be expected. Indeed, by applying these lessons to a cohort of 314 cases with negative clinical exome or genome reports, we could identify the likely causal variant in 54.5%. Our work highlights the need to have a thorough approach to undiagnosed diseases by considering a wide range of challenges rather than a narrow focus on sequencing technologies. It is hoped that by sharing this experience, the yield of undiagnosed disease programs globally can be improved.

Hepatic DKK1-driven steatosis is CD36 dependent

Nonalcoholic fatty liver disease (NAFLD) is prevalent worldwide; about 25% of NAFLD silently progress into steatohepatitis, in which some of them may develop into fibrosis, cirrhosis and liver failure. However, few drugs are available for NAFLD, partly because of an incomplete understanding of its pathogenic mechanisms. Here, using in vivo and in vitro gain- and loss-of-function approaches, we identified up-regulated DKK1 plays a pivotal role in high-fat diet-induced NAFLD and its progression. Mechanistic analysis reveals that DKK1 enhances the capacity of hepatocytes to uptake fatty acids through the ERK-PPARγ-CD36 axis. Moreover, DKK1 increased insulin resistance by activating the JNK signaling, which in turn exacerbates disorders of hepatic lipid metabolism. Our finding suggests that DKK1 may be a potential therapeutic and diagnosis candidate for NAFLD and metabolic disorder progression.

Systems Drug Design for Muscle Invasive Bladder Cancer and Advanced Bladder Cancer by Genome-Wide Microarray Data and Deep Learning Method with Drug Design Specifications

Bladder cancer is the 10th most common cancer worldwide. Due to the lack of understanding of the oncogenic mechanisms between muscle-invasive bladder cancer (MIBC) and advanced bladder cancer (ABC) and the limitations of current treatments, novel therapeutic approaches are urgently needed. In this study, we utilized the systems biology method via genome-wide microarray data to explore the oncogenic mechanisms of MIBC and ABC to identify their respective drug targets for systems drug discovery. First, we constructed the candidate genome-wide genetic and epigenetic networks (GWGEN) through big data mining. Second, we applied the system identification and system order detection method to delete false positives in candidate GWGENs to obtain the real GWGENs of MIBC and ABC from their genome-wide microarray data. Third, we extracted the core GWGENs from the real GWGENs by selecting the significant proteins, genes and epigenetics via the principal network projection (PNP) method. Finally, we obtained the core signaling pathways from the corresponding core GWGEN through the annotations of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway to investigate the carcinogenic mechanisms of MIBC and ABC. Based on the carcinogenic mechanisms, we selected the significant drug targets NFKB1, LEF1 and MYC for MIBC, and LEF1, MYC, NOTCH1 and FOXO1 for ABC. To design molecular drug combinations for MIBC and ABC, we employed a deep neural network (DNN)-based drug-target interaction (DTI) model with drug specifications. The DNN-based DTI model was trained by drug-target interaction databases to predict the candidate drugs for MIBC and ABC, respectively. Subsequently, the drug design specifications based on regulation ability, sensitivity and toxicity were employed as filter criteria for screening the potential drug combinations of Embelin and Obatoclax for MIBC, and Obatoclax, Entinostat and Imiquimod for ABC from their candidate drugs. In conclusion, we not only investigated the oncogenic mechanisms of MIBC and ABC, but also provided promising therapeutic options for MIBC and ABC, respectively.

Gene expression analysis of the Xenopus laevis early limb bud proximodistal axis

BACKGROUND

Limb buds develop as bilateral outgrowths of the lateral plate mesoderm and are patterned along three axes. Current models of proximal to distal patterning of early amniote limb buds suggest that two signals, a distal organizing signal from the apical epithelial ridge (AER, Fgfs) and an opposing proximal (retinoic acid [RA]) act early on pattern this axis.

RESULTS

Transcriptional analysis of stage 51 Xenopus laevis hindlimb buds sectioned along the proximal-distal axis showed that the distal region is distinct from the rest of the limb. Expression of capn8.3, a novel calpain, was located in cells immediately flanking the AER. The Wnt antagonist Dkk1 was AER-specific in Xenopus limbs. Two transcription factors, sall1 and zic5, were expressed in distal mesenchyme. Zic5 has no described association with limb development. We also describe expression of two proximal genes, gata5 and tnn, not previously associated with limb development. Differentially expressed genes were associated with Fgf, Wnt, and RA signaling as well as differential cell adhesion and proliferation.

CONCLUSIONS

We identify new candidate genes for early proximodistal limb patterning. Our analysis of RA-regulated genes supports a role for transient RA gradients in early limb bud in proximal-to-distal patterning in this anamniote model organism.

Puerarin attenuates isoproterenol‑induced myocardial hypertrophy via inhibition of the Wnt/β‑catenin signaling pathway

Myocardial hypertrophy (MH) is an independent risk factor for cardiovascular disease, which in turn lead to arrhythmia or heart failure. Therefore, attention must be paid to formulation of therapeutic strategies for MH. Puerarin is a key bioactive ingredient isolated from Pueraria genera of plants that is beneficial for the treatment of MH. However, its molecular mechanism of action has not been fully determined. In the present study, 40 µM puerarin was demonstrated to be a safe dose for human AC16 cells using Cell Counting Kit‑8 assay. The protective effects of puerarin against MH were demonstrated in AC16 cells stimulated with isoproterenol (ISO). These effects were characterized by a significant decrease in surface area of cells (assessed using fluorescence staining) and mRNA and protein expression levels of MH‑associated biomarkers, including atrial and brain natriuretic peptide, assessed using reverse transcription‑quantitative PCR and western blotting, as well as β‑myosin heavy chain mRNA expression levels. Mechanistically, western blotting demonstrated that puerarin inhibited activation of the Wnt signaling pathway. Puerarin also significantly decreased phosphorylation of p65; this was mediated via crosstalk between the Wnt and NF‑κB signaling pathways. An inhibitor (Dickkopf‑1) and activator (IM‑12) of the Wnt signaling pathway were used to demonstrate that puerarin‑mediated effects alleviated ISO‑induced MH via the Wnt signaling pathway. The results of the present study demonstrated that puerarin pre‑treatment may be a potential therapeutic strategy for preventing ISO‑induced MH and managing MH in the future.

Spatiotemporal map of key signaling factors during early penis development

The formation of the external genitalia is a highly complex developmental process, considering it involves a wide range of cell types and results in sexually dimorphic outcomes. Development is controlled by several secreted signalling factors produced in complex spatiotemporal patterns, including the hedgehog (HH), bone morphogenic protein (BMP), fibroblast growth factor (FGF) and WNT signalling families. Many of these factors act on or are influenced by the actions of the androgen receptor (AR) that is critical to masculinisation. This complexity of expression makes it difficult to conceptualise patterns of potential importance. Mapping expression during key stages of development is needed to develop a comprehensive model of how different cell types interact in formation of external genitalia, and the global regulatory networks at play. This is particularly true in light of the sensitivity of this process to environmental disruption during key stages of development. The goal of this review is to integrate all recent studies on gene expression in early penis development to create a comprehensive spatiotemporal map. This serves as a resource to aid in visualising potentially significant interactions involved in external genital development.

Diagrams of published RNA and protein localisation data for key secreted signalling factors during early penis development.

Unconventional expression patterns are identified that suggest novel signalling axes during development.

Key research gaps and limitations are identified and discussed.

Crosstalk of Brain and Bone-Clinical Observations and Their Molecular Bases

As brain and bone disorders represent major health issues worldwide, substantial clinical investigations demonstrated a bidirectional crosstalk on several levels, mechanistically linking both apparently unrelated organs. While multiple stress, mood and neurodegenerative brain disorders are associated with osteoporosis, rare genetic skeletal diseases display impaired brain development and function. Along with brain and bone pathologies, particularly trauma events highlight the strong interaction of both organs. This review summarizes clinical and experimental observations reported for the crosstalk of brain and bone, followed by a detailed overview of their molecular bases. While brain-derived molecules affecting bone include central regulators, transmitters of the sympathetic, parasympathetic and sensory nervous system, bone-derived mediators altering brain function are released from bone cells and the bone marrow. Although the main pathways of the brain-bone crosstalk remain ‘efferent’, signaling from brain to bone, this review emphasizes the emergence of bone as a crucial ‘afferent’ regulator of cerebral development, function and pathophysiology. Therefore, unraveling the physiological and pathological bases of brain-bone interactions revealed promising pharmacologic targets and novel treatment strategies promoting concurrent brain and bone recovery.

Prenatal low-dose methyltestosterone, but not dihydrotestosterone, treatment induces penile formation in female mice and guinea pigs†

Genital tubercle has bisexual potential before sex differentiation. Females exposed to androgen during sex differentiation show masculinized external genitalia, but the effects of different androgens on tubular urethral and penile formation in females are mostly unknown. In this study, we compared the masculinization effects of commonly used androgens methyltestosterone, dihydrotestosterone, and testosterone on the induction of penile formation in females. Our results suggested that prenatal treatment with low doses of methyltestosterone, but not same doses of dihydrotestosterone or testosterone, could induce penile formation in female mice. The minimum dose of dihydrotestosterone and testosterone for inducing tubular urethral formation in female mice was, respectively, 50 and 20 times higher than that of methyltestosterone. In vivo methyltestosterone treatment induced more nuclear translocation of androgen receptors in genital tubercles of female mice, affected Wnt signaling gene expressions, and then led to similar patterns of cell proliferation and death in developing genital tubercles to those of control males. We further revealed that low-dose methyltestosterone, but not same dose of dihydrotestosterone or testosterone, treatment induced penile formation in female guinea pigs. Exposure of female mouse genital tubercle organ culture to methyltestosterone, dihydrotestosterone, or testosterone could induce nuclear translocation of androgen receptors, suggesting that the differential effect of the three androgens in vivo might be due to the hormonal profile in mother or fetus, rather than the local genital tissue. To understand the differential role of these androgens in masculinization process involved is fundamental to androgen replacement therapy for diseases related to external genital masculinization.

Cross-Species RNA-Seq Study Comparing Transcriptomes of Enriched Osteocyte Populations in the Tibia and Skull

Local site-specific differences between bones in different regions of the skeleton account for their different properties and functions. To identify mechanisms behind these differences, we have performed a cross-species study comparing RNA transcriptomes of cranial and tibial osteocytes, from bones with very different primary functions and physiological responses, collected from the same individual mouse, rat, and rhesus macaque. Bioinformatic analysis was performed to identify 32 genes changed in the same direction between sites and shared across all three species. Several well-established key genes in bone growth and remodeling were upregulated in the tibias of all three species (BMP7, DKK1, FGF1, FRZB, SOST). Many of them associate or crosstalk with the Wnt signaling pathway. These results suggest Wnt signaling-related candidates for different control of regulatory mechanisms in bone homeostasis in the skull and tibia and indicate a different balance between genetically determined structure and feedback mechanisms to strains induced by mechanical loading at the different sites.

Elevated blood levels of Dickkopf-1 are associated with acute infections

INTRODUCTION

Dickkopf-1 (DKK1) is a soluble protein and antagonist of the Wnt/β-catenin signaling pathway. DKK1 is found elevated in serum from patients affected with various types of cancers and in some instances, it is considered a diagnostic and prognostic biomarker. Elevated serum levels of DKK1 have also been detected in animal models of chronic inflammatory diseases. Previous work from our laboratory has demonstrated upregulation of DKK1 in cells and mouse models of the bone marrow failure (BMF) and cancer-prone disease Fanconi anemia (FA). The present study aimed to investigate whether DKK1 blood levels in patients are associated with FA or inflammatory responses to acute infections.

METHODS

Plasma samples were collected from 58 children admitted to the Centre Mère-Enfant Soleil du Centre Hospitalier de Québec-Université Laval with signs of acute infections. Blood plasma specimens were also collected from healthy blood donors at the Héma-Québec blood donor clinic. Plasmas from patients diagnosed with FA were also included in the study. DKK1 levels in blood plasmas were assessed by standard ELISA.

RESULTS

Patients with acute infections showed dramatically high levels of DKK1 (6072 ± 518 pg/ml) in their blood compared to healthy blood donors (1726 ± 95 pg/ml). No correlations were found between DKK1 levels and C reactive protein (CRP) concentration, platelet numbers, or white blood cell counts. Patients with FA showed higher DKK1 plasma levels (3419 ± 147.5 pg/ml) than healthy blood donors (1726 ± 95 pg/ml) but significantly lower than patients with acute infections.

CONCLUSION

These findings suggest that blood DKK1 is elevated in response to infections and perhaps to inflammatory responses.

Presomitic mesoderm-specific expression of the transcriptional repressor Hes7 is controlled by E-box, T-box, and Notch signaling pathways

Somites are a pair of epithelial spheres beside a neural tube and are formed with an accurate periodicity during embryogenesis in vertebrates. It has been known that Hes7 is one of the core clock genes for somitogenesis, and its expression domain is restricted in the presomitic mesoderm (PSM). However, the molecular mechanism of how Hes7 transcription is regulated is not clear. Here, using transgenic mice and luciferase-based reporter assays and in vitro binding assays, we unravel the mechanism by which Hes7 is expressed exclusively in the PSM. We identified a Hes7 essential region residing -1.5 to -1.1 kb from the transcription start site of mouse Hes7, and this region was indispensable for PSM-specific Hes7 expression. We also present detailed analyses of cis-regulatory elements within the Hes7 essential region that directs Hes7 expression in the PSM. Hes7 expression in the PSM was up-regulated through the E-box, T-box, and RBPj-binding element in the Hes7 essential region, presumably through synergistic signaling involving mesogenin1, T-box6 (Tbx6), and Notch. Furthermore, we demonstrate that Tbx18, Ripply2, and Hes7 repress the activation of the Hes7 essential region by the aforementioned transcription factors. Our findings reveal that a unified transcriptional regulatory network involving a Hes7 essential region confers robust PSM-specific Hes7 gene expression.

Small-molecule Wnt agonists correct cleft palates in Pax9 mutant mice in utero

Clefts of the palate and/or lip are among the most common human craniofacial malformations and involve multiple genetic and environmental factors. Defects can only be corrected surgically and require complex life-long treatments. Our studies utilized the well-characterized Pax9^-/- mouse model with a consistent cleft palate phenotype to test small-molecule Wnt agonist therapies. We show that the absence of Pax9 alters the expression of Wnt pathway genes including Dkk1 and Dkk2, proven antagonists of Wnt signaling. The functional interactions between Pax9 and Dkk1 are shown by the genetic rescue of secondary palate clefts in Pax9^-/-Dkk1^f/+;Wnt1Cre embryos. The controlled intravenous delivery of small-molecule Wnt agonists (Dkk inhibitors) into pregnant Pax9^+/- mice restored Wnt signaling and led to the growth and fusion of palatal shelves, as marked by an increase in cell proliferation and osteogenesis in utero, while other organ defects were not corrected. This work underscores the importance of Pax9-dependent Wnt signaling in palatogenesis and suggests that this functional upstream molecular relationship can be exploited for the development of therapies for human cleft palates that arise from single-gene disorders.

Does mouse embryo primordial germ cell activation start before implantation as suggested by single-cell transcriptomics dynamics?

STUDY HYPOTHESIS

Does primordial germ cell (PGC) activation start before mouse embryo implantation, and does the possible regulation of the DNA (cytosine-5-)-methyltransferase 3-like (Dnmt3l) by transcription factor AP-2, gamma (TCFAP2C) have a role in this activation and in the primitive endoderm (PE)-epiblast (EPI) lineage specification?

STUDY FINDING

A burst of expression of PGC markers, such as Dppa3/Stella, Ifitm2/Fragilis, Fkbp6 and Prdm4, is observed from embryonic day (E) 3.25, and some of them, together with the late germ cell markers Zp3, Mcf2 and Morc1, become restricted to the EPI subpopulation at E4.5, while the dynamics analysis of the PE-EPI transitions in the single-cell data suggests that TCFAP2C transitorily represses Dnmt3l in EPI cells at E3.5 and such repression is withdrawn with reactivation of Dnmt3l expression in PE and EPI cells at E4.5.

WHAT IS KNOWN ALREADY

In the mouse preimplantation embryo, cells with the same phenotype take different fates based on the orchestration between topological clues (cell polarity, positional history and division orientation) and gene regulatory rules (at transcriptomics and epigenomics level), prompting the proposal of positional, stochastic and combined models explaining the specification mechanism. PGC specification starts at E6.0-6.5 post-implantation. In view of the important role of DNA methylation in developmental events, the cross-talk between some transcription factors and DNA methyltransferases is of particular relevance. TCFAP2C has a CpG DNA methylation motif that is not methylated in pluripotent cells and that could potentially bind on DNMT3L, the stimulatory DNA methyltransferase co-factor that assists in the process of de novo DNA methylation. Chromatin-immunoprecipitation analysis has demonstrated that Dnmt3l is indeed a target of TCFAP2C.

STUDY DESIGN, SAMPLES/MATERIALS, METHODS

We aimed to assess the timing of early preimplantation events and to understand better the segregation of the inner cell mass (ICM) into PE and EPI. We designed a single-cell transcriptomics dynamics computational study to identify markers of the PE-EPI bifurcation in ICM cells through searching for statistically significant (using the Student's t-test method) differently expressed genes (DEGs) between PE and EPI cells from E3.5 to E4.5. The DEGs common for E3.5 and E4.5 were used as the markers defining the steady states. We collected microarray and next-generation sequencing transcriptomics data from public databases from bulk populations and single cells from mice at E3.25, E3.5 and E4.5. The results are based on three independent single-cell transcriptomics data sets, with a fold change of 3 and P-value <0.01 for the DEG selection.

MAIN RESULTS AND THE ROLE OF CHANCE

The dynamics analysis revealed new transitory E3.5 and steady PE and EPI markers. Among the transitory E3.5 PE markers (Dnmt3l, Dusp4, Cpne8, Akap13, Dcaf12l1, Aaed1, B4galt6, BC100530, Rnpc3, Tfpi, Lgalsl, Ckap4 and Fbxl20), several (Dusp4, Akap13, Cpn8, Dcaf12l1 and Tfpi) are related to the extracellular regulated kinase pathway. We also identified new transitory E3.5 EPI markers (Sgk1, Mal, Ubxn2a, Atg16l2, Gm13102, Tcfap2c, Hexb, Slc1a1, Svip, Liph and Mier3), six new stable PE markers (Sdc4, Cpn1, Dkk1, Havcr1, F2r/Par1 and Slc7a6os) as well as three new stable EPI markers (Zp3, Mcf2 and Hexb), which are known to be late stage germ cell markers. We found that mouse PGC marker activation starts at least at E3.25 preimplantation. The transcriptomics dynamics analyses support the regulation of Dnmt3l expression by TCFAP2C.

LIMITATIONS, REASONS FOR CAUTION

Since the regulation of Dnmt3l by TCFAP2C is based on computational prediction of DNA methylation motifs, Chip-Seq and transcriptomics data, functional studies are required to validate this result.

WIDER IMPLICATIONS OF THE FINDINGS

We identified a collection of previously undescribed E3.5-specific PE and EPI markers, and new steady PE and EPI markers. Identification of these genes, many of which encode cell membrane proteins, will facilitate the isolation and characterization of early PE and EPI populations. Since it is so well established in the literature that mouse PGC specification is a post-implantation event, it was surprising for us to see activation of PGC markers as early as E3.25 preimplantation, and identify the newly found steady EPI markers as late germ cell markers. The discovery of such early activation of PGC markers has important implications in the derivation of germ cells from pluripotent cells (embryonic stem cells or induced pluripotent stem cells), since the initial stages of such derivation resemble early development. The early activation of PGC markers points out the difficulty of separating PGC cells from pluripotent populations. Collectively, our results suggest that the combining of the precision of single-cell omics data with dynamic analysis of time-series data can establish the timing of some developmental stages as earlier than previously thought.

LARGE-SCALE DATA

Not applicable.

STUDY FUNDING AND COMPETING INTERESTS

This work was supported by grants DFG15/14 and DFG15/020 from Diputación Foral de Gipuzkoa (Spain), and grant II14/00016 from I + D + I National Plan 2013-2016 (Spain) and FEDER funds. The authors declare no conflict of interest.

Differential binding of Lef1 and Msx1/2 transcription factors to Dkk1 CNEs correlates with reporter gene expression in vivo

Besides the active Wnt signalling itself, the extracellular inhibition by Dkk1 is important for various embryonic developmental processes, such as optic vesicle differentiation and facial outgrowth. Although a feedback crosstalk of the active Wnt/β-catenin signaling and Dkk1 regulation has been suggested, the control of Dkk1 transcription by the Tcf/Lef1 mediated Wnt signalling and its connection to additional signalling factors has not been elucidated in vivo. Here, we used a combination of transgenic mouse approaches and biochemical analyses to unravel the direct Dkk1 transcriptional regulation via Tcf/Lefs. By using site directed mutagenesis, we tested several conserved Tcf/Lef1 binding sites within Dkk1 conserved non-coding elements (CNEs) and found that these are required for tissue specific reporter expression. In addition a conserved Msx1/2 binding site is required for retinal reporter expression and Msx2 but not Msx1 binds its conserved binding site within CNE195 in the optic cups. Within craniofacial expression domains, Lef1 interferes with Dkk1 directly via two conserved Tcf/Lef1 binding sites in the craniofacial enhancer CNE114, both of which are required for the general craniofacial Dkk1 reporter activation. Furthermore, these Tcf/Lef1 sites are commonly bound in the whisker hair bud mesenchyme but specifically Tcf/Lef1 (no. 2) is required for mandibular activation and repression of maxillar Dkk1 activation. Lastly, we tested the Tcf/Lef1 binding capacities of the Dkk1 promoter and found that although Lef1 binds the Dkk1 promoter, these sites are not sufficient for tissue specific Dkk1 activation. Together, we here present the importance of conserved Tcf/Lef1 and Msx1/2 sites that are required for differential Dkk1 transcriptional reporter activation in vivo. This requirement directly correlates with Lef1 and Msx1/2 interaction with these genomic loci.

Fanconi anemia proteins interact with CtBP1 and modulate the expression of the Wnt antagonist Dickkopf-1

Fanconi anemia (FA) is a genetic disorder characterized by congenital abnormalities, bone marrow failure, and increased susceptibility to cancer. Of the fifteen FA proteins, Fanconi anemia group C (FANCC) is one of eight FA core complex components of the FA pathway. Unlike other FA core complex proteins, FANCC is mainly localized in the cytoplasm, where it is thought to function in apoptosis, redox regulation, cytokine signaling, and other processes. Previously, we showed that regulation of FANCC involved proteolytic processing during apoptosis. To elucidate the biological significance of this proteolytic modification, we searched for molecular interacting partners of proteolytic FANCC fragments. Among the candidates obtained, the transcriptional corepressor protein C-terminal binding protein-1 (CtBP1) interacted directly with FANCC and other FA core complex proteins. Although not required for stability of the FA core complex or ubiquitin ligase activity, CtBP1 is essential for proliferation, cell survival, and maintenance of chromosomal integrity. Expression profiling of CtBP1-depleted and FA-depleted cells revealed that several genes were commonly up- and down-regulated, including the Wnt antagonist Dickkopf-1 (DKK1). These findings suggest that FA and Wnt signaling via CtBP1 could share common effectors.

Effects of lactation and pregnancy on gene expression of endometrium of Holstein cows at day 17 of the estrous cycle or pregnancy

Objectives were to determine effects of lactation and pregnancy on endometrial gene expression on d 17 of the estrous cycle and pregnancy. Heifers (n = 33) were assigned randomly after parturition to lactating (L, n = 17) or nonlactating (NL, n = 16) groups. Cows were subjected to an ovulation synchronization program for a timed artificial insemination (TAI); 10 cows in L and 12 in NL were inseminated. Slaughter occurred 17 d after the day equivalent to TAI, and intercaruncular endometrial tissues were collected. Gene expression was determined by DNA microarray analysis for pregnant (L, n = 8; NL, n = 6) and noninseminated cyclic (L, n = 7; NL, n = 4) cows. Differentially expressed genes were selected with a P-value <0.01 and absolute expression >40. In addition, a fold effect >1.5 was used as a criterion for genes affected by pregnancy. In total, 210 genes were differentially regulated by lactation (136 downregulated and 74 upregulated), and 702 genes were differentially regulated by pregnancy (407 downregulated and 295 upregulated). The interaction effect of pregnancy and lactation affected 61 genes. Genes up- and downregulated in pregnant cows were associated with several gene ontology terms, such as defense response and interferon regulatory factor, cell adhesion, and extracellular matrix. The gene ontology analyses of up- and downregulated genes of lactating cows revealed terms related to immunoglobulin-like fold, immune response, COMM domain, and non-membrane-bounded organelle. Several genes upregulated by lactation, such as IGHG1, IGLL1, IGK, and TRD, were related to immune function, particularly for B cells and γδ T cells. Developmental genes related to limb and neural development and glucose homeostasis (e.g., DKK1, RELN, PDK4) were downregulated by lactation, whereas an interaction was also detected for RELN. The stated genes associated with immune function and developmental genes expressed in the endometrium affected by lactational state are possible candidate genes for interventions to improve fertility of lactating dairy cows.

The Dkk1 dose is critical for eye development

During mammalian ocular development, several signaling pathways control the spatiotemporal highly defined realization of the three-dimensional eye architecture. Given the complexity of these inductive signals, the developing eye is a sensitive organ for several diseases. In this study, we investigated a Dkk1+/- haploinsufficiency during eye development, resulting in coloboma and anterior eye defects, two common developmental eye disorders. Dkk1 impacts eye development from a defined developmental time point on, and is critical for lens separation from the surface ectoderm via β-catenin mediated Pdgfrα and E-cadherin expression. Dkk1 does not impact the dorso ventral retina patterning in general but is critical for Shh dependent Pax2 extension into the midline region. The described results also indicate that the retinal Dkk1 dose is critical for important steps during eye development, such as optic fissure closure and cornea formation. Further analysis of the relationship between Dkk1 and Shh signaling revealed that Dkk1 and Shh coordinatively control anterior head formation and eye induction. During eye development itself, retinal Dkk1 activation is depending on cilia mediated Gli3 regulation. Therefore, our data essentially improve the knowledge of coloboma and anterior eye defects, which are common human eye developmental defects.

Development of the external genitalia: conserved and divergent mechanisms of appendage patterning

Over the past decade, the genetics of external genital development have begun to be understood. Male and female external genitalia develop from the genital tubercle. The early tubercle has a superficial resemblance to the limb bud, but an important distinction is that the limb consists of only mesoderm and ectoderm, whereas the genital tubercle also has an endodermal component, the urethral epithelium. Urethral epithelium, which expresses Sonic hedgehog, acts as a signaling region that controls outgrowth and pattern formation, and ultimately differentiates into the urethral tube. While there are intriguing parallels between limb and genital development, recent studies have identified some key differences, including the role of Fgf signaling. Our understanding of the mechanisms of genital development still lags far behind the limb, and major questions remain to be answered, including the molecular nature of the signals that initiate genital budding, sustain outgrowth, induce tissue polarity and orchestrate urethral tubulogenesis.

Dkk1 regulates ventral midbrain dopaminergic differentiation and morphogenesis

Dickkopf1 (Dkk1) is a Wnt/β-catenin inhibitor that participates in many processes during embryonic development. One of its roles during embryogenesis is to induce head formation, since Dkk1-null mice lack head structures anterior to midbrain. The Wnt/β-catenin pathway is also known to regulate different aspects of ventral midbrain (VM) dopaminergic (DA) neuron development and, in vitro, Dkk1-mediated inhibition of the Wnt/β-catenin pathway improves the DA differentiation in mouse embryonic stem cells (mESC). However, the in vivo function of Dkk1 on the development of midbrain DA neurons remains to be elucidated. Here we examined Dkk1^+/− embryos and found that Dkk1 is required for the differentiation of DA precursors/neuroblasts into DA neurons at E13.5. This deficit persisted until E17.5, when a defect in the number and distribution of VM DA neurons was detected. Furthermore, analysis of the few Dkk1^−/− embryos that survived until E17.5 revealed a more severe loss of midbrain DA neurons and morphogenesis defects. Our results thus show that Dkk1 is required for midbrain DA differentiation and morphogenesis.

Epithelial Wnt/β-catenin signaling regulates palatal shelf fusion through regulation of Tgfβ3 expression

The canonical Wnt/β-catenin signaling plays essential role in development and diseases. Previous studies have implicated the canonical Wnt/β-catenin signaling in the regulation of normal palate development, but functional Wnt/β-catenin signaling and its tissue-specific activities remain to be accurately elucidated. In this study, we show that functional Wnt/β-catenin signaling operates primarily in the palate epithelium, particularly in the medial edge epithelium (MEE) of the developing mouse palatal shelves, consistent with the expression patterns of β-catenin and several Wnt ligands and receptors. Epithelial specific inactivation of β-catenin by the K14-Cre transgenic allele abolishes the canonical Wnt signaling activity in the palatal epithelium and leads to an abnormal persistence of the medial edge seam (MES), ultimately causing a cleft palate formation, a phenotype resembling that in Tgfβ3 mutant mice. Consistent with this phenotype is the down-regulation of Tgfβ3 and suppression of apoptosis in the MEE of the β-catenin mutant palatal shelves. Application of exogenous Tgfβ3 to the mutant palatal shelves in organ culture rescues the midline seam phenotype. On the other hand, expression of stabilized β-catenin in the palatal epithelium also disrupts normal palatogenesis by activating ectopic Tgfβ3 expression in the palatal epithelium and causing an aberrant fusion between the palate shelf and mandible in addition to severely deformed palatal shelves. Collectively, our results demonstrate an essential role for Wnt/β-catenin signaling in the epithelial component at the step of palate fusion during palate development by controlling the expression of Tgfβ3 in the MEE.