Human embryo implantation.

Embryo implantation in humans is interstitial, meaning the entire conceptus embeds in the endometrium before the placental trophoblast invades beyond the uterine mucosa into the underlying inner myometrium. Once implanted, embryo survival pivots on the transformation of the endometrium into an anti-inflammatory placental bed, termed decidua, under homeostatic control of uterine natural killer cells. Here, we examine the evolutionary context of embryo implantation and elaborate on uterine remodelling before and after conception in humans. We also discuss the interactions between the embryo and the decidualising endometrium that regulate interstitial implantation and determine embryo fitness. Together, this Review highlights the precarious but adaptable nature of the implantation process.

Human Leukocyte Antigen F Presents Peptides and Regulates Immunity through Interactions with NK Cell Receptors.

Evidence is mounting that the major histocompatibility complex (MHC) molecule HLA-F (human leukocyte antigen F) regulates the immune system in pregnancy, infection, and autoimmunity by signaling through NK cell receptors (NKRs). We present structural, biochemical, and evolutionary analyses demonstrating that HLA-F presents peptides of unconventional length dictated by a newly arisen mutation (R62W) that has produced an open-ended groove accommodating particularly long peptides. Compared to empty HLA-F open conformers (OCs), HLA-F tetramers bound with human-derived peptides differentially stained leukocytes, suggesting peptide-dependent engagement. Our in vitro studies confirm that NKRs differentiate between peptide-bound and peptide-free HLA-F. The complex structure of peptide-loaded ß2m-HLA-F bound to the inhibitory LIR1 revealed similarities to high-affinity recognition of the viral MHC-I mimic UL18 and a docking strategy that relies on contacts with HLA-F as well as ß2m, thus precluding binding to HLA-F OCs. These findings provide a biochemical framework to understand how HLA-F could regulate immunity via interactions with NKRs.

Interpreting the Evolutionary Echoes of a Protein Complex Essential for Inner-Ear Mechanosensation.

The sensory epithelium of the inner ear, found in all extant lineages of vertebrates, has been subjected to over 500 million years of evolution, resulting in the complex inner ear of modern vertebrates. Inner-ear adaptations are as diverse as the species in which they are found, and such unique anatomical variations have been well studied. However, the evolutionary details of the molecular machinery that is required for hearing are less well known. Two molecules that are essential for hearing in vertebrates are cadherin-23 and protocadherin-15, proteins whose interaction with one another acts as the focal point of force transmission when converting sound waves into electrical signals that the brain can interpret. This "tip-link" interaction exists in every lineage of vertebrates, but little is known about the structure or mechanical properties of these proteins in most non-mammalian lineages. Here, we use various techniques to characterize the evolution of this protein interaction. Results show how evolutionary sequence changes in this complex affect its biophysical properties both in simulations and experiments, with variations in interaction strength and dynamics among extant vertebrate lineages. Evolutionary simulations also characterize how the biophysical properties of the complex in turn constrain its evolution and provide a possible explanation for the increase in deafness-causing mutants observed in cadherin-23 relative to protocadherin-15. Together, these results suggest a general picture of tip-link evolution in which selection acted to modify the tip-link interface, although subsequent neutral evolution combined with varying degrees of purifying selection drove additional diversification in modern tetrapods.

Is there a loophole in Dollo's law? A DevoEvo perspective on irreversibility (of felid dentition).

There is a longstanding interest in whether the loss of complex characters is reversible (so-called "Dollo's law"). Reevolution has been suggested for numerous traits but among the first was Kurtén, who proposed that the presence of the second lower molar (M2 ) of the Eurasian lynx (Lynx lynx) was a violation of Dollo's law because all other Felids lack M2 . While an early and often cited example for the reevolution of a complex trait, Kurtén and Werdelin used an ad hoc parsimony argument to support their. Here I revisit the evidence that M2 reevolved lynx using explicit parsimony and maximum likelihood models of character evolution and find strong evidence that Kurtén and Werdelin were correct-M2 reevolved in E. lynx. Next, I explore the developmental mechanisms which may explain this violation of Dollo's law and suggest that the reevolution of lost complex traits may arise from the reevolution of cis-regulatory elements and protein-protein interactions, which have a longer half-life after silencing that protein coding genes. Finally, I present a developmental model to explain the reevolution M2 in E. lynx, which suggest that the developmental programs required for the establishment of serially homologous characters may never really be lost so long as a single instance of the character remains-thus the gain and loss and regain of serially homologous characters, such mammalian molars, may be developmentally and evolutionarily "simple."

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.

Relaxed constraint and functional divergence of the progesterone receptor (PGR) in the human stem-lineage.

The steroid hormone progesterone, acting through the progesterone receptor (PR), a ligand-activated DNA-binding transcription factor, plays an essential role in regulating nearly every aspect of female reproductive biology. While many reproductive traits regulated by PR are conserved in mammals, Catarrhine primates evolved several derived traits including spontaneous decidualization, menstruation, and a divergent (and unknown) parturition signal, suggesting that PR may also have evolved divergent functions in Catarrhines. There is conflicting evidence, however, whether the progesterone receptor gene (PGR) was positively selected in the human lineage. Here we show that PGR evolved rapidly in the human stem-lineage (as well as other Catarrhine primates), which likely reflects an episode of relaxed selection intensity rather than positive selection. Coincident with the episode of relaxed selection intensity, ancestral sequence resurrection and functional tests indicate that the major human PR isoforms (PR-A and PR-B) evolved divergent functions in the human stem-lineage. These results suggest that the regulation of progesterone signaling by PR-A and PR-B may also have diverged in the human lineage and that non-human animal models of progesterone signaling may not faithfully recapitulate human biology.

Behavior-related gene regulatory networks: A new level of organization in the brain.

Neuronal networks are the standard heuristic model today for describing brain activity associated with animal behavior. Recent studies have revealed an extensive role for a completely distinct layer of networked activities in the brain-the gene regulatory network (GRN)-that orchestrates expression levels of hundreds to thousands of genes in a behavior-related manner. We examine emerging insights into the relationships between these two types of networks and discuss their interplay in spatial as well as temporal dimensions, across multiple scales of organization. We discuss properties expected of behavior-related GRNs by drawing inspiration from the rich literature on GRNs related to animal development, comparing and contrasting these two broad classes of GRNs as they relate to their respective phenotypic manifestations. Developmental GRNs also represent a third layer of network biology, playing out over a third timescale, which is believed to play a crucial mediatory role between neuronal networks and behavioral GRNs. We end with a special emphasis on social behavior, discuss whether unique GRN organization and cis-regulatory architecture underlies this special class of behavior, and review literature that suggests an affirmative answer.

A missense variant in FTCD is associated with arsenic metabolism and toxicity phenotypes in Bangladesh.

Inorganic arsenic (iAs) is a carcinogen, and exposure to iAs via food and water is a global public health problem. iAs-contaminated drinking water alone affects >100 million people worldwide, including ~50 million in Bangladesh. Once absorbed into the blood stream, most iAs is converted to mono-methylated (MMA) and then di-methylated (DMA) forms, facilitating excretion in urine. Arsenic metabolism efficiency varies among individuals, in part due to genetic variation near AS3MT (arsenite methyltransferase; 10q24.32). To identify additional arsenic metabolism loci, we measured protein-coding variants across the human exome for 1,660 Bangladeshi individuals participating in the Health Effects of Arsenic Longitudinal Study (HEALS). Among the 19,992 coding variants analyzed exome-wide, the minor allele (A) of rs61735836 (p.Val101Met) in exon 3 of FTCD (formiminotransferase cyclodeaminase) was associated with increased urinary iAs% (P = 8x10-13), increased MMA% (P = 2x10-16) and decreased DMA% (P = 6x10-23). Among 2,401 individuals with arsenic-induced skin lesions (an indicator of arsenic toxicity and cancer risk) and 2,472 controls, carrying the low-efficiency A allele (frequency = 7%) was associated with increased skin lesion risk (odds ratio = 1.35; P = 1x10-5). rs61735836 is in weak linkage disequilibrium with all nearby variants. The high-efficiency/major allele (G/Valine) is human-specific and eliminates a start codon at the first 5´-proximal Kozak sequence in FTCD, suggesting selection against an alternative translation start site. FTCD is critical for catabolism of histidine, a process that generates one-carbon units that can enter the one-carbon/folate cycle, which provides methyl groups for arsenic metabolism. In our study population, FTCD and AS3MT SNPs together explain ~10% of the variation in DMA% and support a causal effect of arsenic metabolism efficiency on arsenic toxicity (i.e., skin lesions). In summary, this work identifies a coding variant in FTCD associated with arsenic metabolism efficiency, providing new evidence supporting the established link between one-carbon/folate metabolism and arsenic toxicity.

Correction: A missense variant in FTCD is associated with arsenic metabolism and toxicity phenotypes in Bangladesh.

[This corrects the article DOI: 10.1371/journal.pgen.1007984.].

An Ancient Fecundability-Associated Polymorphism Creates a GATA2 Binding Site in a Distal Enhancer of HLA-F.

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, are heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits have been challenging. Here, we explore the functional significance and evolutionary history of a G/A polymorphism at SNP rs2523393, which is an eQTL for HLA-F and is significantly associated with fecundability (the probability of being pregnant within a single menstrual cycle). We replicated the association between the rs2523393 genotype and HLA-F expression by using GTEx data and demonstrate that HLA-F is upregulated in the endometrium during the window of implantation and by progesterone in decidual stromal cells. Next, we show that the rs2523393 A allele creates a GATA2 binding site in a progesterone-responsive distal enhancer that loops to the HLA-F promoter. Remarkably, we found that the A allele is derived in the human lineage and that the G/A polymorphism arose before the divergence of modern and archaic humans and segregates at intermediate to high frequencies across human populations. Remarkably, the derived A allele is has also been identified in a GWAS as a risk allele for multiple sclerosis. These data suggest that the polymorphism is maintained by antagonistic pleiotropy and a reproduction-health tradeoff in human evolution.

Transposable elements are the primary source of novelty in primate gene regulation.

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

An Ancient Fecundability-Associated Polymorphism Switches a Repressor into an Enhancer of Endometrial TAP2 Expression.

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, is heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits has been challenging. Here, we explore the functional significance and evolutionary history of a T/C polymorphism of SNP rs2071473, which we have previously shown is an eQTL for TAP2 and significantly associated with fecundability (time to pregnancy). We replicated the association between the rs2071473 genotype and TAP2 expression by using GTEx data and demonstrated that TAP2 is expressed by decidual stromal cells at the maternal-fetal interface. Next, we showed that rs2071473 is located within a progesterone-responsive cis-regulatory element that functions as a repressor with the T allele and an enhancer with the C allele. Remarkably, we found that this polymorphism arose before the divergence of modern and archaic humans, segregates at intermediate to high frequencies across human populations, and has genetic signatures of long-term balancing selection. This variant has also previously been identified in genome-wide association studies of immune-related disease, suggesting that both alleles are maintained as a result of antagonistic pleiotropy.

Expression Quantitative Trait Locus Mapping Studies in Mid-secretory Phase Endometrial Cells Identifies HLA-F and TAP2 as Fecundability-Associated Genes.

Fertility traits in humans are heritable, however, little is known about the genes that influence reproductive outcomes or the genetic variants that contribute to differences in these traits between individuals, particularly women. To address this gap in knowledge, we performed an unbiased genome-wide expression quantitative trait locus (eQTL) mapping study to identify common regulatory (expression) single nucleotide polymorphisms (eSNPs) in mid-secretory endometrium. We identified 423 cis-eQTLs for 132 genes that were significant at a false discovery rate (FDR) of 1%. After pruning for strong LD (r2 >0.95), we tested for associations between eSNPs and fecundability (the ability to get pregnant), measured as the length of the interval to pregnancy, in 117 women. Two eSNPs were associated with fecundability at a FDR of 5%; both were in the HLA region and were eQTLs for the TAP2 gene (P = 1.3x10-4) and the HLA-F gene (P = 4.0x10-4), respectively. The effects of these SNPs on fecundability were replicated in an independent sample. The two eSNPs reside within or near regulatory elements in decidualized human endometrial stromal cells. Our study integrating eQTL mapping in a primary tissue with association studies of a related phenotype revealed novel genes and associated alleles with independent effects on fecundability, and identified a central role for two HLA region genes in human implantation success.

Regulatory evolution through divergence of a phosphoswitch in the transcription factor CEBPB.

There is an emerging consensus that gene regulation evolves through changes in cis-regulatory elements and transcription factors. Although it is clear how nucleotide substitutions in cis-regulatory elements affect gene expression, it is not clear how amino-acid substitutions in transcription factors influence gene regulation. Here we show that amino-acid changes in the transcription factor CCAAT/enhancer binding protein-ß (CEBPB, also known as C/EBP-ß) in the stem-lineage of placental mammals changed the way it responds to cyclic AMP/protein kinase A (cAMP/PKA) signalling. By functionally analysing resurrected ancestral proteins, we identify three amino-acid substitutions in an internal regulatory domain of CEBPB that are responsible for the novel function. These amino-acid substitutions reorganize the location of key phosphorylation sites, introducing a new site and removing two ancestral sites, reversing the response of CEBPB to GSK-3ß-mediated phosphorylation from repression to activation. We conclude that changing the response of transcription factors to signalling pathways can be an important mechanism of gene regulatory evolution.

Transposable elements are major contributors to the origin, diversification, and regulation of vertebrate long noncoding RNAs.

Advances in vertebrate genomics have uncovered thousands of loci encoding long noncoding RNAs (lncRNAs). While progress has been made in elucidating the regulatory functions of lncRNAs, little is known about their origins and evolution. Here we explore the contribution of transposable elements (TEs) to the makeup and regulation of lncRNAs in human, mouse, and zebrafish. Surprisingly, TEs occur in more than two thirds of mature lncRNA transcripts and account for a substantial portion of total lncRNA sequence (~30% in human), whereas they seldom occur in protein-coding transcripts. While TEs contribute less to lncRNA exons than expected, several TE families are strongly enriched in lncRNAs. There is also substantial interspecific variation in the coverage and types of TEs embedded in lncRNAs, partially reflecting differences in the TE landscapes of the genomes surveyed. In human, TE sequences in lncRNAs evolve under greater evolutionary constraint than their non-TE sequences, than their intronic TEs, or than random DNA. Consistent with functional constraint, we found that TEs contribute signals essential for the biogenesis of many lncRNAs, including ~30,000 unique sites for transcription initiation, splicing, or polyadenylation in human. In addition, we identified ~35,000 TEs marked as open chromatin located within 10 kb upstream of lncRNA genes. The density of these marks in one cell type correlate with elevated expression of the downstream lncRNA in the same cell type, suggesting that these TEs contribute to cis-regulation. These global trends are recapitulated in several lncRNAs with established functions. Finally a subset of TEs embedded in lncRNAs are subject to RNA editing and predicted to form secondary structures likely important for function. In conclusion, TEs are nearly ubiquitous in lncRNAs and have played an important role in the lineage-specific diversification of vertebrate lncRNA repertoires.

Evolutionary forward genomics reveals novel insights into the genes and pathways dysregulated in recurrent early pregnancy loss.

STUDY QUESTION: Are the genes that gained novel expression in the endometria of Eutherian (placental) mammals more likely to be dysregulated in patients with endometrial-associated recurrent early pregnancy loss (REPL)? SUMMARY ANSWER: There was a significant enrichment of genes dysregulated in REPL patients among the Eutherian-specific endometrial genes. WHAT IS KNOWN ALREADY: Pregnancy loss is the most common complication of human pregnancy. REPL has multiple etiologies, including dysregulation of endometrial function, leading to 'suboptimal' implantation. Although the implantation process is tightly regulated in Eutherian (placental) mammals, the molecular factors contributing to dysregulated endometrial gene expression patterns in women with REPL are largely unknown. STUDY DESIGN, SIZE, DURATION: Endometrial biopsies were obtained from 32 REPL patients during the mid-luteal phase, and evaluated for glandular development arrest based on elevated nuclear cyclin E levels in gland cells, and for out-of-phase endometrial development based on histology. Gene expression levels were measured using Illumina Human HT-12v4 BeadChip arrays. PARTICIPANTS/MATERIALS, SETTING, METHODS: Differentially expressed genes were identified between patients with (i) out-of-phase (n = 10) versus normal (n = 22) histological dating and (ii) abnormally elevated (n = 9) versus normal (n = 23) cyclin E levels in the nuclei of endometrial glands, using a likelihood ratio test. Enrichment of dysregulated genes in REPL endometria among Eutherian-specific genes was tested by permutation. Gene ontology and pathway enrichment analyses were carried out for the dysregulated genes. MAIN RESULTS AND THE ROLE OF CHANCE: Fifty-eight and eighty-one genes were identified as differentially expressed at P < 0.001 in women with out-of-phase histological dating and abnormally elevated glandular cyclin E levels, respectively. Genes that were recruited into endometrial expression during the evolution of pregnancy in Eutherian mammals were significantly enriched for dysregulated genes (P = 0.002 for histology, P = 0.021 for cyclin E), as well as for genes involved in immune response and signaling pathways with essential roles in implantation and endometrial biology. LIMITATIONS, REASONS FOR CAUTION: Small sample size limits the statistical power to detect dysregulated genes, and the lack of non-REPL control women does not allow us to test for the contribution of these genes to overall risk of REPL. WIDER IMPLICATIONS OF THE FINDINGS: Enrichment of functional gene categories, as well as genes gained expression in the Eutherian endometria, help to identify molecular etiologies that contribute to normal functioning of the endometrium. These pathways are also strong candidates for successful pregnancy outcomes. Using the evolutionary history of mammalian gene expression in the endometrial tissue may be a promising approach to discover genes involved in female reproductive disorders. STUDY FUNDING/COMPETING INTERESTS: This work is supported by National Institutes of Health (NIH) grant R01 HD21244 to C.O. Authors declare no competing interests.

Evolution of a derived protein-protein interaction between HoxA11 and Foxo1a in mammals caused by changes in intramolecular regulation.

Current models of developmental evolution suggest changes in gene regulation underlie the evolution of morphology. Despite the fact that protein complexes regulate gene expression, the evolution of regulatory protein complexes is rarely studied. Here, we investigate the evolution of a protein-protein interaction (PPI) between Homeobox A11 (HoxA11) and Forkhead box 01A (Foxo1a). Using extant and "resurrected" ancestral proteins, we show that the physical interaction between HoxA11 and Foxo1a originated in the mammalian stem lineage. Functional divergence tests and coimmunoprecipitation with heterologous protein pairs indicate that the evolution of interaction was attributable to changes in HoxA11, and deletion studies demonstrate that the interaction interface is located in the homeodomain region of HoxA11. However, there are no changes in amino acid sequence in the homeodomain region during this time period, indicating that the origin of the derived PPI was attributable to changes outside the binding interface. We infer that the amino acid substitutions in HoxA11 altered Foxo1a's access to the conserved binding interface at the HoxA11 homeodomain. We also found an expansion in the number of paired Hox/Fox binding sites in the genomes of mammalian lineage species suggesting the complex has a biological function. Our data indicate that the physical interaction between HoxA11 and Foxo1a evolved through noninterface changes that facilitate the PPI, which prevents inappropriate interactions, rather than through the evolution of a novel binding interface. We speculate that evolutionary changes of intramolecular regulation have limited pleiotropic effects compared with changes to interaction domains themselves.

Rapid evolution of genes with anti-cancer functions during the origins of large bodies and cancer resistance in elephants.

Elephants have emerged as a model system to study the evolution of body size and cancer resistance because, despite their immense size, they have a very low prevalence of cancer. Previous studies have found that duplication of tumor suppressors at least partly contributes to the evolution of anti-cancer cellular phenotypes in elephants. Still, many other mechanisms must have contributed to their augmented cancer resistance. Here, we use a suite of codon-based maximum-likelihood methods and a dataset of 13,310 protein-coding gene alignments from 261 Eutherian mammals to identify positively selected and rapidly evolving elephant genes. We found 496 genes (3.73% of alignments tested) with statistically significant evidence for positive selection and 660 genes (4.96% of alignments tested) that likely evolved rapidly in elephants. Positively selected and rapidly evolving genes are statistically enriched in gene ontology terms and biological pathways related to regulated cell death mechanisms, DNA damage repair, cell cycle regulation, epidermal growth factor receptor (EGFR) signaling, and immune functions, particularly neutrophil granules and degranulation. All of these biological factors are plausibly related to the evolution of cancer resistance. Thus, these positively selected and rapidly evolving genes are promising candidates for genes contributing to elephant-specific traits, including the evolution of molecular and cellular characteristics that enhance cancer resistance.

Pervasive loss of regulated necrotic cell death genes in elephants, hyraxes, and sea cows (Paenungualta).

Gene loss can promote phenotypic differences between species, for example, if a gene constrains phenotypic variation in a trait, its loss allows for the evolution of a greater range of variation or even new phenotypes. Here, we explore the contribution of gene loss to the evolution of large bodies and augmented cancer resistance in elephants. We used genomes from 17 Afrotherian and Xenarthran species to identify lost genes, i.e., genes that have pseudogenized or been completely lost, and Dollo parsimony to reconstruct the evolutionary history of gene loss across species. We unexpectedly discovered a burst of gene losses in the Afrotherian stem lineage and found that the loss of genes with functions in regulated necrotic cell death modes was pervasive in elephants, hyraxes, and sea cows (Paenungulata). Among the lost genes are MLKL and RIPK3, which mediate necroptosis, and sensors that activate inflammasomes to induce pyroptosis, including AIM2, MEFV, NLRC4, NLRP1, and NLRP6. These data suggest that the mechanisms that regulate necrosis and pyroptosis are either extremely derived or potentially lost in these lineages, which may contribute to the repeated evolution of large bodies and cancer resistance in Paenungulates as well as susceptibility to pathogen infection.

Convergent evolution of endometrial prolactin expression in primates, mice, and elephants through the independent recruitment of transposable elements.

Prolactin (PRL) is a multifunctional signaling molecule best known for its role in regulating lactation in mammals. Systemic PRL is produced by the anterior pituitary, but extrapituitary PRL has also been detected in many tissues including the human endometrium. Prolactin is essential for pregnancy in rodents and one of the most dramatically induced genes in the endometrium during human pregnancy. The promoter for human endometrial Prl is located about 5.8 kb upstream of the pituitary promoter and is derived from a transposable element called MER39. Although it has been shown that prolactin is expressed in the pregnant endometrium of a few mammals other than humans, MER39 has been described as primate specific. Thus, in an effort to understand mechanisms of prolactin regulatory evolution, we sought to determine how uterine prolactin is transcribed in species that lack MER39. Using a variety of complementary strategies, including reverse transcriptase-polymerase chain reaction, 5' rapid amplification of cDNA ends, and whole-transcriptome sequencing, we show that endometrial Prl expression is not a shared character of all placental mammals, as it is not expressed in rabbits, pigs, dogs, or armadillos. We show that in primates, mice, and elephants, prolactin mRNA is transcribed in the pregnant endometrium from alternative promoters, different from the pituitary promoter and different from each other. Moreover, we demonstrate that the spider monkey promoter derives from the long terminal repeat (LTR) element MER39 as in humans, the mouse promoter derives from the LTR element MER77, and the elephant promoter derives from the lineage-specific LINE retrotransposon L1-2_LA. We also find surprising variation of transcriptional start sites within these transposable elements and of Prl splice variants, suggesting a high degree of flexibility in the promoter architecture even among closely related species. Finally, the three groups shown here to express endometrial prolactin-the higher primates, the rodents, and the elephant-represent three of the four lineages that showed adaptive evolution of the Prl gene in an earlier study (Wallis M. 2000. Episodic evolution of protein hormones: molecular evolution of pituitary prolactin. J Mol Evol. 50:465-473), which supports our findings and suggests that the selective forces responsible for accelerated Prl evolution were in the endometrium. This is the first reported case of convergent evolution of gene expression through the independent recruitment of different transposable elements, highlighting the importance of transposable elements in gene regulatory, and potentially adaptive, evolution.