BMP4 triggers regulatory circuits specifying the cardiac mesoderm lineage.

Cardiac lineage specification in the mouse is controlled by TGFß and WNT signaling. From fly to fish, BMP has been identified as an indispensable heart inducer. A detailed analysis of the role of Bmp4 and its effectors Smad1/5, however, was still missing. We show that Bmp4 induces cardiac mesoderm formation in murine embryonic stem cells in vitro. Bmp4 first activates Wnt3 and upregulates Nodal. pSmad1/5 and the WNT effector Tcf3 form a complex, and together with pSmad2/3 activate mesoderm enhancers and Eomes. They then cooperate with Eomes to consolidate the expression of many mesoderm factors, including T. Eomes and T form a positive- feedback loop and open additional enhancers regulating early mesoderm genes, including the transcription factor Mesp1, establishing the cardiac mesoderm lineage. In parallel, the neural fate is suppressed. Our data confirm the pivotal role of Bmp4 in cardiac mesoderm formation in the mouse. We describe in detail the consecutive and cooperative actions of three signaling pathways, BMP, WNT and Nodal, and their effector transcription factors, during cardiac mesoderm specification.

Genome-wide identification of notochord enhancers comprising the regulatory landscape of the brachyury locus in mouse.

The node and notochord are important signaling centers organizing the dorso-ventral patterning of cells arising from neuro-mesodermal progenitors forming the embryonic body anlage. Owing to the scarcity of notochord progenitors and notochord cells, a comprehensive identification of regulatory elements driving notochord-specific gene expression has been lacking. Here, we have used ATAC-seq analysis of FACS-purified notochord cells from Theiler stage 12-13 mouse embryos to identify 8921 putative notochord enhancers. In addition, we established a new model for generating notochord-like cells in culture, and found 3728 of these enhancers occupied by the essential notochord control factors brachyury (T) and/or Foxa2. We describe the regulatory landscape of the T locus, comprising ten putative enhancers occupied by these factors, and confirmed the regulatory activity of three of these elements. Moreover, we characterized seven new elements by knockout analysis in embryos and identified one new notochord enhancer, termed TNE2. TNE2 cooperates with TNE in the trunk notochord, and is essential for notochord differentiation in the tail. Our data reveal an essential role of Foxa2 in directing T-expressing cells towards the notochord lineage.

Fendrr synergizes with Wnt signalling to regulate fibrosis related genes during lung development via its RNA:dsDNA triplex element.

Long non-coding RNAs are a very versatile class of molecules that can have important roles in regulating a cells function, including regulating other genes on the transcriptional level. One of these mechanisms is that RNA can directly interact with DNA thereby recruiting additional components such as proteins to these sites via an RNA:dsDNA triplex formation. We genetically deleted the triplex forming sequence (FendrrBox) from the lncRNA Fendrr in mice and found that this FendrrBox is partially required for Fendrr function in vivo. We found that the loss of the triplex forming site in developing lungs causes a dysregulation of gene programs associated with lung fibrosis. A set of these genes contain a triplex site directly at their promoter and are expressed in lung fibroblasts. We biophysically confirmed the formation of an RNA:dsDNA triplex with target promoters in vitro. We found that Fendrr with the Wnt signalling pathway regulates these genes, implicating that Fendrr synergizes with Wnt signalling in lung fibrosis.

A 37†kb region upstream of brachyury comprising a notochord enhancer is essential for notochord and tail development.

The node-streak border region comprising notochord progenitor cells (NPCs) at the posterior node and neuro-mesodermal progenitor cells (NMPs) in the adjacent epiblast is the prime organizing center for axial elongation in mouse embryos. The T-box transcription factor brachyury (T) is essential for both formation of the notochord and maintenance of NMPs, and thus is a key regulator of trunk and tail development. The T promoter controlling T expression in NMPs and nascent mesoderm has been characterized in detail; however, control elements for T expression in the notochord have not been identified yet. We have generated a series of deletion alleles by CRISPR/Cas9 genome editing in mESCs, and analyzed their effects in mutant mouse embryos. We identified a 37 kb region upstream of T that is essential for notochord function and tailbud outgrowth. Within that region, we discovered a T-binding enhancer required for notochord cell specification and differentiation. Our data reveal a complex regulatory landscape controlling cell type-specific expression and function of T in NMP/nascent mesoderm and node/notochord, allowing proper trunk and tail development.

RAC1 controls progressive movement and competitiveness of mammalian spermatozoa.

Mammalian spermatozoa employ calcium (Ca2+) and cyclic adenosine monophosphate (cAMP) signaling in generating flagellar beat. However, how sperm direct their movement towards the egg cells has remained elusive. Here we show that the Rho small G protein RAC1 plays an important role in controlling progressive motility, in particular average path velocity and linearity. Upon RAC1 inhibition of wild type sperm with the drug NSC23766, progressive movement is impaired. Moreover, sperm from mice homozygous for the genetically variant t-haplotype region (tw5/tw32), which are sterile, show strongly enhanced RAC1 activity in comparison to wild type (+/+) controls, and quickly become immotile in vitro. Sperm from heterozygous (t/+) males, on the other hand, display intermediate RAC1 activity, impaired progressive motility and transmission ratio distortion (TRD) in favor of t-sperm. We show that t/+-derived sperm consist of two subpopulations, highly progressive and less progressive. The majority of highly progressive sperm carry the t-haplotype, while most less progressive sperm contain the wild type (+) chromosome. Dosage-controlled RAC1 inhibition in t/+ sperm by NSC23766 rescues progressive movement of (+)-sperm in vitro, directly demonstrating that impairment of progressive motility in the latter is caused by enhanced RAC1 activity. The combined data show that RAC1 plays a pivotal role in controlling progressive motility in sperm, and that inappropriate, enhanced or reduced RAC1 activity interferes with sperm progressive movement. Differential RAC1 activity within a sperm population impairs the competitiveness of sperm cells expressing suboptimal RAC1 activity and thus their fertilization success, as demonstrated by t/+-derived sperm. In conjunction with t-haplotype triggered TRD, we propose that Rho GTPase signaling is essential for directing sperm towards the egg cells.

Modeling mammalian trunk development in a dish.

Mammalian post-implantation development comprises the coordination of complex lineage decisions and morphogenetic processes shaping the embryo. Despite technological advances, a comprehensive understanding of the dynamics of these processes and of the self-organization capabilities of stem cells and their descendants remains elusive. Building synthetic embryo-like structures from pluripotent embryonic stem cells in vitro promises to fill these knowledge gaps and thereby may prove transformative for developmental biology. Initial efforts to model the post-implantation embryo resulted in structures with compromised morphology (gastruloids). Recent approaches employing modified culture media, an extracellular matrix surrogate or extra-embryonic stem cells, however, succeeded in establishing embryo-like architecture. For example, embedding of gastruloids in Matrigel unlocked self-organization into trunk-like structures with bilateral somites and a neural tube-like structure, together with gut tissue and primordial germ cell-like cells. In this review, we describe the currently available models, discuss how these can be employed to acquire novel biological insights, and detail the imminent challenges for improving current models by in vitro engineering.

Two isoforms of the RAC-specific guanine nucleotide exchange factor TIAM2 act oppositely on transmission ratio distortion by the mouse t-haplotype.

Transmission ratio distortion (TRD) by the mouse t-haplotype, a variant region on chromosome 17, is a well-studied model of non-Mendelian inheritance. It is characterized by the high transmission ratio (up to 99%) of the t-haplotype from t/+ males to their offspring. TRD is achieved by the exquisite ability of the responder (Tcr) to trigger non-Mendelian inheritance of homologous chromosomes. Several distorters (Tcd1-Tcd4), which act cumulatively, together promote the high transmission ratio of Tcr and the t-haplotype. Molecularly, TRD is brought about by deregulation of Rho signaling pathways via the distorter products, which impair sperm motility, and the t-sperm specific rescue of sperm motility by the responder. The t-sperm thus can reach the egg cells faster than +-sperm and fertilize them. Previously we have shown that the responder function is accomplished by a dominant negative form of sperm motility kinase (SMOKTCR), while the distorter functions are accomplished by the Rho G protein regulators TAGAP, FGD2 and NME3 proposed to function in two oppositely acting pathways. Here we identify the RAC1-specific guanine nucleotide exchange factor TIAM2 as modifier of t-haplotype TRD. Tiam2 is expressed in two isoforms, the full-length (Tiam2l) and a short transcript (Tiam2s). Tiam2s expression from the t-allele is strongly increased compared to the wild-type allele. By transgenic approaches we show that Tiam2s enhances t-haplotype transmission, while Tiam2l has the opposite effect. Our data show that a single modifier locus can encode different gene products exerting opposite effects on a trait. They also suggest that the expression ratio of the isoforms determines if the outcome is an enhancing or a suppressive effect on the trait.

BRACHYURY directs histone acetylation to target loci during mesoderm development.

T-box transcription factors play essential roles in multiple aspects of vertebrate development. Here, we show that cooperative function of BRACHYURY (T) with histone-modifying enzymes is essential for mouse embryogenesis. A single point mutation (TY88A) results in decreased histone 3 lysine 27 acetylation (H3K27ac) at T target sites, including the T locus, suggesting that T autoregulates the maintenance of its expression and functions by recruiting permissive chromatin modifications to putative enhancers during mesoderm specification. Our data indicate that T mediates H3K27ac recruitment through a physical interaction with p300. In addition, we determine that T plays a prominent role in the specification of hematopoietic and endothelial cell types. Hematopoietic and endothelial gene expression programs are disrupted in TY88A mutant embryos, leading to a defect in the differentiation of hematopoietic progenitors. We show that this role of T is mediated, at least in part, through activation of a distal Lmo2 enhancer.

Long noncoding RNAs in organogenesis: making the difference.

A large proportion of the cellular transcriptome of higher vertebrates consists of non-protein coding transcripts, among them the long noncoding RNAs (lncRNAs). Although lncRNAs are functionally extremely divergent, many ncRNAs have been shown to interact with chromatin modifying complexes and/or with transcriptional regulators. Via such interactions, many lncRNAs are involved in controlling the activity and expression level of target genes, including important regulators of embryonic processes, and thereby fine-tune gene regulatory networks controlling cell fate, lineage balance, and organogenesis. Intriguingly, an increase in organ complexity during evolution parallels a rise in lncRNA abundance. The current data suggest that lncRNAs support the generation of cell diversity and organ complexity during embryogenesis, and thereby have promoted the evolution of more complex organisms.

Genome-wide association study and meta-analysis identify ISL1 as genome-wide significant susceptibility gene for bladder exstrophy.

The bladder exstrophy-epispadias complex (BEEC) represents the severe end of the uro-rectal malformation spectrum, and is thought to result from aberrant embryonic morphogenesis of the cloacal membrane and the urorectal septum. The most common form of BEEC is isolated classic bladder exstrophy (CBE). To identify susceptibility loci for CBE, we performed a genome-wide association study (GWAS) of 110 CBE patients and 1,177 controls of European origin. Here, an association was found with a region of approximately 220kb on chromosome 5q11.1. This region harbors the ISL1 (ISL LIM homeobox 1) gene. Multiple markers in this region showed evidence for association with CBE, including 84 markers with genome-wide significance. We then performed a meta-analysis using data from a previous GWAS by our group of 98 CBE patients and 526 controls of European origin. This meta-analysis also implicated the 5q11.1 locus in CBE risk. A total of 138 markers at this locus reached genome-wide significance in the meta-analysis, and the most significant marker (rs9291768) achieved a P value of 2.13 × 10-12. No other locus in the meta-analysis achieved genome-wide significance. We then performed murine expression analyses to follow up this finding. Here, Isl1 expression was detected in the genital region within the critical time frame for human CBE development. Genital regions with Isl1 expression included the peri-cloacal mesenchyme and the urorectal septum. The present study identified the first genome-wide significant locus for CBE at chromosomal region 5q11.1, and provides strong evidence for the hypothesis that ISL1 is the responsible candidate gene in this region.

The tissue-specific transcriptomic landscape of the mid-gestational mouse embryo.

Differential gene expression is a prerequisite for the formation of multiple cell types from the fertilized egg during embryogenesis. Understanding the gene regulatory networks controlling cellular differentiation requires the identification of crucial differentially expressed control genes and, ideally, the determination of the complete transcriptomes of each individual cell type. Here, we have analyzed the transcriptomes of six major tissues dissected from mid-gestational (TS12) mouse embryos. Approximately one billion reads derived by RNA-seq analysis provided extended transcript lengths, novel first exons and alternative transcripts of known genes. We have identified 1375 genes showing tissue-specific expression, providing gene signatures for each of the six tissues. In addition, we have identified 1403 novel putative long noncoding RNA gene loci, 439 of which show differential expression. Our analysis provides the first complete transcriptome data for the mouse embryo. It offers a rich data source for the analysis of individual genes and gene regulatory networks controlling mid-gestational development.

Different Concentrations of FGF Ligands, FGF2 or FGF8 Determine Distinct States of WNT-Induced Presomitic Mesoderm.

Presomitic mesoderm (PSM) cells are the precursors of the somites, which flank both sides of the neural tube and give rise to the musculo-skeletal system shaping the vertebrate body. WNT and FGF signaling control the formation of both the PSM and the somites and show a graded distribution with highest levels in the posterior PSM. We have used reporters for the mesoderm/PSM control genes T, Tbx6, and Msgn1 to investigate the differentiation of mouse ESCs from the naïve state via EpiSCs to PSM cells. Here we show that the activation of WNT signaling by CHIR99021 (CH) in combination with FGF ligand induces embryo-like PSM at high efficiency. By varying the FGF ligand concentration, the state of PSM cells formed can be altered. High FGF concentration supports posterior PSM formation, whereas low FGF generates anterior/differentiating PSM, in line with in vivo data. Furthermore, the level of Msgn1 expression depends on the FGF ligand concentration. We also show that Activin/Nodal signaling inhibits CH-mediated PSM induction in EpiSCs, without affecting T-expression. Inversely, Activin/Nodal inhibition enhances PSM induction by WNT/high FGF signaling. The ability to generate PSM cells of either posterior or anterior PSM identity with high efficiency in vitro will promote the investigation of the gene regulatory networks controlling the formation of nascent PSM cells and their switch to differentiating/somitic paraxial mesoderm. Stem Cells 2016;34:1790-1800.

Mechanisms of long noncoding RNA function in development and disease.

Since decades it has been known that non-protein-coding RNAs have important cellular functions. Deep sequencing recently facilitated the discovery of thousands of novel transcripts, now classified as long noncoding RNAs (lncRNAs), in many vertebrate and invertebrate species. LncRNAs are involved in a wide range of cellular mechanisms, from almost all aspects of gene expression to protein translation and stability. Recent findings implicate lncRNAs as key players of cellular differentiation, cell lineage choice, organogenesis and tissue homeostasis. Moreover, lncRNAs are involved in pathological conditions such as cancer and cardiovascular disease, and therefore provide novel biomarkers and pharmaceutical targets. Here we discuss examples illustrating the versatility of lncRNAs in gene control, development and differentiation, as well as in human disease.

Retention of gene products in syncytial spermatids promotes non-Mendelian inheritance as revealed by the t complex responder.

The t complex responder (Tcr) encoded by the mouse t haplotype is able to cause phenotypic differences between t and + sperm derived from t/+ males, leading to non-Mendelian inheritance. This capability of Tcr contradicts the concept of phenotypic equivalence proposed for sperm cells, which develop in a syncytium and actively share gene products. By analyzing a Tcr minigene in hemizygous transgenic mice, we show that Tcr gene products are post-meiotically expressed and are retained in the haploid sperm cells. The wild-type allele of Tcr, sperm motility kinase-1 (Smok1), behaves in the same manner, suggesting that Tcr/Smok reveal a common mechanism prone to evolve non-Mendelian inheritance in mammals.

Genome-wide association study and mouse expression data identify a highly conserved 32 kb intergenic region between WNT3 and WNT9b as possible susceptibility locus for isolated classic exstrophy of the bladder.

Bladder exstrophy-epispadias complex (BEEC), the severe end of the urorectal malformation spectrum, has a profound impact on continence as well as sexual and renal functions. It is widely accepted that for the majority of cases the genetic basis appears to be multifactorial. Here, we report the first study which utilizes genome-wide association methods to analyze a cohort comprising patients presenting the most common BEEC form, classic bladder exstrophy (CBE), to identify common variation associated with risk for isolated CBE. We employed discovery and follow-up samples comprising 218 cases/865 controls and 78 trios in total, all of European descent. Our discovery sample identified a marker near SALL1, showing genome-wide significant association with CBE. However, analyses performed on follow-up samples did not add further support to these findings. We were also able to identify an association with CBE across our study samples (discovery: P = 8.88 × 10(-5); follow-up: P = 0.0025; combined: 1.09 × 10(-6)) in a highly conserved 32 kb intergenic region containing regulatory elements between WNT3 and WNT9B. Subsequent analyses in mice revealed expression for both genes in the genital region during stages relevant to the development of CBE in humans. Unfortunately, we were not able to replicate the suggestive signal for WNT3 and WNT9B in a sample that was enriched for non-CBE BEEC cases (P = 0.51). Our suggestive findings support the hypothesis that larger samples are warranted to identify association of common variation with CBE.

DNA-methylome analysis of mouse intestinal adenoma identifies a tumour-specific signature that is partly conserved in human colon cancer.

Aberrant CpG methylation is a universal epigenetic trait of cancer cell genomes. However, human cancer samples or cell lines preclude the investigation of epigenetic changes occurring early during tumour development. Here, we have used MeDIP-seq to analyse the DNA methylome of APC(Min) adenoma as a model for intestinal cancer initiation, and we present a list of more than 13,000 recurring differentially methylated regions (DMRs) characterizing intestinal adenoma of the mouse. We show that Polycomb Repressive Complex (PRC) targets are strongly enriched among hypermethylated DMRs, and several PRC2 components and DNA methyltransferases were up-regulated in adenoma. We further demonstrate by bisulfite pyrosequencing of purified cell populations that the DMR signature arises de novo in adenoma cells rather than by expansion of a pre-existing pattern in intestinal stem cells or undifferentiated crypt cells. We found that epigenetic silencing of tumour suppressors, which occurs frequently in colon cancer, was rare in adenoma. Quite strikingly, we identified a core set of DMRs, which is conserved between mouse adenoma and human colon cancer, thus possibly revealing a global panel of epigenetically modified genes for intestinal tumours. Our data allow a distinction between early conserved epigenetic alterations occurring in intestinal adenoma and late stochastic events promoting colon cancer progression, and may facilitate the selection of more specific clinical epigenetic biomarkers.

Mutations of the SLIT2-ROBO2 pathway genes SLIT2 and SRGAP1 confer risk for congenital anomalies of the kidney and urinary tract.

Congenital anomalies of the kidney and urinary tract (CAKUT) account for 40-50% of chronic kidney disease that manifests in the first two decades of life. Thus far, 31 monogenic causes of isolated CAKUT have been described, explaining ~12% of cases. To identify additional CAKUT-causing genes, we performed whole-exome sequencing followed by a genetic burden analysis in 26 genetically unsolved families with CAKUT. We identified two heterozygous mutations in SRGAP1 in 2 unrelated families. SRGAP1 is a small GTPase-activating protein in the SLIT2-ROBO2 signaling pathway, which is essential for development of the metanephric kidney. We then examined the pathway-derived candidate gene SLIT2 for mutations in cohort of 749 individuals with CAKUT and we identified 3 unrelated individuals with heterozygous mutations. The clinical phenotypes of individuals with mutations in SLIT2 or SRGAP1 were cystic dysplastic kidneys, unilateral renal agenesis, and duplicated collecting system. We show that SRGAP1 is expressed in early mouse nephrogenic mesenchyme and that it is coexpressed with ROBO2 in SIX2-positive nephron progenitor cells of the cap mesenchyme in developing rat kidney. We demonstrate that the newly identified mutations in SRGAP1 lead to an augmented inhibition of RAC1 in cultured human embryonic kidney cells and that the SLIT2 mutations compromise the ability of the SLIT2 ligand to inhibit cell migration. Thus, we report on two novel candidate genes for causing monogenic isolated CAKUT in humans.

The nucleoside diphosphate kinase gene Nme3 acts as quantitative trait locus promoting non-Mendelian inheritance.

The t-haplotype, a variant form of the t-complex region on mouse chromosome 17, acts as selfish genetic element and is transmitted at high frequencies (> 95%) from heterozygous (t/+) males to their offspring. This phenotype is termed transmission ratio distortion (TRD) and is caused by the interaction of the t-complex responder (Tcr) with several quantitative trait loci (QTL), the t-complex distorters (Tcd1 to Tcd4), all located within the t-haplotype region. Current data suggest that the distorters collectively impair motility of all sperm derived from t/+ males; t-sperm is rescued by the responder, whereas (+)-sperm remains partially dysfunctional. Recently we have identified two distorters as regulators of RHO small G proteins. Here we show that the nucleoside diphosphate kinase gene Nme3 acts as a QTL on TRD. Reduction of the Nme3 dosage by gene targeting of the wild-type allele enhanced the transmission rate of the t-haplotype and phenocopied distorter function. Genetic and biochemical analysis showed that the t-allele of Nme3 harbors a mutation (P89S) that compromises enzymatic activity of the protein and genetically acts as a hypomorph. Transgenic overexpression of the Nme3 t-allele reduced t-haplotype transmission, proving it to be a distorter. We propose that the NME3 protein interacts with RHO signaling cascades to impair sperm motility through hyperactivation of SMOK, the wild-type form of the responder. This deleterious effect of the distorters is counter-balanced by the responder, SMOK(Tcr), a dominant-negative protein kinase exclusively expressed in t-sperm, thus permitting selfish behaviour and preferential transmission of the t-haplotype. In addition, the previously reported association of NME family members with RHO signaling in somatic cell motility and metastasis, in conjunction with our data involving RHO signaling in sperm motility, suggests a functional conservation between mechanisms for motility control in somatic cells and spermatozoa.

Whole-exome resequencing reveals recessive mutations in TRAP1 in individuals with CAKUT and VACTERL association.

Congenital abnormalities of the kidney and urinary tract (CAKUT) account for approximately half of children with chronic kidney disease and they are the most frequent cause of end-stage renal disease in children in the US. However, its genetic etiology remains mostly elusive. VACTERL association is a rare disorder that involves congenital abnormalities in multiple organs including the kidney and urinary tract in up to 60% of the cases. By homozygosity mapping and whole-exome resequencing combined with high-throughput mutation analysis by array-based multiplex PCR and next-generation sequencing, we identified recessive mutations in the gene TNF receptor-associated protein 1 (TRAP1) in two families with isolated CAKUT and three families with VACTERL association. TRAP1 is a heat-shock protein 90-related mitochondrial chaperone possibly involved in antiapoptotic and endoplasmic reticulum stress signaling. Trap1 is expressed in renal epithelia of developing mouse kidney E13.5 and in the kidney of adult rats, most prominently in proximal tubules and in thick medullary ascending limbs of Henle's loop. Thus, we identified mutations in TRAP1 as highly likely causing CAKUT or VACTERL association with CAKUT.

Classic bladder exstrophy: Frequent 22q11.21 duplications and definition of a 414 kb phenocritical region.

BACKGROUND: Classic bladder exstrophy (CBE) is the most common form of the bladder exstrophy and epispadias complex. Previously, we and others have identified four patients with a duplication of 22q11.21 among a total of 96 unrelated CBE patients. METHODS: Here, we investigated whether this chromosomal aberration was commonly associated with CBE/bladder exstrophy and epispadias complex in an extended case-control sample. Multiplex ligation-dependent probe amplification and microarray-based analysis were used to identify 22q11.21 duplications in 244 unrelated bladder exstrophy and epispadias complex patients (including 217 CBE patients) and 665 healthy controls. RESULTS: New duplications of variable size were identified in four CBE patients and one control. Pooling of our previous and present data (eight duplications in 313 CBE patients) yielded a combined odds ratio of 31.86 (95% confidence interval, 4.24-1407.97). Array-based sequence capture and high-throughput targeted re-sequencing established that all breakpoints resided within the low-copy repeats 22A to 22D. Comparison of the eight duplications revealed a 414 kb phenocritical region harboring 12 validated RefSeq genes. Characterization of these 12 candidate genes through whole-mount in situ hybridization of mouse embryos at embryonic day 9.5 suggested that CRKL, THAP7, and LZTR1 are CBE candidate genes. CONCLUSION: Our data suggest that duplication of 22q11.21 increases CBE risk and implicate a phenocritical region in disease formation.