Retroposon analysis and recent geological data suggest near-simultaneous divergence of the three superorders of mammals

Common tenrecs (Tenrec ecaudatus) reduce oxygen consumption in hypoxia and in hypercapnia without concordant changes to body temperature or heart rate

Common tenrecs (Tenrec ecaudatus) are fossorial mammals that use burrows during both active and hibernating seasons in Madagascar and its neighboring islands. Prevailing thought was that tenrecs hibernate for 8-9 months individually, but 13 tenrecs were removed from the same sealed burrow 1 m deep from the surface. Such group hibernation in sealed burrows presumably creates a hypoxic and/or hypercapnic environment and suggests that this placental mammal may have an increased tolerance to hypoxia and hypercapnia. Higher tolerances to hypoxia and hypercapnia have been documented for other mammals capable of hibernation and to determine if this is the case for tenrecs, we exposed them to acute hypoxia (4 h of 16 or 7% O2), progressive hypoxia (2 h of 16, 10 and 4% O2), or progressive hypercapnia (2 h of 2, 5 and 10% CO2) at cold (16 °C) or warm (28 °C) ambient temperatures (Ta). Oxygen equilibrium curves were also constructed on the whole blood of tenrecs at 10, 25, and 37 °C to determine if hemoglobin (Hb)-O2 affinity contributes to hypoxia tolerance. In animals held at 16 °C, normoxic and normocapnic levels of oxygen consumption rate (V ˙ O 2 ), body temperature (Tb), and heart rate (HR) were highly variable between individuals. This inter-individual variation was greatly reduced in animals held at 28 °C for oxygen consumption rate and body temperature. Both hypoxia (acute and progressive) and progressive hypercapnia led to decreases inV ˙ O 2 as well as the variation inV ˙ O 2 between animals held at 16 °C. The fall in oxygen consumption rate in 7% O2 independent of changes in body temperature in tenrecs held at 16 °C is unique and not consistent with the typical hypoxic metabolic response seen in other hibernating species that depends on concomitant falls in Tb. In animals held at 28 °C, exposure to O2 levels as low as 4% and CO2 levels as high as 10% had no significant effect onV ˙ O 2 , HR, or Tb, indicative of high tolerance to both hypoxia and hypercapnia. High variation in heart rate remained between individuals in all gas compositions and at all temperatures. Tenrec Hb-O2 affinity was similar to other homeothermic placental mammals and likely does not contribute to the increased hypoxia tolerance. Ultimately, our results suggest changes in Ta dictate physiological responses to hypoxia or hypercapnia in tenrecs, responses more characteristic of reptiles than of most placental mammals. Given that numerous anatomical and physiological characteristics of tenrecs suggest that they may be representative of an ancestral placental mammal, our findings suggest the typical hypoxic metabolic response evolved later in mammalian evolution.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Phylogenomics reveals an almost perfect polytomy among the almost ungulates (Paenungulata)

Phylogenetic studies have resolved most relationships among Eutherian Orders. However, the branching order of elephants (Proboscidea), hyraxes (Hyracoidea), and sea cows (Sirenia) (i.e., the Paenungulata) has remained uncertain since at least 1758, when Linnaeus grouped elephants and manatees into a single Order (Bruta) to the exclusion of hyraxes. Subsequent morphological, molecular, and large-scale phylogenomic datasets have reached conflicting conclusions on the branching order within Paenungulates. We use a phylogenomic dataset of alignments from 13,388 protein-coding genes across 261 Eutherian mammals to infer phylogenetic relationships within Paenungulates. We find that gene trees almost equally support the three alternative resolutions of Paenungulate relationships and that despite strong support for a Proboscidea+Hyracoidea split in the multispecies coalescent (MSC) tree, there is significant evidence for gene tree uncertainty, incomplete lineage sorting, and introgression among Proboscidea, Hyracoidea, and Sirenia. Indeed, only 8-10% of genes have statistically significant phylogenetic signal to reject the hypothesis of a Paenungulate polytomy. These data indicate little support for any resolution for the branching order Proboscidea, Hyracoidea, and Sirenia within Paenungulata and suggest that Paenungulata may be as close to a real, or at least unresolvable, polytomy as possible.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Phylogenomic Coalescent Analyses of Avian Retroelements Infer Zero-Length Branches at the Base of Neoaves, Emergent Support for Controversial Clades, and Ancient Introgressive Hybridization in Afroaves

Retroelement insertions (RIs) are low-homoplasy characters that are ideal data for addressing deep evolutionary radiations, where gene tree reconstruction errors can severely hinder phylogenetic inference with DNA and protein sequence data. Phylogenomic studies of Neoaves, a large clade of birds (>9000 species) that first diversified near the Cretaceous−Paleogene boundary, have yielded an array of robustly supported, contradictory relationships among deep lineages. Here, we reanalyzed a large RI matrix for birds using recently proposed quartet-based coalescent methods that enable inference of large species trees including branch lengths in coalescent units, clade-support, statistical tests for gene flow, and combined analysis with DNA-sequence-based gene trees. Genome-scale coalescent analyses revealed extremely short branches at the base of Neoaves, meager branch support, and limited congruence with previous work at the most challenging nodes. Despite widespread topological conflicts with DNA-sequence-based trees, combined analyses of RIs with thousands of gene trees show emergent support for multiple higher-level clades (Columbea, Passerea, Columbimorphae, Otidimorphae, Phaethoquornithes). RIs express asymmetrical support for deep relationships within the subclade Afroaves that hints at ancient gene flow involving the owl lineage (Strigiformes). Because DNA-sequence data are challenged by gene tree-reconstruction error, analysis of RIs represents one approach for improving gene tree-based methods when divergences are deep, internodes are short, terminal branches are long, and introgressive hybridization further confounds species−tree inference.

SINEs as Credible Signs to Prove Common Ancestry in the Tree of Life: A Brief Review of Pioneering Case Studies in Retroposon Systematics

Currently, the insertions of SINEs (and other retrotransposed elements) are regarded as one of the most reliable synapomorphies in molecular systematics. The methodological mainstream of molecular systematics is the calculation of nucleotide (or amino acid) sequence divergences under a suitable substitution model. In contrast, SINE insertion analysis does not require any complex model because SINE insertions are unidirectional and irreversible. This straightforward methodology was named the “SINE method,” which resolved various taxonomic issues that could not be settled by sequence comparison alone. The SINE method has challenged several traditional hypotheses proposed based on the fossil record and anatomy, prompting constructive discussions in the Evo/Devo era. Here, we review our pioneering SINE studies on salmon, cichlids, cetaceans, Afrotherian mammals, and birds. We emphasize the power of the SINE method in detecting incomplete lineage sorting by tracing the genealogy of specific genomic loci with minimal noise. Finally, in the context of the whole-genome era, we discuss how the SINE method can be applied to further our understanding of the tree of life.

Methodologies for the De novo Discovery of Transposable Element Families

The discovery and characterization of transposable element (TE) families are crucial tasks in the process of genome annotation. Careful curation of TE libraries for each organism is necessary as each has been exposed to a unique and often complex set of TE families. De novo methods have been developed; however, a fully automated and accurate approach to the development of complete libraries remains elusive. In this review, we cover established methods and recent developments in de novo TE analysis. We also present various methodologies used to assess these tools and discuss opportunities for further advancement of the field.

SINE-Based Phylogenomics Reveal Extensive Introgression and Incomplete Lineage Sorting in Myotis

Using presence/absence data from over 10,000 Ves SINE insertions, we reconstructed a phylogeny for 11 Myotis species. With nearly one-third of individual Ves gene trees discordant with the overall species tree, phylogenetic conflict appears to be rampant in this genus. From the observed conflict, we infer that ILS is likely a major contributor to the discordance. Much of the discordance can be attributed to the hypothesized split between the Old World and New World Myotis clades and with the first radiation of Myotis within the New World. Quartet asymmetry tests reveal signs of introgression between Old and New World taxa that may have persisted until approximately 8 MYA. Our introgression tests also revealed evidence of both historic and more recent, perhaps even contemporary, gene flow among Myotis species of the New World. Our findings suggest that hybridization likely played an important role in the evolutionary history of Myotis and may still be happening in areas of sympatry. Despite limitations arising from extreme discordance, our SINE-based phylogeny better resolved deeper relationships (particularly the positioning of M. brandtii) and was able to identify potential introgression pathways among the Myotis species sampled.

A species-level timeline of mammal evolution integrating phylogenomic data

High-throughput sequencing projects generate genome-scale sequence data for species-level phylogenies1-3. However, state-of-the-art Bayesian methods for inferring timetrees are computationally limited to small datasets and cannot exploit the growing number of available genomes4. In the case of mammals, molecular-clock analyses of limited datasets have produced conflicting estimates of clade ages with large uncertainties5,6, and thus the timescale of placental mammal evolution remains contentious7-10. Here we develop a Bayesian molecular-clock dating approach to estimate a timetree of 4,705 mammal species integrating information from 72 mammal genomes. We show that increasingly larger phylogenomic datasets produce diversification time estimates with progressively smaller uncertainties, facilitating precise tests of macroevolutionary hypotheses. For example, we confidently reject an explosive model of placental mammal origination in the Palaeogene8 and show that crown Placentalia originated in the Late Cretaceous with unambiguous ordinal diversification in the Palaeocene/Eocene. Our Bayesian methodology facilitates analysis of complete genomes and thousands of species within an integrated framework, making it possible to address hitherto intractable research questions on species diversifications. This approach can be used to address other contentious cases of animal and plant diversifications that require analysis of species-level phylogenomic datasets.

Theoretical and practical considerations when using retroelement insertions to estimate species trees in the anomaly zone

A potential shortcoming of concatenation methods for species tree estimation is their failure to account for incomplete lineage sorting. Coalescent methods address this problem but make various assumptions that, if violated, can result in worse performance than concatenation. Given the challenges of analyzing DNA sequences with both concatenation and coalescent methods, retroelement insertions (RIs) have emerged as powerful phylogenomic markers for species tree estimation. Here, we show that two recently proposed quartet-based methods, SDPquartets and ASTRAL_BP, are statistically consistent estimators of the unrooted species tree topology under the coalescent when RIs follow a neutral infinite-sites model of mutation and the expected number of new RIs per generation is constant across the species tree. The accuracy of these (and other) methods for inferring species trees from RIs has yet to be assessed on simulated data sets, where the true species tree topology is known. Therefore, we evaluated eight methods given RIs simulated from four model species trees, all of which have short branches and at least three of which are in the anomaly zone. In our simulation study, ASTRAL_BP and SDPquartets always recovered the correct species tree topology when given a sufficiently large number of RIs, as predicted. A distance-based method (ASTRID_BP) and Dollo parsimony also performed well in recovering the species tree topology. In contrast, unordered, polymorphism, and Camin-Sokal parsimony typically fail to recover the correct species tree topology in anomaly zone situations with more than four ingroup taxa. Of the methods studied, only ASTRAL_BP automatically estimates internal branch lengths (in coalescent units) and support values (i.e. local posterior probabilities). We examined the accuracy of branch length estimation, finding that estimated lengths were accurate for short branches but upwardly biased otherwise. This led us to derive the maximum likelihood (branch length) estimate for when RIs are given as input instead of binary gene trees; this corrected formula produced accurate estimates of branch lengths in our simulation study, provided that a sufficiently large number of RIs were given as input. Lastly, we evaluated the impact of data quantity on species tree estimation by repeating the above experiments with input sizes varying from 100 to 100 000 parsimony-informative RIs. We found that, when given just 1 000 parsimony-informative RIs as input, ASTRAL_BP successfully reconstructed major clades (i.e clades separated by branches > 0.3 CUs) with high support and identified rapid radiations (i.e. shorter connected branches), although not their precise branching order. The local posterior probability was effective for controlling false positive branches in these scenarios.

An annotated checklist of the insectivores (Mammalia, Eulipotyphla) of Laos

An annotated checklist of the insectivorous mammals (order Eulipotyphla) of Laos is presented for the first time. The checklist is based on an extensive literature review and an examination of available museum specimens from world collections. Confirmed fauna of eulipotyphlans includes 11 species from 7 genera and 3 families. Additionally, several species which are known from neighbouring countries and could potentially occur in Laos are included in this checklist. According to these data, the Laotian fauna can include at least 23 species in 11 genera.

Out of Africa: A New Afrotheria Lineage Rises From Extinct South American Mammals

The South American native ungulates (SANUs) are usually overlooked in Eutherian phylogenetic studies. In the rare studies where they were included, the diversity of SANUs was underrated, keeping their evolutionary history poorly known. Some authors recognized the SANUs as a monophyletic lineage and formally named it Meridiungulata. Here, we recognized and defined a new supraordinal lineage of Eutheria, the Sudamericungulata, after performing morphological phylogenetic analyses including all lineages of SANUs and Eutheria. The SANUs resulted as non-monophyletic; thus, Meridiungulata is not a natural group; Litopterna and “Didolodontidae” are Panameriungulata and closer to Laurasiatheria than to other “Meridiungulata” (Astrapotheria, Notoungulata, Pyrotheria, and Xenungulata). The other “Meridiungulata” is grouped in the Sudamericungulata, as a new monophyletic lineage of Afrotheria Paenungulata, and shared a common ancestor with Hyracoidea. The divergence between the African and South American lineages is estimated to Early Paleocene, and their interrelationships support the Atlantogea biogeographic model. Shortly afterward, the Sudamericungulata explosively diversified in its four lineages. Confronting the Sudamericungulata evolutionary patterns and the Cenozoic natural events (such as tectonics and climatic and environmental changes, among others) helps to unveil a new chapter in the evolution of Gondwanan Eutheria, as well as the natural history of South America during the Cenozoic.

100-My history of bornavirus infections hidden in vertebrate genomes

Although viruses have threatened our ancestors for millions of years, prehistoric epidemics of viruses are largely unknown. Endogenous bornavirus-like elements (EBLs) are ancient bornavirus sequences derived from the viral messenger RNAs that were reverse transcribed and inserted into animal genomes, most likely by retrotransposons. These elements can be used as molecular fossil records to trace past bornaviral infections. In this study, we systematically identified EBLs in vertebrate genomes and revealed the history of bornavirus infections over nearly 100 My. We confirmed that ancient bornaviral infections have occurred in diverse vertebrate lineages, especially in primate ancestors. Phylogenetic analyses indicated that primate ancestors were infected with various bornaviral lineages during evolution. EBLs in primate genomes formed clades according to their integration ages, suggesting that bornavirus lineages infected with primate ancestors had changed chronologically. However, some bornaviral lineages may have coexisted with primate ancestors and underwent repeated endogenizations for tens of millions of years. Moreover, a bornaviral lineage that coexisted with primate ancestors also endogenized in the genomes of some ancestral bats. The habitats of these bat ancestors have been reported to overlap with the migration route of primate ancestors. These results suggest that long-term virus-host coexistence expanded the geographic distributions of the bornaviral lineage along with primate migration and may have spread their infections to these bat ancestors. Our findings provide insight into the history of bornavirus infections over geological timescales that cannot be deduced from research using extant viruses alone, thus broadening our perspective on virus-host coevolution.

Hovlinc is a recently evolved class of ribozyme found in human lncRNA

Although naturally occurring catalytic RNA molecules-ribozymes-have attracted a great deal of research interest, very few have been identified in humans. Here, we developed a genome-wide approach to discovering self-cleaving ribozymes and identified a naturally occurring ribozyme in humans. The secondary structure and biochemical properties of this ribozyme indicate that it belongs to an unidentified class of small, self-cleaving ribozymes. The sequence of the ribozyme exhibits a clear evolutionary path, from its appearance between ~130 and ~65 million years ago (Ma), to acquiring self-cleavage activity very recently, ~13-10 Ma, in the common ancestors of humans, chimpanzees and gorillas. The ribozyme appears to be functional in vivo and is embedded within a long noncoding RNA belonging to a class of very long intergenic noncoding RNAs. The presence of a catalytic RNA enzyme in lncRNA creates the possibility that these transcripts could function by carrying catalytic RNA domains.

Phylogenomics and the Genetic Architecture of the Placental Mammal Radiation

The genomes of placental mammals are being sequenced at an unprecedented rate. Alignments of hundreds, and one day thousands, of genomes spanning the rich living and extinct diversity of species offer unparalleled power to resolve phylogenetic controversies, identify genomic innovations of adaptation, and dissect the genetic architecture of reproductive isolation. We highlight outstanding questions about the earliest phases of placental mammal diversification and the promise of newer methods, as well as remaining challenges, toward using whole genome data to resolve placental mammal phylogeny. The next phase of mammalian comparative genomics will see the completion and application of finished-quality, gapless genome assemblies from many ordinal lineages and closely related species. Interspecific comparisons between the most hypervariable genomic loci will likely reveal large, but heretofore mostly underappreciated, effects on population divergence, morphological innovation, and the origin of new species.

Phylogenetic variation in cortical layer II immature neuron reservoir of mammals

The adult mammalian brain is mainly composed of mature neurons. A limited amount of stem cell-driven neurogenesis persists in postnatal life and is reduced in large-brained species. Another source of immature neurons in adult brains is cortical layer II. These cortical immature neurons (cINs) retain developmentally undifferentiated states in adulthood, though they are generated before birth. Here, the occurrence, distribution and cellular features of cINs were systematically studied in 12 diverse mammalian species spanning from small-lissencephalic to large-gyrencephalic brains. In spite of well-preserved morphological and molecular features, the distribution of cINs was highly heterogeneous, particularly in neocortex. While virtually absent in rodents, they are present in the entire neocortex of many other species and their linear density in cortical layer II generally increased with brain size. These findings suggest an evolutionary developmental mechanism for plasticity that varies among mammalian species, granting a reservoir of young cells for the cerebral cortex.

To acquire new skills or recover after injuries, the mammalian brain relies on plasticity, the ability for the brain to change its architecture and its connections during the lifetime of an animal.

Creating new nerve cells is one way to achieve plasticity, but this process is rarer in humans than it is in mammals with smaller brains. In particular, it is absent in the human cortex: this region is enlarged in species with large brains, where it carries out complex tasks such as learning and memory. Producing new cells in the cortex would threaten the stability of the structures that retain long-term memories.

Another route to plasticity is to reshape the connections between existing, mature nerve cells. This process takes place in the human brain during childhood and adolescence, as some connections are strengthened and others pruned away.

An alternative mechanism relies on keeping some nerve cells in an immature, ‘adolescent’ state. When needed, these nerve cells emerge from their state of arrested development and ‘grow up’, connecting with the appropriate brain circuits. This mechanism does not involve producing new nerve cells, and so it would be suitable to maintain plasticity in the cortex. Consistent with this idea, in mice some dormant nerve cells are present in a small, primitive part of the cortex.

La Rosa et al. therefore wanted to determine if the location and number of immature cells in the cortex differed between mammals, and if so, whether these differences depended on brain size. The study spanned 12 mammal species, from small-brained species like mice to larger-brained animals including sheep and non-human primates.

Microscopy imaging was used to identify immature nerve cells in brain samples, which revealed that the cortex in larger-brained species contained more adolescent cells than its mouse counterpart. The difference was greatest in a region called the neocortex, which has evolved most recently. This area is most pronounced in primates – especially humans – where it carries out high-level cognitive tasks.

These results identify immature nerve cells as a potential mechanism for plasticity in the cortex. La Rosa et al. hope that the work will inspire searches for similar reservoirs of young cells in humans, which could perhaps lead to new treatments for brain disorders like dementia.

Retrotransposons spread potential cis-regulatory elements during mammary gland evolution

Acquisition of cis-elements is a major driving force for rewiring a gene regulatory network. Several kinds of transposable elements (TEs), mostly retrotransposons that propagate via a copy-and-paste mechanism, are known to possess transcription factor binding motifs and have provided source sequences for enhancers/promoters. However, it remains largely unknown whether retrotransposons have spread the binding sites of master regulators of morphogenesis and accelerated cis-regulatory expansion involved in common mammalian morphological features during evolution. Here, I demonstrate that thousands of binding sites for estrogen receptor α (ERα) and three related pioneer factors (FoxA1, GATA3 and AP2γ) that are essential regulators of mammary gland development arose from a spreading of the binding motifs by retrotransposons. The TE-derived functional elements serve primarily as distal enhancers and are enriched around genes associated with mammary gland morphogenesis. The source TEs occurred via a two-phased expansion consisting of mainly L2/MIR in a eutherian ancestor and endogenous retrovirus 1 (ERV1) in simian primates and murines. Thus the build-up of potential sources for cis-elements by retrotransposons followed by their frequent utilization by the host (co-option/exaptation) may have a general accelerating effect on both establishing and diversifying a gene regulatory network, leading to morphological innovation.

Polyploidy breaks speciation barriers in Australian burrowing frogs Neobatrachus

Polyploidy has played an important role in evolution across the tree of life but it is still unclear how polyploid lineages may persist after their initial formation. While both common and well-studied in plants, polyploidy is rare in animals and generally less understood. The Australian burrowing frog genus Neobatrachus is comprised of six diploid and three polyploid species and offers a powerful animal polyploid model system. We generated exome-capture sequence data from 87 individuals representing all nine species of Neobatrachus to investigate species-level relationships, the origin and inheritance mode of polyploid species, and the population genomic effects of polyploidy on genus-wide demography. We describe rapid speciation of diploid Neobatrachus species and show that the three independently originated polyploid species have tetrasomic or mixed inheritance. We document higher genetic diversity in tetraploids, resulting from widespread gene flow between the tetraploids, asymmetric inter-ploidy gene flow directed from sympatric diploids to tetraploids, and isolation of diploid species from each other. We also constructed models of ecologically suitable areas for each species to investigate the impact of climate on differing ploidy levels. These models suggest substantial change in suitable areas compared to past climate, which correspond to population genomic estimates of demographic histories. We propose that Neobatrachus diploids may be suffering the early genomic impacts of climate-induced habitat loss, while tetraploids appear to be avoiding this fate, possibly due to widespread gene flow. Finally, we demonstrate that Neobatrachus is an attractive model to study the effects of ploidy on the evolution of adaptation in animals.

Transposable elements as genetic accelerators of evolution: contribution to genome size, gene regulatory network rewiring and morphological innovation

In the current era, as a growing number of genome sequence assemblies have been reported in animals, an in-depth analysis of transposable elements (TEs) is one of the most fundamental and essential studies for evolutionary genomics. Although TEs have, in general, been regarded as non-functional junk/selfish DNA, parasitic elements or harmful mutagens, studies have revealed that TEs have had a substantial and sometimes beneficial impact on host genomes in several ways. First, TEs are themselves diverse and thus provide lineage-specific characteristics to the genomes. Second, because TEs constitute a substantial fraction of animal genomes, they are a major contributing factor to evolutionary changes in genome size and composition. Third, host organisms have co-opted many repetitive sequences as genes, cis-regulatory elements and chromatin domain boundaries, which alter gene regulatory networks and in addition are partly involved in morphological evolution, as has been well documented in mammals. Here, I review the impact of TEs on various aspects of the genome, such as genome size and diversity in animals, as well as the evolution of gene networks and genome architecture in mammals. Given that a number of TE families probably remain to be discovered in many non-model organisms, unknown TEs may have contributed to gene networks in a much wider variety of animals than considered previously.

ILS-Aware Analysis of Low-Homoplasy Retroelement Insertions: Inference of Species Trees and Introgression Using Quartets

DNA sequence alignments have provided the majority of data for inferring phylogenetic relationships with both concatenation and coalescent methods. However, DNA sequences are susceptible to extensive homoplasy, especially for deep divergences in the Tree of Life. Retroelement insertions have emerged as a powerful alternative to sequences for deciphering evolutionary relationships because these data are nearly homoplasy-free. In addition, retroelement insertions satisfy the “no intralocus-recombination” assumption of summary coalescent methods because they are singular events and better approximate neutrality relative to DNA loci commonly sampled in phylogenomic studies. Retroelements have traditionally been analyzed with parsimony, distance, and network methods. Here, we analyze retroelement data sets for vertebrate clades (Placentalia, Laurasiatheria, Balaenopteroidea, Palaeognathae) with 2 ILS-aware methods that operate by extracting, weighting, and then assembling unrooted quartets into a species tree. The first approach constructs a species tree from retroelement bipartitions with ASTRAL, and the second method is based on split-decomposition with parsimony. We also develop a Quartet-Asymmetry test to detect hybridization using retroelements. Both ILS-aware methods recovered the same species-tree topology for each data set. The ASTRAL species trees for Laurasiatheria have consecutive short branch lengths in the anomaly zone whereas Palaeognathae is outside of this zone. For the Balaenopteroidea data set, which includes rorquals (Balaenopteridae) and gray whale (Eschrichtiidae), both ILS-aware methods resolved balaeonopterids as paraphyletic. Application of the Quartet-Asymmetry test to this data set detected 19 different quartets of species for which historical introgression may be inferred. Evidence for introgression was not detected in the other data sets.

Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary

Deciphering the timing of the placental mammal radiation is a longstanding problem in evolutionary biology, but consensus on the tempo and mode of placental diversification remains elusive. Nevertheless, an accurate timetree is essential for understanding the role of important events in Earth history (e.g., Cretaceous Terrestrial Revolution, KPg mass extinction) in promoting the taxonomic and ecomorphological diversification of Placentalia. Archibald and Deutschman described three competing models for the diversification of placental mammals, which are the Explosive, Long Fuse, and Short Fuse Models. More recently, the Soft Explosive Model and Trans-KPg Model have emerged as additional hypotheses for the placental radiation. Here, we review molecular and paleontological evidence for each of these five models including the identification of general problems that can negatively impact divergence time estimates. The Long Fuse Model has received more support from relaxed clock studies than any of the other models, but this model is not supported by morphological cladistic studies that position Cretaceous eutherians outside of crown Placentalia. At the same time, morphological cladistics has a poor track record of reconstructing higher-level relationships among the orders of placental mammals including the results of new pseudoextinction analyses that we performed on the largest available morphological data set for mammals (4,541 characters). We also examine the strengths and weaknesses of different timetree methods (node dating, tip dating, and fossilized birth-death dating) that may now be applied to estimate the timing of the placental radiation. While new methods such as tip dating are promising, they also have problems that must be addressed if these methods are to effectively discriminate among competing hypotheses for placental diversification. Finally, we discuss the complexities of timetree estimation when the signal of speciation times is impacted by incomplete lineage sorting (ILS) and hybridization. Not accounting for ILS results in dates that are older than speciation events. Hybridization, in turn, can result in dates than are younger or older than speciation dates. Disregarding this potential variation in "gene" history across the genome can distort phylogenetic branch lengths and divergence estimates when multiple unlinked genomic loci are combined together in a timetree analysis.

The Melanin-Concentrating Hormone (MCH) System: A Tale of Two Peptides

The melanin-concentrating hormone (MCH) system is a robust integrator of exogenous and endogenous information, modulating arousal and energy balance in mammals. Its predominant function in teleosts, however, is to concentrate melanin in the scales, contributing to the adaptive color change observed in several teleost species. These contrasting functions resulted from a gene duplication that occurred after the teleost divergence, which resulted in the generation of two MCH-coding genes in this clade, which acquired distinctive sequences, distribution, and functions, examined in detail here. We also describe the distribution of MCH immunoreactivity and gene expression in a large number of species, in an attempt to identify its core elements. While initially originated as a periventricular peptide, with an intimate relationship with the third ventricle, multiple events of lateral migration occurred during evolution, making the ventrolateral and dorsolateral hypothalamus the predominant sites of MCH in teleosts and mammals, respectively. Substantial differences between species can be identified, likely reflecting differences in habitat and behavior. This observation aligns well with the idea that MCH is a major integrator of internal and external information, ensuring an appropriate response to ensure the organism’s homeostasis. New studies on the MCH system in species that have not yet been investigated will help us understand more precisely how these habitat changes are connected to the hypothalamic neurochemical circuits, paving the way to new intervention strategies that may be used with pharmacological purposes.

The Coagulation and Immune Systems Are Directly Linked through the Activation of Interleukin-1α by Thrombin

Ancient organisms have a combined coagulation and immune system, and although links between inflammation and hemostasis exist in mammals, they are indirect and slower to act. Here we investigated direct links between mammalian immune and coagulation systems by examining cytokine proproteins for potential thrombin protease consensus sites. We found that interleukin (IL)-1α is directly activated by thrombin. Thrombin cleaved pro-IL-1α at a site perfectly conserved across disparate species, indicating functional importance. Surface pro-IL-1α on macrophages and activated platelets was cleaved and activated by thrombin, while tissue factor, a potent thrombin activator, colocalized with pro-IL-1α in the epidermis. Mice bearing a mutation in the IL-1α thrombin cleavage site (R114Q) exhibited defects in efficient wound healing and rapid thrombopoiesis after acute platelet loss. Thrombin-cleaved IL-1α was detected in humans during sepsis, pointing to the relevance of this pathway for normal physiology and the pathogenesis of inflammatory and thrombotic diseases.

•

Mammalian IL-1α contains a highly conserved thrombin consensus site

•

Thrombin cleavage leads to IL-1α activation and shedding from the cell surface

•

Thrombin activates IL-1α after epidermal wounding and after acute platelet loss

•

Thrombin-cleaved IL-1α is also detected in humans during sepsis

Rapid morphological evolution in placental mammals post-dates the origin of the crown group

Resolving the timing and pattern of early placental mammal evolution has been confounded by conflict among divergence date estimates from interpretation of the fossil record and from molecular-clock dating studies. Despite both fossil occurrences and molecular sequences favouring a Cretaceous origin for Placentalia, no unambiguous Cretaceous placental mammal has been discovered. Investigating the differing patterns of evolution in morphological and molecular data reveals a possible explanation for this conflict. Here, we quantified the relationship between morphological and molecular rates of evolution. We show that, independent of divergence dates, morphological rates of evolution were slow relative to molecular evolution during the initial divergence of Placentalia, but substantially increased during the origination of the extant orders. The rapid radiation of placentals into a highly morphologically disparate Cenozoic fauna is thus not associated with the origin of Placentalia, but post-dates superordinal origins. These findings predict that early members of major placental groups may not be easily distinguishable from one another or from stem eutherians on the basis of skeleto-dental morphology. This result supports a Late Cretaceous origin of crown placentals with an ordinal-level adaptive radiation in the early Paleocene, with the high relative rate permitting rapid anatomical change without requiring unreasonably fast molecular evolutionary rates. The lack of definitive Cretaceous placental mammals may be a result of morphological similarity among stem and early crown eutherians, providing an avenue for reconciling the fossil record with molecular divergence estimates for Placentalia.

The Volcano Rabbit in the Phylogenetic Network of Lagomorphs

The order Lagomorpha unifies pikas (Ochotonidae) and the hares plus rabbits (Leporidae). Phylogenetic reconstructions of the species within Leporidae based on traditional morphological or molecular sequence data provide support for conflicting hypotheses. The retroposon presence/absence patterns analyzed in this study revealed strong support for the broadly accepted splitting of lagomorphs into ochotonids and leporids with Pronolagus as the first divergence in the leporid tree. Furthermore, the retroposon presence/absence patterns nested the rare volcano rabbit, Romerolagus diazi, within an unresolved network of deeper leporid relationships and provide the first homoplasy-free image of incomplete lineage sorting and/or ancestral hybridization/introgression in rapidly radiated Leporidae. At the same time, the strongest retroposon presence/absence signal supports the volcano rabbit as a separate branch between the Pronolagus junction and a unified cluster of the remaining leporids.

Global Rate Variation in Bony Vertebrates

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

Loss of RXFP2 and INSL3 genes in Afrotheria shows that testicular descent is the ancestral condition in placental mammals

Descent of testes from a position near the kidneys into the lower abdomen or into the scrotum is an important developmental process that occurs in all placental mammals, with the exception of five afrotherian lineages. Since soft-tissue structures like testes are not preserved in the fossil record and since key parts of the placental mammal phylogeny remain controversial, it has been debated whether testicular descent is the ancestral or derived condition in placental mammals. To resolve this debate, we used genomic data of 71 mammalian species and analyzed the evolution of two key genes (relaxin/insulin-like family peptide receptor 2 [RXFP2] and insulin-like 3 [INSL3]) that induce the development of the gubernaculum, the ligament that is crucial for testicular descent. We show that both RXFP2 and INSL3 are lost or nonfunctional exclusively in four afrotherians (tenrec, cape elephant shrew, cape golden mole, and manatee) that completely lack testicular descent. The presence of remnants of once functional orthologs of both genes in these afrotherian species shows that these gene losses happened after the split from the placental mammal ancestor. These “molecular vestiges” provide strong evidence that testicular descent is the ancestral condition, irrespective of persisting phylogenetic discrepancies. Furthermore, the absence of shared gene-inactivating mutations and our estimates that the loss of RXFP2 happened at different time points strongly suggest that testicular descent was lost independently in Afrotheria. Our results provide a molecular mechanism that explains the loss of testicular descent in afrotherians and, more generally, highlight how molecular vestiges can provide insights into the evolution of soft-tissue characters.

While fossils of whales with legs demonstrate that these species evolved from legged ancestors, the ancestral state of nonfossilizing soft-tissue structures can only be indirectly inferred. This difficulty is also confounded by uncertainties in the phylogenetic relationships between the animals concerned. A prime example is the case of testicular descent, a developmental process that determines the final position of testes, which occurs in most placental mammals but is absent from several afrotherian lineages. Here, we discovered that afrotherians possess remnants of genes known to be required for testicular descent. These “molecular vestiges” show that testicular descent was already present in the placental ancestor and was subsequently lost in Afrotheria. Our study highlights the potential of molecular vestiges in resolving contradictory ancestral states of soft-tissue characters.

A Late Cretaceous mammal from Brazil and the first radioisotopic age for the Bauru Group

In the last three decades, records of tribosphenidan mammals from India, continental Africa, Madagascar and South America have challenged the notion of a strictly Laurasian distribution of the group during the Cretaceous. Here, we describe a lower premolar from the Late Cretaceous Adamantina Formation, São Paulo State, Brazil. It differs from all known fossil mammals, except for a putative eutherian from the same geologic unity and Deccanolestes hislopi, from the Maastrichtian of India. The incompleteness of the material precludes narrowing down its taxonomic attribution further than Tribosphenida, but it is larger than most coeval mammals and shows a thin layer of parallel crystallite enamel. The new taxon helps filling two major gaps in the fossil record: the paucity of Mesozoic mammals in more northern parts of South America and of tribosphenidans in the Cretaceous of that continent. In addition, high-precision U-Pb geochronology provided a post-Turonian maximal age (≤87.8 Ma) for the type stratum, which is overlain by the dinosaur-bearing Marília Formation, constraining the age of the Adamantina Formation at the site to late Coniacian-late Maastrichtian. This represents the first radioisotopic age for the Bauru Group, a key stratigraphic unit for the study of Cretaceous tetrapods in Gondwana.

A revised phylogeny of macropathine cave crickets (Orthoptera: Rhaphidophoridae) uncovers a paraphyletic Australian fauna

Australian cave crickets are members of the subfamily Macropathinae (Orthoptera: Rhaphidophoridae). The subfamily is thought to have originated prior to the tectonic separation of the supercontinent Gondwana based on distributions of extant lineages and molecular phylogenetic evidence, although the Australian fauna have been underrepresented in previous studies. The current study augments existing multigene data (using 12S, 16S, and 28S rRNA genes) to investigate the placement of the Australian representatives within the Macropathinae and to assess divergence dates of select clades. Results suggest that the endemic Tasmanian genus Parvotettix is the sister lineage to the remaining members of the subfamily, an outcome that presents a paraphyletic Australian fauna in contrast to previous studies. All other Australian taxa represented in this study (Micropathus and Novotettix) emerged as a sister group to the New Zealand and South American macropathine lineages. Estimation of phylogenetic divergence ages among the aforementioned clades were calibrated using two methods, in absence of suitable fossil records: (i) tectonic events depicting the fragmentation of Gondwanan landmasses that invoke vicariant scenarios of present day geographic distributions; and (ii) molecular evolutionary rates. Geological calibrations place the median age of the most recent common ancestor of extant macropathines at ∼125 to ∼165 Ma, whereas analyses derived from molecular substitution rates suggest a considerably younger origin of ∼32 Ma. This phylogenetic study represents the most rigorous taxonomic sampling of the Australian cave cricket fauna to date and stresses the influence of lineage representation on biogeographic inference.

Appropriate fossil calibrations and tree constraints uphold the Mesozoic divergence of solenodons from other extant mammals

The mammalian order Eulipotyphla includes four extant families of insectivorans: Solenodontidae (solenodons); Talpidae (moles); Soricidae (shrews); and Erinaceidae (hedgehogs). Of these, Solenodontidae includes only two extant species, which are endemic to the largest islands of the Greater Antilles: Cuba and Hispaniola. Most molecular studies suggest that eulipotyphlan families diverged from each other across several million years, with the basal split between Solenodontidae and other families occurring in the Late Cretaceous. By contrast, Sato et al. (2016) suggest that eulipotyphlan families diverged from each other in a polytomy ∼58.6 million years ago (Mya). This more recent divergence estimate for Solenodontidae versus other extant eulipotyphlans suggests that solenodons must have arrived in the Greater Antilles via overwater dispersal rather than vicariance. Here, we show that the young timetree estimates for eulipotyphlan families and the polytomy are due to an inverted ingroup-outgroup arrangement of the tree, the result of using Tracer rather than TreeAnnotator to compile interfamilial divergence times, and of not enforcing the monophly of well-established clades such as Laurasiatheria and Eulipotyphla. Finally, Sato et al.'s (2016) timetree includes several zombie lineages where estimated divergence times are much younger than minimum ages that are implied by the fossil record. We reanalyzed Sato et al.'s (2016) original data with enforced monophyly for well-established clades and updated fossil calibrations that eliminate the inference of zombie lineages. Our resulting timetrees, which were compiled with TreeAnnotator rather than Tracer, produce dates that are in good agreement with other recent studies and place the basal split between Solenodontidae and other eulipotyphlans in the Late Cretaceous.

Evolution of the MC5R gene in placental mammals with evidence for its inactivation in multiple lineages that lack sebaceous glands

MC5R is one of five melanocortin receptor genes found in placental mammals. MC5R plays an important role in energy homeostasis and is also expressed in the terminal differentiation of sebaceous glands. Among placental mammals there are multiple lineages that either lack or have degenerative sebaceous glands including Cetacea (whales, dolphins, and porpoises), Hippopotamidae (hippopotamuses), Sirenia (manatees and dugongs), Proboscidea (elephants), Rhinocerotidae (rhinos), and Heterocephalus glaber (naked mole rat). Given the loss or diminution of sebaceous glands in these taxa, we procured MC5R sequences from publicly available genomes and transcriptomes, supplemented by a newly generated sequence for Choeropsis liberiensis (pygmy hippopotamus), to determine if this gene remains intact or is inactivated in association with loss/reduction of sebaceous glands. Our data set includes complete MC5R sequences for 114 placental mammal species including two individuals of Mammuthus primigenius (woolly mammoth) from Oimyakon and Wrangel Island. Complete loss or inactivation of the MC5R gene occurs in multiple placental lineages that have lost sebaceous glands (Cetacea, West Indian manatee, African elephant, white rhinoceros) or are characterized by unusual skin (pangolins, aardvarks). Both M. primigenius individuals share inactivating mutations with the African elephant even though sebaceous glands have been reported in the former. MC5R remains intact in hippopotamuses and the naked mole rat, although slightly elevated dN/dS ratios in these lineages allow for the possibility that the accumulation of inactivating mutations in MC5R may lag behind the relaxation of purifying selection. For Cetacea and Hippopotamidae, the absence of shared inactivating mutations in two different skin genes (MC5R, PSORS1C2) is consistent with the hypothesis that semi-aquatic lifestyles were acquired independently in these clades following divergence from a common ancestor.

Resolving kangaroo phylogeny and overcoming retrotransposon ascertainment bias

Reconstructing phylogeny from retrotransposon insertions is often limited by access to only a single reference genome, whereby support for clades that do not include the reference taxon cannot be directly observed. Here we have developed a new statistical framework that accounts for this ascertainment bias, allowing us to employ phylogenetically powerful retrotransposon markers to explore the radiation of the largest living marsupials, the kangaroos and wallabies of the genera Macropus and Wallabia. An exhaustive in silico screening of the tammar wallaby (Macropus eugenii) reference genome followed by experimental screening revealed 29 phylogenetically informative retrotransposon markers belonging to a family of endogenous retroviruses. We identified robust support for the enigmatic swamp wallaby (Wallabia bicolor) falling within a paraphyletic genus, Macropus. Our statistical approach provides a means to test for incomplete lineage sorting and introgression/hybridization in the presence of the ascertainment bias. Using retrotransposons as “molecular fossils”, we reveal one of the most complex patterns of hemiplasy yet identified, during the rapid diversification of kangaroos and wallabies. Ancestral state reconstruction incorporating the new retrotransposon phylogenetic information reveals multiple independent ecological shifts among kangaroos into more open habitats, coinciding with the Pliocene onset of increased aridification in Australia from ~3.6 million years ago.

Lineage-independent retrotransposition of UTP14 associated with male fertility has occurred multiple times throughout mammalian evolution

In mammals, gamete production is essential for reproductive success. This is particularly true for males where large quantities of sperm are produced to fertilize a limited number of eggs released by the female. Because of this, new genes associated with increased spermatogenic efficiency have been accumulating throughout the evolution of therian mammals. Many of these new genes are testis-specific retrotransposed copies of housekeeping genes located on the X chromosome. Of particular interest are retrotransposed copies of UTP14 that are present in many distantly related eutherian mammals. Analysis of genomic data available in ENSEMBL indicates that these UTP14 retrogenes have arisen independently in the various eutherian clades. They represent an interesting aspect of evolution whereby new homologues of UTP14 have become independently fixed in multiple mammalian lineages due to the reproductive advantage that may be conferred to males. Surprisingly, these genes may also be lost, even after being present within a lineage for millions of years. This phenomenon may potentially be used to delineate evolutionary trees in closely related groups of mammals, particularly in the case of South American primates. Studying these retrogenes will yield new insights into the evolutionary history of male gamete production and the phylogeny of eutherian mammals.

Phylogenetic Conflict in Bears Identified by Automated Discovery of Transposable Element Insertions in Low-Coverage Genomes

Phylogenetic reconstruction from transposable elements (TEs) offers an additional perspective to study evolutionary processes. However, detecting phylogenetically informative TE insertions requires tedious experimental work, limiting the power of phylogenetic inference. Here, we analyzed the genomes of seven bear species using high-throughput sequencing data to detect thousands of TE insertions. The newly developed pipeline for TE detection called TeddyPi (TE detection and discovery for Phylogenetic Inference) identified 150,513 high-quality TE insertions in the genomes of ursine and tremarctine bears. By integrating different TE insertion callers and using a stringent filtering approach, the TeddyPi pipeline produced highly reliable TE insertion calls, which were confirmed by extensive in vitro validation experiments. Analysis of single nucleotide substitutions in the flanking regions of the TEs shows that these substitutions correlate with the phylogenetic signal from the TE insertions. Our phylogenomic analyses show that TEs are a major driver of genomic variation in bears and enabled phylogenetic reconstruction of a well-resolved species tree, despite strong signals for incomplete lineage sorting and introgression. The analyses show that the Asiatic black, sun, and sloth bear form a monophyletic clade, in which phylogenetic incongruence originates from incomplete lineage sorting. TeddyPi is open source and can be adapted to various TE and structural variation callers. The pipeline makes it possible to confidently extract thousands of TE insertions even from low-coverage genomes (∼10×) of nonmodel organisms. This opens new possibilities for biologists to study phylogenies and evolutionary processes as well as rates and patterns of (retro-)transposition and structural variation.

Evolution of brain functions in mammals and LTR retrotransposon-derived genes

In the human genome, there are approximately 30 LTR retrotransposon-derived genes, such as the sushi-ichi retrotransposon homologues (SIRH) and the paraneoplastic Ma antigen (PNMA) family genes. They are derivatives from the original LTR retrotransposons and each gene seems to have its own unique function. PEG10/SIRH1 as well as PEG11/RTL1/SIRH2 and SIRH7/LDOC1 play essential roles in placenta formation, maintenance of fetal capillaries and the differentiation/maturation of a variety of placental cells, respectively. All of this evidence provides strong support for their contribution to the evolution of viviparity in mammals via their eutherian-specific functions. SIRH11/ZCCHC16 is an X-linked gene that encodes a CCHC type of zinc-finger protein that exhibits high sequence identity to the LTR retrotransposon Gag protein and its deletion causes abnormal behavior related to cognition, including attention, impulsivity and working memory, possibly via the locus coeruleus noradrenaergic system in mice. Therefore, we have suggested that the acquisition of SIRH11/ZCCHC16 was involved in eutherian brain evolution. Interestingly, SIRH11/ZCCHC16 displays lineage-specific structural and putative species-specific functional variations in eutherians, suggesting that it contributed to the diversification of eutherians via increasing evolutionary fitness by these changes.

Speciation network in Laurasiatheria: retrophylogenomic signals

Rapid species radiation due to adaptive changes or occupation of new ecospaces challenges our understanding of ancestral speciation and the relationships of modern species. At the molecular level, rapid radiation with successive speciations over short time periods-too short to fix polymorphic alleles-is described as incomplete lineage sorting. Incomplete lineage sorting leads to random fixation of genetic markers and hence, random signals of relationships in phylogenetic reconstructions. The situation is further complicated when you consider that the genome is a mosaic of ancestral and modern incompletely sorted sequence blocks that leads to reconstructed affiliations to one or the other relative, depending on the fixation of their shared ancestral polymorphic alleles. The laurasiatherian relationships among Chiroptera, Perissodactyla, Cetartiodactyla, and Carnivora present a prime example for such enigmatic affiliations. We performed whole-genome screenings for phylogenetically diagnostic retrotransposon insertions involving the representatives bat (Chiroptera), horse (Perissodactyla), cow (Cetartiodactyla), and dog (Carnivora), and extracted among 162,000 preselected cases 102 virtually homoplasy-free, phylogenetically informative retroelements to draw a complete picture of the highly complex evolutionary relations within Laurasiatheria. All possible evolutionary scenarios received considerable retrotransposon support, leaving us with a network of affiliations. However, the Cetartiodactyla-Carnivora relationship as well as the basal position of Chiroptera and an ancestral laurasiatherian hybridization process did exhibit some very clear, distinct signals. The significant accordance of retrotransposon presence/absence patterns and flanking nucleotide changes suggest an important influence of mosaic genome structures in the reconstruction of species histories.

Cretaceous origin of giant rhinoceros beetles (Dynastini; Coleoptera) and correlation of their evolution with the Pangean breakup

The giant rhinoceros beetles (Dynastini, Scarabaeidae, Coleoptera) are distributed in tropical and temperate regions in Asia, America and Africa. Recent molecular phylogenetic studies have revealed that the giant rhinoceros beetles can be divided into three clades representing Asia, America and Africa. Although a correlation between their evolution and the continental drift during the Pangean breakup was suggested, there is no accurate divergence time estimation among the three clades based on molecular data. Moreover, there is a long chronological gap between the timing of the Pangean breakup (Cretaceous: 110-148 Ma) and the emergence of the oldest fossil record (Oligocene: 33 Ma). In this study, we estimated their divergence times based on molecular data, using several combinations of fossil calibration sets, and obtained robust estimates. The inter-continental divergence events among the clades were estimated to have occurred about 99 Ma (Asian clade and others) and 78 Ma (American clade and African clade), both of which are after the Pangean breakup. These estimates suggest their inter-continental divergences occurred by overseas sweepstakes dispersal, rather than by vicariances of the population caused by the Pangean breakup.

Comparative genomics of human stem cell factor (SCF)

Stem cell factor (SCF) is a critical protein with key roles in the cell such as hematopoiesis, gametogenesis and melanogenesis. In the present study a comparative analysis on nucleotide sequences of SCF was performed in Humanoids using bioinformatics tools including NCBI-BLAST, MEGA6, and JBrowse. Our analysis of nucleotide sequences to find closely evolved organisms with high similarity by NCBI-BLAST tools and MEGA6 showed that human and Chimpanzee (Pan troglodytes) were placed into the same cluster. By using JBrowse, we found that SCF in Neanderthal had a single copy number similar to modern human and partly conserved nucleotide sequences. Together, the results approved the gene flow and genetics similarity of SCF among human and P. troglodytes. This may suggest that during evolution, SCF gene transferred partly intact either on the basis of sequence or function from the same ancestors to P. troglodytes, the ancient human like Neanderthal, and then to the modern human.

Resolving the relationships of Paleocene placental mammals

The 'Age of Mammals' began in the Paleocene epoch, the 10 million year interval immediately following the Cretaceous-Palaeogene mass extinction. The apparently rapid shift in mammalian ecomorphs from small, largely insectivorous forms to many small-to-large-bodied, diverse taxa has driven a hypothesis that the end-Cretaceous heralded an adaptive radiation in placental mammal evolution. However, the affinities of most Paleocene mammals have remained unresolved, despite significant advances in understanding the relationships of the extant orders, hindering efforts to reconstruct robustly the origin and early evolution of placental mammals. Here we present the largest cladistic analysis of Paleocene placentals to date, from a data matrix including 177 taxa (130 of which are Palaeogene) and 680 morphological characters. We improve the resolution of the relationships of several enigmatic Paleocene clades, including families of 'condylarths'. Protungulatum is resolved as a stem eutherian, meaning that no crown-placental mammal unambiguously pre-dates the Cretaceous-Palaeogene boundary. Our results support an Atlantogenata-Boreoeutheria split at the root of crown Placentalia, the presence of phenacodontids as closest relatives of Perissodactyla, the validity of Euungulata, and the placement of Arctocyonidae close to Carnivora. Periptychidae and Pantodonta are resolved as sister taxa, Leptictida and Cimolestidae are found to be stem eutherians, and Hyopsodontidae is highly polyphyletic. The inclusion of Paleocene taxa in a placental phylogeny alters interpretations of relationships and key events in mammalian evolutionary history. Paleocene mammals are an essential source of data for understanding fully the biotic dynamics associated with the end-Cretaceous mass extinction. The relationships presented here mark a critical first step towards accurate reconstruction of this important interval in the evolution of the modern fauna.

Waking the undead: Implications of a soft explosive model for the timing of placental mammal diversification

The explosive, long fuse, and short fuse models represent competing hypotheses for the timing of placental mammal diversification. Support for the explosive model, which posits both interordinal and intraordinal diversification after the KPg mass extinction, derives from morphological cladistic studies that place Cretaceous eutherians outside of crown Placentalia. By contrast, most molecular studies favor the long fuse model wherein interordinal cladogenesis occurred in the Cretaceous followed by intraordinal cladogenesis after the KPg boundary. Phillips (2016) proposed a soft explosive model that allows for the emergence of a few lineages (Xenarthra, Afrotheria, Euarchontoglires, Laurasiatheria) in the Cretaceous, but otherwise agrees with the explosive model in positing the majority of interordinal diversification after the KPg mass extinction. Phillips (2016) argues that rate transference errors associated with large body size and long lifespan have inflated previous estimates of interordinal divergence times, and further suggests that most interordinal divergences are positioned after the KPg boundary when rate transference errors are avoided through the elimination of calibrations in large-bodied and/or long lifespan clades. Here, we show that rate transference errors can also occur in the opposite direction and drag forward estimated divergence dates when calibrations in large-bodied/long lifespan clades are omitted. This dragging forward effect results in the occurrence of more than half a billion years of 'zombie lineages' on Phillips' preferred timetree. By contrast with ghost lineages, which are a logical byproduct of an incomplete fossil record, zombie lineages occur when estimated divergence dates are younger than the minimum age of the oldest crown fossils. We also present the results of new timetree analyses that address the rate transference problem highlighted by Phillips (2016) by deleting taxa that exceed thresholds for body size and lifespan. These analyses recover all interordinal divergence times in the Cretaceous and are consistent with the long fuse model of placental diversification. Finally, we outline potential problems with morphological cladistic analyses of higher-level relationships among placental mammals that may account for the perceived discrepancies between molecular and paleontological estimates of placental divergence times.

Mammal madness: is the mammal tree of life not yet resolved?

Most molecular phylogenetic studies place all placental mammals into four superordinal groups, Laurasiatheria (e.g. dogs, bats, whales), Euarchontoglires (e.g. humans, rodents, colugos), Xenarthra (e.g. armadillos, anteaters) and Afrotheria (e.g. elephants, sea cows, tenrecs), and estimate that these clades last shared a common ancestor 90-110 million years ago. This phylogeny has provided a framework for numerous functional and comparative studies. Despite the high level of congruence among most molecular studies, questions still remain regarding the position and divergence time of the root of placental mammals, and certain 'hard nodes' such as the Laurasiatheria polytomy and Paenungulata that seem impossible to resolve. Here, we explore recent consensus and conflict among mammalian phylogenetic studies and explore the reasons for the remaining conflicts. The question of whether the mammal tree of life is or can be ever resolved is also addressed.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism

The notion of species as reproductively isolated units related through a bifurcating tree implies that gene trees should generally agree with the species tree and that sister taxa should not share polymorphisms unless they diverged recently and should be equally closely related to outgroups. It is now possible to evaluate this model systematically. We sequenced multiple individuals from 27 described taxa representing the entire Arabidopsis genus. Cluster analysis identified seven groups, corresponding to described species that capture the structure of the genus. However, at the level of gene trees, only the separation of Arabidopsis thaliana from the remaining species was universally supported, and, overall, the amount of shared polymorphism demonstrated that reproductive isolation was considerably more recent than the estimated divergence times. We uncovered multiple cases of past gene flow that contradict a bifurcating species tree. Finally, we showed that the pattern of divergence differs between gene ontologies, suggesting a role for selection.

An LTR Retrotransposon-Derived Gene Displays Lineage-Specific Structural and Putative Species-Specific Functional Variations in Eutherians

Amongst the 11 eutherian-specific genes acquired from a sushi-ichi retrotransposon is the CCHC type zinc-finger protein-encoding gene SIRH11/ZCCHC16. Its contribution to eutherian brain evolution is implied because of its involvement in cognitive function in mice, possibly via the noradrenergic system. Although, the possibility that Sirh11/Zcchc16 functions as a non-coding RNA still remains, dN/dS ratios in pairwise comparisons between its orthologs have provided supportive evidence that it acts as a protein. It became a pseudogene in armadillos (Cingulata) and sloths (Pilosa), the only two extant orders of xenarthra, which prompted us to examine the lineage-specific variations of SIRH11/ZCCHC16 in eutherians. We examined the predicted SIRH11/ZCCHC16 open reading frame (ORF) in 95 eutherian species based on the genomic DNA information in GenBank. A large variation in the SIRH11/ZCCHC16 ORF was detected in several lineages. These include a lack of a CCHC RNA-binding domain in its C-terminus, observed in gibbons (Hylobatidae: Primates) and megabats (Megachiroptera: Chiroptera). A lack of the N-terminal half, on the other hand, was observed in New World monkeys (Platyrrhini: Primates) and species belonging to New World and African Hystricognaths (Caviomorpha and Bathyergidae: Rodents) along with Cetacea and Ruminantia (Cetartiodactyla). Among the hominoids, interestingly, three out of four genera of gibbons have lost normal SIRH11/ZCCHC16 function by deletion or the lack of the CCHC RNA-binding domain. Our extensive dN/dS analysis suggests that such truncated SIRH11/ZCCHC16 ORFs are functionally diversified even within lineages. Combined, our results show that SIRH11/ZCCHC16 may contribute to the diversification of eutherians by lineage-specific structural changes after its domestication in the common eutherian ancestor, followed by putative species-specific functional changes that enhanced fitness and occurred as a consequence of complex natural selection events.

Alu SINE analyses of 3,000-year-old human skeletal remains: a pilot study

Background

As Short Interspersed Elements (SINEs), human-specific Alu elements can be used for population genetic studies. Very recent inserts are polymorphic within and between human populations. In a sample of 30 elements originating from three different Alu subfamilies, we investigated whether they are preserved in prehistorical skeletal human remains from the Bronze Age Lichtenstein cave in Lower Saxony, Germany. In the present study, we examined a prehistoric triad of father, mother and daughter.

Results

For 26 of the 30 Alu loci investigated, definite results were obtained. We were able to demonstrate that presence/absence analyses of Alu elements can be conducted on individuals who lived 3,000 years ago. The preservation of the ancient DNA (aDNA) is good enough in two out of three ancient individuals to routinely allow the amplification of 500 bp fragments. The third individual revealed less well-preserved DNA, which results in allelic dropout or complete amplification failures. We here present an alternative molecular approach to deal with these degradation phenomena by using internal Alu subfamily specific primers producing short fragments of approximately 150 bp.

Conclusions

Our data clearly show the possibility of presence/absence analyses of Alu elements in individuals from the Lichtenstein cave. Thus, we demonstrate that our method is reliably applicable for aDNA samples with good or moderate DNA preservation. This method will be very useful for further investigations with more Alu loci and larger datasets. Human population genetic studies and other large-scale investigations would provide insight into Alu SINE-based microevolutionary processes in humans during the last few thousand years and help us comprehend the evolutionary dynamics of our genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0063-y) contains supplementary material, which is available to authorized users.

Incomplete Lineage Sorting and Hybridization Statistics for Large-Scale Retroposon Insertion Data

Ancient retroposon insertions can be used as virtually homoplasy-free markers to reconstruct the phylogenetic history of species. Inherited, orthologous insertions in related species offer reliable signals of a common origin of the given species. One prerequisite for such a phylogenetically informative insertion is that the inserted element was fixed in the ancestral population before speciation; if not, polymorphically inserted elements may lead to random distributions of presence/absence states during speciation and possibly to apparently conflicting reconstructions of their ancestry. Fortunately, such misleading fixed cases are relatively rare but nevertheless, need to be considered. Here, we present novel, comprehensive statistical models applicable for (1) analyzing any pattern of rare genomic changes, (2) testing and differentiating conflicting phylogenetic reconstructions based on rare genomic changes caused by incomplete lineage sorting or/and ancestral hybridization, and (3) differentiating between search strategies involving genome information from one or several lineages. When the new statistics are applied, in non-conflicting cases a minimum of three elements present in both of two species and absent in a third group are considered significant support (p<0.05) for the branching of the third from the other two, if all three of the given species are screened equally for genome or experimental data. Five elements are necessary for significant support (p<0.05) if a diagnostic locus derived from only one of three species is screened, and no conflicting markers are detected. Most potentially conflicting patterns can be evaluated for their significance and ancestral hybridization can be distinguished from incomplete lineage sorting by considering symmetric or asymmetric distribution of rare genomic changes among possible tree configurations. Additionally, we provide an R-application to make the new KKSC insertion significance test available for the scientific community at http://retrogenomics.uni-muenster.de:3838/KKSC_significance_test/.

The presence/absence patterns of transposed elements, so called jumping genes, provide invaluable information about evolution. Unfortunately, there is still no clear all-encompassing analysis of the statistical significance of insertion patterns, and the single existing model of insertion data is no longer sufficient for the emerging genomic era. Here, we have provided a comprehensive statistical framework for testing the significance of support for phylogenetic hypotheses derived from genome-level presence/absence data such as retroposon insertions and for evaluating such data for different evolutionary scenarios, including polytomy, incomplete lineage sorting, and ancestral hybridization. This statistical framework is especially important for high-throughput applications of current and upcoming genome projects due to its treatment of unlimited numbers of testable markers, and is embedded in a user-friendly R-application available to the scientific community online. Finally, a reliable, adaptable calculation for the significance of support for phylogenetic trees derived from retroposon presence/absence data is now available.