Comparative single cell analysis reveals complex patterns of cell type and cell signaling innovations at the fetal-maternal interface.

The decidual-placental interface is one of the most diverse and rapidly evolving tissues in mammals. Its origin as a chimeric fetal-maternal tissue poses a unique evolutionary puzzle. We present single-cell RNA sequencing atlases from the fetal-maternal interfaces of the opossum, a marsupial, the Malagasy common tenrec, an afrotherian with primitive reproductive features, and mouse, guinea pig, and human. Invasive trophoblast shares a common transcriptomic signature across eutherians, which we argue represents a cell type family that radiated following the evolution of hemochorial placentation. We find evidence that the eutherian decidual stromal cell evolved stepwise from a predecidual state retained in Tenrec , followed by a second decidual cell type originating in Boreoeutheria with endocrine characteristics. We reconstruct ligand-receptor signaling to test evolutionary hypotheses at scale. Novel trophoblast and decidual cell types display strong integration into signaling networks compared to other cells. Additionally, we find consistent disambiguation between fetal and maternal signaling. Using phylogenetic analysis, we infer the cell-cell signaling network of the Placental common ancestor, and identify increased rates of signaling evolution in Euarchontoglires. Together, our findings reveal novel cell type identities and cell signaling dynamics at the mammalian fetal-maternal interface.

Reading the palimpsest of cell interactions: What questions may we ask of the data?

Biological function depends on the composition and structure of the organism, the latter describing the organization of interactions between parts. While cells in multicellular organisms are capable of a remarkable degree of autonomy, most functions do require cell communication: the coordination of functions (growth, differentiation, and apoptosis), the compartmentalization of cellular processes, and the integration of cells into higher levels of structural organization. A wealth of data on putative cell interactions has become available, yet its biological interpretation depends on our expectations about the structure of interaction networks. Here, we attempt to formulate basic questions to ask when interpreting cell interaction data. We build on the understanding that cells fulfill two general functions: the integrity-maintaining and the organismal service function. We derive the expected patterns of cell interactions considering two intertwined aspects: the functional and the evolutionary. Based on these, we propose guidelines for analysis and interpretation of transcriptional cell-interactome data.

A common allele increases endometrial Wnt4 expression, with antagonistic implications for pregnancy, reproductive cancers, and endometriosis.

The common human SNP rs3820282 is associated with multiple phenotypes including gestational length and likelihood of endometriosis and cancer, presenting a paradigmatic pleiotropic variant. Deleterious pleiotropic mutations cause the co-occurrence of disorders either within individuals, or across population. When adverse and advantageous effects are combined, pleiotropy can maintain high population frequencies of deleterious alleles. To reveal the causal molecular mechanisms of this pleiotropic SNP, we introduced this substitution into the mouse genome by CRISPR/Cas 9. Previous work showed that rs3820282 introduces a high-affinity estrogen receptor alpha-binding site at the Wnt4 locus. Here, we show that this mutation upregulates Wnt4 transcription in endometrial stroma, following the preovulatory estrogen peak. Effects on uterine transcription include downregulation of epithelial proliferation and induction of progesterone-regulated pro-implantation genes. We propose that these changes increase uterine permissiveness to embryo invasion, whereas they decrease resistance to invasion by cancer and endometriotic foci in other estrogen-responsive tissues.

Experimental and phylogenetic evidence for correlated gene expression evolution in endometrial and skin fibroblasts.

Gene expression change is a dominant mode of evolution. Mutations, however, can affect gene expression in multiple cell types. Therefore, gene expression evolution in one cell type can lead to similar gene expression changes in another cell type. Here, we test this hypothesis by investigating dermal skin fibroblasts (SFs) and uterine endometrial stromal fibroblasts (ESFs). The comparative dataset consists of transcriptomes from cultured SF and ESF of nine mammalian species. We find that evolutionary changes in gene expression in SF and ESF are highly correlated. The experimental dataset derives from a SCID mouse strain selected for slow cancer growth leading to substantial gene expression changes in SFs. We compared the gene expression profiles of SF with that of ESF and found a significant correlation between them. We discuss the implications of these findings for the evolutionary correlation between placental invasiveness and vulnerability to metastatic cancer.

Normalizing need not be the norm: count-based math for analyzing single-cell data.

Counting transcripts of mRNA are a key method of observation in modern biology. With advances in counting transcripts in single cells (single-cell RNA sequencing or scRNA-seq), these data are routinely used to identify cells by their transcriptional profile, and to identify genes with differential cellular expression. Because the total number of transcripts counted per cell can vary for technical reasons, the first step of many commonly used scRNA-seq workflows is to normalize by sequencing depth, transforming counts into proportional abundances. The primary objective of this step is to reshape the data such that cells with similar biological proportions of transcripts end up with similar transformed measurements. But there is growing concern that normalization and other transformations result in unintended distortions that hinder both analyses and the interpretation of results. This has led to an intense focus on optimizing methods for normalization and transformation of scRNA-seq data. Here, we take an alternative approach, by avoiding normalization and transformation altogether. We abandon the use of distances to compare cells, and instead use a restricted algebra, motivated by measurement theory and abstract algebra, that preserves the count nature of the data. We demonstrate that this restricted algebra is sufficient to draw meaningful and practical comparisons of gene expression through the use of the dot product and other elementary operations. This approach sidesteps many of the problems with common transformations, and has the added benefit of being simpler and more intuitive. We implement our approach in the package countland, available in python and R.

Reframing research on evolutionary novelty and co-option: Character identity mechanisms versus deep homology.

A central topic in research at the intersection of development and evolution is the origin of novel traits. Despite progress on understanding how developmental mechanisms underlie patterns of diversity in the history of life, the problem of novelty continues to challenge researchers. Here we argue that research on evolutionary novelty and the closely associated phenomenon of co-option can be reframed fruitfully by: (1) specifying a conceptual model of mechanisms that underwrite character identity, (2) providing a richer and more empirically precise notion of co-option that goes beyond common appeals to "deep homology", and (3) attending to the nature of experimental interventions that can determine whether and how the co-option of identity mechanisms can help to explain novel character origins. This reframing has the potential to channel future investigation to make substantive progress on the problem of evolutionary novelty. To illustrate this potential, we apply our reframing to two case studies: treehopper helmets and beetle horns.

The hierarchical basis of serial homology and evolutionary novelty.

Given the pervasiveness of gene sharing in evolution and the extent of homology across the tree of life, why is everything not homologous with everything else? The continuity and overlapping genetic contributions to diverse traits across lineages seem to imply that no discrete determination of homology is possible. Although some argue that the widespread overlap in parts and processes should be acknowledged as "partial" homology, this threatens a broad base of presumed comparative morphological knowledge accepted by most biologists. Following a long scientific tradition, we advocate a strategy of "theoretical articulation" that introduces further distinctions to existing concepts to produce increased contrastive resolution among the labels used to represent biological phenomena. We pursue this strategy by drawing on successful patterns of reasoning from serial homology at the level of gene sequences to generate an enriched characterization of serial homology as a hierarchical, phylogenetic concept. Specifically, we propose that the concept of serial homology should be applied primarily to repeated but developmentally individualized body parts, such as cell types, differentiated body segments, or epidermal appendages. For these characters, a phylogenetic history can be reconstructed, similar to families of paralogous genes, endowing the notion of serial homology with a hierarchical, phylogenetic interpretation. On this basis, we propose a five-fold theoretical classification that permits a more fine-grained mapping of diverse trait-types. This facilitates answering the question of why everything is not homologous with everything else, as well as how novelty is possible given that any new character possesses evolutionary precursors. We illustrate the fecundity of our account by reference to debates over insect wing serial homologs and vertebrate paired appendages.

Cooption of polyalanine tract into a repressor domain in the mammalian transcription factor HoxA11.

An enduring problem in biology is explaining how novel functions of genes originated and how those functions diverge between species. Despite detailed studies on the functional evolution of a few proteins, the molecular mechanisms by which protein functions have evolved are almost entirely unknown. Here, we show that a polyalanine tract in the homeodomain transcription factor HoxA11 arose in the stem-lineage of mammals and functions as an autonomous repressor module by physically interacting with the PAH domains of SIN3 proteins. These results suggest that long polyalanine tracts, which are common in transcription factors and often associated with disease, may tend to function as repressor domains and can contribute to the diversification of transcription factor functions despite the deleterious consequences of polyalanine tract expansion.

The value of broad taxonomic comparisons in evolutionary medicine: Disease is not a trait but a state of a trait!

In this short paper, we argue that there is a fundamental connection between the medical sciences and evolutionary biology as both are sciences of biological variation. Medicine studies pathological variation among humans (and domestic animals in veterinary medicine) and evolutionary biology studies variation within and among species in general. A key principle of evolutionary biology is that genetic differences among species have arisen first from mutations originating within populations. This implies a mechanistic continuity between variation among individuals within a species and variation between species. This fact motivates research that seeks to leverage comparisons among species to unravel the genetic basis of human disease vulnerabilities. This view also implies that genetically caused diseases can be understood as extreme states of an underlying trait, that is, an axis of variation, rather than distinct traits, as often assumed in GWAS studies. We illustrate these points with a number of examples as diverse as anatomical birth defects, cranio-facial variation, preeclampsia and vulnerability to metastatic cancer.

Paracrine HB-EGF signaling reduce enhanced contractile and energetic state of activated decidual fibroblasts by rebalancing SRF-MRTF-TCF transcriptional axis.

Multiple parallels exist between placentation and cancer dissemination at molecular, cellular, and anatomical levels, presenting placentation as a unique model to mechanistically understand the onset of cancer metastasis. In humans, interaction of placenta and the endometrium results eventually in deep invasion of placental extravillous trophoblasts (EVTs) into the maternal stroma, a process similar to stromal trespass by disseminating carcinoma cells. In anticipation of implantation, endometrial fibroblasts (ESFs) undergo a process called decidualization during the secretory phase of the menstrual cycle. Decidualization, among other substantial changes associated with ESF differentiation, also involves a component of fibroblast activation, and myofibroblast transformation. Here, using traction force microscopy, we show that increased cellular contractility in decidualized ESFs is reversed after interaction with EVTs. We also report here the large changes in energetic state of ESFs upon decidualization, showing increased oxidative phosphorylation, mitochondrial competency and ATP generation, as well as enhanced aerobic glycolysis, presenting mechanical contractility and energetic state as new functional hallmarks for decidualization. These energetic changes accompanying the marked increase in contractile force generation in decidualization were reduced in the presence of EVTs. We also show that increase in decidual contractility and mechanical resistance to invasion is achieved by SRF-MRTF transcriptional activation, achieved via increased phosphorylation of fibroblast-specific myosin light chain 9 (MYL9). EVT induced paracrine secretion of Heparin Binding Epidermal Growth Factor (HBEGF), a potent MAPK activator, which shifts the balance of SRF association away from MRTF based transcription, reducing decidual ESF contractility and mechanical resistance to placental invasion. Our results identify a new axis of intercellular communication in the placental bed modulating stromal force generation and resistance to invasion with concurrent downregulation of cellular energetics. These findings have important implications for implantation related disorders, as well as stromal control of cancer dissemination.

Evolution of higher mesenchymal CD44 expression in the human lineage: A gene linked to cancer malignancy.

CD44 is an extracellular matrix receptor implicated in cancer progression. CD44 increases the invasibility of skin (SF) and endometrial stromal fibroblasts (ESF) by cancer and trophoblast cells. We reasoned that the evolution of CD44 expression can affect both, the fetal-maternal interaction through CD44 in ESF as well as vulnerability to malignant cancer through expression in SF. We studied the evolution of CD44 expression in mammalian SF and ESF and demonstrate that in the human lineage evolved higher CD44 expression. Isoform expression in cattle and human is very similar suggesting that differences in invasibility are not due to the nature of expressed isoforms. We then asked whether the concerted gene expression increase in both cell types is due to shared regulatory mechanisms or due to cell type-specific factors. Reporter gene experiments with cells and cis-regulatory elements from human and cattle show that the difference of CD44 expression is due to cis effects as well as cell type-specific trans effects. These results suggest that the concerted expression increase is likely due to selection acting on both cell types because the evolutionary change in cell type-specific factors requires selection on cell type-specific functions. This scenario implies that the malignancy enhancing effects of elevated CD44 expression in humans likely evolved as a side-effect of positive selection on a yet unidentified other function of CD44. A possible candidate is the anti-fibrotic effect of CD44 but there are no reliable data showing that humans and primates are less fibrotic than other mammals.

Proteotype coevolution and quantitative diversity across 11 mammalian species.

Evolutionary profiling has been largely limited to the nucleotide level. Using consistent proteomic methods, we quantified proteomic and phosphoproteomic layers in fibroblasts from 11 common mammalian species, with transcriptomes as reference. Covariation analysis indicates that transcript and protein expression levels and variabilities across mammals remarkably follow functional role, with extracellular matrix-associated expression being the most variable, demonstrating strong transcriptome-proteome coevolution. The biological variability of gene expression is universal at both interindividual and interspecies scales but to a different extent. RNA metabolic processes particularly show higher interspecies versus interindividual variation. Our results further indicate that while the ubiquitin-proteasome system is strongly conserved in mammals, lysosome-mediated protein degradation exhibits remarkable variation between mammalian lineages. In addition, the phosphosite profiles reveal a phosphorylation coevolution network independent of protein abundance.

Female Genital Variation Far Exceeds that of Male Genitalia: A Review of Comparative Anatomy of Clitoris and the Female Lower Reproductive Tract in Theria.

A review of the literature on the anatomy of the lower female genital tract in therian mammals reveals, contrary to the general perception, a large amount of inter-specific variation. Variation in female external genitalia is anatomically more radical than that in the male genitalia. It includes the absence of whole anatomical units, like the cervix in many Xenarthra, or the absence of the urogenital sinus (UGS), as well as the complete spatial separation of the external clitoral parts from the genital canal (either vagina or UGS). A preliminary phylogenetic analysis shows two patterns. Some morphs are unique to early branching clades, like the absence of the cervix, while others arose multiple times independently, like the flattening out or loss of the UGS, or the extreme elongation of the clitoris. Based on available information, the ancestral eutherian configuration of the external female genitalia included a cervix, a single vaginal segment, a tubular UGS, and an unperforated clitoris close to the entrance of the genital canal. The evidence for either bilobed or unitary glandes clitorides is ambivalent. Despite the wealth of information available, many gaps in knowledge remain and will require a community-wide effort to come to a more robust model of female genital evolutionary patterns.

Histone H3K4me3 breadth in hypoxia reveals endometrial core functions and stress adaptation linked to endometriosis.

Trimethylation of histone H3 at lysine 4 (H3K4me3) is a marker of active promoters. Broad H3K4me3 promoter domains have been associated with cell type identity, but H3K4me3 dynamics upon cellular stress have not been well characterized. We assessed this by exposing endometrial stromal cells to hypoxia, which is a major cellular stress condition. We observed that hypoxia modifies the existing H3K4me3 marks and that promoter H3K4me3 breadth rather than height correlates with transcription. Broad H3K4me3 domains mark genes for endometrial core functions and are maintained or selectively extended upon hypoxia. Hypoxic extension of H3K4me3 breadth associates with stress adaptation genes relevant for the survival of endometrial cells including transcription factor KLF4, for which we found increased protein expression in the stroma of endometriosis lesions. These results substantiate the view on broad H3K4me3 as a marker of cell identity genes and reveal participation of H3K4me3 extension in cellular stress adaptation.

Tracing the cis-regulatory changes underlying the endometrial control of placental invasion.

Among eutherian (placental) mammals, placental embedding into the maternal endometrium exhibits great differences, from being deeply invasive (e.g., humans) to noninvasive (e.g., cattle). The degree of invasion of placental trophoblasts is positively correlated with the rate of cancer malignancy. Previously, we have shown that fibroblasts from different species offer different levels of resistance to the invading trophoblasts as well as to cancer cell invasion. Here we present a comparative genomic investigation revealing cis-regulatory elements underlying these interspecies differences in invasibility. We identify transcription factors that regulate proinvasibility and antiinvasibility genes in stromal cells. Using an in vitro invasibility assay combined with CRISPR-Cas9 gene knockout, we found that the transcription factors GATA2 and TFDP1 strongly influence the invasibility of endometrial and skin fibroblasts. This work identifies genomic mechanisms explaining species differences in stromal invasibility, paving the way to therapies targeting stromal characteristics to regulate placental invasion, wound healing, and cancer dissemination.

Co-option of stress mechanisms in the origin of evolutionary novelties.

It is widely accepted that stressful conditions can facilitate evolutionary change. The mechanisms elucidated thus far accomplish this with a generic increase in heritable variation that facilitates more rapid adaptive evolution, often via plastic modifications of existing characters. Through scrutiny of different meanings of stress in biological research, and an explicit recognition that stressors must be characterized relative to their effect on capacities for maintaining functional integrity, we distinguish between: (1) previously identified stress-responsive mechanisms that facilitate evolution by maintaining an adaptive fit with the environment, and (2) the co-option of stress-responsive mechanisms that are specific to stressors leading to the origin of novelties via compensation. Unlike standard accounts of gene co-option that identify component sources of evolutionary change, our model documents the cost-benefit trade-offs and thereby explains how one mechanism-an immediate response to acute stress-is transformed evolutionarily into another-routine protection from recurring stressors. We illustrate our argument with examples from cell type origination as well as processes and structures at higher levels of organization. These examples suggest a general principle of evolutionary origination based on the capacity to switch between regulatory states related to reproduction and proliferation versus survival and differentiation.

Single-cell analysis of prostaglandin E2-induced human decidual cell in vitro differentiation: a minimal ancestral deciduogenic signal†.

The decidua is a hallmark of reproduction in many placental mammals. Differentiation of decidual stromal cells is known to be induced by progesterone and the cyclic AMP/protein kinase A (cAMP/PKA) pathway. Several candidates have been identified as the physiological stimulus for adenylyl cyclase activation, but their relative importance remains unclear. To bypass this uncertainty, the standard approach for in vitro experiments uses membrane-permeable cAMP and progestin. We phylogenetically infer that prostaglandin E2 (PGE2) likely was the signal that ancestrally induced decidualization in conjunction with progesterone. This suggests that PGE2 and progestin should be able to activate the core gene regulatory network of decidual cells. To test this prediction, we performed a genome-wide study of gene expression in human endometrial fibroblasts decidualized with PGE2 and progestin. Comparison to a cAMP-based protocol revealed shared activation of core decidual genes and decreased induction of senescence-associated genes. Single-cell transcriptomics of PGE2-mediated decidualization revealed a distinct, early-activated state transitioning to a differentiated decidual state. PGE2-mediated decidualization was found to depend upon progestin-dependent induction of PGE2 receptor 2 (PTGER2) which in turn leads to PKA activation upon PGE2 stimulation. Progesterone-dependent induction of PTGER2 is absent in opossum, an outgroup taxon of placental mammals which is incapable of decidualization. Together, these findings suggest that the origin of decidualization involved the evolution of progesterone-dependent activation of the PGE2/PTGER2/PKA axis, facilitating entry into a PKA-dominant rather than AKT-dominant cellular state. We propose the use of PGE2 for in vitro decidualization as an alternative to 8-Br-cAMP.

The Coevolution of Placentation and Cancer.

Analogies between placentation, in particular the behavior of trophoblast cells, and cancer have been noted since the beginning of the twentieth century. To what degree these can be explained as a consequence of the evolution of placentation has been unclear. In this review, we conclude that many similarities between trophoblast and cancer cells are shared with other, phylogenetically older processes than placentation. The best candidates for cancer hallmarks that can be explained by the evolution of eutherian placenta are mechanisms of immune evasion. Another dimension of the maternal accommodation of the placenta with an impact on cancer malignancy is the evolution of endometrial invasibility. Species with lower degrees of placental invasion tend to have lower vulnerability to cancer malignancy. We finally identify several areas in which one could expect to see coevolutionary changes in placental and cancer biology but that, to our knowledge, have not been explored.

Molecular Evolution of CatSper in Mammals and Function of Sperm Hyperactivation in Gray Short-Tailed Opossum.

Males have evolved species-specifical sperm morphology and swimming patterns to adapt to different fertilization environments. In eutherians, only a small fraction of the sperm overcome the diverse obstacles in the female reproductive tract and successfully migrate to the fertilizing site. Sperm arriving at the fertilizing site show hyperactivated motility, a unique motility pattern displaying asymmetric beating of sperm flagella with increased amplitude. This motility change is triggered by Ca2+ influx through the sperm-specific ion channel, CatSper. However, the current understanding of the CatSper function and its molecular regulation is limited in eutherians. Here, we report molecular evolution and conservation of the CatSper channel in the genome throughout eutherians and marsupials. Sequence analyses reveal that CatSper proteins are slowly evolved in marsupials. Using an American marsupial, gray short-tailed opossum (Monodelphis domestica), we demonstrate the expression of CatSper in testes and its function in hyperactivation and unpairing of sperm. We demonstrate that a conserved IQ-like motif in CatSperζ is required for CatSperζ interaction with the pH-tuned Ca2+ sensor, EFCAB9, for regulating CatSper activity. Recombinant opossum EFCAB9 can interact with mouse CatSperζ despite high sequence divergence of CatSperζ among CatSper subunits in therians. Our finding suggests that molecular characteristics and functions of CatSper are evolutionarily conserved in gray short-tailed opossum, unraveling the significance of sperm hyperactivation and fertilization in marsupials for the first time.