The evolution of multiple colour mechanisms is correlated with diversification in sunbirds (Nectariniidae).

How and why certain groups become speciose is a key question in evolutionary biology. Novel traits that enable diversification by opening new ecological niches are likely important mechanisms. However, ornamental traits can also promote diversification by opening up novel sensory niches and thereby creating novel inter-specific interactions. More specifically, ornamental colours may enable more precise and/or easier species recognition, and may act as key innovations by increasing the number of species-specific patterns and promoting diversification. While the influence of colouration on diversification is well-studied, the influence of the mechanisms that produce those colours (e.g. pigmentary, nanostructural) is less so, even though the ontogeny and evolution of these mechanisms differ. We estimated a new phylogenetic tree for 121 sunbird species and combined colour data of 106 species with a range of phylogenetic tools to test the hypothesis that the evolution of novel colour mechanisms increases diversification in sunbirds, one of the most colourful bird clades. Results suggest that (1) the evolution of novel colour mechanisms expands the visual sensory niche, increasing the number of achievable colours. (2) Structural colouration diverges more readily across the body than pigment-based colouration, enabling an increase in colour complexity. (3) Novel colour mechanisms might minimize trade-offs between natural and sexual selection such that colour can function both as camouflage and conspicuous signal. (4) Despite structural colours being more colourful and mobile, only melanin-based colouration is positively correlated with net diversification. Together, these findings explain why colour distances increase with increasing number of sympatric species, even though packing of colour space would predict otherwise.

Ancestral reconstruction of the MotA stator subunit reveals that conserved residues far from the pore are required to drive flagellar motility.

The bacterial flagellar motor (BFM) is a rotary nanomachine powered by the translocation of ions across the inner membrane through the stator complex. The stator complex consists of two membrane proteins: MotA and MotB (in H+-powered motors), or PomA and PomB (in Na+-powered motors). In this study, we used ancestral sequence reconstruction (ASR) to probe which residues of MotA correlate with function and may have been conserved to preserve motor function. We reconstructed 10 ancestral sequences of MotA and found four of them were motile in combination with contemporary Escherichia coli MotB and in combination with our previously published functional ancestral MotBs. Sequence comparison between wild-type (WT) E. coli MotA and MotA-ASRs revealed 30 critical residues across multiple domains of MotA that were conserved among all motile stator units. These conserved residues included pore-facing, cytoplasm-facing, and MotA-MotA intermolecular facing sites. Overall, this work demonstrates the role of ASR in assessing conserved variable residues in a subunit of a molecular complex.

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations.

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale.

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNAIle, ND1, tRNALeu2, 16S, tRNAVal, 12S, tRNAPhe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae).

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4-10x, Veronica: 12-18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.

Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales.

Evolutionary models account for either population- or species-level processes but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth-death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology.[Caninae; fossilized birth-death; molecular clock; multispecies coalescent; phylogenetics; species trees.].

Flagellar export apparatus and ATP synthetase: Homology evidenced by synteny predating the Last Universal Common Ancestor.

We report evidence further supporting homology between proteins in the F1 FO -ATP synthetase and the bacterial flagellar motor (BFM). BFM proteins FliH, FliI, and FliJ have been hypothesized to be homologous to FO -b + F1 -δ, F1 -α/ß, and F1 -γ, with similar structure and interactions. We conduct a further test by constructing a gene order dataset, examining the order of fliH, fliI, and fliJ genes across the phylogenetic breadth of flagellar and nonflagellar type 3 secretion systems, and comparing this to published surveys of gene order in the F1 FO -ATP synthetase, its N-ATPase relatives, and the bacterial/archaeal V- and A-type ATPases. Strikingly, the fliHIJ gene order was deeply conserved, with the few exceptions appearing derived, and exactly matching the widely conserved F-ATPase gene order atpFHAG, coding for subunits b-δ-α-γ. The V/A-type ATPases have a similar conserved gene order. Our results confirm homology between these systems, and suggest a rare case of synteny conserved over billions of years, predating the Last Universal Common Ancestor (LUCA).

Corrigendum: Ancestral Sequence Reconstructions of MotB Are Proton-Motile and Require MotA for Motility.

[This corrects the article DOI: 10.3389/fmicb.2020.625837.].

Trait-dependent dispersal in rails (Aves: Rallidae): Historical biogeography of a cosmopolitan bird clade.

The ability of lineages to disperse over evolutionary timescales may be influenced by the gain or loss of traits after adaptation to new ecological conditions. For example, rails (Aves: Rallidae) have many cases of flightless insular endemic species that presumably evolved after flying ancestors dispersed over large ocean barriers and became isolated. Nonetheless, the details of how flying and its loss have influenced the clade's historical biogeography are unknown, as is the importance of other predictors of dispersal such as the geographic distance between regions. Here, we used a dated phylogeny of 158 species of rails to compare trait-dependent and trait-independent biogeography models in BioGeoBEARS. We evaluated a probabilistic historical biogeographical model that allows geographic range and flight to co-evolve and influence dispersal ability on a phylogeny. The best-fitting dispersal model was a trait-dependent dispersal (DEC + j + x + t21 + m1) that accrued 85.2% of the corrected Akaike Information Criterion (AICc) model weight. The distance-dependence parameter, x was estimated at -0.54, ranging from -0.49 to -0.65 across models, suggesting that a doubling of dispersal distance results in an approximately 31% decrease in dispersal rate (2-0.54 = 0.69). The estimated rate of loss of flight (t21) was similar across all models (~0.029 loss events per lineage per million years). The multiplier on dispersal rate when a lineage is non-flying, m1, is estimated to be 0.38 under this model. Surprisingly, the estimate of m1 was not 0.0, probably because the loss of flight is so common in the rails that entire clades of flightless species are found in the data, forcing the model to attribute some dispersal to flightless lineages. These results indicate that long-distance dispersal over macroevolutionary timespans can be modelled, rather than simply attributed to chance, allowing support for different hypotheses to be quantified and limitations to be identified. Overall, by combining new analytical methods with a comprehensive phylogeny, we use a quantitative framework to show how traits influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Ancestral Sequence Reconstructions of MotB Are Proton-Motile and Require MotA for Motility.

The bacterial flagellar motor (BFM) is a nanomachine that rotates the flagellum to propel many known bacteria. The BFM is powered by ion transit across the cell membrane through the stator complex, a membrane protein. Different bacteria use various ions to run their BFM, but the majority of BFMs are powered by either proton (H+) or sodium (Na+) ions. The transmembrane (TM) domain of the B-subunit of the stator complex is crucial for ion selectivity, as it forms the ion channel in complex with TM3 and TM4 of the A-subunit. In this study, we reconstructed and engineered thirteen ancestral sequences of the stator B-subunit to evaluate the functional properties and ionic power source of the stator proteins at reconstruction nodes to evaluate the potential of ancestral sequence reconstruction (ASR) methods for stator engineering and to test specific motifs previously hypothesized to be involved in ion-selectivity. We found that all thirteen of our reconstructed ancient B-subunit proteins could assemble into functional stator complexes in combination with the contemporary Escherichia coli MotA-subunit to restore motility in stator deleted E. coli strains. The flagellar rotation of the thirteen ancestral MotBs was found to be Na+ independent which suggested that the F30/Y30 residue was not significantly correlated with sodium/proton phenotype, in contrast to what we had reported previously. Additionally, four among the thirteen reconstructed B-subunits were compatible with the A-subunit of Aquifex aeolicus and able to function in a sodium-independent manner. Overall, this work demonstrates the use of ancestral reconstruction to generate novel stators and quantify which residues are correlated with which ionic power source.

Statistical Comparison of Trait-Dependent Biogeographical Models Indicates That Podocarpaceae Dispersal Is Influenced by Both Seed Cone Traits and Geographical Distance.

The ability of lineages to disperse long distances over evolutionary timescales may be influenced by the gain or loss of traits adapted to enhance local, ecological dispersal. For example, some species in the southern conifer family Podocarpaceae have fleshy cones that encourage bird dispersal, but it is unknown how this trait has influenced the clade's historical biogeography, or its importance compared with other predictors of dispersal such as the geographic distance between regions. We answer these questions quantitatively by using a dated phylogeny of 197 species of southern conifers (Podocarpaceae and their sister family Araucariaceae) to statistically compare standard, trait-independent biogeography models with new BioGeoBEARS models where an evolving trait can influence dispersal probability, and trait history, biogeographical history, and model parameters are jointly inferred. We validate the method with simulation-inference experiments. Comparing all models, those that include trait-dependent dispersal accrue 87.5% of the corrected Akaike Information Criterion (AICc) model weight. Averaged across all models, lineages with nonfleshy cones had a dispersal probability multiplier of 0.49 compared with lineages with fleshy cones. Distance is included as a predictor of dispersal in all credible models (100% model weight). However, models with changing geography earned only 22.0% of the model weight, and models submerging New Caledonia/New Zealand earned only 0.01%. The importance of traits and distance suggests that long-distance dispersal over macroevolutionary timespans should not be thought of as a highly unpredictable chance event. Instead, long-distance dispersal can be modeled, allowing statistical model comparison to quantify support for different hypotheses.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate.

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Sodium-powered stators of the bacterial flagellar motor can generate torque in the presence of phenamil with mutations near the peptidoglycan-binding region.

The bacterial flagellar motor powers the rotation that propels the swimming bacteria. Rotational torque is generated by harnessing the flow of ions through ion channels known as stators which couple the energy from the ion gradient across the inner membrane to rotation of the rotor. Here, we used error-prone PCR to introduce single point mutations into the sodium-powered Vibrio alginolyticus/Escherichia coli chimeric stator PotB and selected for motors that exhibited motility in the presence of the sodium-channel inhibitor phenamil. We found single mutations that enable motility under phenamil occurred at two sites: (i) the transmembrane domain of PotB, corresponding to the TM region of the PomB stator from V. alginolyticus and (ii) near the peptidoglycan binding region that corresponds to the C-terminal region of the MotB stator from E. coli. Single cell rotation assays confirmed that individual flagellar motors could rotate in up to 100 µM phenamil. Using phylogenetic logistic regression, we found correlation between natural residue variation and ion source at positions corresponding to PotB F22Y, but not at other sites. Our results demonstrate that it is not only the pore region of the stator that moderates motility in the presence of ion-channel blockers.

Taxon cycle predictions supported by model-based inference in Indo-Pacific trap-jaw ants (Hymenoptera: Formicidae: Odontomachus).

Nonequilibrium dynamics and non-neutral processes, such as trait-dependent dispersal, are often missing from quantitative island biogeography models despite their potential explanatory value. One of the most influential nonequilibrium models is the taxon cycle, but it has been difficult to test its validity as a general biogeographical framework. Here, we test predictions of the taxon cycle model using six expected phylogenetic patterns and a time-calibrated phylogeny of Indo-Pacific Odontomachus (Hymenoptera: Formicidae: Ponerinae), one of the ant genera that E.O. Wilson used when first proposing the hypothesis. We used model-based inference and a newly developed trait-dependent dispersal model to jointly estimate ancestral biogeography, ecology (habitat preferences for forest interiors, vs. "marginal" habitats, such as savannahs, shorelines, disturbed areas) and the linkage between ecology and dispersal rates. We found strong evidence that habitat shifts from forest interior to open and disturbed habitats increased macroevolutionary dispersal rate. In addition, lineages occupying open and disturbed habitats can give rise to both island endemics re-occupying only forest interiors and taxa that re-expand geographical ranges. The phylogenetic predictions outlined in this study can be used in future work to evaluate the relative weights of neutral (e.g., geographical distance and area) and non-neutral (e.g., trait-dependent dispersal) processes in historical biogeography and community ecology.

Including autapomorphies is important for paleontological tip-dating with clocklike data, but not with non-clock data.

Tip-dating, where fossils are included as dated terminal taxa in Bayesian dating inference, is an increasingly popular method. Data for these studies often come from morphological character matrices originally developed for non-dated, and usually parsimony, analyses. In parsimony, only shared derived characters (synapomorphies) provide grouping information, so many character matrices have an ascertainment bias: they omit autapomorphies (unique derived character states), which are considered uninformative. There has been no study of the effect of this ascertainment bias in tip-dating, but autapomorphies can be informative in model-based inference. We expected that excluding autapomorphies would shorten the morphological branchlengths of terminal branches, and thus bias downwards the time branchlengths inferred in tip-dating. We tested for this effect using a matrix for Carboniferous-Permian eureptiles where all autapomorphies had been deliberately coded. Surprisingly, date estimates are virtually unchanged when autapomorphies are excluded, although we find large changes in morphological rate estimates and small effects on topological and dating confidence. We hypothesized that the puzzling lack of effect on dating was caused by the non-clock nature of the eureptile data. We confirm this explanation by simulating strict clock and non-clock datasets, showing that autapomorphy exclusion biases dating only for the clocklike case. A theoretical solution to ascertainment bias is computing the ascertainment bias correction (Mkparsinf), but we explore this correction in detail, and show that it is computationally impractical for typical datasets with many character states and taxa. Therefore we recommend that palaeontologists collect autapomorphies whenever possible when assembling character matrices.

Dating the Species Network: Allopolyploidy and Repetitive DNA Evolution in American Daisies (Melampodium sect. Melampodium, Asteraceae).

Allopolyploidy has played an important role in the evolution of the flowering plants. Genome mergers are often accompanied by significant and rapid alterations of genome size and structure via chromosomal rearrangements and altered dynamics of tandem and dispersed repetitive DNA families. Recent developments in sequencing technologies and bioinformatic methods allow for a comprehensive investigation of the repetitive component of plant genomes. Interpretation of evolutionary dynamics following allopolyploidization requires both the knowledge of parentage and the age of origin of an allopolyploid. Whereas parentage is typically inferred from cytogenetic and phylogenetic data, age inference is hampered by the reticulate nature of the phylogenetic relationships. Treating subgenomes of allopolyploids as if they belonged to different species (i.e., no recombination among subgenomes) and applying cross-bracing (i.e., putting a constraint on the age difference of nodes pertaining to the same event), we can infer the age of allopolyploids within the framework of the multispecies coalescent within BEAST2. Together with a comprehensive characterization of the repetitive DNA fraction using the RepeatExplorer pipeline, we apply the dating approach in a group of closely related allopolyploids and their progenitor species in the plant genus Melampodium (Asteraceae). We dated the origin of both the allotetraploid, Melampodium strigosum, and its two allohexaploid derivatives, Melampodium pringlei and Melampodium sericeum, which share both parentage and the direction of the cross, to the Pleistocene ($<$1.4 Ma). Thus, Pleistocene climatic fluctuations may have triggered formation of allopolyploids possibly in short intervals, contributing to difficulties in inferring the precise temporal order of allopolyploid species divergence of M. sericeum and M. pringlei. The relatively recent origin of the allopolyploids likely played a role in the near-absence of major changes in the repetitive fraction of the polyploids' genomes. The repetitive elements most affected by the postpolyploidization changes represented retrotransposons of the Ty1-copia lineage Maximus and, to a lesser extent, also Athila elements of Ty3-gypsy family.

Recent origin and rapid speciation of Neotropical orchids in the world's richest plant biodiversity hotspot.

The Andean mountains of South America are the most species-rich biodiversity hotspot worldwide with c. 15% of the world's plant species, in only 1% of the world's land surface. Orchids are a key element of the Andean flora, and one of the most prominent components of the Neotropical epiphyte diversity, yet very little is known about their origin and diversification. We address this knowledge gap by inferring the biogeographical history and diversification dynamics of the two largest Neotropical orchid groups (Cymbidieae and Pleurothallidinae), using two unparalleled, densely sampled orchid phylogenies (including more than 400 newly generated DNA sequences), comparative phylogenetic methods, geological and biological datasets. We find that the majority of Andean orchid lineages only originated in the last 20-15 million yr. Andean lineages are derived from lowland Amazonian ancestors, with additional contributions from Central America and the Antilles. Species diversification is correlated with Andean orogeny, and multiple migrations and recolonizations across the Andes indicate that mountains do not constrain orchid dispersal over long timescales. Our study sheds new light on the timing and geography of a major Neotropical diversification, and suggests that mountain uplift promotes species diversification across all elevational zones.

Empirical and Bayesian approaches to fossil-only divergence times: A study across three reptile clades.

Estimating divergence times on phylogenies is critical in paleontological and neontological studies. Chronostratigraphically-constrained fossils are the only direct evidence of absolute timing of species divergence. Strict temporal calibration of fossil-only phylogenies provides minimum divergence estimates, and various methods have been proposed to estimate divergences beyond these minimum values. We explore the utility of simultaneous estimation of tree topology and divergence times using BEAST tip-dating on datasets consisting only of fossils by using relaxed morphological clocks and birth-death tree priors that include serial sampling (BDSS) at a constant rate through time. We compare BEAST results to those from the traditional maximum parsimony (MP) and undated Bayesian inference (BI) methods. Three overlapping datasets were used that span 250 million years of archosauromorph evolution leading to crocodylians. The first dataset focuses on early Sauria (31 taxa, 240 chars.), the second on early Archosauria (76 taxa, 400 chars.) and the third on Crocodyliformes (101 taxa, 340 chars.). For each dataset three time-calibrated trees (timetrees) were calculated: a minimum-age timetree with node ages based on earliest occurrences in the fossil record; a 'smoothed' timetree using a range of time added to the root that is then averaged over zero-length internodes; and a tip-dated timetree. Comparisons within datasets show that the smoothed and tip-dated timetrees provide similar estimates. Only near the root node do BEAST estimates fall outside the smoothed timetree range. The BEAST model is not able to overcome limited sampling to correctly estimate divergences considerably older than sampled fossil occurrence dates. Conversely, the smoothed timetrees consistently provide node-ages far older than the strict dates or BEAST estimates for morphologically conservative sister-taxa when they sit on long ghost lineages. In this latter case, the relaxed-clock model appears to be correctly moderating the node-age estimate based on the limited morphological divergence. Topologies are generally similar across analyses, but BEAST trees for crocodyliforms differ when clades are deeply nested but contain very old taxa. It appears that the constant-rate sampling assumption of the BDSS tree prior influences topology inference by disfavoring long, unsampled branches.

Model selection in statistical historical biogeography of Neotropical insects-The Exophthalmus genus complex (Curculionidae: Entiminae).

Statistical historical biogeographic methods rely on models that represent various biogeographic processes. Until recently model selection in this domain was not widely used, and the impact of differential model selection on inferring biogeographic scenarios was not well understood. Focusing on Neotropical weevils in the Exophthalmus genus complex (EGC) (Insecta: Curculionidae: Entiminae), we compare three commonly used biogeographic models - DIVA (Dispersal-Vicariance Analysis), DEC (Dispersal-Extinction-Cladogenesis) and BayArea (Bayesian Analysis of Biogeography), and examine the impact of modeling founder-event jump dispersal on historical biogeographic reconstructions. We also investigate the biogeographic events that have shaped patterns of distribution, diversification, and endemism in this weevil lineage. We sample representatives of 65 species of the EGC and 26 outgroup terminals from the Neotropics, including Caribbean islands and the mainland. We reconstruct a molecular phylogeny based on six genes and apply molecular dating using a relaxed clock with three fossil calibration points. Historical biogeographic estimations and alternative biogeographic models are computed and compared with the R package BioGeoBEARS. Model selection strongly favors biogeographic models that include founder-event jump dispersal. Without modeling jump dispersal, estimations based on the three biogeographic models are dramatically different, especially for early-diverging nodes. When jump dispersal is included, the three biogeographic models perform similarly. Accordingly, we show that the Neotropical mainland was colonized by Caribbean species in the early Miocene, and that in situ diversification accounts for a majority (∼75%) of the biogeographic events in the EGC. Our study highlights the need to assess wide-ranging historical biogeographic processes - including founder-event jump dispersal - for best-fitting statistical Caribbean biogeographic reconstructions. Moreover, colonization of the Neotropical mainland from the Caribbean reinforces the notion that islands can be an important source of continental diversity.