H3K4me3 breadth is linked to cell identity and transcriptional consistency.

Trimethylation of histone H3 at lysine 4 (H3K4me3) is a chromatin modification known to mark the transcription start sites of active genes. Here, we show that H3K4me3 domains that spread more broadly over genes in a given cell type preferentially mark genes that are essential for the identity and function of that cell type. Using the broadest H3K4me3 domains as a discovery tool in neural progenitor cells, we identify novel regulators of these cells. Machine learning models reveal that the broadest H3K4me3 domains represent a distinct entity, characterized by increased marks of elongation. The broadest H3K4me3 domains also have more paused polymerase at their promoters, suggesting a unique transcriptional output. Indeed, genes marked by the broadest H3K4me3 domains exhibit enhanced transcriptional consistency and [corrected] increased transcriptional levels, and perturbation of H3K4me3 breadth leads to changes in transcriptional consistency. Thus, H3K4me3 breadth contains information that could ensure transcriptional precision at key cell identity/function genes.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes.

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Uterine injury during diestrus leads to placental and embryonic defects in future pregnancies in mice†.

Uterine injury from procedures such as Cesarean sections (C-sections) often have severe consequences on subsequent pregnancy outcomes, leading to disorders such as placenta previa, placenta accreta, and infertility. With rates of C-section at ~30% of deliveries in the USA and projected to continue to climb, a deeper understanding of the mechanisms by which these pregnancy disorders arise and opportunities for intervention are needed. Here we describe a rodent model of uterine injury on subsequent in utero outcomes. We observed three distinct phenotypes: increased rates of resorption and death, embryo spacing defects, and placenta accreta-like features of reduced decidua and expansion of invasive trophoblasts. We show that the appearance of embryo spacing defects depends entirely on the phase of estrous cycle at the time of injury. Using RNA-seq, we identified perturbations in the expression of components of the COX/prostaglandin pathway after recovery from injury, a pathway that has previously been demonstrated to play an important role in embryo spacing. Therefore, we demonstrate that uterine damage in this mouse model causes morphological and molecular changes that ultimately lead to placental and embryonic developmental defects.

Molecular Signatures of Placentation and Secretion Uncovered in Poeciliopsis Maternal Follicles.

Placentation evolved many times independently in vertebrates. Although the core functions of all placentas are similar, we know less about how this similarity extends to the molecular level. Here, we study Poeciliopsis, a unique genus of live-bearing fish that have independently evolved complex placental structures at least three times. The maternal follicle is a key component of these structures. It envelops yolk-rich eggs and is morphologically simple in lecithotrophic species but has elaborate villous structures in matrotrophic species. Through sequencing, the follicle transcriptome of a matrotrophic, Poeciliopsis retropinna, and lecithotrophic, P. turrubarensis, species we found genes known to be critical for placenta function expressed in both species despite their difference in complexity. Additionally, when we compare the transcriptome of different river populations of P. retropinna, known to vary in maternal provisioning, we find differential expression of secretory genes expressed specifically in the top layer of villi cells in the maternal follicle. This provides some of the first evidence that the placental structures of Poeciliopsis function using a secretory mechanism rather than direct contact with maternal circulation. Finally, when we look at the expression of placenta proteins at the maternal-fetal interface of a larger sampling of Poeciliopsis species, we find expression of key maternal and fetal placenta proteins in their cognate tissue types of all species, but follicle expression of prolactin is restricted to only matrotrophic species. Taken together, we suggest that all Poeciliopsis follicles are poised for placenta function but require expression of key genes to form secretory villi.

The Genomes of the Livebearing Fish Species Poeciliopsis retropinna and Poeciliopsis turrubarensis Reflect Their Different Reproductive Strategies.

The evolution of a placenta is predicted to be accompanied by rapid evolution of genes involved in processes that regulate mother-offspring interactions during pregnancy, such as placenta formation, embryonic development, and nutrient transfer to offspring. However, these predictions have only been tested in mammalian species, where only a single instance of placenta evolution has occurred. In this light, the genus Poeciliopsis is a particularly interesting model for placenta evolution, because in this genus a placenta has evolved independently from the mammalian placenta. Here, we present and compare genome assemblies of two species of the livebearing fish genus Poeciliopsis (family Poeciliidae) that differ in their reproductive strategy: Poeciliopsis retropinna which has a well-developed complex placenta and P. turrubarensis which lacks a placenta. We applied different assembly strategies for each species: PacBio sequencing for P. retropinna (622-Mb assembly, scaffold N50 of 21.6 Mb) and 10× Genomics Chromium technology for P. turrubarensis (597-Mb assembly, scaffold N50 of 4.2 Mb). Using the high contiguity of these genome assemblies and near-completeness of gene annotations to our advantage, we searched for gene duplications and performed a genome-wide scan for genes evolving under positive selection. We find rapid evolution in major parts of several molecular pathways involved in parent-offspring interaction in P. retropinna, both in the form of gene duplications as well as positive selection. We conclude that the evolution of the placenta in the genus Poeciliopsis is accompanied by rapid evolution of genes involved in similar genomic pathways as found in mammals.

PRG2 and AQPEP are misexpressed in fetal membranes in placenta previa and percreta†.

The obstetrical conditions placenta accreta spectrum (PAS) and placenta previa are a significant source of pregnancy-associated morbidity and mortality, yet the specific molecular and cellular underpinnings of these conditions are not known. In this study, we identified misregulated gene expression patterns in tissues from placenta previa and percreta (the most extreme form of PAS) compared with control cases. By comparing this gene set with existing placental single-cell and bulk RNA-Seq datasets, we show that the upregulated genes predominantly mark extravillous trophoblasts. We performed immunofluorescence on several candidate molecules and found that PRG2 and AQPEP protein levels are upregulated in both the fetal membranes and the placental disk in both conditions. While this increased AQPEP expression remains restricted to trophoblasts, PRG2 is mislocalized and is found throughout the fetal membranes. Using a larger patient cohort with a diverse set of gestationally aged-matched controls, we validated PRG2 as a marker for both previa and PAS and AQPEP as a marker for only previa in the fetal membranes. Our findings suggest that the extraembryonic tissues surrounding the conceptus, including both the fetal membranes and the placental disk, harbor a signature of previa and PAS that is characteristic of EVTs and that may reflect increased trophoblast invasiveness.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

HEB and E2A function as SMAD/FOXH1 cofactors.

Nodal signaling, mediated through SMAD transcription factors, is necessary for pluripotency maintenance and endoderm commitment. We identified a new motif, termed SMAD complex-associated (SCA), that is bound by SMAD2/3/4 and FOXH1 in human embryonic stem cells (hESCs) and derived endoderm. We demonstrate that two basic helix-loop-helix (bHLH) proteins-HEB and E2A-bind the SCA motif at regions overlapping SMAD2/3 and FOXH1. Furthermore, we show that HEB and E2A associate with SMAD2/3 and FOXH1, suggesting they form a complex at critical target regions. This association is biologically important, as E2A is critical for mesendoderm specification, gastrulation, and Nodal signal transduction in Xenopus tropicalis embryos. Taken together, E proteins are novel Nodal signaling cofactors that associate with SMAD2/3 and FOXH1 and are necessary for mesendoderm differentiation.

XenMine: A genomic interaction tool for the Xenopus community.

The Xenopus community has embraced recent advances in sequencing technology, resulting in the accumulation of numerous RNA-Seq and ChIP-Seq datasets. However, easily accessing and comparing datasets generated by multiple laboratories is challenging. Thus, we have created a central space to view, search and analyze data, providing essential information on gene expression changes and regulatory elements present in the genome. XenMine (www.xenmine.org) is a user-friendly website containing published genomic datasets from both Xenopus tropicalis and Xenopus laevis. We have established an analysis pipeline where all published datasets are uniformly processed with the latest genome releases. Information from these datasets can be extracted and compared using an array of pre-built or custom templates. With these search tools, users can easily extract sequences for all putative regulatory domains surrounding a gene of interest, identify the expression values of a gene of interest over developmental time, and analyze lists of genes for gene ontology terms and publications. Additionally, XenMine hosts an in-house genome browser that allows users to visualize all available ChIP-Seq data, extract specifically marked sequences, and aid in identifying important regulatory elements within the genome. Altogether, XenMine is an excellent tool for visualizing, accessing and querying analyzed datasets rapidly and efficiently.

Copy number variation is a fundamental aspect of the placental genome.

Discovery of lineage-specific somatic copy number variation (CNV) in mammals has led to debate over whether CNVs are mutations that propagate disease or whether they are a normal, and even essential, aspect of cell biology. We show that 1,000 N polyploid trophoblast giant cells (TGCs) of the mouse placenta contain 47 regions, totaling 138 Megabases, where genomic copies are underrepresented (UR). UR domains originate from a subset of late-replicating heterochromatic regions containing gene deserts and genes involved in cell adhesion and neurogenesis. While lineage-specific CNVs have been identified in mammalian cells, classically in the immune system where V(D)J recombination occurs, we demonstrate that CNVs form during gestation in the placenta by an underreplication mechanism, not by recombination nor deletion. Our results reveal that large scale CNVs are a normal feature of the mammalian placental genome, which are regulated systematically during embryogenesis and are propagated by a mechanism of underreplication.

RNA sequencing reveals a diverse and dynamic repertoire of the Xenopus tropicalis transcriptome over development.

The Xenopus embryo has provided key insights into fate specification, the cell cycle, and other fundamental developmental and cellular processes, yet a comprehensive understanding of its transcriptome is lacking. Here, we used paired end RNA sequencing (RNA-seq) to explore the transcriptome of Xenopus tropicalis in 23 distinct developmental stages. We determined expression levels of all genes annotated in RefSeq and Ensembl and showed for the first time on a genome-wide scale that, despite a general state of transcriptional silence in the earliest stages of development, approximately 150 genes are transcribed prior to the midblastula transition. In addition, our splicing analysis uncovered more than 10,000 novel splice junctions at each stage and revealed that many known genes have additional unannotated isoforms. Furthermore, we used Cufflinks to reconstruct transcripts from our RNA-seq data and found that ∼13.5% of the final contigs are derived from novel transcribed regions, both within introns and in intergenic regions. We then developed a filtering pipeline to separate protein-coding transcripts from noncoding RNAs and identified a confident set of 6686 noncoding transcripts in 3859 genomic loci. Since the current reference genome, XenTro3, consists of hundreds of scaffolds instead of full chromosomes, we also performed de novo reconstruction of the transcriptome using Trinity and uncovered hundreds of transcripts that are missing from the genome. Collectively, our data will not only aid in completing the assembly of the Xenopus tropicalis genome but will also serve as a valuable resource for gene discovery and for unraveling the fundamental mechanisms of vertebrate embryogenesis.

Maternal bias and escape from X chromosome imprinting in the midgestation mouse placenta.

To investigate the epigenetic landscape at the interface between mother and fetus, we provide a comprehensive analysis of parent-of-origin bias in the mouse placenta. Using F1 interspecies hybrids between mus musculus (C57BL/6J) and mus musculus castaneus, we sequenced RNA from 23 individual midgestation placentas, five late stage placentas, and two yolk sac samples and then used SNPs to determine whether transcripts were preferentially generated from the maternal or paternal allele. In the placenta, we find 103 genes that show significant and reproducible parent-of-origin bias, of which 78 are novel candidates. Most (96%) show a strong maternal bias which we demonstrate, via multiple mathematical models, pyrosequencing, and FISH, is not due to maternal decidual contamination. Analysis of the X chromosome also reveals paternal expression of Xist and several genes that escape inactivation, most significantly Alas2, Fhl1, and Slc38a5. Finally, sequencing individual placentas allowed us to reveal notable expression similarity between littermates. In all, we observe a striking preference for maternal transcription in the midgestation mouse placenta and a dynamic imprinting landscape in extraembryonic tissues, reflecting the complex nature of epigenetic pathways in the placenta.

Chromatin immunoprecipitation and deep sequencing in Xenopus tropicalis and Xenopus laevis.

Chromatin immunoprecipitation and deep sequencing (ChIP-SEQ) represents a powerful tool for identifying the genomic targets of transcription factors, chromatin remodeling factors, and histone modifications. The frogs Xenopus laevis and Xenopus tropicalis have historically been outstanding model systems for embryology and cell biology, with emerging utility as highly accessible embryos for genome-wide studies. Here we focus on the particular strengths and limitations of Xenopus cell biology and genomics as they apply to ChIP-SEQ, and outline a methodology for ChIP-SEQ in both species, providing detailed strategies for sample preparation, antibody selection, quality control, sequencing library preparation, and basic analysis.

Genomic analysis of gastrulation and organogenesis in the mouse.

We developed a comprehensive dataset that samples the mouse transcriptome every 6 hr, from gastrulation through organogenesis. We observe an abrupt increase in overall transcript diversity at the onset of organogenesis (e8.0); the genes that comprise these changes are preferentially clustered along chromosome 7 and contain a significant enrichment of Gli binding sites. Furthermore, we identify seven dominant patterns of gene expression during gastrulation and organogenesis. Genes clustered according to these seven patterns constitute distinct functional classes, including a cluster enriched for gastrulation and pluripotency genes, two clusters differentially regulating localization and ion metabolism, and three clusters involved in discrete aspects of organogenesis. The last cluster is defined by a dramatic transient decrease in the expression of genes that regulate RNA processing and the cell cycle. Drosophila homologs of these genes are also coordinately downregulated following gastrulation, suggesting that the combined function of these genes has been conserved during metazoan evolution.

PTK7/CCK-4 is a novel regulator of planar cell polarity in vertebrates.

In addition to the apical-basal polarity pathway operating in epithelial cells, a planar cell polarity (PCP) pathway establishes polarity within the plane of epithelial tissues and is conserved from Drosophila to mammals. In Drosophila, a 'core' group of PCP genes including frizzled (fz), flamingo/starry night, dishevelled (dsh), Van Gogh/strabismus and prickle, function to regulate wing hair, bristle and ommatidial polarity. In vertebrates, the PCP pathway regulates convergent extension movements and neural tube closure, as well as the orientation of stereociliary bundles of sensory hair cells in the inner ear. Here we show that a mutation in the mouse protein tyrosine kinase 7 (PTK7) gene, which encodes an evolutionarily conserved transmembrane protein with tyrosine kinase homology, disrupts neural tube closure and stereociliary bundle orientation, and shows genetic interactions with a mutation in the mouse Van Gogh homologue vangl2. We also show that PTK7 is dynamically localized during hair cell polarization, and that the Xenopus homologue of PTK7 is required for neural convergent extension and neural tube closure. These results identify PTK7 as a novel regulator of PCP in vertebrates.

Siva plays a critical role in mouse embryonic development.

The Siva protein, named after the Hindu God of Destruction, plays important roles in apoptosis in various contexts, including downstream of death receptor activation or p53 tumor suppressor engagement. The function of Siva in organismal development and homeostasis, however, has remained uncharacterized. Here, we generate Siva knockout mice to characterize the physiological function of Siva in vivo. Interestingly, we find that Siva deficiency causes early embryonic lethality accompanied by multiple phenotypes, including developmental delay, abnormal neural tube closure, and defective placenta and yolk sac formation. Examination of Siva expression during embryogenesis shows that Siva is expressed in both embryonic and extra-embryonic tissues, including within the mesoderm, which may explain the vascular defects observed in the placenta and yolk sac. The embryonic phenotypes caused by Siva loss are not rescued by p53 deficiency, nor do they resemble those of p53 null embryos, suggesting that the embryonic function of Siva is not related to the p53 pathway. Moreover, loss of the Ripk3 necroptosis protein does not rescue the observed lethality or developmental defects, suggesting that Siva may play a non-apoptotic role in development. Collectively, these studies reveal a key role for Siva in proper embryonic development.

Investigating human placentation and pregnancy using first trimester chorionic villi.

Chorionic villus sampling (CVS), routinely used for prenatal diagnosis of cytogenetic disorders, also possesses great potential for the study of placentation. To better understand villus biology, human placentation, and how these relate to pregnancy outcomes, we examined the morphology and transcriptomes of villi obtained via CVS from 10 to 14 weeks of pregnancy and correlated these with pregnancy attributes and clinical outcomes. First, we established a morphological scoring system based on three main villus features: branching, budding and vascularization. We then tested whether morphology scores were predictive of pregnancy attributes and clinical outcomes. Finally, we used RNA sequencing to assess the transcriptional basis of villus morphology and tested the hypothesis that gene expression may predict pregnancy outcomes. We demonstrate that villus morphology varies tremendously between patients, irrespective of gestational age, and that transcriptional differences are highly predictive of villus morphology. We show that pre-eclampsia markers are associated with villi with low morphology scores. Additionally, we identify SVEP1 as a possible biomarker for defining gestational age. Overall, chorionic villi in the first trimester remain one of the few means to correlate placental function with pregnancy outcome and these samples are a valuable and increasingly rare resource.

Characterization of mouse Rsk4 as an inhibitor of fibroblast growth factor-RAS-extracellular signal-regulated kinase signaling.

Receptor tyrosine kinase (RTK) signals regulate the specification of a varied array of tissue types by utilizing distinct modules of proteins to elicit diverse effects. The RSK proteins are part of the RTK signal transduction pathway and are thought to relay these signals by acting downstream of extracellular signal-regulated kinase (ERK). In this study we report the identification of ribosomal S6 kinase 4 (Rsk4) as an inhibitor of RTK signals. Among the RSK proteins, RTK inhibition is specific to RSK4 and, in accordance, is dependent upon a region of the RSK4 protein that is divergent from other RSK family members. We demonstrate that Rsk4 inhibits the transcriptional activation of specific targets of RTK signaling as well as the activation of ERK. Developmentally, Rsk4 is expressed in extraembryonic tissue, where RTK signals are known to have critical roles. Further examination of Rsk4 expression in the extraembryonic tissues demonstrates that its expression is inversely correlated with the presence of activated ERK 1/2. These studies demonstrate a new and divergent function for RSK4 and support a role for RSK proteins in the specification of RTK signals during early mouse development.

Desynchronizing Embryonic Cell Division Waves Reveals the Robustness of Xenopus laevis Development.

The early Xenopus laevis embryo is replete with dynamic spatial waves. One such wave, the cell division wave, emerges from the collective cell division timing of first tens and later hundreds of cells throughout the embryo. Here, we show that cell division waves do not propagate between neighboring cells and do not rely on cell-to-cell coupling to maintain their division timing. Instead, intrinsic variation in division period autonomously and gradually builds these striking patterns of cell division. Disrupting this pattern of division by placing embryos in a temperature gradient resulted in highly asynchronous entry to the midblastula transition and misexpression of the mesodermal marker Xbra. Remarkably, this gene expression defect is corrected during involution, resulting in delayed yet normal Xbra expression and viable embryos. This implies the existence of a previously unknown mechanism for normalizing mesodermal gene expression during involution.