Likelihood-based tests of topologies in phylogenetics

Cryogenian Origins of Multicellularity in Archaeplastida

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

Performance of tree-building methods using a morphological dataset and a well-supported Hexapoda phylogeny

Recently, many studies have addressed the performance of phylogenetic tree-building methods (maximum parsimony, maximum likelihood, and Bayesian inference), focusing primarily on simulated data. However, for discrete morphological data, there is no consensus yet on which methods recover the phylogeny with better performance. To address this lack of consensus, we investigate the performance of different methods using an empirical dataset for hexapods as a model. As an empirical test of performance, we applied normalized indices to effectively measure accuracy (normalized Robinson-Foulds metric, nRF) and precision, which are measured via resolution, one minus Colless' consensus fork index (1-CFI). Additionally, to further explore phylogenetic accuracy and support measures, we calculated other statistics, such as the true positive rate (statistical power) and the false positive rate (type I error), and constructed receiver operating characteristic plots to visualize the relationship between these statistics. We applied the normalized indices to the reconstructed trees from the reanalyses of an empirical discrete morphological dataset from extant Hexapoda using a well-supported phylogenomic tree as a reference. Maximum likelihood and Bayesian inference applying the k-state Markov (Mk) model (without or with a discrete gamma distribution) performed better, showing higher precision (resolution). Additionally, our results suggest that most available tree topology tests are reliable estimators of the performance measures applied in this study. Thus, we suggest that likelihood-based methods and tree topology tests should be used more often in phylogenetic tree studies based on discrete morphological characters. Our study provides a fair indication that morphological datasets have robust phylogenetic signal.

Performance of Topology Tests under Extreme Selection Bias

Tree tests like the Kishino-Hasegawa (KH) test and chi-square test suffer a selection bias that tests like the Shimodaira-Hasegawa (SH) test and approximately unbiased test were intended to correct. We investigate tree-testing performance in the presence of severe selection bias. The SH test is found to be very conservative and, surprisingly, its uncorrected analog, the KH test has low Type I error even in the presence of extreme selection bias, leading to a recommendation that the SH test be abandoned. A chi-square test is found to usually behave well and but to require correction in extreme cases. We show how topology testing procedures can be used to get support values for splits and compare the likelihood-based support values to the approximate likelihood ratio test (aLRT) support values. We find that the aLRT support values are reasonable even in settings with severe selection bias that they were not designed for. We also show how they can be used to construct tests of topologies and, in doing so, point out a multiple comparisons issue that should be considered when looking at support values for splits.

A 50-million-year-old, three-dimensionally preserved bat skull supports an early origin for modern echolocation

Bats are among the most recognizable, numerous, and widespread of all mammals. But much of their fossil record is missing, and bat origins remain poorly understood, as do the relationships of early to modern bats. Here, we describe a new early Eocene bat that helps bridge the gap between archaic stem bats and the hyperdiverse modern bat radiation of more than 1,460 living species. Recovered from ∼50 million-year-old cave sediments in the Quercy Phosphorites of southwestern France, Vielasia sigei's remains include a near-complete, three-dimensionally preserved skull-the oldest uncrushed bat cranium yet found. Phylogenetic analyses of a 2,665 craniodental character matrix, with and without 36.8 kb of DNA sequence data, place Vielasia outside modern bats, with total evidence tip-dating placing it sister to the crown clade. Vielasia retains the archaic dentition and skeletal features typical of early Eocene bats, but its inner ear shows specializations found in modern echolocating bats. These features, which include a petrosal only loosely attached to the basicranium, an expanded cochlea representing ∼25% basicranial width, and a long basilar membrane, collectively suggest that the kind of laryngeal echolocation used by most modern bats predates the crown radiation. At least 23 individuals of V. sigei are preserved together in a limestone cave deposit, indicating that cave roosting behavior had evolved in bats by the end of the early Eocene; this period saw the beginning of significant global climate cooling that may have been an evolutionary driver for bats to first congregate in caves.

A molecular phylogeny for all 21 families within Chiroptera (bats)

Bats, members of the Chiroptera order, rank as the second most diverse group among mammals. Recent molecular systematic studies on bats have successfully classified 21 families within two suborders: Yinpterochiroptera and Yangochiroptera. Nevertheless, the phylogeny within these 21 families has remained a subject of controversy. In this study, we have employed a balanced approach to establish a robust family-level phylogenetic hypothesis for bats, utilizing a more comprehensive molecular dataset. This dataset includes representative species from all 21 bat families, resulting in a reduced level of missing genetic information. The resulting phylogenetic tree comprises 21 lineages that are strongly supported, each corresponding to one of the bat families. Our findings support to place the Emballonuroidea superfamily as the basal lineage of Yangochiroptera, and that Myzopodidae should be situated as a basal lineage of Emballonuroidea, forming a sister relationship with the clade consisting of Nycteridae and Emballonuridae. Finally, we have conducted dating analyses on this newly resolved phylogenetic tree, providing divergence times for each bat family. Collectively, our study has employed a relatively comprehensive molecular dataset to establish a more robust phylogeny encompassing all 21 bat families. This improved phylogenetic framework will significantly contribute to our understanding of evolutionary processes, ecological roles, disease dynamics, and biodiversity conservation in the realm of bats.

A probable koala from the Oligocene of central Australia provides insights into early diprotodontian evolution

Diprotodontians are the morphologically and ecologically most diverse order of marsupials. However, an approximately 30-million-year gap in the Australian terrestrial vertebrate fossil record means that the first half of diprotodontian evolution is unknown. Fossil taxa from immediately either side of this gap are therefore critical for reconstructing the early evolution of the order. Here we report the likely oldest-known koala relatives (Phascolarctidae), from the late Oligocene Pwerte Marnte Marnte Local Fauna (central Australia). These include coeval species of Madakoala and Nimiokoala, as well as a new probable koala (?Phascolarctidae). The new taxon, Lumakoala blackae gen. et sp. nov., was comparable in size to the smallest-known phascolarctids, with body-mass estimates of 2.2-2.6 kg. Its bunoselenodont upper molars retain the primitive metatherian condition of a continuous centrocrista, and distinct stylar cusps B and D which lacked occlusion with the hypoconid. This structural arrangement: (1) suggests a morphocline within Phascolarctidae from bunoselenodonty to selenodonty; and (2) better clarifies the evolutionary transitions between molar morphologies within Vombatomorphia. We hypothesize that the molar form of Lumakoala blackae approximates the ancestral condition of the suborder Vombatiformes. Furthermore, it provides a plausible link between diprotodontians and the putative polydolopimorphians Chulpasia jimthorselli and Thylacotinga bartholomaii from the early Eocene Tingamarra Local Fauna (eastern Australia), which we infer as having molar morphologies consistent with stem diprotodontians.

CRISPR Comparison Toolkit: Rapid Identification, Visualization, and Analysis of CRISPR Array Diversity

CRISPR-Cas systems provide immunity against mobile genetic elements (MGEs) through sequence-specific targeting by spacer sequences encoded in CRISPR arrays. Spacers are highly variable between microbial strains and can be acquired rapidly, making them well suited for use in strain typing of closely related organisms. However, no tools are currently available to automate the process of reconstructing strain histories using CRISPR spacers. We therefore developed the CRISPR Comparison Toolkit (CCTK) to enable analyses of array relationships. The CCTK includes tools to identify arrays, analyze relationships between arrays using CRISPRdiff and CRISPRtree, and predict targets of spacers. CRISPRdiff visualizes arrays and highlights the similarities between them. CRISPRtree infers a phylogenetic tree from array relationships and presents a hypothesis of the evolutionary history of the arrays. The CCTK unifies several CRISPR analysis tools into a single command line application, including the first tool to infer phylogenies from array relationships.

Basal Anseriformes from the Early Paleogene of North America and Europe

We describe nearly complete skeletons of basal Anseriformes from the Latest Paleocene to the early Eocene of North America and Europe. Collectively, these birds appear to be representative of anseriforms near the divergence of Anhimae and Anseres, but their exact positions relative to these clades remains uncertain. A new family, Anachronornithidae nov. fam., is erected on the basis of one of these, Anachronornis anhimops nov. gen., nov. gen. et sp., to which the others cannot be confidently assigned. The new fossils augment a growing collection of early Pan-Anseriformes, which in their diversity do not paint an unambiguous picture of phylogeny or character state evolution on the path to or within crown-Anseriformes. Anachronornis nov. gen. is similar in some aspects of both cranial and postcranial anatomy to other well-represented early Paleogene Anseriformes and members of Anseres, such as Presbyornis Wetmore, 1926. However, it exhibits a more landfowl-like bill, like that of Anhimae and unlike the spatulate bill of Anseres. Additional specimens of similar basal Anseriformes of uncertain affinities from the early Eocene of North America and Europe further complicate interpretation of character state polarity due to the mosaicism of primitive and derived characters they exhibit.

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods

Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

Distance-Based Phylogenetic Placement with Statistical Support

Simple Summary

Phylogenetic placement seeks to find the optimal position for a new query species on an existing backbone tree. Fast and accurate distance-based phylogenetic placement methods lack the crucial feature of estimating the support values for various placements of a query sequence. This study presents both parametric and nonparametric methods for measuring the support values of distance-based phylogenetic placements.

Abstract

Phylogenetic identification of unknown sequences by placing them on a tree is routinely attempted in modern ecological studies. Such placements are often obtained from incomplete and noisy data, making it essential to augment the results with some notion of uncertainty. While the standard likelihood-based methods designed for placement naturally provide such measures of uncertainty, the newer and more scalable distance-based methods lack this crucial feature. Here, we adopt several parametric and nonparametric sampling methods for measuring the support of phylogenetic placements that have been obtained with the use of distances. Comparing the alternative strategies, we conclude that nonparametric bootstrapping is more accurate than the alternatives. We go on to show how bootstrapping can be performed efficiently using a linear algebraic formulation that makes it up to 30 times faster and implement this optimized version as part of the distance-based placement software APPLES. By examining a wide range of applications, we show that the relative accuracy of maximum likelihood (ML) support values as compared to distance-based methods depends on the application and the dataset. ML is advantageous for fragmentary queries, while distance-based support values are more accurate for full-length and multi-gene datasets. With the quantification of uncertainty, our work fills a crucial gap that prevents the broader adoption of distance-based placement tools.

Gene expression phylogenies and ancestral transcriptome reconstruction resolves major transitions in the origins of pregnancy

Structural and physiological changes in the female reproductive system underlie the origins of pregnancy in multiple vertebrate lineages. In mammals, the glandular portion of the lower reproductive tract has transformed into a structure specialized for supporting fetal development. These specializations range from relatively simple maternal nutrient provisioning in egg-laying monotremes to an elaborate suite of traits that support intimate maternal-fetal interactions in Eutherians. Among these traits are the maternal decidua and fetal component of the placenta, but there is considerable uncertainty about how these structures evolved. Previously, we showed that changes in uterine gene expression contributes to several evolutionary innovations during the origins of pregnancy (Mika et al., 2021b). Here, we reconstruct the evolution of entire transcriptomes ('ancestral transcriptome reconstruction') and show that maternal gene expression profiles are correlated with degree of placental invasion. These results indicate that an epitheliochorial-like placenta evolved early in the mammalian stem-lineage and that the ancestor of Eutherians had a hemochorial placenta, and suggest maternal control of placental invasiveness. These data resolve major transitions in the evolution of pregnancy and indicate that ancestral transcriptome reconstruction can be used to study the function of ancestral cell, tissue, and organ systems.

Widespread Interspecific Phylogenetic Tree Incongruence Between Mosquito-Borne and Insect-Specific Flaviviruses at Hotspots Originally Identified in Zika Virus

Intraspecies (homologous) phylogenetic incongruence, or ‘tree conflict’ between different loci within the same genome of mosquito-borne flaviviruses (MBFV), was first identified in dengue virus (DENV) and subsequently in Japanese encephalitis virus (JEV), St Louis encephalitis virus, and Zika virus (ZIKV). Recently, the first evidence of phylogenetic incongruence between interspecific members of the MBFV was reported in ZIKV and its close relative, Spondweni virus. Uniquely, these hybrid proteomes were derived from four incongruent trees involving an Aedes-associated DENV node (1 tree) and three different Culex-associated flavivirus nodes (3 trees). This analysis has now been extended across a wider spectrum of viruses within the MBFV lineage targeting the breakpoints between phylogenetic incongruent loci originally identified in ZIKV. Interspecies phylogenetic incongruence at these breakpoints was identified in 10 of 50 viruses within the MBFV lineage, representing emergent Aedes and Culex-associated viruses including JEV, West Nile virus, yellow fever virus, and insect-specific viruses. Thus, interspecies phylogenetic incongruence is widespread amongst the flaviviruses and is robustly associated with the specific breakpoints that coincide with the interspecific phylogenetic incongruence previously identified, inferring they are ‘hotspots’. The incongruence amongst the emergent MBFV group was restricted to viruses within their respective associated epidemiological boundaries. This MBFV group was RY-coded at the third codon position (‘wobble codon’) to remove transition saturation. The resulting ‘wobble codon’ trees presented a single topology for the entire genome that lacked any robust evidence of phylogenetic incongruence between loci. Phylogenetic interspecific incongruence was therefore observed for exactly the same loci between amino acid and the RY-coded ‘wobble codon’ alignments and this incongruence represented either a major part, or the entire genomes. Maximum likelihood codon analysis revealed positive selection for the incongruent lineages. Positive selection could result in the same locus producing two opposing trees. These analyses for the clinically important MBFV suggest that robust interspecific phylogenetic incongruence resulted from amino acid selection. Convergent or parallel evolutions are evolutionary processes that would explain the observation, whilst interspecific recombination is unlikely.

Phylogenomic analysis of evolutionary relationships in Ranitomeya poison frogs (Family Dendrobatidae) using ultraconserved elements

The use of genome-scale data in phylogenetics has enabled recent strides in determining the relationships between taxa that are taxonomically problematic because of extensive morphological variation. Here, we employ a phylogenomic approach to infer evolutionary relationships within Ranitomeya (Anura: Dendrobatidae), an Amazonian lineage of poison frogs consisting of 16 species with remarkable diversity in color pattern, range size, and parental care behavior. We infer phylogenies with all described species of Ranitomeya from ultraconserved nuclear genomic elements (UCEs) and also estimate divergence times. Our results differ from previous analyses regarding interspecific relationships. Notably, we find that R. toraro and R. defleri are not sister species but rather distantly related, contrary to previous analyses based on smaller genetic datasets. We recover R. uakarii as paraphyletic, designate certain populations formerly assigned to R. fantastica from Peru as R. summersi, and transfer the French Guianan and eastern Brazilian R. amazonica populations to R. variabilis. By clarifying both inter- and intraspecific relationships within Ranitomeya, our study paves the way for future tests of hypotheses on color pattern evolution and historical biogeography.

Phylogenetic Analysis of SARS-CoV-2 Data Is Difficult

Numerous studies covering some aspects of SARS-CoV-2 data analyses are being published on a daily basis, including a regularly updated phylogeny on nextstrain.org. Here, we review the difficulties of inferring reliable phylogenies by example of a data snapshot comprising a quality-filtered subset of 8,736 out of all 16,453 virus sequences available on May 5, 2020 from gisaid.org. We find that it is difficult to infer a reliable phylogeny on these data due to the large number of sequences in conjunction with the low number of mutations. We further find that rooting the inferred phylogeny with some degree of confidence either via the bat and pangolin outgroups or by applying novel computational methods on the ingroup phylogeny does not appear to be credible. Finally, an automatic classification of the current sequences into subclasses using the mPTP tool for molecular species delimitation is also, as might be expected, not possible, as the sequences are too closely related. We conclude that, although the application of phylogenetic methods to disentangle the evolution and spread of COVID-19 provides some insight, results of phylogenetic analyses, in particular those conducted under the default settings of current phylogenetic inference tools, as well as downstream analyses on the inferred phylogenies, should be considered and interpreted with extreme caution.

Genome-based local dynamics of canine rabies virus epidemiology, transmission, and evolution in Davao City, Philippines, 2018-2019

Rabies is a fatal zoonotic and neglected tropical disease caused by the rabies virus (RABV) and is associated with neuronal dysfunction and death, with dogs as the predominant carrier. The Philippines plans to eradicate rabies by 2022, but this is challenged with sub-optimal coverage of vaccination programs coupled with sustained transmission chains, making it unable to eradicate the disease. We investigated the dynamics of canine rabies in the highly urbanized Davao City of the Philippines and its neighboring localities by assessing genetic relationships, transmission patterns, selection pressure, and recombination events using the whole genome sequence of 49 RABV cases from June 2018 to May 2019, majority of which (46%) were from the district of Talomo, Davao City. Although phylogeographic clustering was observed, local variants also exhibited genetic sub-lineages. Phylogenetic and spatial transmission analysis provided evidence for intra- and inter-city transmission predominantly through the Talomo district of Davao City. Around 84% of the cases were owned dogs, but the genetic similiarity of RABVs from stray and owned dogs further alluded to the role of the former as transmission vectors. The high rate of improper vaccination among the affected dogs (80%) was also a likely contributor to transmission. The RABV population under Investigation is generally under strong purifying selection with no evidence of vaccine evasion due to the genetic homogeneity of viruses from vaccinated and improperly vaccinated dogs. However, some homologous recombination (HR) events were identified along the G and L genes, also predominantly associated with viruses from Talomo. The complementary findings on epidemiology, transmission, and recombination for Talomo suggest that high incidence areas can be seeds for virus dispersal and evolution. We recommend further Investigations on the possibility of HR in future large-scale genome studies. Finally, districts associated with these phenomena can be targeted for evidence-based local strategies that can help break RABV transmission chains and prevent emergence of novel strains in Davao City.

Felsenstein Phylogenetic Likelihood

In 1981, the Journal of Molecular Evolution (JME) published an article entitled “Evolutionary trees from DNA sequences: A maximum likelihood approach” by Joseph (Joe) Felsenstein (J Mol Evol 17:368–376, 1981). This groundbreaking work laid the foundation for the emerging field of statistical phylogenetics, providing a tractable way of finding maximum likelihood (ML) estimates of evolutionary trees from DNA sequence data. This paper is the second most cited (more than 9000 citations) in JME after Kimura’s (J Mol Evol 16:111–120, 1980) seminal paper on a model of nucleotide substitution (with nearly 20,000 citations). On the occasion of the 50th anniversary of JME, we elaborate on the significance of Felsenstein’s ML approach to estimating phylogenetic trees.

Genetic diversity of North American vertebrates in protected areas

Protected areas play a crucial role in the conservation of biodiversity; however, it is unclear if these areas have an influence on genetic diversity. As a first step towards addressing this issue, we compare the genetic diversity inside and outside of protected areas. We tested the null hypothesis that there is no difference between genetic diversity inside compared to outside of protected areas in 44 vertebrate species. By automatically skimming the Global Biodiversity Information Facility (GBIF) and the National Center for Biotechnology Information (NCBI) GenBank we obtained genetic and geographical data to be repurposed and reanalysed. Novel pipelines were used to automate the process of assigning individuals to inside or outside of protected areas and then used to calculate different measures of intraspecific diversity. Forty-eight percent of examined species showed a significant difference in the amount of nucleotide diversity they contained inside compared to outside of protected areas, with similar numbers of species containing more or less genetic diversity inside compared to outside. Although our simulation testing suggests that this result is not an artefact of sampling, it is unclear what factors influence the relative amount of genetic diversity in protected areas across species.

The Location of Substitutions and Bacterial Genome Arrangements

Increasing evidence supports the notion that different regions of a genome have unique rates of molecular change. This variation is particularly evident in bacterial genomes where previous studies have reported gene expression and essentiality tend to decrease, whereas substitution rates usually increase with increasing distance from the origin of replication. Genomic reorganization such as rearrangements occur frequently in bacteria and allow for the introduction and restructuring of genetic content, creating gradients of molecular traits along genomes. Here, we explore the interplay of these phenomena by mapping substitutions to the genomes of Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti, quantifying how many substitutions have occurred at each position in the genome. Preceding work indicates that substitution rate significantly increases with distance from the origin. Using a larger sample size and accounting for genome rearrangements through ancestral reconstruction, our analysis demonstrates that the correlation between the number of substitutions and the distance from the origin of replication is significant but small and inconsistent in direction. Some replicons had a significantly decreasing trend (E. coli and the chromosome of S. meliloti), whereas others showed the opposite significant trend (B. subtilis, Streptomyces, pSymA and pSymB in S. meliloti). dN, dS, and ω were examined across all genes and there was no significant correlation between those values and distance from the origin. This study highlights the impact that genomic rearrangements and location have on molecular trends in some bacteria, illustrating the importance of considering spatial trends in molecular evolutionary analysis. Assuming that molecular trends are exclusively in one direction can be problematic.

Early vertebrate origin of CTCFL, a CTCF paralog, revealed by proximity-guided shark genome scaffolding

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in diverse animals. The gene encoding this protein has a paralog which was first identified to be expressed exclusively in the testis in mammals and designated as CTCFL (also called BORIS). CTCFL orthologs were reported only among amniotes, and thus CTCFL was once thought to have arisen in the amniote lineage. In this study, we identified elasmobranch CTCFL orthologs, and investigated its origin with the aid of a shark genome assembly improved by proximity-guided scaffolding. Our analysis employing evolutionary interpretation of syntenic gene location suggested an earlier timing of the gene duplication between CTCF and CTCFL than previously thought, that is, around the common ancestor of extant vertebrates. Also, our transcriptomic sequencing revealed a biased expression of the catshark CTCFL in the testis, suggesting the origin of the tissue-specific localization in mammals more than 400 million years ago. To understand the historical process of the functional consolidation of the long-standing chromatin regulator CTCF, its additional paralogs remaining in some of the descendant lineages for spatially restricted transcript distribution should be taken into consideration.

The influence of a priori grouping on inference of genetic clusters: simulation study and literature review of the DAPC method

Inference of genetic clusters is a key aim of population genetics, sparking development of numerous analytical methods. Within these, there is a conceptual divide between finding de novo structure versus assessment of a priori groups. Recently developed, Discriminant Analysis of Principal Components (DAPC), combines discriminant analysis (DA) with principal component (PC) analysis. When applying DAPC, the groups used in the DA (specified a priori or described de novo) need to be carefully assessed. While DAPC has rapidly become a core technique, the sensitivity of the method to misspecification of groups and how it is being empirically applied, are unknown. To address this, we conducted a simulation study examining the influence of a priori versus de novo group designations, and a literature review of how DAPC is being applied. We found that with a priori groupings, distance between genetic clusters reflected underlying F_ST. However, when migration rates were high and groups were described de novo there was considerable inaccuracy, both in terms of the number of genetic clusters suggested and placement of individuals into those clusters. Nearly all (90.1%) of 224 studies surveyed used DAPC to find de novo clusters, and for the majority (62.5%) the stated goal matched the results. However, most studies (52.3%) omit key run parameters, preventing repeatability and transparency. Therefore, we present recommendations for standard reporting of parameters used in DAPC analyses. The influence of groupings in genetic clustering is not unique to DAPC, and researchers need to consider their goal and which methods will be most appropriate.

Nonbifurcating Phylogenetic Tree Inference via the Adaptive LASSO

Phylogenetic tree inference using deep DNA sequencing is reshaping our understanding of rapidly evolving systems, such as the within-host battle between viruses and the immune system. Densely sampled phylogenetic trees can contain special features, including sampled ancestors in which we sequence a genotype along with its direct descendants, and polytomies in which multiple descendants arise simultaneously. These features are apparent after identifying zero-length branches in the tree. However, current maximum-likelihood based approaches are not capable of revealing such zero-length branches. In this paper, we find these zero-length branches by introducing adaptive-LASSO-type regularization estimators for the branch lengths of phylogenetic trees, deriving their properties, and showing regularization to be a practically useful approach for phylogenetics.

Nanopore sequencing-based genome assembly and evolutionary genomics of circum-basmati rice

Background

The circum-basmati group of cultivated Asian rice (Oryza sativa) contains many iconic varieties and is widespread in the Indian subcontinent. Despite its economic and cultural importance, a high-quality reference genome is currently lacking, and the group’s evolutionary history is not fully resolved. To address these gaps, we use long-read nanopore sequencing and assemble the genomes of two circum-basmati rice varieties.

Results

We generate two high-quality, chromosome-level reference genomes that represent the 12 chromosomes of Oryza. The assemblies show a contig N50 of 6.32 Mb and 10.53 Mb for Basmati 334 and Dom Sufid, respectively. Using our highly contiguous assemblies, we characterize structural variations segregating across circum-basmati genomes. We discover repeat expansions not observed in japonica—the rice group most closely related to circum-basmati—as well as the presence and absence variants of over 20 Mb, one of which is a circum-basmati-specific deletion of a gene regulating awn length. We further detect strong evidence of admixture between the circum-basmati and circum-aus groups. This gene flow has its greatest effect on chromosome 10, causing both structural variation and single-nucleotide polymorphism to deviate from genome-wide history. Lastly, population genomic analysis of 78 circum-basmati varieties shows three major geographically structured genetic groups: Bhutan/Nepal, India/Bangladesh/Myanmar, and Iran/Pakistan.

Conclusion

The availability of high-quality reference genomes allows functional and evolutionary genomic analyses providing genome-wide evidence for gene flow between circum-aus and circum-basmati, describes the nature of circum-basmati structural variation, and reveals the presence/absence variation in this important and iconic rice variety group.

Comparative Genomics and Phylogenomic Analysis of the Genus Salinivibrio

In the genomic era phylogenetic relationship among prokaryotes can be inferred from the core orthologous genes (OGs) or proteins in order to elucidate their evolutionary history and current taxonomy should benefits of that. The genus Salinivibrio belongs to the family Vibrionaceae and currently includes only five halophilic species, in spite the fact that new strains are very frequently isolated from hypersaline environments. Species belonging to this genus have undergone several reclassifications and, moreover, there are many strains of Salinivibrio with available genomes which have not been affiliated to the existing species or have been wrongly designated. Therefore, a phylogenetic study using the available genomic information is necessary to clarify the relationships of existing strains within this genus and to review their taxonomic affiliation. For that purpose, we have also sequenced the first complete genome of a Salinivibrio species, Salinivibrio kushneri AL184T, which was employed as a reference to order the contigs of the draft genomes of the type strains of the current species of this genus, as well as to perform a comparative analysis with all the other available Salinivibrio sp. genomes. The genome of S. kushneri AL184T was assembled in two circular chromosomes (with sizes of 2.84 Mb and 0.60 Mb, respectively), as typically occurs in members of the family Vibrionaceae, with nine complete ribosomal operons, which might explain the fast growing rate of salinivibrios cultured under laboratory conditions. Synteny analysis among the type strains of the genus revealed a high level of genomic conservation in both chromosomes, which allow us to hypothesize a slow speciation process or homogenization events taking place in this group of microorganisms to be tested experimentally in the future. Phylogenomic and orthologous average nucleotide identity (OrthoANI)/average amino acid identity (AAI) analyses also evidenced the elevated level of genetic relatedness within members of this genus and allowed to group all the Salinivibrio strains with available genomes in seven separated species. Genome-scale attribute study of the salinivibrios identified traits related to polar flagellum, facultatively anaerobic growth and osmotic response, in accordance to the phenotypic features described for species of this genus.

Multi-gene incongruence consistent with hybridisation in Cladocopium (Symbiodiniaceae), an ecologically important genus of coral reef symbionts

Coral reefs rely on their intracellular dinoflagellate symbionts (family Symbiodiniaceae) for nutritional provision in nutrient-poor waters, yet this association is threatened by thermally stressful conditions. Despite this, the evolutionary potential of these symbionts remains poorly characterised. In this study, we tested the potential for divergent Symbiodiniaceae types to sexually reproduce (i.e. hybridise) within Cladocopium, the most ecologically prevalent genus in this family. With sequence data from three organelles (cob gene, mitochondrion; psbA^ncr region, chloroplast; and ITS2 region, nucleus), we utilised the Incongruence Length Difference test, Approximately Unbiased test, tree hybridisation analyses and visual inspection of raw data in stepwise fashion to highlight incongruences between organelles, and thus provide evidence of reticulate evolution. Using this approach, we identified three putative hybrid Cladocopium samples among the 158 analysed, at two of the seven sites sampled. These samples were identified as the common Cladocopium types C40 or C1 with respect to the mitochondria and chloroplasts, but the rarer types C3z, C3u and C1# with respect to their nuclear identity. These five Cladocopium types have previously been confirmed as evolutionarily distinct and were also recovered in non-incongruent samples multiple times, which is strongly suggestive that they sexually reproduced to produce the incongruent samples. A concomitant inspection of next generation sequencing data for these samples suggests that other plausible explanations, such as incomplete lineage sorting or the presence of co-dominance, are much less likely. The approach taken in this study allows incongruences between gene regions to be identified with confidence, and brings new light to the evolutionary potential within Symbiodiniaceae.

Sponges Lack ParaHox Genes

Addressing the origin of axial-patterning machinery is essential for understanding the evolution of animal form. Historically, sponges, a lineage that branched off early in animal evolution, were thought to lack Hox and ParaHox genes, suggesting that these critical axial-patterning genes arose after sponges diverged. However, a recent study has challenged this long-held doctrine by claiming to identify ParaHox genes (Cdx family) in two calcareous sponge species, Sycon ciliatum and Leucosolenia complicata. We reanalyzed the main data sets in this paper and analyzed an additional data set that expanded the number of bilaterians represented and removed outgroup homeodomains. As in the previous study, our Neighbor-Joining analyses of the original data sets recovered a clade that included sponge and Cdx genes, whereas Bayesian analyses placed these sponge genes within the NKL subclass of homeodomains. Unlike the original study, only one of our two maximum-likelihood analyses was congruent with Cdx genes in sponges. Our analyses of our additional data set led to the sponge genes consistently being placed within the NKL subclass of homeodomains regardless of method or model. Our results show more support for these sponge genes belonging to the NKL subclass, and therefore imply that Hox and ParaHox genes arose after Porifera diverged from the rest of animals.

Phylogeny of Paleozoic limbed vertebrates reassessed through revision and expansion of the largest published relevant data matrix

The largest published phylogenetic analysis of early limbed vertebrates (Ruta M, Coates MI. 2007. Journal of Systematic Palaeontology 5:69-122) recovered, for example, Seymouriamorpha, Diadectomorpha and (in some trees) Caudata as paraphyletic and found the "temnospondyl hypothesis" on the origin of Lissamphibia (TH) to be more parsimonious than the "lepospondyl hypothesis" (LH)-though only, as we show, by one step. We report 4,200 misscored cells, over half of them due to typographic and similar accidental errors. Further, some characters were duplicated; some had only one described state; for one, most taxa were scored after presumed relatives. Even potentially continuous characters were unordered, the effects of ontogeny were not sufficiently taken into account, and data published after 2001 were mostly excluded. After these issues are improved-we document and justify all changes to the matrix-but no characters are added, we find (Analysis R1) much longer trees with, for example, monophyletic Caudata, Diadectomorpha and (in some trees) Seymouriamorpha; Ichthyostega either crownward or rootward of Acanthostega; and Anthracosauria either crownward or rootward of Temnospondyli. The LH is nine steps shorter than the TH (R2; constrained) and 12 steps shorter than the "polyphyly hypothesis" (PH-R3; constrained). Brachydectes (Lysorophia) is not found next to Lissamphibia; instead, a large clade that includes the adelogyrinids, urocordylid "nectrideans" and aïstopods occupies that position. As expected from the taxon/character ratio, most bootstrap values are low. Adding 56 terminal taxa to the original 102 increases the resolution (and decreases most bootstrap values). The added taxa range in completeness from complete articulated skeletons to an incomplete lower jaw. Even though the lissamphibian-like temnospondyls Gerobatrachus, Micropholis and Tungussogyrinus and the extremely peramorphic salamander Chelotriton are added, the difference between LH (R4; unconstrained) and TH (R5) rises to 10 steps, that between LH and PH (R6) to 15; the TH also requires several more regains of lost bones than the LH. Casineria, in which we tentatively identify a postbranchial lamina, emerges rather far from amniote origins in a gephyrostegid-chroniosuchian grade. Bayesian inference (Analysis EB, settings as in R4) mostly agrees with R4. High posterior probabilities are found for Lissamphibia (1.00) and the LH (0.92); however, many branches remain weakly supported, and most are short, as expected from the small character sample. We discuss phylogeny, approaches to coding, methods of phylogenetics (Bayesian inference vs. equally weighted vs. reweighted parsimony), some character complexes (e.g. preaxial/postaxial polarity in limb development), and prospects for further improvement of this matrix. Even in its revised state, the matrix cannot provide a robust assessment of the phylogeny of early limbed vertebrates. Sufficient improvement will be laborious-but not difficult.

Phylogenomic incongruence, hypothesis testing, and taxonomic sampling: The monophyly of characiform fishes

Phylogenomic studies using genome-wide datasets are quickly becoming the state of the art for systematics and comparative studies, but in many cases, they result in strongly supported incongruent results. The extent to which this conflict is real depends on different sources of error potentially affecting big datasets (assembly, stochastic, and systematic error). Here, we apply a recently developed methodology (GGI or gene genealogy interrogation) and data curation to new and published datasets with more than 1000 exons, 500 ultraconserved element (UCE) loci, and transcriptomic sequences that support incongruent hypotheses. The contentious non-monophyly of the order Characiformes proposed by two studies is shown to be a spurious outcome induced by sample contamination in the transcriptomic dataset and an ambiguous result due to poor taxonomic sampling in the UCE dataset. By exploring the effects of number of taxa and loci used for analysis, we show that the power of GGI to discriminate among competing hypotheses is diminished by limited taxonomic sampling, but not equally sensitive to gene sampling. Taken together, our results reinforce the notion that merely increasing the number of genetic loci for a few representative taxa is not a robust strategy to advance phylogenetic knowledge of recalcitrant groups. We leverage the expanded exon capture dataset generated here for Characiformes (206 species in 23 out of 24 families) to produce a comprehensive phylogeny and a revised classification of the order.

Functional and evolutionary characterization of a secondary metabolite gene cluster in budding yeasts

Secondary metabolites are key in how organisms from all domains of life interact with their environment and each other. The iron-binding molecule pulcherrimin was described a century ago, but the genes responsible for its production in budding yeasts have remained uncharacterized. Here, we used phylogenomic footprinting on 90 genomes across the budding yeast subphylum Saccharomycotina to identify the gene cluster associated with pulcherrimin production. Using targeted gene replacements in Kluyveromyces lactis, we characterized the four genes that make up the cluster, which likely encode two pulcherriminic acid biosynthesis enzymes, a pulcherrimin transporter, and a transcription factor involved in both biosynthesis and transport. The requirement of a functional putative transporter to utilize extracellular pulcherrimin-complexed iron demonstrates that pulcherriminic acid is a siderophore, a chelator that binds iron outside the cell for subsequent uptake. Surprisingly, we identified homologs of the putative transporter and transcription factor genes in multiple yeast genera that lacked the biosynthesis genes and could not make pulcherrimin, including the model yeast Saccharomyces cerevisiae We deleted these previously uncharacterized genes and showed they are also required for pulcherrimin utilization in S. cerevisiae, raising the possibility that other genes of unknown function are linked to secondary metabolism. Phylogenetic analyses of this gene cluster suggest that pulcherrimin biosynthesis and utilization were ancestral to budding yeasts, but the biosynthesis genes and, subsequently, the utilization genes, were lost in many lineages, mirroring other microbial public goods systems that lead to the rise of cheater organisms.

The Interrelationships of Land Plants and the Nature of the Ancestral Embryophyte

The evolutionary emergence of land plant body plans transformed the planet. However, our understanding of this formative episode is mired in the uncertainty associated with the phylogenetic relationships among bryophytes (hornworts, liverworts, and mosses) and tracheophytes (vascular plants). Here we attempt to clarify this problem by analyzing a large transcriptomic dataset with models that allow for compositional heterogeneity between sites. Zygnematophyceae is resolved as sister to land plants, but we obtain several distinct relationships between bryophytes and tracheophytes. Concatenated sequence analyses that can explicitly accommodate site-specific compositional heterogeneity give more support for a mosses-liverworts clade, "Setaphyta," as the sister to all other land plants, and weak support for hornworts as the sister to all other land plants. Bryophyte monophyly is supported by gene concatenation analyses using models explicitly accommodating lineage-specific compositional heterogeneity and analyses of gene trees. Both maximum-likelihood analyses that compare the fit of each gene tree to proposed species trees and Bayesian supertree estimation based on gene trees support bryophyte monophyly. Of the 15 distinct rooted relationships for embryophytes, we reject all but three hypotheses, which differ only in the position of hornworts. Our results imply that the ancestral embryophyte was more complex than has been envisaged based on topologies recognizing liverworts as the sister lineage to all other embryophytes. This requires many phenotypic character losses and transformations in the liverwort lineage, diminishes inconsistency between phylogeny and the fossil record, and prompts re-evaluation of the phylogenetic affinity of early land plant fossils, the majority of which are considered stem tracheophytes.

Distribution of endangered Italian gudgeon Romanogobio benacensis (Cypriniformes, Cyprinidae, Gobioninae) with remarks on distinguishing morphological characters

Distribution data on many freshwater fish species in Croatia are scarce and species identifications are difficult, requiring further detailed studies. This paper presents a report of the Italian gudgeon Romanogobiobenacensis from the Mirna River in the Istra Peninsula in Croatia, in the south-east from its previously known distribution range. The identification of R.benacensis in Croatia was supported by a morphological comparison with R.benacensis from Italy and Slovenia, the common gudgeon Gobiogobio, and the Danubian gudgeon Gobioobtusirostris from geographically close locations. A combination of character states (number of scales between anus and anal-fin origin, branched pectoral-fin rays, lateral-line scales, total, abdominal, and caudal vertebrae, and the size and number of lateral blotches) distinguishes R.benacensis from both G.gobio and G.obtusirostris. The phylogenetic analyses using mitochondrial sequences of cytochrome b gene confirmed that specimens from the Mirna River belong to R.benacensis. Also, Reka River system (Adriatic Sea basin) in Slovenia is inhabited by a possibly introduced Danubian gudgeon, G.obtusirostris, and not by R.benacensis.

Bayesian selection of misspecified models is overconfident and may cause spurious posterior probabilities for phylogenetic trees

The Bayesian method is noted to produce spuriously high posterior probabilities for phylogenetic trees in analysis of large datasets, but the precise reasons for this overconfidence are unknown. In general, the performance of Bayesian selection of misspecified models is poorly understood, even though this is of great scientific interest since models are never true in real data analysis. Here we characterize the asymptotic behavior of Bayesian model selection and show that when the competing models are equally wrong, Bayesian model selection exhibits surprising and polarized behaviors in large datasets, supporting one model with full force while rejecting the others. If one model is slightly less wrong than the other, the less wrong model will eventually win when the amount of data increases, but the method may become overconfident before it becomes reliable. We suggest that this extreme behavior may be a major factor for the spuriously high posterior probabilities for evolutionary trees. The philosophical implications of our results to the application of Bayesian model selection to evaluate opposing scientific hypotheses are yet to be explored, as are the behaviors of non-Bayesian methods in similar situations.

Phylogeny of Vibrio vulnificus from the Analysis of the Core-Genome: Implications for Intra-Species Taxonomy

Vibrio vulnificus (Vv) is a multi-host pathogenic species currently subdivided into three biotypes (Bts). The three Bts are human-pathogens, but only Bt2 is also a fish-pathogen, an ability that is conferred by a transferable virulence-plasmid (pVvbt2). Here we present a phylogenomic analysis from the core genome of 80 Vv strains belonging to the three Bts recovered from a wide range of geographical and ecological sources. We have identified five well-supported phylogenetic groups or lineages (L). L1 comprises a mixture of clinical and environmental Bt1 strains, most of them involved in human clinical cases related to raw seafood ingestion. L2 is formed by a mixture of Bt1 and Bt2 strains from various sources, including diseased fish, and is related to the aquaculture industry. L3 is also linked to the aquaculture industry and includes Bt3 strains exclusively, mostly related to wound infections or secondary septicemia after farmed-fish handling. Lastly, L4 and L5 include a few strains of Bt1 associated with specific geographical areas. The phylogenetic trees for ChrI and II are not congruent to one another, which suggests that inter- and/or intra-chromosomal rearrangements have been produced along Vv evolution. Further, the phylogenetic trees for each chromosome and the virulence plasmid were also not congruent, which also suggests that pVvbt2 has been acquired independently by different clones, probably in fish farms. From all these clones, the one with zoonotic capabilities (Bt2-Serovar E) has successfully spread worldwide. Based on these results, we propose a new updated classification of the species based on phylogenetic lineages rather than on Bts, as well as the inclusion of all Bt2 strains in a pathovar with the particular ability to cause fish vibriosis, for which we suggest the name "piscis."

Studying AMPK in an Evolutionary Context

The AMPK protein kinase forms the heart of a complex network controlling the metabolic activities in a eukaryotic cell. Unraveling the steps by which this pathway evolved from its primordial roots in the last eukaryotic common ancestor to its present status in contemporary species has the potential to shed light on the evolution of eukaryotes. A homolog search for the proteins interacting in this pathway is considerably straightforward. However, interpreting the results, when reconstructing the evolutionary history of the pathway over larger evolutionary distances, bears a number of pitfalls. With this in mind, we present a protocol to trace a metabolic pathway across contemporary species and backward in evolutionary time. Alongside the individual analysis steps, we provide guidelines for data interpretation generalizing beyond the analysis of AMPK.

Phylogenomics

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. The abundance of genomic data for an enormous variety of organisms has enabled phylogenomic inference of many groups, and this has motivated the development of many computer programs implementing the associated methods. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Using Minimum Bootstrap support for Splits to Construct Confidence Regions for Trees

Many of the estimated topologies in phylogenetic studies are presented with the bootstrap support for each of the splits in the topology indicated. If phylogenetic estimation is unbiased, high bootstrap support for a split suggests that there is a good deal of certainty that the split actually is present in the tree and low bootstrap support suggests that one or more of the taxa on one side of the estimated split might in reality be located with taxa on the other side. In the latter case the follow-up questions about how many and which of the taxa could reasonably be incorrectly placed as well as where they might alternatively be placed are not addressed through the presented bootstrap support. We present here an algorithm that finds the set of all trees with minimum bootstrap support for their splits greater than some given value. The output is a ranked list of trees, ranked according to the minimum bootstrap supports for splits in the trees. The number of such trees and their topologies provides useful supplementary information in bootstrap analyses about the reasons for low bootstrap support for splits. We also present ways of quantifying low bootstrap support by considering the set of all topologies with minimum bootstrap greater than some quantity as providing a confidence region of topologies. Using a double bootstrap we are able to choose a cutoff so that the set of topologies with minimum bootstrap support for a split greater than that cutoff gives an approximate 95% confidence region. As with bootstrap support one advantage of the methods is that they are generally applicable to the wide variety of phylogenetic estimation methods.

An African Origin of the Eurylaimides (Passeriformes) and the Successful Diversification of the Ground-Foraging Pittas (Pittidae)

The Eurylaimides is one of the few passerine groups with a pantropical distribution. In this study, we generated a multi-calibrated tree with 83% of eurylaimid species diversity based on 30 molecular loci. Particular attention was given to the monotypic Sapayoidae to reconstruct the biogeography of this radiation. We conducted several topological tests including nonoverlapping subsampling of the concatenated alignment and coalescent species tree reconstruction. These tests firmly placed the South American Sapayoidae as the sister group to all other Eurylaimides families (split at ∼28 Ma), with increasing branch support as highly variable sites were removed. This topology is consistent with the breakup of the insular connection between Africa and South America (Atlantogea) that took place between the middle Eocene and the early Oligocene. We recovered Africa as the cradle of the core Eurylaimides, and this result is supported by all African lineages corresponding to the oldest splits within each family in this group. Our timescale suggests that desertification and the uplift of the Tibetan Plateau caused a parallel divergence between African and Asian lineages in all major clades in the core Eurylaimides at 22-9 Ma. We also propose that the ground-foraging behavior in the Pittidae ancestor allowed the pitta lineage to thrive and coexist with the older arboreal lineages of the core Eurylaimides. In contrast, the diversification of pittas in Australia was likely hindered by direct competition with the endemic ground-foraging oscines that had been well established in that continent since the Eocene.

Which molecular markers for assessing which taxonomic level? The case study of the mite family Phytoseiidae (Acari: Mesostigmata)

The use of molecular markers for resolving systematics issues has improved our knowledge of life history. However, for the taxa studied herein-the predatory mite family Phytoseiidae-molecular phylogeny is impeded by a lack of suitable markers for deeper taxonomic levels. This study aims (i) to establish DNA amplification protocols for molecular markers known to resolve supraspecific nodes in other taxa, (ii) to determine their individual performance in assessing the clustering of species, genera, tribes and subfamilies, and (iii) to characterize the additional information provided when markers are concatenated. A new phylogenetic index is proposed based on ecological concepts, considering trees as a community of nodes. New and efficient protocols for DNA amplification of six molecular markers are provided. The concatenated tree globally provides more robust and reliable information, especially for deeper nodes. However, for assessing species identification and within-genera phylogenies, the combined use of six markers does not seem necessary, underlining the need to resize experiments depending on their taxonomic objectives. Finally, this study lays the methodological foundations with which to test the present Phytoseiidae classification as the first phylogeny obtained shows incongruence with the present morphological classification.

Phylogenomic analyses of Crassiclitellata support major Northern and Southern Hemisphere clades and a Pangaean origin for earthworms

BACKGROUND

Earthworms (Crassiclitellata) are a diverse group of annelids of substantial ecological and economic importance. Earthworms are primarily terrestrial infaunal animals, and as such are probably rather poor natural dispersers. Therefore, the near global distribution of earthworms reflects an old and likely complex evolutionary history. Despite a long-standing interest in Crassiclitellata, relationships among and within major clades remain unresolved.

METHODS

In this study, we evaluate crassiclitellate phylogenetic relationships using 38 new transcriptomes in combination with publicly available transcriptome data. Our data include representatives of nearly all extant earthworm families and a representative of Moniligastridae, another terrestrial annelid group thought to be closely related to Crassiclitellata. We use a series of differentially filtered data matrices and analyses to examine the effects of data partitioning, missing data, compositional and branch-length heterogeneity, and outgroup inclusion.

RESULTS AND DISCUSSION

We recover a consistent, strongly supported ingroup topology irrespective of differences in methodology. The topology supports two major earthworm clades, each of which consists of a Northern Hemisphere subclade and a Southern Hemisphere subclade. Divergence time analysis results are concordant with the hypothesis that these north-south splits are the result of the breakup of the supercontinent Pangaea.

CONCLUSIONS

These results support several recently proposed revisions to the classical understanding of earthworm phylogeny, reveal two major clades that seem to reflect Pangaean distributions, and raise new questions about earthworm evolutionary relationships.

Unique Phylogenetic Distributions of the Ska and Dam1 Complexes Support Functional Analogy and Suggest Multiple Parallel Displacements of Ska by Dam1

Faithful chromosome segregation relies on kinetochores, the large protein complexes that connect chromatin to spindle microtubules. Although human and yeast kinetochores are largely homologous, they track microtubules with the unrelated protein complexes Ska (Ska-C, human) and Dam1 (Dam1-C, yeast). We here uncovered that Ska-C and Dam1-C are both widespread among eukaryotes, but in an exceptionally inverse manner, supporting their functional analogy. Within the complexes, all Ska-C and various Dam1-C subunits are ancient paralogs, showing that gene duplication shaped these complexes. We examined various evolutionary scenarios to explain the nearly mutually exclusive patterns of Ska-C and Dam1-C in present-day species. We propose that Ska-C was present in the last eukaryotic common ancestor, that subsequently Dam1-C displaced Ska-C in an early fungus and was horizontally transferred to diverse non-fungal lineages, displacing Ska-C in these lineages too.