Comparative analysis of 12 water lily plastid genomes reveals genomic divergence and evolutionary relationships in early flowering plants

The aquatic plant Nymphaea, a model genus of the early flowering plant lineage Nymphaeales and family Nymphaeaceae, has been extensively studied. However, the availability of chloroplast genome data for this genus is incomplete, and phylogenetic relationships within the order Nymphaeales remain controversial. In this study, 12 chloroplast genomes of Nymphaea were assembled and analyzed for the first time. These genomes were 158,290-160,042 bp in size and contained 113 non-repeat genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. We also report on codon usage, RNA editing sites, microsatellite structures, and new repetitive sequences in this genus. Comparative genomics revealed that expansion and contraction of IR regions can lead to changes in the gene numbers. Additionally, it was observed that the highly variable regions of the chloroplast genome were mainly located in intergenic regions. Furthermore, the phylogenetic tree showed the order Nymphaeales was divided into three families, and the genus Nymphaea can be divided into five (or three) subgenera, with the subgenus Nymphaea being the oldest. The divergence times of nymphaealean taxa were analyzed, with origins of the order Nymphaeales and family Nymphaeaceae being about 194 and 131 million years, respectively. The results of the phylogenetic analysis and estimated divergence times will be useful for future evolutionary studies of basal angiosperm lineages.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00242-0.

Hidden in the riffles: A new suckermouth catfish (Mochokidae, Chiloglanis) from the middle Zambezi River system, Zimbabwe

The recent surge in the discovery of hidden diversity within rheophilic taxa, particularly in West and East Africa, prompted a closer examination of the extent to which the current taxonomy may obscure the diversity of riffle-dwelling suckermouth catfishes in the genus Chiloglanis in southern Africa. Currently, the region comprises eight valid species within this genus. Seven of them have relatively narrow geographic distribution ranges except for C.neumanni, which is considered to be widely distributed, occurring from the Buzi River system in the south, and its northern limit being the eastward draining river systems in Tanzania. Recent surveys of the middle Zambezi River system revealed Chiloglanis specimens that were distinguishable from the known species of the genus from southern Africa. Integration of molecular and morphological data indicated that these specimens from the Mukwadzi River represent a new species to science, herein described as Chiloglaniscarnatus Mutizwa, Bragança & Chakona, sp. nov. This species is readily distinguished from its southern African congeners by the possession of a distinctive extended dermal tissue covering the base of the dorsal fin and the possession of ten mandibular teeth (vs 8, 12, or 14 in the other taxa). Results from this study add to the growing evidence of a high level of undocumented diversity within riffle-dwelling taxa in southern Africa.

Nucleotide Substitution Model Selection Is Not Necessary for Bayesian Inference of Phylogeny With Well-Behaved Priors

Model selection aims to choose the most adequate model for the statistical analysis at hand. The model must be complex enough to capture the complexity of the data but should be simple enough not to overfit. In phylogenetics, the most common model selection scenario concerns selecting an adequate substitution and partition model for sequence evolution to infer a phylogenetic tree. Previously, several studies showed that substitution model under-parameterization can bias phylogenetic studies. Here, we explored the impact of substitution model over-parameterization in a Bayesian statistical framework. We performed simulations under the simplest substitution model, the Jukes-Cantor model, and compare posterior estimates of phylogenetic tree topologies and tree length under the true model to the most complex model, the $\text{GTR}+\Gamma+\text{I}$ substitution model, including over-splitting the data into additional subsets (i.e., applying partitioned models). We explored 4 choices of prior distributions: the default substitution model priors of MrBayes, BEAST2, and RevBayes and a newly devised prior choice (Tame). Our results show that Bayesian inference of phylogeny is robust to substitution model over-parameterization and over-partitioning but only under our new prior settings. All 3 current default priors introduced biases for the estimated tree length. We conclude that substitution and partition model selection are superfluous steps in Bayesian phylogenetic inference pipelines if well-behaved prior distributions are applied and more effort should focus on more complex and biologically realistic substitution models.

Phylogenetic Relationships of the Mutualistic Fungi Associated with Macrotermes subhyalinus in Oman

The symbiotic association between fungus-gardening termites Macrotermes and its fungal symbiont has a moderate degree of specificity-although the symbiotic fungi (Termitomyces) form a monophyletic clade, there is not a one-to-one association between termite species and their fungus-garden associates. Here, we aim to determine the origin and phylogenetic relationships of Termitomyces in Oman. We used sequences of the internal transcribed spacer region (ITS) and the nuclear large subunit ribosomal RNA (LSU rRNA, 25S) gene and analyzed these with sequences of Termitomyces from other geographic areas. We find no evidence for more than a single colonization of Oman by Termitomyces. Unexpectedly, we find Termitomyces in Oman is most closely related to the symbiont of M. subhyalinus in West Africa rather than to those of geographically closer populations in East Africa.

Phylodynamic approaches to studying avian influenza virus

Avian influenza viruses can cause severe disease in domestic and wild birds and are a pandemic threat. Phylodynamics is the study of how epidemiological, evolutionary, and immunological processes can interact to shape viral phylogenies. This review summarizes how phylodynamic methods have and could contribute to the study of avian influenza viruses. Specifically, we assess how phylodynamics can be used to examine viral spread within and between wild or domestic bird populations at various geographical scales, identify factors associated with virus dispersal, and determine the order and timing of virus lineage movement between geographic regions or poultry production systems. We discuss factors that can complicate the interpretation of phylodynamic results and identify how future methodological developments could contribute to improved control of the virus.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Variational Phylodynamic Inference Using Pandemic-scale Data

The ongoing global pandemic has sharply increased the amount of data available to researchers in epidemiology and public health. Unfortunately, few existing analysis tools are capable of exploiting all of the information contained in a pandemic-scale data set, resulting in missed opportunities for improved surveillance and contact tracing. In this paper, we develop the variational Bayesian skyline (VBSKY), a method for fitting Bayesian phylodynamic models to very large pathogen genetic data sets. By combining recent advances in phylodynamic modeling, scalable Bayesian inference and differentiable programming, along with a few tailored heuristics, VBSKY is capable of analyzing thousands of genomes in a few minutes, providing accurate estimates of epidemiologically relevant quantities such as the effective reproduction number and overall sampling effort through time. We illustrate the utility of our method by performing a rapid analysis of a large number of SARS-CoV-2 genomes, and demonstrate that the resulting estimates closely track those derived from alternative sources of public health data.

Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

Turcichondrostoma, a new genus for the Leuciscidae (Teleostei: Cypriniformes) from Southwestern Anatolia

Turcichondrostoma, a new genus, from the Southwestern Anatolia is distinguished by having fewer gill rakers on first gill arch and morphologies of premaxilla and dentary bones. In addition, as a result of the phylogenetic analyses based on combine data set (mtDNA COI + Cytb) sequences (1706 bp.), Turcichondrostoma genus was recovered with high posterior probability value (BI PP:1.0) and strong-supported bootstrap value (ML BP: 100%) among the former Chondrostoma groups. Also, high K2P mean genetic distance values (more than 7.84%) differentiated genus Turcichondrostoma from the other genera of former Chondrostoma group. The results of both morphological-osteological and molecular analyses are congruent with each other. The results of this study revealed that the genus Turcichondrostoma is easily distinguished from the genera in Chondrostoma group.

Accuracy in near-perfect virus phylogenies

Phylogenetic trees from real-world data often include short edges with very few substitutions per site, which can lead to partially resolved trees and poor accuracy. Theory indicates that the number of sites needed to accurately reconstruct a fully resolved tree grows at a rate proportional to the inverse square of the length of the shortest edge. However, when inferred trees are partially resolved due to short edges, “accuracy” should be defined as the rate of discovering false splits (clades on a rooted tree) relative to the actual number found. Thus, accuracy can be high even if short edges are common. Specifically, in a “near-perfect” parameter space in which trees are large, the tree length ξ (the sum of all edge lengths) is small, and rate variation is minimal, the expected false positive rate is less than ξ∕3; the exact value depends on tree shape and sequence length. This expected false positive rate is far below the false negative rate for small ξ and often well below 5% even when some assumptions are relaxed. We show this result analytically for maximum parsimony and explore its extension to maximum likelihood using theory and simulations. For hypothesis testing, we show that measures of split “support” that rely on bootstrap resampling consistently imply weaker support than that implied by the false positive rates in near-perfect trees. The near-perfect parameter space closely fits several empirical studies of human virus diversification during outbreaks and epidemics, including Ebolavirus, Zika virus, and SARS-CoV-2, reflecting low substitution rates relative to high transmission/sampling rates in these viruses.

Relative Efficiencies of Simple and Complex Substitution Models in Estimating Divergence Times in Phylogenomics

The conventional wisdom in molecular evolution is to apply parameter-rich models of nucleotide and amino acid substitutions for estimating divergence times. However, the actual extent of the difference between time estimates produced by highly complex models compared with those from simple models is yet to be quantified for contemporary data sets that frequently contain sequences from many species and genes. In a reanalysis of many large multispecies alignments from diverse groups of taxa, we found that the use of the simplest models can produce divergence time estimates and credibility intervals similar to those obtained from the complex models applied in the original studies. This result is surprising because the use of simple models underestimates sequence divergence for all the data sets analyzed. We found three fundamental reasons for the observed robustness of time estimates to model complexity in many practical data sets. First, the estimates of branch lengths and node-to-tip distances under the simplest model show an approximately linear relationship with those produced by using the most complex models applied on data sets with many sequences. Second, relaxed clock methods automatically adjust rates on branches that experience considerable underestimation of sequence divergences, resulting in time estimates that are similar to those from complex models. And, third, the inclusion of even a few good calibrations in an analysis can reduce the difference in time estimates from simple and complex models. The robustness of time estimates to model complexity in these empirical data analyses is encouraging, because all phylogenomics studies use statistical models that are oversimplified descriptions of actual evolutionary substitution processes.

Comparative chloroplast genomes: insights into the evolution of the chloroplast genome of Camellia sinensis and the phylogeny of Camellia

BACKGROUND

Chloroplast genome resources can provide useful information for the evolution of plant species. Tea plant (Camellia sinensis) is among the most economically valuable member of Camellia. Here, we determined the chloroplast genome of the first natural triploid Chinary type tea ('Wuyi narcissus' cultivar of Camellia sinensis var. sinensis, CWN) and conducted the genome comparison with the diploid Chinary type tea (Camellia sinensis var. sinensis, CSS) and two types of diploid Assamica type teas (Camellia sinensis var. assamica: Chinese Assamica type tea, CSA and Indian Assamica type tea, CIA). Further, the evolutionary mechanism of the chloroplast genome of Camellia sinensis and the relationships of Camellia species based on chloroplast genome were discussed.

RESULTS

Comparative analysis showed the evolutionary dynamics of chloroplast genome of Camellia sinensis were the repeats and insertion-deletions (indels), and distribution of the repeats, indels and substitutions were significantly correlated. Chinese tea and Indian tea had significant differences in the structural characteristic and the codon usage of the chloroplast genome. Analysis of sequence characterized amplified region (SCAR) using sequences of the intergenic spacers (trnE/trnT) showed none of 292 different Camellia sinensis cultivars had similar sequence characteristic to triploid CWN, but the other four Camellia species did. Estimations of the divergence time showed that CIA diverged from the common ancestor of two Assamica type teas about 6.2 Mya (CI: 4.4-8.1 Mya). CSS and CSA diverged to each other about 0.8 Mya (CI: 0.4-1.5 Mya). Moreover, phylogenetic clustering was not exactly consistent with the current taxonomy of Camellia.

CONCLUSIONS

The repeat-induced and indel-induced mutations were two important dynamics contributed to the diversification of the chloroplast genome in Camellia sinensis, which were not mutually exclusive. Chinese tea and Indian tea might have undergone different selection pressures. Chloroplast transfer occurred during the polyploid evolution in Camellia sinensis. In addition, our results supported the three different domestication origins of Chinary type tea, Chinese Assamica type tea and Indian Assamica type tea. And, the current classification of some Camellia species might need to be further discussed.

Cryptic speciation in populations of the genus Aphyocharax (Characiformes: Characidae) from eastern Amazon coastal river drainages and surroundings revealed by single locus species delimitation methods

Abstract Recent studies in eastern Amazon coastal drainages and their surroundings have revealed new fish species that sometimes exhibit little morphological differentiation (cryptic species). Thus, we used a DNA-based species delimitation approach to test if populations showing the morphotype and typical character states of the Aphyocharax avary holotype correspond either to A. avary or A. brevicaudatus, two known species from the region, or if they form independent lineages, indicating cryptic speciation. WP and GMYC analyses recovered five lineages (species) in the ingroup, while a bPTP analysis delimited three lineages. ABGD analyses produced two possible results: one corroborating the WP and GMYC methods and another corroborating the bPTP method. All methods indicate undescribed cryptic species in the region and show variation from at least 1 to 4 species in the ingroup, depending on the approach, corroborating previous studies, and revealing this region as a possible hotspot for discovering undescribed fish species.

An integrative approach to the systematics of the Berthella californica species complex (Heterobranchia: Pleurobranchidae)

Berthella californica (W. H. Dall, 1900) is a widespread species of heterobranch sea slug distributed across the North Pacific Ocean, from Korea and Japan to the Galapagos Islands. Two distinct morphotypes are observed in B. californica, which differ in external coloration, egg-mass morphology and geographic distribution (with the exception of a small range overlap in Southern California). Molecular and morphological data obtained in this study reveals that these two morphotypes constitute distinct species. The name B. californica (type locality: San Pedro, California) is retained for the southern morphotype, whereas the name Berthella chacei (J. Q. Burch, 1944) (type locality: Crescent City, California) is resurrected for the northern morphotype. Moreover, molecular phylogenetic analyses recovered B. californica as sister to Berthellina, in a well-supported clade separate from Berthella, suggesting that the classification of B. californica may need additional revision.

Candidate cases of poecilogony in Neogastropoda: implications for the systematics of the genus Raphitoma Bellardi, 1847

Poecilogony is the intraspecific variation in developmental mode, with larvae of different types produced by the same individual, population or species. It is very rare among marine invertebrates, and in gastropods has long been described only in a few opisthobranchs. The physiological and regulatory mechanisms underlying larval evolutionary transitions, such as loss of planktotrophy that occurred repeatedly in many caenogastropod lineages, are still largely unknown. We have studied the inter- v. intraspecific variation in larval development in the north-east Atlantic neogastropod genus Raphitoma Bellardi, 1847, starting with an iterative taxonomy approach: 17 morphology-based Preliminary Species Hypotheses were tested against a COI molecular-distance-based method (ABGD), and the retained species hypotheses were eventually inspected for reciprocal monophyly on a multilocus dataset. We subsequently performed an ancestral state reconstruction on an ultrametric tree of the 10 Final Species Hypotheses, time-calibrated by fossils, revealing that the interspecific changes were planktotrophy > lecithotrophy, and all have occurred in the Pleistocene, after 2.5 million years ago. This is suggestive of a major role played by Pleistocene Mediterranean oceanographic conditions – enhanced oligotrophy, unpredictable availability of water column resources – likely to favour loss of planktotrophy. Within this group of species, which has diversified after the Miocene, we identified one pair of sibling species differing in their larval development, Raphitoma cordieri (Payraudeau, 1826) and R. horrida (Monterosato, 1884). However, we also identified two Final Species Hypotheses, each comprising individuals with both larval developmental types. Our working hypothesis is that they correspond to one or two poecilogonous species. If confirmed by other nuclear markers, this would be the first documentation of poecilogony in the Neogastropoda, and the second in the whole Caenogastropoda. Although sibling species with different developmental strategies may offer good models to study some evolutionary aspects, poecilogonous taxa are optimally suited for identifying regulatory and developmental mechanisms underlying evolutionary transitions.

A new species of Aphyocharax Günther, 1868 (Characiformes, Characidae) from the Maracaçumé river basin, eastern Amazon

A new species of Aphyocharax is described from the Maracaçumé river basin, eastern Amazon, based on morphological and molecular data. The new species differs from all its congeners, mainly by possessing the upper caudal-fin lobe longer than the lower one in mature males, and other characters related to teeth counts, colour pattern, and body depth at dorsal-fin origin. In addition, the new species is corroborated by a haplotype phylogenetic analyses based on the Cytochrome B (Cytb) mitochondrial gene, where its haplotypes are grouped into an exclusive lineage, supported by maximum posterior probability value, a species delimitation method termed the Wiens and Penkrot analysis (WP).

A new cryptic species of Hyphessobrycon Durbin, 1908 (Characiformes, Characidae) from the Eastern Amazon, revealed by integrative taxonomy

Hyphessobryconcarusp. nov. is described based on five different and independent methods of species delimitation, making the hypothesis of this new species supported by an integrative taxonomy perspective. This new species has a restricted distribution, occurring just in the upper Pindaré river drainage, Mearim river basin, Brazil. It is a member of the rosy tetra clade, which is characterized mainly by the presence of a dark brown or black blotch on dorsal fin and absence of a midlateral stripe on the body. Hyphessobryconcarusp. nov. is distinguished from the members of this clade mainly by the shape of its humeral spot, possessing few irregular inconspicuous vertically arranged chromatophores in the humeral region, or sometimes a very thin and inconspicuous humeral spot, and other characters related to teeth count, and color pattern. The phylogenetic position of the new species within the rosy tetra clade was based on molecular phylogenetic analysis using sequences of the mitochondrial gene cytochrome oxidase subunit 1. In addition, a new clade (here termed Hyphessobryconmicropterus clade) within the rosy tetra clade is proposed based on molecular data, comprising H.carusp. nov., H.micropterus, H.piorskii, and H.simulatus, and with H.carusp. nov. and H.piorskii recovered as sister species. Our results suggest cryptic speciation in the rosy tetra clade and, more specifically, in the H.micropterus clade. We recommend the use of integrative taxonomy for future taxonomic revisions and species descriptions when dealing with species complexes and groups containing possible cryptic species.

A Systematist's Guide to Estimating Bayesian Phylogenies From Morphological Data

Phylogenetic trees are crucial to many aspects of taxonomic and comparative biology. Many researchers have adopted Bayesian methods to estimate their phylogenetic trees. In this family of methods, a model of morphological evolution is assumed to have generated the data observed by the researcher. These models make a variety of assumptions about the evolution of morphological characters, and these assumptions are translated into mathematics as parameters. The incorporation of prior distributions further allows researchers to quantify their prior beliefs about the value any one parameter can take. How to translate biological knowledge into mathematical language is difficult, and can be confusing to many biologists. This review aims to help systematics researchers understand the biological meaning of common models and assumptions. Using examples from the insect fossil record, I will demonstrate empirically what assumptions mean in concrete terms, and discuss how researchers can use and understand Bayesian methods for phylogenetic estimation.

Bayesian and likelihood phylogenetic reconstructions of morphological traits are not discordant when taking uncertainty into consideration: a comment on Puttick et al

Puttick et al. (2017 Proc. R. Soc. B284, 20162290 (doi:10.1098/rspb.2016.2290)) performed a simulation study to compare accuracy among methods of inferring phylogeny from discrete morphological characters. They report that a Bayesian implementation of the Mk model (Lewis 2001 Syst. Biol.50, 913-925 (doi:10.1080/106351501753462876)) was most accurate (but with low resolution), while a maximum-likelihood (ML) implementation of the same model was least accurate. They conclude by strongly advocating that Bayesian implementations of the Mk model should be the default method of analysis for such data. While we appreciate the authors' attempt to investigate the accuracy of alternative methods of analysis, their conclusion is based on an inappropriate comparison of the ML point estimate, which does not consider confidence, with the Bayesian consensus, which incorporates estimation credibility into the summary tree. Using simulation, we demonstrate that ML and Bayesian estimates are concordant when confidence and credibility are comparably reflected in summary trees, a result expected from statistical theory. We therefore disagree with the conclusions of Puttick et al. and consider their prescription of any default method to be poorly founded. Instead, we recommend caution and thoughtful consideration of the model or method being applied to a morphological dataset.

Phylogeny and species delimitation based on molecular approaches on the species of the Australoheros autrani group (Teleostei, Cichlidae), with biogeographic comments

Three distinct and independent molecular-based species delimitation analyses were performed among the species and populations included within theAustraloherosautranigroup, based on sequences of the mitochondrial gene Cytochrome b: a tree-based method proposed by Wiens and Penkrot (WP), a Character-based DNA Barcoding (CBB) and coalescent species delimitation method termed the Bayesian Implementation of the Poisson tree processes (bPTP). The congruence of WP and CBB delimited 11 independent lineages (species), while the bPTP delimited just nine lineages. We did not favour any of the methods, and we considered the possibility of two slightly variant scenarios. A time-calibrated phylogenetic analysis is proposed based on the predominant congruence of the results of these three species delimitation methods herein applied. The monophyly of theA.autranispecies group was highly supported with maximum node support value and diagnosed by 11 nucleotide substitutions. The sister clade of theA.autranispecies group is the clade comprisingA.sp. Timbé do Sul andA.minuano. The phylogenetic analysis supports three main clades within theA.autranispecies group, supported by maximum node support value, with the Southern Mata Atlântica clade as the most basal clade. Divergence time estimates indicate that the diversification of theAustraloherosoriginated during the early Neogene, but only in the late Neogene did the processes of diversification in the southeast and north regions occur. Diversification within theAustraloherosautranispecies group occurred synchronically for the three main clades during the beginning of the Quaternary. It is demonstrated that molecular characters are valuable tools for species recognition, particularly in speciose groups with inconspicuous or difficult to record morphological characters. The resulting phylogeny of theAustraloherosautranigroup is highly compatible with the geological and biogeographic scenarios proposed for the Neogene and Quarternary shaping of the extant river basins of eastern Brazil. Despite the origin of theA.autranigroup being dated to the late Miocene, species level diversification occurred in the Pleistocene and was probably driven by headwater capture events and sea-level fluctuations.

Three new species of the ‘ Geophagus’ brasiliensis species group from the northeast Brazil (Cichlidae, Geophagini)

Morphological characters and phylogenetic trees generated by analyses of segments of two mitochondrial genes cytochrome b and cytochrome c oxidase I support recognition of three new species of the ‘Geophagus’ brasiliensis species group from coastal basins of northeast Brazil. All new species were diagnosed by exclusive morphological characters and exclusive nucleotide transformations. Geophagusrufomarginatus sp. n., from the Rio Buranhém Basin, is distinguished from all other species of the group by dorsal-fin lappets with red edges, the presence of longitudinal series of small light blue spots between the anal-fin spines and rays, and non-denticulated gill-rakers; it is closely related to G.brasiliensis and G.iporangensis. Geophagusmultiocellus sp. n., from the Rio de Contas Basin, is distinguished from all other species of the group by having small pale blue spots with minute bright blue dots at their centres, that are often vertically coalesced to form short bars on the caudal fin. Geophagussantosi sp. n., from the Rio Mariana Basin, is distinguished from all other species of the group by having blue stripes parallel to the dorsal and anal fin rays on their longest portions. Geophagusmultiocellus and G.santosi belong to the same clade of G.itapicuruensis. The clade composed by the Rio Paraguaçu Basin species was recovered as the sister group of the other species of the ‘G.’ brasiliensis species group.

A new species of Hyphessobrycon Durbin from northeastern Brazil: evidence from morphological data and DNA barcoding (Characiformes, Characidae)

A new species of Hyphessobrycon is described for the upper Munim and Preguiças river basins, northeastern Brazil, supported by morphological and molecular species delimitation methods. This new species belongs to the Hyphessobryconsensu stricto group, as it has the three main diagnostic character states of this assemblage: presence of a dark brown or black blotch on the dorsal fin, absence of a black midlateral stripe on its flank and the position of Weberian apparatus upward horizontal through dorsal margin of operculum. Our phylogenetic analysis also supported the allocation of the new species in this group; however, it was not possible to recover the species sister-group. Pristellamaxillaris and Moenkhausiahemigrammoides were recovered as the sister-clade of the Hyphessobryconsensu stricto group.

Hypothesis-Testing Demands Trustworthy Data-A Simulation Approach to Inferential Statistics Advocating the Research Program Strategy

In psychology as elsewhere, the main statistical inference strategy to establish empirical effects is null-hypothesis significance testing (NHST). The recent failure to replicate allegedly well-established NHST-results, however, implies that such results lack sufficient statistical power, and thus feature unacceptably high error-rates. Using data-simulation to estimate the error-rates of NHST-results, we advocate the research program strategy (RPS) as a superior methodology. RPS integrates Frequentist with Bayesian inference elements, and leads from a preliminary discovery against a (random) H₀-hypothesis to a statistical H₁-verification. Not only do RPS-results feature significantly lower error-rates than NHST-results, RPS also addresses key-deficits of a “pure” Frequentist and a standard Bayesian approach. In particular, RPS aggregates underpowered results safely. RPS therefore provides a tool to regain the trust the discipline had lost during the ongoing replicability-crisis.

The Effect of Nonreversibility on Inferring Rooted Phylogenies

Most phylogenetic models assume that the evolutionary process is stationary and reversible. In addition to being biologically improbable, these assumptions also impair inference by generating models under which the likelihood does not depend on the position of the root. Consequently, the root of the tree cannot be inferred as part of the analysis. Yet identifying the root position is a key component of phylogenetic inference because it provides a point of reference for polarizing ancestor-descendant relationships and therefore interpreting the tree. In this paper, we investigate the effect of relaxing the unrealistic reversibility assumption and allowing the position of the root to be another unknown. We propose two hierarchical models that are centered on a reversible model but perturbed to allow nonreversibility. The models differ in the degree of structure imposed on the perturbations. The analysis is performed in the Bayesian framework using Markov chain Monte Carlo methods for which software is provided. We illustrate the performance of the two nonreversible models in analyses of simulated data using two types of topological priors. We then apply the models to a real biological data set, the radiation of polyploid yeasts, for which there is robust biological opinion about the root position. Finally, we apply the models to a second biological alignment for which the rooted tree is controversial: the ribosomal tree of life. We compare the two nonreversible models and conclude that both are useful in inferring the position of the root from real biological data.

Tanglegrams: A Reduction Tool for Mathematical Phylogenetics

Many discrete mathematics problems in phylogenetics are defined in terms of the relative labeling of pairs of leaf-labeled trees. These relative labelings are naturally formalized as tanglegrams, which have previously been an object of study in coevolutionary analysis. Although there has been considerable work on planar drawings of tanglegrams, they have not been fully explored as combinatorial objects until recently. In this paper, we describe how many discrete mathematical questions on trees "factor" through a problem on tanglegrams, and how understanding that factoring can simplify analysis. Depending on the problem, it may be useful to consider a unordered version of tanglegrams, and/or their unrooted counterparts. For all of these definitions, we show how the isomorphism types of tanglegrams can be understood in terms of double cosets of the symmetric group, and we investigate their automorphisms. Understanding tanglegrams better will isolate the distinct problems on leaf-labeled pairs of trees and reveal natural symmetries of spaces associated with such problems.

Genome-wide analysis of the Brachypodium distachyon (L.) P. Beauv. Hsp90 gene family reveals molecular evolution and expression profiling under drought and salt stresses

The structure, evolution, and function of heat shock proteins 90 (Hsp90s) have been investigated in great detail in fungi and animals. However, studies on the Hsp90 genes in plants are generally limited. Brachypodium distachyon (L.) P. Beauv., as a model plant for cereal crops, has become a potential biofuel grass. During its long evolution, the Hsp90 gene family in Brachypodium has developed some strategies to cope with adverse environments. How the Hsp90 gene family in Brachypodium evolved in different plant lineages and what its role is in plant responses to drought and salt stresses remains to be elucidated. We used a set of different bioinformatics tools to identify 94 Hsp90 genes from 10 species representing four plant lineages and classified into three subgroups. Eight BdHsp90 genes were detected from B. distachyon. The number of exon-intron structures differed in each subgroup, and the motif analysis revealed that these genes were relatively conservative in each group. The fragments duplication and tandem duplication, which are the prime powers for functional diversity, generally occurred during the duplication of the whole plant genome. Transcriptional analysis of the BdHsp90 genes under salt and drought stress conditions indicated that the expression of these genes was delayed or increased at different stress time points; The expression was more affected in that of Bradi3g39630, Bradi4g06370, and Bradi1g30130. Our findings suggest the involvement of BdHsp90s in plant abiotic stress response, and further consolidate our views on the stress response mechanism of Hsp90 in general.

Phylogeography of the termite Macrotermes gilvus and insight into ancient dispersal corridors in Pleistocene Southeast Asia

Dispersal of soil-dwelling organisms via the repeatedly exposed Sunda shelf through much of the Pleistocene in Southeast Asia has not been studied extensively, especially for invertebrates. Here we investigated the phylogeography of an endemic termite species, Macrotermes gilvus (Hagen), to elucidate the spatiotemporal dynamics of dispersal routes of terrestrial fauna in Pleistocene Southeast Asia. We sampled 213 termite colonies from 66 localities throughout the region. Independently inherited microsatellites and mtDNA markers were used to infer the phylogeographic framework of M. gilvus. Discrete phylogeographic analysis and molecular dating based on fossil calibration were used to infer the dynamics of M. gilvus dispersal in time and space across Southeast Asia. We found that the termite dispersal events were consistently dated within the Pleistocene time frame. The dispersal pattern was multidirectional, radiating eastwards and southwards out of Indochina, which was identified as the origin for dispersal events. We found no direct dispersal events between Sumatra and Borneo despite the presence of a terrestrial connection between them during the Pleistocene. Instead, central Java served as an important link allowing termite colonies to be established in Borneo and Sumatra. Our findings support the hypothesis of a north-south dispersal corridor in Southeast Asia and suggest the presence of alternative dispersal routes across Sundaland during the Pleistocene. For the first time, we also propose that a west-east dispersal through over-water rafting likely occurred across the Pleistocene South China Sea. We found at least two independent entry routes for terrestrial species to infiltrate Sumatra and Borneo at different times.

Early-branching euteleost relationships: areas of congruence between concatenation and coalescent model inferences

Phylogenetic inference based on evidence from DNA sequences has led to significant strides in the development of a stable and robustly supported framework for the vertebrate tree of life. To date, the bulk of those advances have relied on sequence data from a small number of genome regions that have proven unable to produce satisfactory answers to consistently recalcitrant phylogenetic questions. Here, we re-examine phylogenetic relationships among early-branching euteleostean fish lineages classically grouped in the Protacanthopterygii using DNA sequence data surrounding ultraconserved elements. We report and examine a dataset of thirty-four OTUs with 17,957 aligned characters from fifty-three nuclear loci. Phylogenetic analysis is conducted in concatenated, joint gene trees and species tree estimation and summary coalescent frameworks. All analytical frameworks yield supporting evidence for existing hypotheses of relationship for the placement of Lepidogalaxias salamandroides, monophyly of the Stomiatii and the presence of an esociform + salmonid clade. Lepidogalaxias salamandroides and the Esociformes + Salmoniformes are successive sister lineages to all other euteleosts in the majority of analyses. The concatenated and joint gene trees and species tree analysis types produce high support values for this arrangement. However, inter-relationships of Argentiniformes, Stomiatii and Neoteleostei remain uncertain as they varied by analysis type while receiving strong and contradictory indices of support. Topological differences between analysis types are also apparent within the otomorph and the percomorph taxa in the data set. Our results identify concordant areas with strong support for relationships within and between early-branching euteleost lineages but they also reveal limitations in the ability of larger datasets to conclusively resolve other aspects of that phylogeny.

Diversity of culturable filamentous Ascomycetes in the eastern South Pacific Ocean off Chile

Our study reports the diversity of culturable mycoplankton in the eastern South Pacific Ocean off Chile to contribute with novel knowledge on taxonomy of filamentous fungi isolated from distinct physicochemical and biological marine environments. We characterized spatial distribution of isolates, evaluated their viability and assessed the influence of organic substrate availability on fungal development. Thirty-nine Operational Taxonomic Units were identified from 99 fungal strains isolated from coastal and oceanic waters by using Automatic Barcode Gap Discovery. All Operational Taxonomic Units belonged to phylum Ascomycota and orders Eurotiales, Dothideales, Sordariales and Hypocreales, mainly Penicillium sp. (82%); 11 sequences did not match existing species in GenBank, suggesting occurrence of novel fungal taxa. Our results suggest that fungal communities in the South Pacific Ocean off Chile appear to thrive in a wide range of environmental conditions in the ocean and that substrate availability may be a factor influencing fungal viability in the ocean.

Peste des petits ruminants (PPR): A neglected tropical disease in Maghreb region of North Africa and its threat to Europe

BACKGROUND

Peste des petits ruminants (PPR) is a contagious disease listed by the World Organisation for Animal health (OIE) as being a specific hazard. It affects sheep, goats, and wild ungulates, and is prevalent throughout the developing world particularly Asia, the Middle East, and Africa. PPR has been targeted for eradication by 2030 by the Food and Agriculture Organization of the United Nations (FAO) and the OIE, after the successful eradication of the related disease, rinderpest in cattle. PPR was first reported in 1942 in the Ivory Coast in Western Africa and has since extended its range in Asia, the Middle East, and Africa posing an immediate threat of incursion into Europe, South East Asia and South Africa. Although robust vaccines are available, the use of these vaccines in a systematic and rational manner is not widespread, resulting in this devastating disease becoming an important neglected tropical disease in the developing world.

METHODOLOGY

We isolated and characterized the PPR virus from an outbreak in Cheraga, northern Algeria, during October 2015 by analyzing the partial N-gene sequence in comparison with other viruses from the Maghreb region. As well as sequencing the full length viral genome and performing real-time RT-PCR on clinical samples. Maximum-likelihood and Bayesian temporal and phylogeographic analyses were performed to assess the persistence and spread of PPRV circulation from Eastern Africa in the Maghreb region of North Africa.

CONCLUSIONS

Recent PPR outbreaks in Cheraga, in the northern part of Algiers (October 2015) and North-West Morocco (June, 2015) highlight that PPRV has spread to the northern border of North Africa and may pose a threat of introduction to Europe. Phylogeographic analysis suggests that lineage IV PPRV has spread from Eastern Africa, most likely from the Sudan 2000 outbreak, into Northern Africa resulting in the 2008 Moroccan outbreak. Maximum-likelihood and Bayesian analysis shows that these North African viruses cluster closely together suggesting the existence of continual regional circulation. Considering the same virus is circulating in Algeria, Morocco and Tunisia, implementation of a common Maghreb PPR eradication strategy would be beneficial for the region.

Inferring Trees

Molecular evolution can reveal the relationship between sets of homologous sequences and the patterns of change that occur during their evolution. An important aspect of these studies is the inference of a phylogenetic tree, which explicitly describes evolutionary relationships between homologous sequences. This chapter provides an introduction to evolutionary trees and how to infer them from sequence data using some commonly used inferential methodology. It focuses on statistical methods for inferring trees and how to assess the confidence one should have in any resulting tree, with a particular emphasis on the underlying assumptions of the methods and how they might affect the tree estimate. There is also some discussion of the underlying algorithms used to perform tree search and recommendations regarding the performance of different algorithms. Finally, there are a few practical guidelines, including how to combine multiple software packages to improve inference, and a comparison between Bayesian and Maximum likelihood phylogenetics.

Out of southern Africa: Origin, biogeography and age of the Aizooideae (Aizoaceae)

The Aizooideae is an early-diverging lineage within the Aizoaceae. It is most diverse in southern Africa, but also has endemic species in Australasia, Eurasia and South America. We derived a phylogenetic hypothesis from Bayesian and Maximum Likelihood analyses of plastid DNA-sequences. We find that one of the seven genera, the fynbos-endemic Acrosanthes, does not belong to the Aizooideae, but is an ancient sister-lineage to the subfamilies Mesembryanthemoideae & Ruschioideae. Galenia and Plinthus are embedded inside Aizoon and Aizoanthemum is polyphyletic. The Namibian endemic Tetragonia schenckii is sister to Tribulocarpus of the Sesuvioideae. For the Aizooideae, we explored their possible age by means of relaxed Bayesian dating and used Bayesian Binary MCMC reconstruction of ancestral areas to investigate their area of origin. Early diversification occurred in southern Africa in the Eocene-Oligocene, with a split into a mainly African lineage and an Eurasian-Australasian-African-South American lineage. These subsequently radiated in the early Miocene. For Tetragonia, colonisation of Australasia via long-distance dispersal from Eurasia gave rise to the Australasian lineage from which there were subsequent dispersals to South America and Southern Africa. Despite the relatively old age of the Aizooideae, more than half the species have radiated since the Pleiocene, coinciding with the large and rapid diversification of the Ruschioideae. The lineage made up of Tetragonia schenckii &Tribulocarpus split from the remainder of the Sesuvioideae already in the mid Oligocene and its disjunct distribution between Namibia and north-east Africa may be the result of a previously wider distribution within an early Arid African flora. Our reconstruction of ancestral character-states indicates that the expanding keels giving rise to hygrochastic fruits originated only once, i.e. after the split of the Sesuvioideae from the remainder of the Aizoaceae and that they were subsequently lost many times. Variously winged and spiky fruits, adapted to dispersal by wind and animals, have evolved independently in the Aizooideae and the Sesuvioideae. There is then a greater diversity of dispersal systems in the earlier lineages than in the Mesembryanthemoideae and Ruschioideae, where dispersal is mainly achieved by rain.

Resolving Cypriniformes relationships using an anchored enrichment approach

BACKGROUND

Cypriniformes (minnows, carps, loaches, and suckers) is the largest group of freshwater fishes in the world (~4300 described species). Despite much attention, previous attempts to elucidate relationships using molecular and morphological characters have been incongruent. In this study we present the first phylogenomic analysis using anchored hybrid enrichment for 172 taxa to represent the order (plus three out-group taxa), which is the largest dataset for the order to date (219 loci, 315,288 bp, average locus length of 1011 bp).

RESULTS

Concatenation analysis establishes a robust tree with 97 % of nodes at 100 % bootstrap support. Species tree analysis was highly congruent with the concatenation analysis with only two major differences: monophyly of Cobitoidei and placement of Danionidae.

CONCLUSIONS

Most major clades obtained in prior molecular studies were validated as monophyletic, and we provide robust resolution for the relationships among these clades for the first time. These relationships can be used as a framework for addressing a variety of evolutionary questions (e.g. phylogeography, polyploidization, diversification, trait evolution, comparative genomics) for which Cypriniformes is ideally suited.

Origin and diversification of living cycads: a cautionary tale on the impact of the branching process prior in Bayesian molecular dating

Background

Bayesian relaxed-clock dating has significantly influenced our understanding of the timeline of biotic evolution. This approach requires the use of priors on the branching process, yet little is known about their impact on divergence time estimates. We investigated the effect of branching priors using the iconic cycads. We conducted phylogenetic estimations for 237 cycad species using three genes and two calibration strategies incorporating up to six fossil constraints to (i) test the impact of two different branching process priors on age estimates, (ii) assess which branching prior better fits the data, (iii) investigate branching prior impacts on diversification analyses, and (iv) provide insights into the diversification history of cycads.

Results

Using Bayes factors, we compared divergence time estimates and the inferred dynamics of diversification when using Yule versus birth-death priors. Bayes factors were calculated with marginal likelihood estimated with stepping-stone sampling. We found striking differences in age estimates and diversification dynamics depending on prior choice. Dating with the Yule prior suggested that extant cycad genera diversified in the Paleogene and with two diversification rate shifts. In contrast, dating with the birth-death prior yielded Neogene diversifications, and four rate shifts, one for each of the four richest genera. Nonetheless, dating with the two priors provided similar age estimates for the divergence of cycads from Ginkgo (Carboniferous) and their crown age (Permian). Of these, Bayes factors clearly supported the birth-death prior.

Conclusions

These results suggest the choice of the branching process prior can have a drastic influence on our understanding of evolutionary radiations. Therefore, all dating analyses must involve a model selection process using Bayes factors to select between a Yule or birth-death prior, in particular on ancient clades with a potential pattern of high extinction. We also provide new insights into the history of cycad diversification because we found (i) periods of extinction along the long branches of the genera consistent with fossil data, and (ii) high diversification rates within the Miocene genus radiations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0347-8) contains supplementary material, which is available to authorized users.

Diversification and the evolution of dispersal ability in the tribe Brassiceae (Brassicaceae)

BACKGROUND AND AIMS

Dispersal and establishment ability can influence evolutionary processes such as geographic isolation, adaptive divergence and extinction probability. Through these population-level dynamics, dispersal ability may also influence macro-evolutionary processes such as species distributions and diversification. This study examined patterns of evolution of dispersal-related fruit traits, and how the evolution of these traits is correlated with shifts in geographic range size, habitat and diversification rates in the tribe Brassiceae (Brassicaceae).

METHODS

The phylogenetic analysis included 72 taxa sampled from across the Brassiceae and included both nuclear and chloroplast markers. Dispersal-related fruit characters were scored and climate information for each taxon was retrieved from a database. Correlations between fruit traits, seed characters, habitat, range and climate were determined, together with trait-dependent diversification rates.

KEY RESULTS

It was found that the evolution of traits associated with limited dispersal evolved only in association with compensatory traits that increase dispersal ability. The evolution of increased dispersal ability occurred in multiple ways through the correlated evolution of different combinations of fruit traits. The evolution of traits that increase dispersal ability was in turn associated with larger seed size, increased geographic range size and higher diversification rates.

CONCLUSIONS

This study provides evidence that the evolution of increased dispersal ability and larger seed size, which may increase establishment ability, can also influence macro-evolutionary processes, possibly by increasing the propensity for long-distance dispersal. In particular, it may increase speciation and consequent diversification rates by increasing the likelihood of geographic and thereby reproductive isolation.

Genetic and epidemiological insights into the emergence of peste des petits ruminants virus (PPRV) across Asia and Africa

Small ruminants are important components in the livelihood of millions of households in many parts of the world. The spread of the highly contagious peste des petits ruminants (PPR) disease, which is caused by an RNA virus, PPRV, across Asia and Africa remains a major concern. The present study explored the evolutionary and epidemiological dynamics of PPRV through the analyses of partial N-gene and F-gene sequences of the virus. All the four previously described PPRV lineages (I-IV) diverged from their common ancestor during the late-19(th) to early-20(th) century. Among the four lineages, PPRV-IV showed pronounced genetic structuring across the region; however, haplotype sharing among the geographic regions, together with the presence of multiple genetic clusters within a country, indicates the possibility of frequent mobility of the diseased individuals across the region. The gradual decline in the effective number of infections suggests a limited genetic variation, which could be attributed to the effective vaccination that has been practiced since 1990s. However, the movement of infected animals across the region likely contributes to the spread of PPRV-IV. No evidence of positive selection was identified from this study.

Erosive processes after tectonic uplift stimulate vicariant and adaptive speciation: evolution in an Afrotemperate-endemic paper daisy genus

BACKGROUND

The role of tectonic uplift in stimulating speciation in South Africa's only alpine zone, the Drakensberg, has not been explicitly examined. Tectonic processes may influence speciation both through the creation of novel habitats and by physically isolating plant populations. We use the Afrotemperate endemic daisy genus Macowania to explore the timing and mode (geographic versus adaptive) of speciation in this region. Between sister species pairs we expect high morphological divergence where speciation has happened in sympatry (adaptive) while with geographic (vicariant) speciation we may expect to find less morphological divergence and a greater degree of allopatry. A dated molecular phylogenetic hypothesis for Macowania elucidates species' relationships and is used to address the potential impact of uplift on diversification. Morphological divergence of a small sample of reproductive and vegetative characters, used as a proxy for adaptive divergence, is measured against species' range distributions to estimate mode of speciation across two subclades in the genus.

RESULTS

The Macowania crown age is consistent with the hypothesis of post-uplift diversification, and we find evidence for both vicariant and adaptive speciation between the two subclades within Macowania. Both subclades exhibit strong signals of range allopatry, suggesting that geographic isolation was important in speciation. One subclade, associated with dry, rocky environments at high altitudes, shows very little morphological and ecological differentiation but high range allopatry. The other subclade occupies a greater variety of habitats and exhibits far greater morphological differentiation, but contains species with overlapping distribution ranges.

CONCLUSIONS

Species in Macowania are likely to have diversified in response to tectonic uplift, and we invoke uplift and uplift-mediated erosion as the main drivers of speciation. The greater relative morphological divergence in sympatric species of Macowania indicates that speciation in the non-sympatric taxa may not have required obvious adaptive differences, implying that simple geographic isolation was the driving force for speciation ('neutral speciation').