Heterogeneous models place the root of the placental mammal phylogeny

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Current genome sequencing initiatives across a wide range of life forms offer significant potential to enhance our understanding of evolutionary relationships and support transformative biological and medical applications. Species trees play a central role in many of these applications; however, despite the widespread availability of genome assemblies, accurate inference of species trees remains challenging due to the limited automation, substantial domain expertise, and computational resources required by conventional methods. To address this limitation, we present ROADIES, a fully automated pipeline to infer species trees starting from raw genome assemblies. In contrast to the prominent approach, ROADIES incorporates a unique strategy of randomly sampling segments of the input genomes to generate gene trees. This eliminates the need for predefining a set of loci, limiting the analyses to a fixed number of genes, and performing the cumbersome gene annotation and/or whole genome alignment steps. ROADIES also eliminates the need to infer orthology by leveraging existing discordance-aware methods that allow multicopy genes. Using the genomic datasets from large-scale sequencing efforts across four diverse life forms (placental mammals, pomace flies, birds, and budding yeasts), we show that ROADIES infers species trees that are comparable in quality to the state-of-the-art studies but in a fraction of the time and effort, including on challenging datasets with rampant gene tree discordance and complex polyploidy. With its speed, accuracy, and automation, ROADIES has the potential to vastly simplify species tree inference, making it accessible to a broader range of scientists and applications.

Conservation and Variability in Mitochondrial Genomes of Perlodidae: Insights from Comparative Mitogenomics

The mitochondrial genomes of three stoneflies, e.g., Tibetisoperla wangluyui Huo and Du, 2021, Perlodinella kozlovi Klapálek, 1912 and Perlodinella epiproctalis (Zwick, 1997), were sequenced in this study, with lengths 16,043 bp, 16,024 bp, and 16,071 bp, respectively. Each mitogenome contained 37 genes including 22 tRNAs, two ribosomal RNAs, 13 protein-coding genes (PCGs), and a noncoding control region (CR). In general, standard ATN start and TAN termination codons were evident in the PCGs. Meanwhile, in this paper, three newly published mitochondrial genomes and 11 existing mitochondrial genomes of the Perlodidae from NCBI were analyzed. Among the 13 PCGs in the mitochondrial genome of Perlodidae, the lengths of atp6, atp8, cox2, cox3, cytb, nad1, nad2, nad3, and nad4 are exactly the same, and the length of cox1 is 1536-1569 bp. The length of nad4L is 297, but the length of Arcynopteryx dichroa is 300. The length of nad5 ranges from 1732 bp to 1752 bp, while that of nad6 ranges from 525 bp to 534 bp. The length of rrnL is between 1292 and 391 bp, and the length of rrnS is between 793 and 869 bp. In addition, we found that atp8 in Isoperlinae started with GTG as a start codon but in Perlodinae, it started with ATG. Despite these advances, mitochondrial genome data from the Perlodidae are still needed.

CAT-Posterior Mean Site Frequencies Improves Phylogenetic Modeling Under Maximum Likelihood and Resolves Tardigrada as the Sister of Arthropoda Plus Onychophora

Tardigrada, the water bears, are microscopic animals with walking appendages that are members of Ecdysozoa, the clade of molting animals that also includes Nematoda (round worms), Nematomorpha (horsehair worms), Priapulida (penis worms), Kinorhyncha (mud dragons), Loricifera (loricated animals), Arthropoda (insects, spiders, centipedes, crustaceans, and their allies), and Onychophora (velvet worms). The phylogenetic relationships within Ecdysozoa are still unclear, with analyses of molecular and morphological data yielding incongruent results. Accounting for across-site compositional heterogeneity using mixture models that partition sites in frequency categories, CATegories (CAT)-based models, has been shown to improve fit in Bayesian analyses. However, CAT-based models such as CAT-Poisson or CAT-GTR (where CAT is combined with a General Time Reversible matrix to account for replacement rate heterogeneity) have proven difficult to implement in maximum likelihood. Here, we use CAT-posterior mean site frequencies (CAT-PMSF), a new method to export dataset-specific mixture models (CAT-Poisson and CAT-GTR) parameterized using Bayesian methods to maximum likelihood software. We developed new maximum likelihood-based model adequacy tests using parametric bootstrap and show that CAT-PMSF describes across-site compositional heterogeneity better than other across-site compositionally heterogeneous models currently implemented in maximum likelihood software. CAT-PMSF suggests that tardigrades are members of Panarthropoda, a lineage also including Arthropoda and Onychophora. Within Panarthropoda, our results favor Tardigrada as sister to Onychophora plus Arthropoda (the Lobopodia hypothesis). Our results illustrate the power of CAT-PMSF to model across-site compositionally heterogeneous datasets in the maximum likelihood framework and clarify the relationships between the Tardigrada and the Ecdysozoa.

Saltational Episodes of Reticulate Evolution in the Drosophila saltans Species Group

Phylogenomics reveals reticulate evolution to be widespread across taxa, but whether reticulation is due to low statistical power or it is a true evolutionary pattern remains a field of study. Here, we investigate the phylogeny and quantify reticulation in the Drosophila saltans species group, a Neotropical clade of the subgenus Sophophora comprising 23 species whose relationships have long been problematic. Phylogenetic analyses revealed conflicting topologies between the X chromosome, autosomes and the mitochondria. We extended the ABBA-BABA test of asymmetry in phylogenetic discordance to cases where no "true" species tree could be inferred, and applied our new test (called 2A2B) to whole genome data and to individual loci. We used four strategies, two based on our new assemblies using either conserved genes or ≥50 kb-long syntenic blocks with conserved collinearity across Neotropical Sophophora, and two consisted of windows from pseudo-reference genomes aligned to either an ingroup or outgroup species. Evidence for reticulation varied among the strategies, being lowest in the synteny-based approach, where it did not exceed ∼7% of the blocks in the most conflicting species quartets. High incidences of reticulation were restricted to three nodes on the tree that coincided with major paleogeographical events in South America. Our results identify possible technical biases in quantifying reticulate evolution and indicate that episodic rapid radiations have played a major role in the evolution of a largely understudied Neotropical clade.

Whole-genome-based phylogenetic analyses provide new insights into the evolution of springtails (Hexapoda: Collembola)

Springtails (Collembola) stand as one of the most abundant, widespread, and ancient terrestrial arthropods on earth. However, their evolutionary history and deep phylogenetic relationships remain elusive. In this study, we employed phylogenomic approaches to elucidate the basal relationships among Collembola. We sampled whole-genome data representing all major collembolan lineages in proportion to their known diversity. To account for potential phylogenomic biases, we implemented various data extraction, locus sampling, and signal filtering strategies to generate matrices. Subsequently, we applied a diverse array of tree-searching and rate-modelling methods to reconstruct the phylogeny. Our analyses, utilizing different matrices and methods, converged on the same unrooted relationships among collembolan ingroups, supporting the current ordinal classification and challenging the monophyly of Arthropleona and Symphypleona s.l. However, discrepancies across analyses existed in the root of Collembola. Among various root positions, those based on more informative matrices and biologically realistic models, favoring a basal topology of Entomobryomorpha + (Symphypleona s.s. + (Neelipleona + Poduromorpha)), were supported by subsequent methodological assessment, topology tests, and rooting analyses. This optimal topology suggests multiple independent reduction of the pronotum in non-poduromorph orders and aligns with the plesiomorphic status of neuroendocrine organs and epicuticular structure of Entomobryomorpha. Fossil-calibrated dating analyses based on the optimal topology indicated late-Paleozoic to mid-Mesozoic origins of the crown Collembola and four orders. In addition, our results questioned the monophyly of Isotomidae and Neanuridae, underscoring the need for further attention to the systematics of these families. Overall, this study provides novel insights into the phylogenetic backbone of Collembola, which will inform future studies on the systematics, ecology, and evolution of this significant arthropod lineage.

Femora nutrient foramina and aerobic capacity in giant extinct xenarthrans

Nutrient foramina are small openings in the periosteal surface of the mid-shaft region of long bones that traverse the cortical layer and reach the medullary cavity. They are important for the delivery of nutrients and oxygen to bone tissue and are crucial for the repair and remodeling of bones over time. The nutrient foramina in the femur's diaphysis are related to the energetic needs of the femur and have been shown to be related to the maximum metabolic rate (MMR) of taxa. Here, we investigate the relationship between nutrient foramen size and body mass as a proxy to the aerobic capacity of taxa in living and extinct xenarthrans, including living sloths, anteaters, and armadillos, as well as extinct xenarthrans such as glyptodonts, pampatheres, and ground sloths. Seventy femora were sampled, including 20 from extant taxa and 50 from extinct taxa. We obtained the blood flow rate (Q̇) based on foramina area and performed PGLS and phylogenetic ANCOVA in order to explore differences among mammalian groups. Our results show that, among mammals, taxa commonly associated with lower metabolism like living xenarthrans showed relatively smaller foramina, while the foramina of giant extinct xenarthrans like ground sloths and glyptodonts overlapped with non-xenarthran placentals. Consequently, Q̇ estimations indicated aerobic capacities comparable to other placental giant taxa like elephants or some ungulates. Furthermore, the estimation of the MMR for fossil giant taxa showed similar results, with almost all taxa showing high values except for those for which strong semi-arboreal or fossorial habits have been proposed. Moreover, the results are compatible with the diets predicted for extinct taxa, which indicate a strong consumption of grass similar to ungulates and in contrast to the folivorous or insectivorous diets of extant xenarthrans. The ancestral reconstruction of the MMR values indicated a lack of a common pattern for all xenarthrans, strongly supporting the occurrence of low metabolic rates in extant forms due to their particular dietary preferences and arboreal or fossorial habits. Our results highlight the importance of considering different evidence beyond the phylogenetic position of extinct taxa, especially when extinct forms are exceptionally different from their extant relatives. Future studies evaluating the energetic needs of giant extinct xenarthrans should not assume lower metabolic rates for these extinct animals based solely on their phylogenetic position and the observations on their extant relatives.

Genome and life-history evolution link bird diversification to the end-Cretaceous mass extinction

Complex patterns of genome evolution associated with the end-Cretaceous [Cretaceous-Paleogene (K-Pg)] mass extinction limit our understanding of the early evolutionary history of modern birds. Here, we analyzed patterns of avian molecular evolution and identified distinct macroevolutionary regimes across exons, introns, untranslated regions, and mitochondrial genomes. Bird clades originating near the K-Pg boundary exhibited numerous shifts in the mode of molecular evolution, suggesting a burst of genomic heterogeneity at this point in Earth's history. These inferred shifts in substitution patterns were closely related to evolutionary shifts in developmental mode, adult body mass, and patterns of metabolic scaling. Our results suggest that the end-Cretaceous mass extinction triggered integrated patterns of evolution across avian genomes, physiology, and life history near the dawn of the modern bird radiation.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Codon usage bias analysis in the mitochondrial genomes of five Rhingia Scopoli (Diptera, Syrphidae, Eristalinae) species

This study presents the sequencing and annotation of mitochondrial genomes from five Rhingia species of the family Syrphidae, focusing on codon bias. Each species possessed 22 tRNAs genes, 13 protein-coding genes, 2 rRNAs genes, and a control region, without any observed gene rearrangements. Nucleotide composition analysis revealed a higher AT content compared with GC content, indicating AT enrichment. Neutrality plot, Parity rule 2 bias, and effective number of codons plot analyses collectively indicated that natural selection primarily influences the codon usage bias in the five Rhingia species. Relative synonymous codon usage analysis identified the optimal codons for Rhingia binotata, R. fromosana, R. campestris, R. louguanensis, and R. xanthopoda as 10, 14, 10, 11, and 12, respectively, all ending with A/U and exhibiting AT preference. Phylogenetic analysis, based on maximum likelihood and Bayesian inference methods applied to three datasets, confirmed the monophyly of Rhingia. In conclusion, this research establishes a foundation for understanding the phylogenetic evolution and codon usage patterns in Rhingia, offering valuable for future studies.

Auditory brainstem responses in the nine-banded armadillo (Dasypus novemcinctus)

The auditory brainstem response (ABR) to tone burst stimuli of thirteen frequencies ranging from 0.5 to 48 kHz was recorded in the nine-banded armadillo (Dasypus novemcinctus), the only extant member of the placental mammal superorder Xenarthra in North America. The armadillo ABR consisted of five main peaks that were visible within the first 10 ms when stimuli were presented at high intensities. The latency of peak I of the armadillo ABR increased as stimulus intensity decreased by an average of 20 μs/dB. Estimated frequency-specific thresholds identified by the ABR were used to construct an estimate of the armadillo audiogram describing the mean thresholds of the eight animals tested. The majority of animals tested (six out of eight) exhibited clear responses to stimuli from 0.5 to 38 kHz, and two animals exhibited responses to stimuli of 48 kHz. Across all cases, the lowest thresholds were observed for frequencies from 8 to 12 kHz. Overall, we observed that the armadillo estimated audiogram bears a similar pattern as those observed using ABR in members of other mammalian clades, including marsupials and later-derived placental mammals.

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.

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

Multiple Outgroups Can Cause Random Rooting in Phylogenomics

Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.

Mammalian dental diversity: an evolutionary template for regenerative dentistry

The discovery of odontogenic mechanisms essential for regenerating dental tissues and eventually developing a biomimetic artificial whole tooth for replacement is an ongoing aspiration for dental clinicians and researchers. Studying the diversity, development and evolutionary changes of mammalian dentitions can provide key insights into the mechanisms of odontogenesis that can be harnessed for regenerative dental medicine. A myriad of influences is expected to have shaped the dentitions of mammals and our objective is to highlight the contributions of phylogeny, functional adaptation, and development to tooth shape. Innovations in tooth shape analysis will be discussed, such as in imaging methodologies and quantitative comparisons, molecular biology approaches to phylogeny and the ontogenetic basis of tooth form. Study of the inter- and intra-species differences in tooth form as well as dental anomalies has provided clues toward the mechanisms of evolutionary change in dental form. Thus, phenotypic variation in tooth shape will also be discussed, including the role of development in creating tooth shape differences that evolutionary selection pressures may act upon. Functional adaptations have occurred in the context of the phylogenetic signal of primitive mammals, and predecessors to each phylogenetic branch, and examples will be discussed within members of the Order Carnivora, the Superfamily Suoidea and the Order Primates. The comparative study of mammalian tooth shape holds the potential to inform dental research areas, such as etiopathogeneses of dental variation and tooth shape anomalies, molecular mechanisms of tooth development and functional issues. Ultimately, insights from these research areas can be potentially translated for futuristic clinical applications like regeneration of various tooth tissue layers and eventually full tooth replacement.

A burst of genomic innovation at the origin of placental mammals mediated embryo implantation

The origin of embryo implantation in mammals ~148 million years ago was a dramatic shift in reproductive strategy, yet the molecular changes that established mammal implantation are largely unknown. Although progesterone receptor signalling predates the origin of mammals and is highly conserved in, and critical for, successful mammal pregnancy, it alone cannot explain the origin and subsequent diversity of implantation strategies throughout the placental mammal radiation. MiRNAs are known to be flexible and dynamic regulators with a well-established role in the pathophysiology of mammal placenta. We propose that a dynamic core microRNA (miRNA) network originated early in placental mammal evolution, responds to conserved mammal pregnancy cues (e.g. progesterone), and facilitates species-specific responses. Here we identify 13 miRNA gene families that arose at the origin of placental mammals and were subsequently retained in all descendent lineages. The expression of these miRNAs in response to early pregnancy molecules is regulated in a species-specific manner in endometrial epithelia of species with extreme implantation strategies (i.e. bovine and human). Furthermore, this set of miRNAs preferentially target proteins under positive selective pressure on the ancestral eutherian lineage. Discovery of this core embryo implantation toolkit and specifically adapted proteins helps explain the origin and evolution of implantation in mammals.

Analysis of the Mitogenomes of Two Helotid Species Provides New Insights into the Phylogenetic Relationship of the Basal Cucujoidea (Insecta: Coleoptera)

Helotid beetles are commonly found in places where sap flows from tree trunks and in crevices in bark. The Helotidae family is a rare and primitive group of Cucujoidea. To date, no complete mitochondrial (mt) genome has been sequenced for this family. To better understand the characteristics of the mt genome and the evolution of Cucujoidea, we sequenced and annotated the complete mt genomes of Helota thoracica (Ritsema, 1895) and Helota yehi Lee, 2017 using next-generation sequencing. These are the first record of Helotidae mt genomes. The RNA secondary structures of both species were also predicted in this study. The mt genomes of H. thoracica and H. yehi are circular, with total lengths of 16,112 bp and 16,401 bp, respectively. After comparing the mt genomes of H. thoracica and H. yehi, we observed the gene arrangement, codon usage patterns, base content, and RNA secondary structures of both species to be similar, which has also been noted in other Coleoptera insects. The nucleotide sequence of the coding regions and the control region has small differences. The phylogenetic analysis indicated that Helotidae and Protocucujidae are sister groups and revealed the relationship between seven families; however, the validity of the two series (Erotylid series and Nitidulid series) as larger groups in the superfamily was not supported. The mt phylogenomic relationships have strong statistical support. Therefore, the division of Cucujoidea into series should be re-examined. Our results will provide a better understanding of the mt genome and phylogeny of Helotidae and Cucujoidea and will provide valuable molecular markers for further genetic studies.

Eighteen mitochondrial genomes of Syrphidae (Insecta: Diptera: Brachycera) with a phylogenetic analysis of Muscomorpha

In this study, 18 mitochondrial genomes (mitogenomes) of Syrphidae were sequenced. These mitogenomes ranged from 15,648 to 16,405 bp and contained 37 genes that were similar to those from other Syrphidae species. Most protein-coding genes (PCGs) started with a standard ATN codon and ended with TAA/G. All transfer RNAs (tRNAs) could be folded into the cloverleaf secondary structure except tRNA-Ser (AGN), which lacks a dihydrouridine arm. The secondary structures of ribosomal RNAs (rRNAs) were predicted. Six domains (III is absent in arthropods) and 44 helices were included in the 16S rRNA, and three domains and 24 helices were included in the 12S rRNA. We found three conserved fragments in all syrphid mitogenomes. Phylogenetic analyses were performed based on the nucleotide data of 13 PCGs and two rRNAs from 76 Muscomorpha and three outgroup species. In results the paraphyly of Aschiza and Schizophora were supported, the Acalyptratae was also paraphyletic but the relationships of its superfamilies were difficult to determine, the monophyly of Calyptratea was supported with the relationships of Oestroidea and Muscoidea need to be further reconsidered. Within Syrphidae the monophyly of family level was supported, the Syrphinae were clustered into one branch, while the paraphyly of Eristalinae was still well supported.

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GC_CV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GC_CV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GC_CV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family

BACKGROUND

Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes.

RESULTS

We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae.

CONCLUSIONS

This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.

Phylogenomic discordance suggests polytomies along the backbone of the large genus Solanum

PREMISE

Evolutionary studies require solid phylogenetic frameworks, but increased volumes of phylogenomic data have revealed incongruent topologies among gene trees in many organisms both between and within genomes. Some of these incongruences indicate polytomies that may remain impossible to resolve. Here we investigate the degree of gene-tree discordance in Solanum, one of the largest flowering plant genera that includes the cultivated potato, tomato, and eggplant, as well as 24 minor crop plants.

METHODS

A densely sampled species-level phylogeny of Solanum is built using unpublished and publicly available Sanger sequences comprising 60% of all accepted species (742 spp.) and nine regions (ITS, waxy, and seven plastid markers). The robustness of this topology is tested by examining a full plastome dataset with 140 species and a nuclear target-capture dataset with 39 species of Solanum (Angiosperms353 probe set).

RESULTS

While the taxonomic framework of Solanum remained stable, gene tree conflicts and discordance between phylogenetic trees generated from the target-capture and plastome datasets were observed. The latter correspond to regions with short internodal branches, and network analysis and polytomy tests suggest the backbone is composed of three polytomies found at different evolutionary depths. The strongest area of discordance, near the crown node of Solanum, could potentially represent a hard polytomy.

CONCLUSIONS

We argue that incomplete lineage sorting due to rapid diversification is the most likely cause for these polytomies, and that embracing the uncertainty that underlies them is crucial to understand the evolution of large and rapidly radiating lineages.

Adaptive Evolution in TRIF Leads to Discordance between Human and Mouse Innate Immune Signaling

The TIR domain-containing adapter inducing IFN-β (TRIF) protein is an innate immune system protein that mediates the MyD88-independent toll-like receptor response pathway in mice and humans. Previously, we identified positive selection at seven distinct residues in mouse TRIF (mTRIF), as compared with human and other mammalian orthologs, thus predicting protein functional shift in mTRIF. We reconstructed TRIF for the most recent common ancestor of mouse and human, and mutated this at the seven sites to their extant mouse/human states. We overexpressed these TRIF mutants in immortalized human and mouse cell lines and monitored TRIF-dependent cytokine production and gene expression induction. We show that optimal TRIF function in human and mouse is dependent on the identity of the positively selected sites. These data provide us with molecular data relating observed differences in response between mouse and human MyD88-independent signaling in the innate immune system with protein functional change.

The evolution of insect biodiversity

Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.

Comparative mitogenome phylogeography of two anteater genera (Tamandua and Myrmecophaga; Myrmecophagidae, Xenarthra): Evidence of discrepant evolutionary traits

The species within Xenarthra (sloths, anteaters, and armadillos) are quintessential South American mammals. Of the three groups, Vermilingua (anteaters) contains the fewest extant and paleontological species. Here, we sampled and sequenced the entire mitochondrial genomes (mitogenomes) of two Tamandua species (Tamandua tetradactyla and Tamandua mexicana) (n=74) from Central and South America, as well as Myrmecophaga tridactyla (n=41) from South America. Within Tamandua, we detected three different haplogroups. The oldest (THI) contained many specimens with the T. tetradactyla morphotype (but also several with the T. mexicana morphotype) and originated in southeastern South America (currently Uruguay) before moving towards northern South America, where the THII haplogroup originated. THII primarily contained specimens with the T. mexicana morphotype (but also several with the T. tetradactyla morphotype) and was distributed in Central America, Colombia, and Ecuador. THI and THII yielded a genetic distance of 4%. THII originated in either northern South America or "in situ" in Central America with haplogroup THIII, which consisted of ~50% T. mexicana and 50% T. tetradactyla phenotypes. THIII was mostly located in the same areas as THII, i.e., Central America, Ecuador, and Colombia, though mainly in the latter. The three haplogroups overlapped in Colombia and Ecuador. Thus, T. tetradactyla and T. mexicana were not reciprocally monophyletic. For this reason, we considered that a unique species of Tamandua likely exists, i.e., T. tetradactyla. In contrast to Tamandua, M. tridactyla did not show different morphotypes throughout its geographical range in the Neotropics. However, two very divergent genetic haplogroups (MHI and MHII), with a genetic distance of ~10%, were detected. The basal haplogroup, MHI, originated in northwestern South America, whereas the more geographically derived haplogroup, MHII, overlapped with MHI, but also expanded into central and southern South America. Thus, Tamandua migrated from south to north whereas Myrmecophaga migrated from north to south. Our results also showed that temporal mitochondrial diversification for Tamandua began during the Late Pliocene and Upper Pleistocene, but for Myrmecophaga began during the Late Miocene. Furthermore, both taxa showed elevated levels of mitochondrial genetic diversity. Tamandua showed more evidence of female population expansion than Myrmecophaga. Tamandua experienced population expansion ~0.6-0.17 million years ago (Mya), whereas Myrmecophaga showed possible population expansion ~0.3-0.2 Mya. However, both taxa experienced a conspicuous female decline in the last 10 000-20 000 years. Our results also showed little spatial genetic structure for both taxa. However, several analyses revealed higher spatial structure in Tamandua than in Myrmecophaga. Therefore, Tamandua and Myrmecophaga were not subjected to the same biogeographical, geological, or climatological events in shaping their genetic structures.

Structural features of the mitogenome of the leafhopper genus Cladolidia (Hemiptera: Cicadellidae: Coelidiinae) and phylogenetic implications in Cicadellidae

The first two complete mitogenomes of the leafhopper genus Cladolidia (C. biungulata and C. robusta) were sequenced and annotated to further explore the phylogeny of Cladolidia. Both the newly sequenced mitogenomes have a typical circular structure, with lengths of 15,247 and 15,376 bp and A + T contents of 78.2% and 78%, respectively. We identified a highly conserved genome organization in the two Cladolidia spp. through comparative analysis that included the following assessments: genome content, gene order, nucleotide composition, codon usage, amino acid composition, and tRNA secondary structure. Moreover, we detected the base heterogeneity of Cicadellidae mitogenomic data and constructed phylogenetic trees using the nucleotide alignments of 12 subfamilies of 58 leafhopper species. We noted a weak heterogeneity in the base composition among the Cicadellidae mitogenomes. Phylogenetic analyses showed that the monophyly of each subfamily was generally well supported in the family Cicadellidae; the main topology was as follows: (Deltocephalinae + (Treehoppers + ((Megophthalminae + (Macropsinae + (Hylicinae + (Coelidiinae +Iassinae)) + (Idiocerinae + (Cicadellinae + (Typhlocybinae + (Mileewinae + (Evacanthinae +Ledrinae)))))))))). Within Coelidiinae, phylogenetic analyses revealed that C. biungulata and C. robusta belong to Coelidiinae and the monophyly of Cladolidia is well supported. In addition, on the basis of complete mitogenome phylogenetic analysis and the comparison of morphological characteristics, we further confirm the genus Olidiana as a paraphyletic group, suggesting that the genus may need taxonomic revisions.

Is Phylotranscriptomics as Reliable as Phylogenomics?

Phylogenomics, the study of phylogenetic relationships among taxa based on their genome sequences, has emerged as the preferred phylogenetic method because of the wealth of phylogenetic information contained in genome sequences. Genome sequencing, however, can be prohibitively expensive, especially for taxa with huge genomes and when many taxa need sequencing. Consequently, the less costly phylotranscriptomics has seen an increased use in recent years. Phylotranscriptomics reconstructs phylogenies using DNA sequences derived from transcriptomes, which are often orders of magnitude smaller than genomes. However, in the absence of corresponding genome sequences, comparative analyses of transcriptomes can be challenging and it is unclear whether phylotranscriptomics is as reliable as phylogenomics. Here, we respectively compare the phylogenomic and phylotranscriptomic trees of 22 mammals and 15 plants that have both sequenced nuclear genomes and publicly available RNA sequencing data from multiple tissues. We found that phylotranscriptomic analysis can be sensitive to orthologous gene identification. When a rigorous method for identifying orthologs is employed, phylogenomic and phylotranscriptomic trees are virtually identical to each other, regardless of the tissue of origin of the transcriptomes and whether the same tissue is used across species. These findings validate phylotranscriptomics, brighten its prospect, and illustrate the criticality of reliable ortholog detection in such practices.

Phylogenomic conflict coincides with rapid morphological innovation

Evolutionary biologists have long been fascinated with the episodes of rapid phenotypic innovation that underlie the emergence of major lineages. Although our understanding of the environmental and ecological contexts of such episodes has steadily increased, it has remained unclear how population processes contribute to emergent macroevolutionary patterns. One insight gleaned from phylogenomics is that gene-tree conflict, frequently caused by population-level processes, is often rampant during the origin of major lineages. With the understanding that phylogenomic conflict is often driven by complex population processes, we hypothesized that there may be a direct correspondence between instances of high conflict and elevated rates of phenotypic innovation if both patterns result from the same processes. We evaluated this hypothesis in six clades spanning vertebrates and plants. We found that the most conflict-rich regions of these six clades also tended to experience the highest rates of phenotypic innovation, suggesting that population processes shaping both phenotypic and genomic evolution may leave signatures at deep timescales. Closer examination of the biological significance of phylogenomic conflict may yield improved connections between micro- and macroevolution and increase our understanding of the processes that shape the origin of major lineages across the Tree of Life.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Interrogating Genomic-Scale Data to Resolve Recalcitrant Nodes in the Spider Tree of Life

Genome-scale data sets are converging on robust, stable phylogenetic hypotheses for many lineages; however, some nodes have shown disagreement across classes of data. We use spiders (Araneae) as a system to identify the causes of incongruence in phylogenetic signal between three classes of data: exons (as in phylotranscriptomics), noncoding regions (included in ultraconserved elements [UCE] analyses), and a combination of both (as in UCE analyses). Gene orthologs, coded as amino acids and nucleotides (with and without third codon positions), were generated by querying published transcriptomes for UCEs, recovering 1,931 UCE loci (codingUCEs). We expected that congeners represented in the codingUCE and UCEs data would form clades in the presence of phylogenetic signal. Noncoding regions derived from UCE sequences were recovered to test the stability of relationships. Phylogenetic relationships resulting from all analyses were largely congruent. All nucleotide data sets from transcriptomes, UCEs, or a combination of both recovered similar topologies in contrast with results from transcriptomes analyzed as amino acids. Most relationships inferred from low-occupancy data sets, containing several hundreds of loci, were congruent across Araneae, as opposed to high occupancy data matrices with fewer loci, which showed more variation. Furthermore, we found that low-occupancy data sets analyzed as nucleotides (as is typical of UCE data sets) can result in more congruent relationships than high occupancy data sets analyzed as amino acids (as in phylotranscriptomics). Thus, omitting data, through amino acid translation or via retention of only high occupancy loci, may have a deleterious effect in phylogenetic reconstruction.

Phylogenetic tree building in the genomic age

Knowing phylogenetic relationships among species is fundamental for many studies in biology. An accurate phylogenetic tree underpins our understanding of the major transitions in evolution, such as the emergence of new body plans or metabolism, and is key to inferring the origin of new genes, detecting molecular adaptation, understanding morphological character evolution and reconstructing demographic changes in recently diverged species. Although data are ever more plentiful and powerful analysis methods are available, there remain many challenges to reliable tree building. Here, we discuss the major steps of phylogenetic analysis, including identification of orthologous genes or proteins, multiple sequence alignment, and choice of substitution models and inference methodologies. Understanding the different sources of errors and the strategies to mitigate them is essential for assembling an accurate tree of life.

Gene Fusions Derived by Transcriptional Readthrough are Driven by Segmental Duplication in Human

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products.

A Novel Test for Absolute Fit of Evolutionary Models Provides a Means to Correctly Identify the Substitution Model and the Model Tree

A novel test is described that visualizes the absolute model-data fit of the substitution and tree components of an evolutionary model. The test utilizes statistics based on counts of character state matches and mismatches in alignments of observed and simulated sequences. This comparison is used to assess model-data fit. In simulations conducted to evaluate the performance of the test, the test estimator was able to identify both the correct tree topology and substitution model under conditions where the Goldman-Cox test-which tests the fit of a substitution model to sequence data and is also based on comparing simulated replicates with observed data-showed high error rates. The novel test was found to identify the correct tree topology within a wide range of DNA substitution model misspecifications, indicating the high discriminatory power of the test. Use of this test provides a practical approach for assessing absolute model-data fit when testing phylogenetic hypotheses.

Large-scale genomic sequence data resolve the deepest divergences in the legume phylogeny and support a near-simultaneous evolutionary origin of all six subfamilies

Phylogenomics is increasingly used to infer deep-branching relationships while revealing the complexity of evolutionary processes such as incomplete lineage sorting, hybridization/introgression and polyploidization. We investigate the deep-branching relationships among subfamilies of the Leguminosae (or Fabaceae), the third largest angiosperm family. Despite their ecological and economic importance, a robust phylogenetic framework for legumes based on genome-scale sequence data is lacking. We generated alignments of 72 chloroplast genes and 7621 homologous nuclear-encoded proteins, for 157 and 76 taxa, respectively. We analysed these with maximum likelihood, Bayesian inference, and a multispecies coalescent summary method, and evaluated support for alternative topologies across gene trees. We resolve the deepest divergences in the legume phylogeny despite lack of phylogenetic signal across all chloroplast genes and the majority of nuclear genes. Strongly supported conflict in the remainder of nuclear genes is suggestive of incomplete lineage sorting. All six subfamilies originated nearly simultaneously, suggesting that the prevailing view of some subfamilies as 'basal' or 'early-diverging' with respect to others should be abandoned, which has important implications for understanding the evolution of legume diversity and traits. Our study highlights the limits of phylogenetic resolution in relation to rapid successive speciation.

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.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae)

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

Rare coral under the genomic microscope: timing and relationships among Hawaiian Montipora

Background

Evolutionary patterns of scleractinian (stony) corals are difficult to infer given the existence of few diagnostic characters and pervasive phenotypic plasticity. A previous study of Hawaiian Montipora (Scleractinia: Acroporidae) based on five partial mitochondrial and two nuclear genes revealed the existence of a species complex, grouping one of the rarest known species (M. dilatata, which is listed as Endangered by the International Union for Conservation of Nature - IUCN) with widespread corals of very different colony growth forms (M. flabellata and M. cf. turgescens). These previous results could result from a lack of resolution due to a limited number of markers, compositional heterogeneity or reflect biological processes such as incomplete lineage sorting (ILS) or introgression.

Results

All 13 mitochondrial protein-coding genes from 55 scleractinians (14 lineages from this study) were used to evaluate if a recent origin of the M. dilatata species complex or rate heterogeneity could be compromising phylogenetic inference. Rate heterogeneity detected in the mitochondrial data set seems to have no significant impacts on the phylogenies but clearly affects age estimates. Dating analyses show different estimations for the speciation of M. dilatata species complex depending on whether taking compositional heterogeneity into account (0.8 [0.05–2.6] Myr) or assuming rate homogeneity (0.4 [0.14–0.75] Myr). Genomic data also provided evidence of introgression among all analysed samples of the complex. RADseq data indicated that M. capitata colour morphs may have a genetic basis.

Conclusions

Despite the volume of data (over 60,000 SNPs), phylogenetic relationships within the M. dilatata species complex remain unresolved most likely due to a recent origin and ongoing introgression. Species delimitation with genomic data is not concordant with the current taxonomy, which does not reflect the true diversity of this group. Nominal species within the complex are either undergoing a speciation process or represent ecomorphs exhibiting phenotypic polymorphisms.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1476-2) contains supplementary material, which is available to authorized users.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

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.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

Analyzing Contentious Relationships and Outlier Genes in Phylogenomics

Recent studies have demonstrated that conflict is common among gene trees in phylogenomic studies, and that less than one percent of genes may ultimately drive species tree inference in supermatrix analyses. Herein, we examined two data sets where supermatrix and coalescent-based species trees conflict. We identified two highly influential "outlier" genes in each data set. When removed from each data set, the inferred supermatrix trees matched the topologies obtained from coalescent analyses. We also demonstrate that, while the outlier genes in the vertebrate data set have been shown in a previous study to be the result of errors in orthology detection, the outlier genes from a plant data set did not exhibit any obvious systematic error, and therefore, may be the result of some biological process yet to be determined. While topological comparisons among a small set of alternate topologies can be helpful in discovering outlier genes, they can be limited in several ways, such as assuming all genes share the same topology. Coalescent species tree methods relax this assumption but do not explicitly facilitate the examination of specific edges. Coalescent methods often also assume that conflict is the result of incomplete lineage sorting. Herein, we explored a framework that allows for quickly examining alternative edges and support for large phylogenomic data sets that does not assume a single topology for all genes. For both data sets, these analyses provided detailed results confirming the support for coalescent-based topologies. This framework suggests that we can improve our understanding of the underlying signal in phylogenomic data sets by asking more targeted edge-based questions.

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.

Improving phylogenetic inference of core Chlorophyta using chloroplast sequences with strong phylogenetic signals and heterogeneous models

Phylogenetic relationships within the green algal phylum Chlorophyta have proven difficult to resolve. The core Chlorophyta include Chlorophyceae, Ulvophyceae, Trebouxiophyceae, Pedinophyceae and Chlorodendrophyceae, but the relationships among these classes remain unresolved and the monophyly of Ulvophyceae and Trebouxiophyceae are highly controversial. We analyzed a dataset of 101 green algal species and 73 protein-coding genes sampled from complete and partial chloroplast genomes, including six newly sequenced ulvophyte genomes (Blidingia minima NIES-1837, Ulothrix zonata, Halochlorococcum sp. NIES-1838, Scotinosphaera sp. NIES-154, Caulerpa brownii and Cephaleuros sp. HZ-2017). We applied the Tree Certainty (TC) score to quantify the level of incongruence between phylogenetic trees in chloroplast genomic datasets, and show that the conflicting phylogenetic trees of core Chlorophyta stem from the most GC-heterogeneous sites. With removing the most GC-heterogeneous sites, our chloroplast phylogenomic analyses using heterogeneous models consistently support monophyly of the Chlorophyceae and of the Trebouxiophyceae, but the Ulvophyceae was resolved as polyphyletic. Our analytical framework provides an efficient approach to reconstruct the optimal phylogenetic relationships by minimizing conflicting signals.

Sequencing of complete mitochondrial genomes confirms synonymization of Hyalomma asiaticum asiaticum and kozlovi, and advances phylogenetic hypotheses for the Ixodidae

Phylogeny of hard ticks (Ixodidae) remains unresolved. Mitochondrial genomes (mitogenomes) are increasingly used to resolve phylogenetic controversies, but remain unavailable for the entire large Hyalomma genus. Hyalomma asiaticum is a parasitic tick distributed throughout the Asia. As a result of great morphological variability, two subspecies have been recognised historically; until a morphological data-based synonymization was proposed. However, this hypothesis was never tested using molecular data. Therefore, objectives of this study were to: 1. sequence the first Hyalomma mitogenome; 2. scrutinise the proposed synonymization using molecular data, i.e. complete mitogenomes of both subspecies: H. a. asiaticum and kozlovi; 3. conduct phylogenomic and comparative analyses of all available Ixodidae mitogenomes. Results corroborate the proposed synonymization: the two mitogenomes are almost identical (99.6%). Genomic features of both mitogenomes are standard for Metastriata; which includes the presence of two control regions and all three "Tick-Box" motifs. Gene order and strand distribution are perfectly conserved for the entire Metastriata group. Suspecting compositional biases, we conducted phylogenetic analyses (29 almost complete mitogenomes) using homogeneous and heterogeneous (CAT) models of substitution. The results were congruent, apart from the deep-level topology of prostriate ticks (Ixodes): the homogeneous model produced a monophyletic Ixodes, but the CAT model produced a paraphyletic Ixodes (and thereby Prostriata), divided into Australasian and non-Australasian clades. This topology implies that all metastriate ticks have evolved from the ancestor of the non-Australian branch of prostriate ticks. Metastriata was divided into three clades: 1. Amblyomminae and Rhipicephalinae (Rhipicephalus, Hyalomma, Dermacentor); 2. Haemaphysalinae and Bothriocrotoninae, plus Amblyomma sphenodonti; 3. Amblyomma elaphense, basal to all Metastriata. We conclude that mitogenomes have the potential to resolve the long-standing debate about the evolutionary history of ticks, but heterogeneous evolutionary models should be used to alleviate the effects of compositional heterogeneity on deep-level relationships.

Phylogenomic Resolution of the Phylogeny of Laurasiatherian Mammals: Exploring Phylogenetic Signals within Coding and Noncoding Sequences

The interordinal relationships of Laurasiatherian mammals are currently one of the most controversial questions in mammalian phylogenetics. Previous studies mainly relied on coding sequences (CDS) and seldom used noncoding sequences. Here, by data mining public genome data, we compiled an intron data set of 3,638 genes (all introns from a protein-coding gene are considered as a gene) (19,055,073 bp) and a CDS data set of 10,259 genes (20,994,285 bp), covering all major lineages of Laurasiatheria (except Pholidota). We found that the intron data contained stronger and more congruent phylogenetic signals than the CDS data. In agreement with this observation, concatenation and species-tree analyses of the intron data set yielded well-resolved and identical phylogenies, whereas the CDS data set produced weakly supported and incongruent results. Further analyses showed that the phylogeny inferred from the intron data is highly robust to data subsampling and change in outgroup, but the CDS data produced unstable results under the same conditions. Interestingly, gene tree statistical results showed that the most frequently observed gene tree topologies for the CDS and intron data are identical, suggesting that the major phylogenetic signal within the CDS data is actually congruent with that within the intron data. Our final result of Laurasiatheria phylogeny is (Eulipotyphla,((Chiroptera, Perissodactyla),(Carnivora, Cetartiodactyla))), favoring a close relationship between Chiroptera and Perissodactyla. Our study 1) provides a well-supported phylogenetic framework for Laurasiatheria, representing a step towards ending the long-standing "hard" polytomy and 2) argues that intron within genome data is a promising data resource for resolving rapid radiation events across the tree of life.

Novel and divergent genes in the evolution of placental mammals

Analysis of genome sequences within a phylogenetic context can give insight into the mode and tempo of gene and protein evolution, including inference of gene ages. This can reveal whether new genes arose on particular evolutionary lineages and were recruited for new functional roles. Here, we apply MCL clustering with all-versus-all reciprocal BLASTP to identify and phylogenetically date 'Homology Groups' among vertebrate proteins. Homology Groups include new genes and highly divergent duplicate genes. Focusing on the origin of the placental mammals within the Eutheria, we identify 357 novel Homology Groups that arose on the stem lineage of Placentalia, 87 of which are deduced to play core roles in mammalian biology as judged by extensive retention in evolution. We find the human homologues of novel eutherian genes are enriched for expression in preimplantation embryo, brain, and testes, and enriched for functions in keratinization, reproductive development, and the immune system.