Marsupials have monoallelic MEST expression with a conserved antisense lncRNA but MEST is not imprinted.

The imprinted isoform of the Mest gene in mice is involved in key mammalian traits such as placental and fetal growth, maternal care and mammary gland maturation. The imprinted isoform has a distinct differentially methylated region (DMR) at its promoter in eutherian mammals but in marsupials, there are no differentially methylated CpG islands between the parental alleles. Here, we examined similarities and differences in the MEST gene locus across mammals using a marsupial, the tammar wallaby, a monotreme, the platypus, and a eutherian, the mouse, to investigate how imprinting of this gene evolved in mammals. By confirming the presence of the short isoform in all mammalian groups (which is imprinted in eutherians), this study suggests that an alternative promoter for the short isoform evolved at the MEST gene locus in the common ancestor of mammals. In the tammar, the short isoform of MEST shared the putative promoter CpG island with an antisense lncRNA previously identified in humans and an isoform of a neighbouring gene CEP41. The antisense lncRNA was expressed in tammar sperm, as seen in humans. This suggested that the conserved lncRNA might be important in the establishment of MEST imprinting in therian mammals, but it was not imprinted in the tammar. In contrast to previous studies, this study shows that MEST is not imprinted in marsupials. MEST imprinting in eutherians, therefore must have occurred after the marsupial-eutherian split with the acquisition of a key epigenetic imprinting control region, the differentially methylated CpG islands between the parental alleles.

The evolutionary loss of the Eh1 motif in FoxE1 in the lineage of placental mammals.

Forkhead box E1 (FoxE1) protein is a transcriptional regulator known to play a major role in the development of the thyroid gland. By performing sequence alignments, we detected a deletion in FoxE1, which occurred in the evolution of mammals, near the point of divergence of placental mammals. This deletion led to the loss of the majority of the Eh1 motif, which was important for interactions with transcriptional corepressors. To investigate a potential mechanism for this deletion, we analyzed replication through the deletion area in mammalian cells with two-dimensional gel electrophoresis, and in vitro, using a primer extension reaction. We demonstrated that the area of the deletion presented an obstacle for replication in both assays. The exact position of polymerization arrest in primer extension indicated that it was most likely caused by a quadruplex DNA structure. The quadruplex structure hypothesis is also consistent with the exact borders of the deletion. The exact roles of these evolutionary changes in FoxE1 family proteins are still to be determined.

Patterns of Evolution of TRIM Genes Highlight the Evolutionary Plasticity of Antiviral Effectors in Mammals.

The innate immune system of mammals is formed by a complex web of interacting proteins, which together constitute the first barrier of entry for infectious pathogens. Genes from the E3-ubiquitin ligase tripartite motif (TRIM) family have been shown to play an important role in the innate immune system by restricting the activity of different retrovirus species. For example, TRIM5 and TRIM22 have both been associated with HIV restriction and are regarded as crucial parts of the antiretroviral machinery of mammals. Our analyses of positive selection corroborate the great significance of these genes for some groups of mammals. However, we also show that many species lack TRIM5 and TRIM22 altogether. By analyzing a large number of mammalian genomes, here we provide the first comprehensive view of the evolution of these genes in eutherians, showcasing that the pattern of accumulation of TRIM genes has been dissimilar across mammalian orders. Our data suggest that these differences are caused by the evolutionary plasticity of the immune system of eutherians, which have adapted to use different strategies to combat retrovirus infections. Altogether, our results provide insights into the dissimilar evolution of a representative family of restriction factors, highlighting an example of adaptive and idiosyncratic evolution in the innate immune system.

Lack of evidence for coevolution between oxytocin receptor N-terminal variants and monogamy in placental mammals.

Oxytocin (OXT) is a neurohypophyseal hormone that influences a wide range of affiliative behaviors, such as pair-bonding and infant care, across mammals. The effects of OXT depend significantly on an adequate interaction with its receptor, OXTR. OXTR belongs to the G-protein coupled receptor family. The extracellular N-terminal domain of OXTR interacts with the linear C-terminal tail of OXT and is required for OXT binding. Across mammalian species there is a genetic diversity in OXTR terminal sequence. Previous work on primates has shown an association between OXTR phylogeny and monogamy. However, it is not clear whether this variation coevolved with either mating system (monogamy) or infant care behaviors (such as allomaternal care). Here, we take a phylogenetic comparative and evolutionary modeling approach across a wide range of placental mammals (n = 60) to test whether OXTR N-terminal variants co-evolved with either monogamy or allomaternal care behaviors. Our results indicate that the diversity in OXTR N-terminal region is unlikely to provide the underlying genetic bases for variation in mating system and/or allomaternal behavior as we find no evidence for co-evolution between protein sequence and affiliative behaviors. Hence, the role played by OXT in influencing affiliative behaviors is unlikely to be mediated by the genetic diversity of its receptor.

Placental Mammals Acquired Functional Sequences in NRK for Regulating the CK2-PTEN-AKT Pathway and Placental Cell Proliferation.

The molecular evolution processes underlying the acquisition of the placenta in eutherian ancestors are not fully understood. Mouse NCK-interacting kinase (NIK)-related kinase (NRK) is expressed highly in the placenta and plays a role in preventing placental hyperplasia. Here, we show the molecular evolution of NRK, which confers its function for inhibiting placental cell proliferation. Comparative genome analysis identified NRK orthologs across vertebrates, which share the kinase and citron homology (CNH) domains. Evolutionary analysis revealed that NRK underwent extensive amino acid substitutions in the ancestor of placental mammals and has been since conserved. Biochemical analysis of mouse NRK revealed that the CNH domain binds to phospholipids, and a region in NRK binds to and inhibits casein kinase-2 (CK2), which we named the CK2-inhibitory region (CIR). Cell culture experiments suggest the following: 1) Mouse NRK is localized at the plasma membrane via the CNH domain, where the CIR inhibits CK2. 2) This mitigates CK2-dependent phosphorylation and inhibition of PTEN and 3) leads to the inhibition of AKT signaling and cell proliferation. Nrk deficiency increased phosphorylation levels of PTEN and AKT in mouse placenta, supporting our hypothesis. Unlike mouse NRK, chicken NRK did not bind to phospholipids and CK2, decrease phosphorylation of AKT, or inhibit cell proliferation. Both the CNH domain and CIR have evolved under purifying selection in placental mammals. Taken together, our study suggests that placental mammals acquired the phospholipid-binding CNH domain and CIR in NRK for regulating the CK2-PTEN-AKT pathway and placental cell proliferation.

Whole embryonic detection of maternal microchimeric cells highlights significant differences in their numbers among individuals.

During pregnancy in placental mammals, small numbers of maternal cells (maternal microchimeric cells, or MMc cells) migrate into the fetus and persist decades, or perhaps for the rest of their lives, and higher frequencies of MMc cells are reported to correlate with variety of phenomena, such as immune tolerance, tissue repair, and autoimmune diseases. While detection of these MMc cells is considered in all pregnancies, their frequency differs largely according to tissue type and disease cases, and it remains unclear whether the number of MMc cells differs significantly among embryos in normal pregnancies. Here, for the first time, we developed a whole embryonic detection method for MMc cells using transgenic mice and counted live MMc cells in each individual embryo. Using this technique, we found that the number of MMc cells was comparable in most of the analyzed embryos; however, around 500 times higher number of MMc cells was detected in one embryo at the latest stage. This result suggests that the number of MMc cells could largely differ in rare cases with unknown underlying mechanisms. Our methodology provides a basis for testing differences in the numbers of MMc cells among individual embryos and for analyzing differences in MMc cell type repertoires in future studies. These data could provide a hint toward understanding the mechanisms underlying the variety of apparently inconsistent MMc-related phenomena.

A eutherian-specific microRNA controls the translation of Satb2 in a model of cortical differentiation.

Cerebral cortical development is controlled by key transcription factors that specify the neuronal identities in the different layers. The mechanisms controlling their expression in distinct cells are only partially known. We investigated the expression and stability of Tbr1, Bcl11b, Fezf2, Satb2, and Cux1 mRNAs in single developing mouse cortical cells. We observe that Satb2 mRNA appears much earlier than its protein and in a set of cells broader than expected, suggesting an initial inhibition of its translation, subsequently released during development. Mechanistically, Satb2 3'UTR modulates protein translation of GFP reporters during mouse corticogenesis. We select miR-541, a eutherian-specific miRNA, and miR-92a/b as the best candidates responsible for SATB2 inhibition, being strongly expressed in early and reduced in late progenitor cells. Their inactivation triggers robust and premature SATB2 translation in both mouse and human cortical cells. Our findings indicate RNA interference as a major mechanism in timing cortical cell identities.

The role of eutherian-specific RTL1 in the nervous system and its implications for the Kagami-Ogata and Temple syndromes.

RTL1 (also termed paternal expressed 11 (PEG11)) is considered the major imprinted gene responsible for the placental and fetal/neonatal muscle defects that occur in the Kagami-Ogata and Temple syndromes (KOS14 and TS14, respectively). However, it remains elusive whether RTL1 is also involved in their neurological symptoms, such as behavioral and developmental delay/intellectual disability, feeding difficulties, motor delay, and delayed speech. Here, we demonstrate that the mouse RTL1 protein is widely expressed in the central nervous system (CNS), including the limbic system. Importantly, two disease model mice with over- and under-expression of Rtl1 exhibited reduced locomotor activity, increased anxiety, and impaired amygdala-dependent cued fear, demonstrating that Rtl1 also plays an important role in the CNS. These results indicate that the KOS14 and TS14 are neuromuscular as well as neuropsychiatric diseases caused by irregular CNS RTL1 expression, presumably leading to impaired innervation of motor neurons to skeletal muscles as well as malfunction of the hippocampus-amygdala complex. It is of considerable interest that eutherian-specific RTL1 is expressed in mammalian- and eutherian-specific brain structures, that is, the corticospinal tract and corpus callosum, respectively, suggesting that RTL1 might have contributed to the acquisition of both these structures themselves and fine motor skill in eutherian brain evolution.

Signatures of conserved and unique molecular features in Afrotheria.

Afrotheria is a clade of African-origin species with striking dissimilarities in appearance and habitat. In this study, we compared whole proteome sequences of six Afrotherian species to obtain a broad viewpoint of their underlying molecular make-up, to recognize potentially unique proteomic signatures. We find that 62% of the proteomes studied here, predominantly involved in metabolism, are orthologous, while the number of homologous proteins between individual species is as high as 99.5%. Further, we find that among Afrotheria, L. africana has several orphan proteins with 112 proteins showing < 30% sequence identity with their homologues. Rigorous sequence searches and complementary approaches were employed to annotate 156 uncharacterized protein sequences and 28 species-specific proteins. For 122 proteins we predicted potential functional roles, 43 of which we associated with protein- and nucleic-acid binding roles. Further, we analysed domain content and variations in their combinations within Afrotheria and identified 141 unique functional domain architectures, highlighting proteins with potential for specialized functions. Finally, we discuss the potential relevance of highly represented protein families such as MAGE-B2, olfactory receptor and ribosomal proteins in L. africana and E. edwardii, respectively. Taken together, our study reports the first comparative study of the Afrotherian proteomes and highlights salient molecular features.

Folding Stability of Pax9 Intronic G-Quadruplex Correlates with Relative Molar Size in Eutherians.

Eutherian dentition has been the focus of a great deal of studies in the areas of evolution, development, and genomics. The development of molar teeth is regulated by an antero-to-posterior cascade mechanism of activators and inhibitors molecules, where the relative sizes of the second (M2) and third (M3) molars are dependent of the inhibitory influence of the first molar (M1). Higher activator/inhibitor ratios will result in higher M2/M1 or M3/M1. Pax9 has been shown to play a key role in tooth development. We have previously shown that a G-quadruplex in the first intron of Pax9 can modulate the splicing efficiency. Using a sliding window approach with we analyzed the association of the folding energy (Mfe) of the Pax9 first intron with the relative molar sizes in 42 mammalian species, representing 9 orders. The Mfe of two regions located in the first intron of Pax9 were shown to be significantly associated with the M2/M1 and M3/M1 areas and mesiodistal lengths. The first region is located at the intron beginning and can fold into a stable G4 structure, whereas the second is downstream the G4 and 265 bp from intron start. Across species, the first intron of Pax9 varied in G-quadruplex structural stability. The correlations were further increased when the Mfe of the two sequences were added. Our results indicate that this region has a role in the evolution of the mammalian dental pattern by influencing the relative size of the molars.

Emergent Coordination of the CHKB and CPT1B Genes in Eutherian Mammals: Implications for the Origin of Brown Adipose Tissue.

Mitochondrial fatty acid oxidation (FAO) contributes to the proton motive force that drives ATP synthesis in many mammalian tissues. In eutherian (placental) mammals, brown adipose tissue (BAT) can also dissipate this proton gradient through uncoupling protein 1 (UCP1) to generate heat, but the evolutionary events underlying the emergence of BAT are unknown. An essential step in FAO is the transport of cytoplasmic long chain acyl-coenzyme A (acyl-CoA) into the mitochondrial matrix, which requires the action of carnitine palmitoyltransferase 1B (CPT1B) in striated muscle and BAT. In eutherians, the CPT1B gene is closely linked to the choline kinase beta (CHKB) gene, which is transcribed from the same DNA strand and terminates just upstream of CPT1B. CHKB is a rate-limiting enzyme in the synthesis of phosphatidylcholine (PC), a predominant mitochondrial membrane phospholipid, suggesting that the coordinated expression of CHKB and CPT1B may cooperatively enhance mitochondrial FAO. The present findings show that transcription of the eutherian CHKB and CPT1B genes is linked within a unitary epigenetic domain targeted to the CHKB gene, and that that this regulatory linkage appears to have resulted from an intergenic deletion in eutherians that significantly altered the distribution of CHKB and CPT1B expression. Informed by the timing of this event relative to the emergence of BAT, the phylogeny of CHKB-CPT1B synteny, and the insufficiency of UCP1 to account for eutherian BAT, these data support a mechanism for the emergence of BAT based on the acquisition of a novel capacity for adipocyte FAO in a background of extant UCP1.

The hypercholinergic brain of the Cape golden mole (Chrysochloris asiatica).

Studies detailing the anatomy of the brain of the golden moles are few. A recent study indicated that in the Hottentot golden mole (a member of the Amblysominae clade), there was a broad, atypical, distribution of cholinergic interneurons in the olfactory bulb, cerebral cortex, hippocampus and amygdala. To determine whether this broad distribution of cholinergic neurons is shared by other species of golden mole, we here examine the brain of the Cape golden mole (a member of the Chrysochlorinae clade, representing the second major clade within the family Chrysochloridae). Our analyses indicates the presence of a similar widespread distribution of cholinergic interneurons in the Cape golden mole. Thus, we conclude that these features are derived morphological traits in the brains of golden moles. In addition, we describe the nuclei generally considered to be part of the typical cholinergic system in mammals. Whereas the vast majority of these generally reported cholinergic nuclei were the same as recorded in other Eutherian mammals, it was noted that the cholinergic nuclei involved in oculomotion were substantially reduced in size, or absent in the case of the abducens nucleus. In addition, there was an absence of the cholinergic medial septal nucleus, but the presence of a cholinergic lateral septal nucleus. The laterodorsal and pedunculopontine tegmental nuclei evince regions where the cholinergic neurons are densely packed. These are atypical features of the mammalian cholinergic system, which when combined with the widespread atypical distribution of cholinergic interneurons, reveals a family-specific complement of cholinergic nuclei in the Chrysochloridae.

Molecular simulation of SARS-CoV-2 spike protein binding to pangolin ACE2 or human ACE2 natural variants reveals altered susceptibility to infection.

We constructed complex models of SARS-CoV-2 spike protein binding to pangolin or human ACE2, the receptor for virus transmission, and estimated the binding free energy changes using molecular dynamics simulation. SARS-CoV-2 can bind to both pangolin and human ACE2, but has a significantly lower binding affinity for pangolin ACE2 due to the increased binding free energy (9.5 kcal mol-1). Human ACE2 is among the most polymorphous genes, for which we identified 317 missense single-nucleotide variations (SNVs) from the dbSNP database. Three SNVs, E329G (rs143936283), M82I (rs267606406) and K26R (rs4646116), had a significant reduction in binding free energy, which indicated higher binding affinity than wild-type ACE2 and greater susceptibility to SARS-CoV-2 infection for people with them. Three other SNVs, D355N (rs961360700), E37K (rs146676783) and I21T (rs1244687367), had a significant increase in binding free energy, which indicated lower binding affinity and reduced susceptibility to SARS-CoV-2 infection.

Molecular evolution of autophagy rate-limiting factor LAMP2 in placental mammals.

Autophagy is the cellular process of removal of misfolded or damaged macromolecules and organelles. Experimental studies have demonstrated autophagy as a major mechanism of lifespan extension in long-lived mammals such as bats and mole rat rodents. Moreover, the role of this biological process has been well documented in protection against age-associated diseases and viral infection. However, studies on the molecular adaptive changes of autophagy factors during evolution are scarce. Here, we conducted a bioinformatics study of the molecular evolution of the Lysosomal Associated Membrane Protein 2 (LAMP2), as a rate-limiting factor in the lysosomal degradation stage of autophagy (the communal step of two of autophagy types: macroautophagy and chaperone-mediated). Analyzing LAMP2 across placental mammals, our phylogenetic-based maximum likelihood analyses indicate that the majority of the coding sites undergo purifying selection. However, around 27% of sites display a relaxation of purifying constraints (average ω = 0.42128), among which, 14 particular sites undergo positive selection (ω > 1). These sites are mostly located in the first luminal domain of LAMP2 (N-domain), with a hotspot region in the 135-144 codons interval. Therefore, the N-domain may account for the functional diversity and regulation of LAMP2. In addition, the identified positive selection sites could act as key regulatory sites in the LAMP2 function. On the other hand, testing the rate of evolution in LAMP2 along different clades of placental mammals revealed a relatively relaxed evolution in LAMP2 along megabats' clade. It is not clear yet whether an expedited evolution of LAMP2 in megabats has contributed to their reported up-regulation of autophagy. Finally, our data indicate positive selection sites along the ancestral branch of the clades of rodents, mouse-related rodents, and mole-rats; and suggest the potentially important regulatory role of these sites in LAMP2. Identifying the residues under positive selection, our findings pave the way for future experimental investigations to define how these selective substitutions have functionally affected autophagy.

Evolutionary Analysis of Bile Acid-Conjugating Enzymes Reveals a Complex Duplication and Reciprocal Loss History.

To fulfill their physiological functions, bile acids are conjugated with amino acids. In humans, conjugation is catalyzed by bile acid coenzyme A: amino acid N-acyltransferase (BAAT), an enzyme with a highly conserved catalytic triad in its active site. Interestingly, the conjugated amino acids are highly variable among mammals, with some species conjugating bile acids with both glycine and taurine, whereas others conjugate only taurine. The genetic origin of these bile acid conjugation differences is unknown. Here, we tested whether mutations in BAAT's catalytic triad could explain bile acid conjugation differences. Our comparative analysis of 118 mammals first revealed that the ancestor of placental mammals and marsupials possessed two genes, BAAT and BAATP1, that arose by a tandem duplication. This duplication was followed by numerous gene losses, including BAATP1 in humans. Losses of either BAAT or BAATP1 largely happened in a reciprocal fashion, suggesting that a single conjugating enzyme is generally sufficient for mammals. In intact BAAT and BAATP1 genes, we observed multiple changes in the catalytic triad between Cys and Ser residues. Surprisingly, although mutagenesis experiments with the human enzyme have shown that replacing Cys for Ser greatly diminishes the glycine-conjugating ability, across mammals we found that this residue provides little power in predicting the experimentally measured amino acids that are conjugated with bile acids. This suggests that the mechanism of BAAT's enzymatic function is incompletely understood, despite relying on a classic catalytic triad. More generally, our evolutionary analysis indicates that results of mutagenesis experiments may not easily be extrapolatable to other species.

Mercury, Lead and Cadmium Concentrations in Talpa occidentalis and in Their Digeneans of the Genus Ityogonimus.

INTRODUCTION: Many parasites living in aquatic ecosystems are useful indicators of environmental health. On the other hand, information is scarcer with respect to the use of helminth parasites of vertebrates living in terrestrial ecosystems as monitoring tools for toxic element environmental pollution. The present study evaluates the suitability of the model Talpa occidentalis/Ityogonimus spp. as a bioindicator system for mercury (Hg), lead (Pb) and cadmium (Cd) contamination in agricultural soils from Asturias (Spain). METHODS: Kidney and liver samples collected from T. occidentalis specimens (n = 36) and Ityogonimus spp. samples collected from 14 infected hosts were analyzed by ICP-MS. RESULTS: The highest mean levels of Hg and Pb were found in Ityogonimus individuals (20.9 and 12.4 µg g-1 wet weight, respectively). Considering renal and hepatic concentrations in T. occidentalis, bioaccumulation factors of Ityogonimus for Hg were 83.7 and 58.6, respectively, whereas concerning Pb bioaccumulation factors were 38.2 and 82.9, respectively. No bioaccumulation was detected in Ityogonimus in the case of Cd. CONCLUSIONS: More studies involving digenean parasites of small mammals are needed, especially when biomonitoring environmental toxic element pollution in terrestrial ecosystems. The present results support the above-mentioned model as a suitable biomonitoring system to evaluate environmental Hg and Pb contamination in terrestrial non-urban Iberian habitats. Similar models involving other species (Talpa spp./Ityogonimus spp.) might be used in a much wider geographical range.

A placental mammal-specific microRNA cluster acts as a natural brake for sociability in mice.

Aberrant synaptic function is thought to underlie social deficits in neurodevelopmental disorders such as autism and schizophrenia. Although microRNAs have been shown to regulate synapse development and plasticity, their potential involvement in the control of social behaviour in mammals remains unexplored. Here, we show that deletion of the placental mammal-specific miR379-410 cluster in mice leads to hypersocial behaviour, which is accompanied by increased excitatory synaptic transmission, and exaggerated expression of ionotropic glutamate receptor complexes in the hippocampus. Bioinformatic analyses further allowed us to identify five "hub" microRNAs whose deletion accounts largely for the upregulation of excitatory synaptic genes observed, including Cnih2, Dlgap3, Prr7 and Src. Thus, the miR379-410 cluster acts a natural brake for sociability, and interfering with specific members of this cluster could represent a therapeutic strategy for the treatment of social deficits in neurodevelopmental disorders.

A new research paradigm for bivariate allometry: combining ANOVA and non-linear regression.

A novel statistical routine is presented here for exploring and comparing patterns of allometric variation in two or more groups of subjects. The routine combines elements of the analysis of variance (ANOVA) with non-linear regression to achieve the equivalent of an analysis of covariance (ANCOVA) on curvilinear data. The starting point is a three-parameter power equation to which a categorical variable has been added to identify membership by each subject in a specific group or treatment. The protocol differs from earlier ones in that different assumptions can be made about the form for random error in the full statistical model (i.e. normal and homoscedastic, normal and heteroscedastic, lognormal and heteroscedastic). The general equation and several modifications thereof were used to study allometric variation in field metabolic rates of marsupial and placental mammals. The allometric equations for both marsupials and placentals have an explicit, non-zero intercept, but the allometric exponent is higher in the equation for placentals than in that for marsupials. The approach followed here is extraordinarily versatile, and it has wider application in allometry than standard ANCOVA performed on logarithmic transformations.