Gene trees and species trees are not the same

Museum specimens shedding light on the evolutionary history and cryptic diversity of the hedgehog family Erinaceidae

The family Erinaceidae encompasses 27 extant species in two subfamilies: Erinaceinae, which includes spiny hedgehogs, and Galericinae, which comprises silky-furred gymnures and moonrats. Although they are commonly recognized by the general public, their phylogenetic history remains incompletely understood, and several species have never been included in any molecular analyses. Additionally, previous research suggested that the species diversity of Erinaceidae might be underestimated. In this study, we sequenced the mitochondrial genomes of 29 individuals representing 18 erinaceid species using 18 freshly collected tissue and 11 historical museum specimens. We also integrated previously published data for a concatenated analysis. We aimed to elucidate the evolutionary relationships within Erinaceidae, estimate divergence times, and uncover potential underestimated species diversity. Our data finely resolved intergeneric and interspecific relationships and presented the first molecular evidence for the phylogenetic position of Mesechinus wangi, Paraechinus micropus, and P. nudiventris. Our results revealed a sister relationship between Neotetracus and Neohylomys gymnures, as well as a sister relationship between Hemiechinus and Mesechinus, supporting previous hypotheses. Additionally, our findings provided a novel phylogenetic position for Paraechinus aethiopicus, placing it in a basal position within the genus. Furthermore, our study uncovered cryptic species diversity within Hylomys suillus as well as in Neotetracus sinensis, Atelerix albiventris, P. aethiopicus, and Hemiechinus auratus, most of which have been previously overlooked.

Phylogenomic analyses of ochrophytes (stramenopiles) with an emphasis on neglected lineages

Ochrophyta is a photosynthetic lineage that crowns the phylogenetic tree of stramenopiles, one of the major eukaryotic supergroups. Due to their ecological impact as a major primary producer, ochrophytes are relatively well-studied compared to the rest of the stramenopiles, yet their evolutionary relationships remain poorly understood. This is in part due to a number of missing lineages in large-scale multigene analyses, and an apparently rapid radiation leading to many short internodes between ochrophyte subgroups in the tree. These short internodes are also found across deep-branching lineages of stramenopiles with limited phylogenetic signal, leaving many relationships controversial overall. We have addressed this issue with other deep-branching stramenopiles recently, and now examine whether contentious relationships within the ochrophytes may be resolved with the help of filling in missing lineages in an updated phylogenomic dataset of ochrophytes, along with exploring various gene filtering criteria to identify the most phylogenetically informative genes. We generated ten new transcriptomes from various culture collections and a single-cell isolation from an environmental sample, added these to an existing phylogenomic dataset, and examined the effects of selecting genes with high phylogenetic signal or low phylogenetic noise. For some previously contentious relationships, we find a variety of analyses and gene filtering criteria consistently unite previously unstable groupings with strong statistical support. For example, we recovered a robust grouping of Eustigmatophyceae with Raphidophyceae-Phaeophyceae-Xanthophyceae while Olisthodiscophyceae formed a sister-lineage to Pinguiophyceae. Selecting genes with high phylogenetic signal or data quality recovered more stable topologies. Overall, we find that adding under-represented groups across different lineages is still crucial in resolving phylogenetic relationships, and discrete gene properties affect lineages of stramenopiles differently. This is something which may be explored to further our understanding of the molecular evolution of stramenopiles.

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting (ILS), introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call mixtures across sites and trees (MAST). This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of ILS in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of 4 Platyrrhine species for which standard concatenated maximum likelihood (ML) and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e., the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyze a concatenated alignment using ML while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae)

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

Deep-time phylogenetic inference by paleoproteomic analysis of dental enamel

In temperate and subtropical regions, ancient proteins are reported to survive up to about 2 million years, far beyond the known limits of ancient DNA preservation in the same areas. Accordingly, their amino acid sequences currently represent the only source of genetic information available to pursue phylogenetic inference involving species that went extinct too long ago to be amenable for ancient DNA analysis. Here we present a complete workflow, including sample preparation, mass spectrometric data acquisition and computational analysis, to recover and interpret million-year-old dental enamel protein sequences. During sample preparation, the proteolytic digestion step, usually an integral part of conventional bottom-up proteomics, is omitted to increase the recovery of the randomly degraded peptides spontaneously generated by extensive diagenetic hydrolysis of ancient proteins over geological time. Similarly, we describe other solutions we have adopted to (1) authenticate the endogenous origin of the protein traces we identify, (2) detect and validate amino acid variation in the ancient protein sequences and (3) attempt phylogenetic inference. Sample preparation and data acquisition can be completed in 3-4 working days, while subsequent data analysis usually takes 2-5 days. The workflow described requires basic expertise in ancient biomolecules analysis, mass spectrometry-based proteomics and molecular phylogeny. Finally, we describe the limits of this approach and its potential for the reconstruction of evolutionary relationships in paleontology and paleoanthropology.

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Testing plastomes and nuclear ribosomal DNA sequences as the next-generation DNA barcodes for species identification and phylogenetic analysis in Acer

BACKGROUND

Acer is a taxonomically intractable and speciose genus that contains over 150 species. It is challenging to distinguish Acer species only by morphological method due to their abundant variations. Plastome and nuclear ribosomal DNA (nrDNA) sequences are recommended as powerful next-generation DNA barcodes for species discrimination. However, their efficacies were still poorly studied. The current study will evaluate the application of plastome and nrDNA in species identification and perform phylogenetic analyses for Acer.

RESULT

Based on a collection of 83 individuals representing 55 species (c. 55% of Chinese species) from 13 sections, our barcoding analyses demonstrated that plastomes exhibited the highest (90.47%) species discriminatory power among all plastid DNA markers, such as the standard plastid barcodes matK + rbcL + trnH-psbA (61.90%) and ycf1 (76.19%). And the nrDNA (80.95%) revealed higher species resolution than ITS (71.43%). Acer plastomes show abundant interspecific variations, however, species identification failure may be due to the incomplete lineage sorting (ILS) and chloroplast capture resulting from hybridization. We found that the usage of nrDNA contributed to identifying those species that were unidentified by plastomes, implying its capability to some extent to mitigate the impact of hybridization and ILS on species discrimination. However, combining plastome and nrDNA is not recommended given the cytonuclear conflict caused by potential hybridization. Our phylogenetic analysis covering 19 sections (95% sections of Acer) and 128 species (over 80% species of this genus) revealed pervasive inter- and intra-section cytonuclear discordances, hinting that hybridization has played an important role in the evolution of Acer.

CONCLUSION

Plastomes and nrDNA can significantly improve the species resolution in Acer. Our phylogenetic analysis uncovered the scope and depth of cytonuclear conflict in Acer, providing important insights into its evolution.

Implications of gene tree heterogeneity on downstream phylogenetic analyses: A case study employing the Fair Proportion index

Many questions in evolutionary biology require the specification of a phylogeny for downstream phylogenetic analyses. However, with the increasingly widespread availability of genomic data, phylogenetic studies are often confronted with conflicting signal in the form of genomic heterogeneity and incongruence between gene trees and the species tree. This raises the question of determining what data and phylogeny should be used in downstream analyses, and to what extent the choice of phylogeny (e.g., gene trees versus species trees) impacts the analyses and their outcomes. In this paper, we study this question in the realm of phylogenetic diversity indices, which provide ways to prioritize species for conservation based on their relative evolutionary isolation on a phylogeny, and are thus one example of downstream phylogenetic analyses. We use the Fair Proportion (FP) index, also known as the evolutionary distinctiveness score, and explore the variability in species rankings based on gene trees as compared to the species tree for several empirical data sets. Our results indicate that prioritization rankings among species vary greatly depending on the underlying phylogeny, suggesting that the choice of phylogeny is a major influence in assessing phylogenetic diversity in a conservation setting. While we use phylogenetic diversity conservation as an example, we suspect that other types of downstream phylogenetic analyses such as ancestral state reconstruction are similarly affected by genomic heterogeneity and incongruence. Our aim is thus to raise awareness of this issue and inspire new research on which evolutionary information (species trees, gene trees, or a combination of both) should form the basis for analyses in these settings.

Molecular and morphometric study of Brazilian populations of Psychodopygus davisi

In this study, we analysed the molecular and morphometric differences of several populations of the putative sand fly vector Psychodopygus davisi (Root, 1934) (Diptera, Psychodidae, Phlebotominae) in Brazil. We amplified the 658 base pair fragments of the DNA barcoding region-cytochrome c oxidase subunit 1 (COI) gene-for 57 specimens of P. davisi and three specimens of Psychodopygus claustrei (Abonnenc, Léger & Fauran, 1979). We merged our data with public sequences of the same species available from GenBank. Then, the combined dataset-87 sequences and 20 localities-was analysed using population structure analysis and different species delimitation approaches. Geometric morphometry of wings was performed for 155 specimens of P. davisi populations from the North, Midwest and Southeast Brazilian regions, analysing the differences in centroid sizes and canonical variates. Molecular analysis indicated high intraspecific genetic distance values for P. davisi (maximum p distance = 5.52%). All algorithms identified P. davisi and P. claustrei as distinct molecular taxonomic units, despite the low interspecific distance (p distance to the nearest neighbour = 4.79%). P. davisi sequences were split into four genetic clusters by population structure analysis and at least five genetic lineages using intermediate scenarios of the species delimitation algorithms. The species validation analysis of BPP strongly supported the five-species model in our dataset. We found high genetic diversity in this taxon, which is in agreement with its wide geographic distribution in Brazil. Furthermore, the wing analysis showed that specimens from the Southeast Region of Brazil are different from those in the North and the Midwest. The evolutionary patterns of P. davisi populations in Brazil suggest the presence of candidate species, which need to be validated in future studies using a more comprehensive approach with both genomic data and morphological characters.

Likelihood-Based Tests of Species Tree Hypotheses

Likelihood-based tests of phylogenetic trees are a foundation of modern systematics. Over the past decade, an enormous wealth and diversity of model-based approaches have been developed for phylogenetic inference of both gene trees and species trees. However, while many techniques exist for conducting formal likelihood-based tests of gene trees, such frameworks are comparatively underdeveloped and underutilized for testing species tree hypotheses. To date, widely used tests of tree topology are designed to assess the fit of classical models of molecular sequence data and individual gene trees and thus are not readily applicable to the problem of species tree inference. To address this issue, we derive several analogous likelihood-based approaches for testing topologies using modern species tree models and heuristic algorithms that use gene tree topologies as input for maximum likelihood estimation under the multispecies coalescent. For the purpose of comparing support for species trees, these tests leverage the statistical procedures of their original gene tree-based counterparts that have an extended history for testing phylogenetic hypotheses at a single locus. We discuss and demonstrate a number of applications, limitations, and important considerations of these tests using simulated and empirical phylogenomic data sets that include both bifurcating topologies and reticulate network models of species relationships. Finally, we introduce the open-source R package SpeciesTopoTestR (SpeciesTopology Tests in R) that includes a suite of functions for conducting formal likelihood-based tests of species topologies given a set of input gene tree topologies.

A Comprehensive Overview of Baboon Phylogenetic History

Baboons (genus Papio) are an intriguing study system to investigate complex evolutionary processes and the evolution of social systems. An increasing number of studies over the last 20 years has shown that considerable incongruences exist between phylogenies based on morphology, mitochondrial, and nuclear sequence data of modern baboons, and hybridization and introgression have been suggested as the main drivers of these patterns. Baboons, therefore, present an excellent opportunity to study these phenomena and their impact on speciation. Advances both in geographic and genomic coverage provide increasing details on the complexity of the phylogeography of baboons. Here, we compile the georeferenced genetic data of baboons and review the current knowledge on baboon phylogeny, discuss the evolutionary processes that may have shaped the patterns that we observe today, and propose future avenues for research.

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Subspecific Nomenclature of Giardia duodenalis in the Light of a Compared Population Genomics of Pathogens

Genetic and genomic data have long recognized that the species Giardia duodenalis is subdivided into at least eight genetic clusters that have been named "assemblages" by specialists in the field. Some of these assemblages have been given the status of species, with Linnean binames. In the framework of the predominant clonal evolution model (PCE), we have shown that, from an evolutionary point of view, G. duodenalis assemblages are equatable to "near-clades", that is to say: clades whose discreteness is somewhat clouded by occasional genetic exchange, but remain discrete and stable in space and time. The implications of this evolutionary status for the species described within G. duodenalis are discussed in light of the most recent genetic and genomic studies. The pattern of this species' subspecific genetic variability and genetic clustering appears to be very similar to the ones of various parasitic, fungal and bacteria species. This underlines the relevance of a compared population genomics of pathogenic species allowed by the broad framework of the PCE model.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.

Effects of discordance between species and gene trees on phylogenetic diversity conservation

Phylogenetic diversity indices such as the Fair Proportion (FP) index are frequently discussed as prioritization criteria in biodiversity conservation. They rank species according to their contribution to overall diversity by taking into account the unique and shared evolutionary history of each species as indicated by its placement in an underlying phylogenetic tree. Traditionally, phylogenetic trees were inferred from single genes and the resulting gene trees were assumed to be a valid estimate for the species tree, i.e., the "true" evolutionary history of the species under consideration. However, nowadays it is common to sequence whole genomes of hundreds or thousands of genes, and it is often the case that conflicting genealogical histories exist in different genes throughout the genome, resulting in discordance between individual gene trees and the species tree. Here, we analyze the effects of gene and species tree discordance on prioritization decisions based on the FP index. In particular, we consider the ranking order of taxa induced by (i) The FP index on a species tree, and (ii) The expected FP index across all gene tree histories associated with the species tree. On the one hand, we show that for particular tree shapes, the two rankings always coincide. On the other hand, we show that for all leaf numbers greater than or equal to five, there exist species trees for which the two rankings differ. Finally, we illustrate the variability in the rankings obtained from the FP index across different gene tree and species tree estimates for an empirical multilocus mammal data set.

Different Traits, Different Evolutionary Pathways: Insights from Salamandrina (Amphibia, Caudata)

Species delimitation is often based on a single or very few genetic or phenetic traits, something which leads to misinterpretations and often does not provide information about evolutionary processes. Here, we investigated the diversity pattern of multiple phenetic traits of the two extant species of Salamandrina, a genus split only after molecular traits had been studied but the two species of which are phenetically very similar. The phenetic traits we studied are size, external body shape and head colour pattern, in a model comparison framework using non-linear mixed models and unsupervised and supervised clustering. Overall, we found high levels of intra-specific variability for body size and shape, depending on population belonging and habitat, while differences between species were generally lower. The habitat the salamanders dwell in also seems important for colour pattern. Basing on our findings, from the methodological point of view, we suggest (i) to take into account the variability at population level when testing for higher level variability, and (ii) a semi-supervised learning approach to high dimensional data. We also showed that different phenotypic traits of the same organism could result from different evolutionary routes. Local adaptation is likely responsible for body size and shape variability, with selective pressures more similar across species than within them. Head colour pattern also depends on habitat, differently from ventral colour pattern (not studied in this paper) which likely evolved under genetic drift.

Inference of Gene Flow between Species under Misspecified Models

Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species.

Microevolution and subspecific taxonomy of Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, is a highly polymorphic species, subdivided into 6 main evolutionary lineages or near-clades (formerly discrete typing units or DTUs). An additional near-clade (TC-bat) has recently been evidenced. This pattern is considered to be the result of predominant clonal evolution (PCE). PCE is compatible with occasional mating/hybridization, which do not break the prevalent pattern of clonal evolution, the main trait of it being the presence of Multigene Bifurcating Trees (MGBTs) at all evolutionary levels ("clonal frame"). The development of highly resolutive genetic (microsatellites*) and genomic (sequencing and multi-single nucleotide polymorphism {SNP}* typing) markers shows that PCE also operates at a microevolutionary* level within each of the near-clades ("Russian doll pattern"), in spite of occasional meiosis and hybridization events. Within each near-clade, one can evidence widespread clonal multilocus genotypes*, linkage disequilibrium*, Multigene Bifurcating Trees and lesser near-clades. The within near-clade population structure is like a miniature picture of that of the whole species, suggesting gradual rather than saltatory evolution. Additional data are required to evaluate the stability of these lesser near-clades in the long run and to evaluate the need for an adequate nomenclature for this microevolutionary level.

Unveiling new perspective of phylogeography, genetic diversity, and population dynamics of Southeast Asian and Pacific chickens

The complex geographic and temporal origins of chicken domestication have attracted wide interest in molecular phylogeny and phylogeographic studies as they continue to be debated up to this day. In particular, the population dynamics and lineage-specific divergence time estimates of chickens in Southeast Asia (SEA) and the Pacific region are not well studied. Here, we analyzed 519 complete mitochondrial DNA control region sequences and identified 133 haplotypes with 70 variable sites. We documented 82.7% geographically unique haplotypes distributed across major haplogroups except for haplogroup C, suggesting high polymorphism among studied individuals. Mainland SEA (MSEA) chickens have higher overall genetic diversity than island SEA (ISEA) chickens. Phylogenetic trees and median-joining network revealed evidence of a new divergent matrilineage (i.e., haplogroup V) as a sister-clade of haplogroup C. The maximum clade credibility tree estimated the earlier coalescence age of ancestral D-lineage (i.e., sub-haplogroup D2) of continental chickens (3.7 kya; 95% HPD 1985-4835 years) while island populations diverged later at 2.1 kya (95% HPD 1467-2815 years). This evidence of earlier coalescence age of haplogroup D ancestral matriline exemplified dispersal patterns to the ISEA, and thereafter the island clade diversified as a distinct group.

Bayesian phylogenetic inference using relaxed-clocks and the multispecies coalescent

The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.

Impact of Ghost Introgression on Coalescent-based Species Tree Inference and Estimation of Divergence Time

The species studied in any evolutionary investigation generally constitute a small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," i.e., unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has rarely been studied and is poorly understood. Here, we use mathematical analysis and simulations to examine the robustness of species tree methods based on the multispecies coalescent model to introgression from a ghost or extant lineage. We found that many results originally obtained for introgression between extant species can easily be extended to ghost introgression, such as the strongly interactive effects of incomplete lineage sorting (ILS) and introgression on the occurrence of anomalous gene trees (AGTs). The relative performance of the summary species tree method (ASTRAL) and the full-likelihood method (*BEAST) varies under different introgression scenarios, with the former being more robust to gene flow between non-sister species whereas the latter performing better under certain conditions of ghost introgression. When an outgroup ghost (defined as a lineage that diverged before the most basal species under investigation) acts as the donor of the introgressed genes, the time of root divergence among the investigated species generally was overestimated, whereas ingroup introgression, as commonly perceived, can only lead to underestimation. In many cases of ingroup introgression that may or may not involve ghost lineages, the stronger the ILS, the higher the accuracy achieved in estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression.

Conservation genomics of urban populations of Streamside Salamander (Ambystoma barbouri)

In Tennessee, populations of the state endangered Streamside Salamander (Ambystoma barbouri) are in decline as their distribution lies mostly within rapidly developing areas in the Nashville Basin. Information regarding the partitioning of genetic variation among populations of A. barbouri and the taxonomic status of these populations relative to northern populations and their congener, the Small-mouthed Salamander (A. texanum), have important implications for management and conservation of this species. Here we combined mitochondrial sequencing and genome-wide single nucleotide polymorphism (SNP) data generated using Genotyping-by-Sequencing (GBS) to investigate patterns of genetic variation within Tennessee populations of A. barbouri, to assess their relationship to populations in Kentucky, and to examine their phylogenetic relationship to the closely related A. texanum. Results from phylogenetic reconstructions reveal a complex history of Tennessee A. barbouri populations with regards to northern populations, unisexual A. barbouri, and A. texanum. Patterns of mitochondrial sequence variation suggest that A. barbouri may have originated within Tennessee and expanded north multiple times into Kentucky, Ohio, Indiana, and West Virginia. Phylogenetic reconstructions based on genome-wide SNP data contradict results based on mitochondrial DNA and correspond to geographic and taxonomic boundaries. Variation in allele frequencies at SNP genotypes, as identified by multivariate analyses and Bayesian assignment tests, identified three evolutionary significant units (ESUs) for A. barbouri within Tennessee. Collectively, these results emphasize the need for prioritizing conservation needs for Tennessee populations of A. barbouri to ensure the long-term persistence of this species.

Genomic insights into zokors' phylogeny and speciation in China

SignificancePhylogenies are the basis of many ecological and evolutionary studies. However, zokor phylogeny and speciation patterns are heavily debated. This study disentangled the phylogeny and speciation of zokors genomically. Six species of the Eospalax were separated into high-altitude E. baileyi and E. smithi and the rest four low-altitude species by recent Qinghai-Tibet Plateau uplift 3.6 million y ago. E. rothschildi and E. smithi speciated south of the Qinling-Huaihe Line, where refuges were supplied during glaciation. Introgression and incomplete lineage sorting led to the complex phylogeny of zokors. Genomic islands were formed due to ancient polymorphisms and divergence hitchhiking. This study concluded that climatic, geological, and tectonic events shaped the phylogeny and speciation of zokors in China.

Two new glassfrogs (Centrolenidae: Hyalinobatrachium) from Ecuador, with comments on the endangered biodiversity of the Andes

Background

The Tropical Andes is the world's most biodiverse hotspot. This region contains >1,000 amphibian species, more than half of which are endemic. Herein we describe two new glassfrog species (Centrolenidae: Hyalinobatrachium) that we discovered within relatively unexplored and isolated localities of the Ecuadorian Andes.

Methods

We employed morphological, acoustic, and molecular methods to test the hypothesis that Hyalinobatrachium mashpi sp. nov and H. nouns sp. nov. are species new to science. Following standard methods, we generated mitochondrial sequences (16S) of 37 individuals in the genus Hyalinobatrachium. We inferred the phylogenetic relationships of the two new species in comparison to all other glassfrogs using Maximum Likelihood. In addition to describing the call of H. mashpi sp. nov., we performed a discriminant analysis of principal components (DAPC) with the advertisement call characteristics of several congeners.

Results

Based on an integrative taxonomy approach, we describe two new species. Morphological traits and the inferred phylogeny unambiguously place the new taxa in the genus Hyalinobatrachium. Both species are distinguished from other glassfrogs mainly by their dorsal coloration (i.e., dorsum lime green with small light yellow spots, head usually with interorbital bar) and transparent pericardium (i.e., the heart is visible through the ventral skin). The new species exhibit a high morphological similarity (i.e., cryptic) and occur within relatively close geographical proximity (closest aerial distance = 18.9 km); however, their uncorrected p distance for the mitochondrial gene 16S is 4.6-4.7%, a value that greatly exceeds the genetic distance between closely related species of centrolenid frogs. The DAPC revealed that the advertisement call of H. mashpi sp. nov. is acoustically distinct.

Discussion

Our findings are congruent with several previous studies that report a high degree of endemism in the Toisán mountain range, which appears to be isolated from the main Andean cordillera for some amphibian groups. We recommend that both H. mashpi sp. nov. and H. nouns sp. nov. be listed as Endangered, following IUCN criteria. These new species provide another example of cryptic diversity in the Andes-further evidence that the region fosters much more biodiversity than we have the resources to catalog. Threatened by mining and other exploitative industries, these glassfrogs and many other yet-to-be-discovered Andean species highlight the dire need for effective conservation measures-especially in northwestern Ecuador.

Multispecies coalescent and its applications to infer species phylogenies and cross-species gene flow

Multispecies coalescent (MSC) is the extension of the single-population coalescent model to multiple species. It integrates the phylogenetic process of species divergences and the population genetic process of coalescent, and provides a powerful framework for a number of inference problems using genomic sequence data from multiple species, including estimation of species divergence times and population sizes, estimation of species trees accommodating discordant gene trees, inference of cross-species gene flow and species delimitation. In this review, we introduce the major features of the MSC model, discuss full-likelihood and heuristic methods of species tree estimation and summarize recent methodological advances in inference of cross-species gene flow. We discuss the statistical and computational challenges in the field and research directions where breakthroughs may be likely in the next few years.

Species limits and phylogeographic structure in two genera of solitary African mole-rats Georychus and Heliophobius

African mole-rats (Bathyergidae) are an intensively studied family of subterranean rodents including three highly social and three solitary genera. Although their phylogenetic interrelations are clear, genetic diversity and the number of species within each genus is much less certain. Among the solitary genera, Heliophobius and Georychus were for a long time considered as monotypic, but molecular studies demonstrated strong phylogeographic structure within each genus and proposed that they represent complexes of cryptic species. The present study re-evaluates their internal genetic/phylogenetic structure using a combination of methodological approaches. We generated datasets of one mitochondrial and six specifically selected nuclear markers as well as of a large number of double digest restriction site associated (ddRAD) loci and then applied species delimitation analyses based on the multispecies coalescent model or clustering on co-ancestry matrices. The population structure was largely congruent across all analyses, but the methods differed in their resolution scale when determining distinct gene pools. While the multispecies coalescent model distinguished five Georychus and between eleven to thirteen Heliophobius gene pools in both Sanger sequenced and ddRAD loci, two clustering algorithms revealed significantly finer or coarser structure in ddRAD based co-ancestry matrices. Tens of clusters were distinguished by fineRADstructure and one (in Georychus) or two clusters (in Heliophobius) by Infomap. The divergence dating of the bathyergid phylogeny estimated that diversification within both genera coincided with the onset of the Pleistocene and was likely driven by repeated large-scale climatic changes. Based on this updated genetic evidence, we suggest recognizing one species of Georychus and two species of Heliophobius, corresponding to a northern and southern major lineage, separated by the Eastern Arc Mountains. Yet, the final taxonomic revision should await integrated evidence stemming from e.g.. morphological, ecological, or behavioral datasets.

Phylogenetic signal in the vocalizations of vocal learning and vocal non-learning birds

Some animal vocalizations develop reliably in the absence of relevant experience, but an intriguing subset of animal vocalizations is learned: they require acoustic models during ontogeny in order to develop, and the learner's vocal output reflects those models. To what extent do such learned vocalizations reflect phylogeny? We compared the degree to which phylogenetic signal is present in vocal signals from a wide taxonomic range of birds, including both vocal learners (songbirds) and vocal non-learners. We used publically available molecular phylogenies and developed methods to analyse spectral and temporal features in a carefully curated collection of high-quality recordings of bird songs and bird calls, to yield acoustic distance measures. Our methods were initially developed using pairs of closely related North American and European bird species, and then applied to a non-overlapping random stratified sample of European birds. We found strong similarity in acoustic and genetic distances, which manifested itself as a significant phylogenetic signal, in both samples. In songbirds, both learned song and (mostly) unlearned calls allowed reconstruction of phylogenetic trees nearly isomorphic to the phylogenetic trees derived from genetic analysis. We conclude that phylogeny and inheritance constrain vocal structure to a surprising degree, even in learned birdsong. This article is part of the theme issue 'Vocal learning in animals and humans'.

Homoploid hybrids, allopolyploids, and high ploidy levels characterize the evolutionary history of a western North American quillwort (Isoëtes) complex

Polyploidy and hybridization are important processes in seed-free plant evolution. However, a major gap lies in our understanding of how these processes affect the evolutionary history of high-ploidy systems. The heterosporous lycophyte genus Isoëtes is a lineage with many putative hybrids and high-level polyploid taxa (ranging from tetraploid to dodecaploid). Here, we use a complex of western North American Isoëtes, to understand the role of hybridization and high-level polyploidy in generating and maintaining novel diversity. To uncover these processes, we use restriction-site associated DNA sequencing (RADseq), multiple alleles of a single low-copy nuclear marker, whole plastomes, cytology (genome size estimates and chromosome counts), and reproductive status (fertile or sterile). With this dataset, we show that hybridization occurs easily between species in this complex and is bidirectional between identical, but not different, cytotypes. Furthermore, we show that fertile allopolyploids appear to have formed repeatedly from sterile homoploid and interploid hybrids. We propose that low prezygotic reproductive barriers and a high frequency of whole-genome duplication allow for high-level polyploid systems to generate novel lineages, and that these mechanisms may be important in shaping extant Isoëtes diversity.

Most Genomic Loci Misrepresent the Phylogeny of an Avian Radiation Because of Ancient Gene Flow

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis, and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the "winner-takes-all democratic majority tree" is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger data sets, especially whole-genome data sets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.[Bifurcation; introgression; recombination; reticulation; Phylloscopus.].

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

DNA barcodes combined with multi-locus data of representative taxa can generate reliable higher-level phylogenies

Taxa are frequently labeled incertae sedis when their placement is debated at ranks above the species level, such as their subgeneric, generic, or subtribal placement. This is a pervasive problem in groups with complex systematics due to difficulties in identifying suitable synapomorphies. In this study, we propose combining DNA barcodes with a multilocus backbone phylogeny in order to assign taxa to genus or other higher-level categories. This sampling strategy generates molecular matrices containing large amounts of missing data that are not distributed randomly: barcodes are sampled for all representatives, and additional markers are sampled only for a small percentage. We investigate the effects of the degree and randomness of missing data on phylogenetic accuracy using simulations for up to 100 markers in 1000-tips trees, as well as a real case: the subtribe Polyommatina (Lepidoptera: Lycaenidae), a large group including numerous species with unresolved taxonomy. Our simulation tests show that when a strategic and representative selection of species for higher-level categories has been made for multigene sequencing (approximately one per simulated genus), the addition of this multigene backbone DNA data for as few as 5–10% of the specimens in the total data set can produce high-quality phylogenies, comparable to those resulting from 100% multigene sampling. In contrast, trees based exclusively on barcodes performed poorly. This approach was applied to a 1365-specimen data set of Polyommatina (including ca. 80% of described species), with nearly 8% of representative species included in the multigene backbone and the remaining 92% included only by mitochondrial COI barcodes, a phylogeny was generated that highlighted potential misplacements, unrecognized major clades, and placement for incertae sedis taxa. We use this information to make systematic rearrangements within Polyommatina, and to describe two new genera. Finally, we propose a systematic workflow to assess higher-level taxonomy in hyperdiverse groups. This research identifies an additional, enhanced value of DNA barcodes for improvements in higher-level systematics using large data sets. [Birabiro; DNA barcoding; incertae sedis; Kipepeo; Lycaenidae; missing data; phylogenomic; phylogeny; Polyommatina; supermatrix; systematics; taxonomy]

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

A refined proposal for the origin of dogs: the case study of Gnirshöhle, a Magdalenian cave site

Dogs are known to be the oldest animals domesticated by humans. Although many studies have examined wolf domestication, the geographic and temporal origin of this process is still being debated. To address this issue, our study sheds new light on the early stages of wolf domestication during the Magdalenian period (16–14 ka cal BP) in the Hegau Jura region (Southwestern Germany and Switzerland). By combining morphology, genetics, and isotopes, our multidisciplinary approach helps to evaluate alternate processes driving the early phases of domestication. The isotope analysis uncovered a restricted, low δ¹⁵N protein diet for all analyzed Gnirshöhle specimens, while morphological examinations and phylogenetic relationships did not unequivocally assign them to one or the other canid lineage. Intriguingly, the newly generated mitochondrial canid genomes span the entire genetic diversity of modern dogs and wolves. Such high mitochondrial diversity could imply that Magdalenian people tamed and reared animals originating from different wolf lineages. We discuss our results in light of three ecological hypotheses and conclude that both domestication and the existence of a specialized wolf ecomorph are highly probable. However, due to their proximity to humans and a restricted diet, we propose domestication as the most likely scenario explaining the patterns observed herein.

Ancient Antibiotics, Ancient Resistance

As the spread of antibiotic resistance threatens our ability to treat infections, avoiding the return of a preantibiotic era requires the discovery of new drugs. While therapeutic use of antibiotics followed by the inevitable selection of resistance is a modern phenomenon, these molecules and the genetic determinants of resistance were in use by environmental microbes long before humans discovered them. In this review, we discuss evidence that antibiotics and resistance were present in the environment before anthropogenic use, describing techniques including direct sampling of ancient DNA and phylogenetic analyses that are used to reconstruct the past. We also pay special attention to the ecological and evolutionary forces that have shaped the natural history of antibiotic biosynthesis, including a discussion of competitive versus signaling roles for antibiotics, proto-resistance, and substrate promiscuity of biosynthetic and resistance enzymes. Finally, by applying an evolutionary lens, we describe concepts governing the origins and evolution of biosynthetic gene clusters and cluster-associated resistance determinants. These insights into microbes' use of antibiotics in nature, a game they have been playing for millennia, can provide inspiration for discovery technologies and management strategies to combat the growing resistance crisis.

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe)

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson's D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.

Complexity of the simplest species tree problem

The multispecies coalescent model provides a natural framework for species tree estimation accounting for gene-tree conflicts. Although a number of species tree methods under the multispecies coalescent have been suggested and evaluated using simulation, their statistical properties remain poorly understood. Here, we use mathematical analysis aided by computer simulation to examine the identifiability, consistency, and efficiency of different species tree methods in the case of three species and three sequences under the molecular clock. We consider four major species-tree methods including concatenation, two-step, independent-sites maximum likelihood, and maximum likelihood. We develop approximations that predict that the probit transform of the species tree estimation error decreases linearly with the square root of the number of loci. Even in this simplest case, major differences exist among the methods. Full-likelihood methods are considerably more efficient than summary methods such as concatenation and two-step. They also provide estimates of important parameters such as species divergence times and ancestral population sizes,whereas these parameters are not identifiable by summary methods. Our results highlight the need to improve the statistical efficiency of summary methods and the computational efficiency of full likelihood methods of species tree estimation.

Molecular systematics of Thouarella (Octocorallia:Primnoidae) with the description of three new species from the Southern Ocean based on combined molecular and morphological evidence

Thouarella Gray, 1870, is one of the most speciose genera among gorgonians of the family Primnoidae (Cnidaria:Octocorallia:Anthozoa), being remarkably diverse in the Antarctic and sub-Antarctic seafloor. However, their diversity in the Southern Ocean is likely underestimated. Phylogenetic analyses of mitochondrial and nuclear DNA markers were integrated with species delimitation approaches as well as morphological colonial and polyps features and skeletal SEM examinations to describe and illustrate three new species within Thouarella, from the Weddell Sea, Southern Ocean: T. amundseni sp. nov., T. dolichoespinosa sp. nov. and T. pseudoislai sp. nov. Our species delimitation results suggest, for the first time, the potential presence of Antarctic and sub-Antarctic cryptic species of primnoids, based on the likely presence of sibling species within T. undulata and T. crenelata. With the three new species here described, the global diversity of Thouarella has increased to 41 species, 15 of which are endemic to the Antarctic and sub-Antarctic waters. Consequently, our results provide new steps for uncovering the shelf benthonic macrofauna’s hidden diversity in the Southern Ocean. Finally, we recommend using an integrative taxonomic framework in this group of organisms and species delimitation approaches because the distinctions between some Thouarella species based only on a superficial examination of their macro- and micromorphological features is, in many cases, limited.

Response to Comment on "Ancient origins of allosteric activation in a Ser-Thr kinase"

Park et al question one out of seven findings from Hadzipasic et al: whether TPX2 allosterically regulates the oldest Aurora. We had already addressed the two concerns raised-sparse sequence sampling and not forcing the gene to the species tree-before publication. Moreover, we believe their ancestral sequence reconstruction would be consistent with a nonallosteric common ancestor, and we show large sequence differences caused by species tree-enforced gene trees.

Molecular systematics and phylogeography of a widespread Neotropical avian lineage: evidence for cryptic speciation with protracted gene flow throughout the Late Quaternary

Here we use an integrative approach, including coalescent-based methods, isolation–migration and species distribution models, to infer population structure, divergence times and diversification in the two species of the genus Cymbilaimus (Aves, Thamnophilidae). Our results support a recent and rapid diversification with both incomplete lineage sorting and gene flow shaping the evolutionary history of Cymbilaimus. The spatio-temporal pattern of cladogenesis suggests that Cymbilaimus originated in the north/western portion of cis-Andean South America and then diversified into the Brazilian Shield and Central America after consolidation of the modern Amazonian drainage and the Andean range. This evolutionary scenario is explained by cycles of range expansion and dispersal, followed by isolation, and recurrent gene flow, during the last 1.2 Myr. Our results agree with those recently reported for other closely related suboscine lineages, whereby the window of introgression between closely related taxa remains open for up to a few million years after their original split. In Cymbilaimus, introgression was recurrent between C. lineatus and C. sanctaemariae, even after they acquired vocal and ecological differentiation, supporting the claim that at least in Neotropical suboscines, full reproductive compatibility may take millions of years to evolve and cannot be interpreted as synonymous with a lack of speciation.

ASAP: assemble species by automatic partitioning

Here, we describe Assemble Species by Automatic Partitioning (ASAP), a new method to build species partitions from single locus sequence alignments (i.e., barcode data sets). ASAP is efficient enough to split data sets as large 10⁴ sequences into putative species in several minutes. Although grounded in evolutionary theory, ASAP is the implementation of a hierarchical clustering algorithm that only uses pairwise genetic distances, avoiding the computational burden of phylogenetic reconstruction. Importantly, ASAP proposes species partitions ranked by a new scoring system that uses no biological prior insight of intraspecific diversity. ASAP is a stand-alone program that can be used either through a graphical web-interface or that can be downloaded and compiled for local usage. We have assessed its power along with three others programs (ABGD, PTP and GMYC) on 10 real COI barcode data sets representing various degrees of challenge (from small and easy cases to large and complicated data sets). We also used Monte-Carlo simulations of a multispecies coalescent framework to assess the strengths and weaknesses of ASAP and the other programs. Through these analyses, we demonstrate that ASAP has the potential to become a major tool for taxonomists as it proposes rapidly in a full graphical exploratory interface relevant species hypothesis as a first step of the integrative taxonomy process.

Identification of the taxonomic status of Scinax nebulosus and Scinax constrictus (Scinaxinae, Anura) based on molecular markers

The validation of many anuran species is based on a strictly descriptive, morphological analysis of a small number of specimens with a limited geographic distribution. The Scinax Wagler, 1830 genus is a controversial group with many doubtful taxa and taxonomic uncertainties, due a high number of cryptic species. One example is the pair of species Scinax constrictus and Scinax nebulosus, which share a similar morphology. Scinax constrictus is restricted to the Brazilian Cerrado savanna, while S. nebulosus is widely distributed throughout northern South America. Despite the validation of many anuran species, discriminations based only on morphological traits is quite difficult due to the high conservative morphology of some groups. In this context, the present study uses mitochondrial and nuclear genes to provide a more consistent diagnosis and test the validity of S. constrictus as a distinct species from S. nebulosus, as well as evaluate the position of these taxa within the Scinax genus. The topologies obtained herein uphold the monophyletic status of Scinax based on all molecular markers assessed in this study, in all analytical approaches, with high levels of statistical support.

Evolutionary Study of the Crassphage Virus at Gene Level

crAss-like viruses are a putative family of bacteriophages recently discovered. The eponym of the clade, crAssphage, is an enteric bacteriophage estimated to be present in at least half of the human population and it constitutes up to 90% of the sequences in some human fecal viral metagenomic datasets. We focused on the evolutionary dynamics of the genes encoded on the crAssphage genome. By investigating the conservation of the genes, a consistent variation in the evolutionary rates across the different functional groups was found. Gene duplications in crAss-like genomes were detected. By exploring the differences among the functional categories of the genes, we confirmed that the genes encoding capsid proteins were the most ubiquitous, despite their overall low sequence conservation. It was possible to identify a core of proteins whose evolutionary trees strongly correlate with each other, suggesting their genetic interaction. This group includes the capsid proteins, which are thus established as extremely suitable for rebuilding the phylogenetic tree of this viral clade. A negative correlation between the ubiquity and the conservation of viral protein sequences was shown. Together, this study provides an in-depth picture of the evolution of different genes in crAss-like viruses.