Py PHLAWD : A python tool for phylogenetic dataset construction

The diversity of social complexity in termites

Sociality underpins major evolutionary transitions and significantly influences the structure and function of complex ecosystems. Social insects, seen as the pinnacle of sociality, have traits like obligate sterility that are considered 'master traits', used as single phenotypic measures of this complexity. However, evidence is mounting that completely aligning both phenotypic and evolutionary social complexity, and having obligate sterility central to both, is erroneous. We hypothesize that obligate and functional sterility are insufficient in explaining the diversity of phenotypic social complexity in social insects. To test this, we explore the relative importance of these sterility traits in an understudied but diverse taxon: the termites. We compile the largest termite social complexity dataset to date, using specimen and literature data. We find that although functional and obligate sterility explain a significant proportion of variance, neither trait is an adequate singular proxy for the phenotypic social complexity of termites. Further, we show both traits have only a weak association with the other social complexity traits within termites. These findings have ramifications for our general comprehension of the frameworks of phenotypic and evolutionary social complexity, and their relationship with sterility.

Towards a species-level phylogeny for Neotropical Myrtaceae: Notes on topology and resources for future studies

PREMISE

Increasingly complete phylogenies underpin studies in systematics, ecology, and evolution. Myrteae (Myrtaceae), with ~2700 species, is a key component of the exceptionally diverse Neotropical flora, but given its complicated taxonomy, automated assembling of molecular supermatrices from public databases often lead to unreliable topologies due to poor species identification.

METHODS

Here, we build a taxonomically verified molecular supermatrix of Neotropical Myrteae by assembling 3909 published and 1004 unpublished sequences from two nuclear and seven plastid molecular markers. We infer a time-calibrated phylogenetic tree that covers 712 species of Myrteae (~28% of the total diversity in the clade) and evaluate geographic and taxonomic gaps in sampling.

RESULTS

The tree inferred from the fully concatenated matrix mostly reflects the topology of the plastid data set and there is a moderate to strong incongruence between trees inferred from nuclear and plastid partitions. Large, species-rich genera are still the poorest sampled within the group. Eastern South America is the best-represented area in proportion to its species diversity, while Western Amazon, Mesoamerica, and the Caribbean are the least represented.

CONCLUSIONS

We provide a time-calibrated tree that can be more reliably used to address finer-scale eco-evolutionary questions that involve this group in the Neotropics. Gaps to be filled by future studies include improving representation of taxa and areas that remain poorly sampled, investigating causes of conflict between nuclear and plastid partitions, and the role of hybridization and incomplete lineage sorting in relationships that are poorly supported.

Short-wavelength-sensitive 1 (SWS1) opsin gene duplications and parallel visual pigment tuning support ultraviolet communication in damselfishes (Pomacentridae)

Damselfishes (Pomacentridae) are one of the most behaviourally diverse, colourful and species-rich reef fish families. One remarkable characteristic of damselfishes is their communication in ultraviolet (UV) light. Not only are they sensitive to UV, they are also prone to have UV-reflective colours and patterns enabling social signalling. Using more than 50 species, we aimed to uncover the evolutionary history of UV colour and UV vision in damselfishes. All damselfishes had UV-transmitting lenses, expressed the UV-sensitive SWS1 opsin gene, and most displayed UV-reflective patterns and colours. We find evidence for several tuning events across the radiation, and while SWS1 gene duplications are generally very rare among teleosts, our phylogenetic reconstructions uncovered two independent duplication events: one close to the base of the most species-rich clade in the subfamily Pomacentrinae, and one in a single Chromis species. Using amino acid comparisons, we found that known spectral tuning sites were altered several times in parallel across the damselfish radiation (through sequence change and duplication followed by sequence change), causing repeated shifts in peak spectral absorbance of around 10 nm. Pomacentrinae damselfishes expressed either one or both copies of SWS1, likely to further finetune UV-signal detection and differentiation. This highly advanced and modified UV vision among damselfishes, in particular the duplication of SWS1 among Pomacentrinae, might be seen as a key evolutionary innovation that facilitated the evolution of the exuberant variety of UV-reflectance traits and the diversification of this coral reef fish lineage.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Host phylogeny and environment shape the diversity of salamander skin bacterial communities

The composition and diversity of animal-associated microbial communities are shaped by multiple ecological and evolutionary processes acting at different spatial and temporal scales. Skin microbiomes are thought to be strongly influenced by the environment due to the direct interaction of the host's skin with the external media. As expected, the diversity of amphibian skin microbiomes is shaped by climate and host sampling habitats, whereas phylogenetic effects appear to be weak. However, the relative strength of phylogenetic and environmental effects on salamander skin microbiomes remains poorly understood. Here, we analysed sequence data from 1164 adult salamanders of 44 species to characterise and compare the diversity and composition of skin bacteria. We assessed the relative contribution of climate, host sampling habitat, and host phylogeny to the observed patterns of bacterial diversity. We found that bacterial alpha diversity was mainly associated with host sampling habitat and climate, but that bacterial beta diversity was more strongly associated with host taxonomy and phylogeny. This phylogenetic effect predominantly occurred at intermediate levels of host divergence (0-50 Mya). Our results support the importance of environmental factors shaping the diversity of salamander skin microbiota, but also support host phylogenetic history as a major factor shaping these bacterial communities.

Fossil berries reveal global radiation of the nightshade family by the early Cenozoic

Fossil discoveries can transform our understanding of plant diversification over time and space. Recently described fossils in many plant families have pushed their known records farther back in time, pointing to alternative scenarios for their origin and spread. Here, we describe two new Eocene fossil berries of the nightshade family (Solanaceae) from the Esmeraldas Formation in Colombia and the Green River Formation in Colorado (USA). The placement of the fossils was assessed using clustering and parsimony analyses based on 10 discrete and five continuous characters, which were also scored in 291 extant taxa. The Colombian fossil grouped with members of the tomatillo subtribe, and the Coloradan fossil aligned with the chili pepper tribe. Along with two previously reported early Eocene fossils from the tomatillo genus, these findings indicate that Solanaceae were distributed at least from southern South America to northwestern North America by the early Eocene. Together with two other recently discovered Eocene berries, these fossils demonstrate that the diverse berry clade and, in turn, the entire nightshade family, is much older and was much more widespread in the past than previously thought.

Recent progress on methods for estimating and updating large phylogenies

With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

OpenTree: A Python Package for Accessing and Analyzing Data from the Open Tree of Life

The Open Tree of Life project constructs a comprehensive, dynamic, and digitally available tree of life by synthesizing published phylogenetic trees along with taxonomic data. Open Tree of Life provides web-service application programming interfaces (APIs) to make the tree estimate, unified taxonomy, and input phylogenetic data available to anyone. Here, we describe the Python package opentree, which provides a user friendly Python wrapper for these APIs and a set of scripts and tutorials for straightforward downstream data analyses. We demonstrate the utility of these tools by generating an estimate of the phylogenetic relationships of all bird families, and by capturing a phylogenetic estimate for all taxa observed at the University of California Merced Vernal Pools and Grassland Reserve.[Evolution; open science; phylogenetics; Python; taxonomy.].

Combining molecular and geographical data to infer the phylogeny of Lamiales and its dispersal patterns in and out of the tropics

Lamiales is one of the most intractable orders of flowering plants, with several changes in family composition, and circumscription throughout history. The order is worldwide distributed, occurring in tropical forests and frozen habitats. In this study, a comprehensive phylogeny of Lamiales was reconstructed using DNA sequences. The tree was used to infer dispersal patterns, focusing on the tropics and extratropics. Molecular and species geographic data available from public repositories were combined to address both objectives. A total of 6,910 species, and 842 genera of Lamiales were sampled using the Python tool PyPHLAWD. The tree was inferred using RAxML, and recovered a monophyletic Lamiales. All 26 families were recovered as monophyletic with high support. The families Bignoniaceae, and Plantaginaceae are remarkable examples. The first emerged as monophyletic and included tribe Jacarandeae, while the later emerged as monophyletic in its sensu lato and included both the tribes Angelonieae, and Gratioleae. Distribution points for all species were retrieved from GBIF. After filtering, 1,136,425 records were retained. Species were coded as present in extratropical or tropical environments. The in and out of the tropics dispersal patterns were inferred using a maximum likelihood approach that identifies hidden rate changes. The model recovered higher rates of transition from extratropics to tropics, estimating two rates of state transitions. When ancestral states are considered, more discrete transitions from extratropics to tropics were observed. The extratropical state was also inferred for the crown node of Lamiales and old nested nodes, revealing a rare pattern of transitions to the tropics throughout the upper Cretaceous and Tertiary. A significant phylogenetic signal was recovered for the in and out of the tropics dispersal patterns, showing that state transitions are not frequent enough to erase the effect of tree structure on the data.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

The growing and vital role of botanical gardens in climate change research

Botanical gardens make unique contributions to climate change research, conservation, and public engagement. They host unique resources, including diverse collections of plant species growing in natural conditions, historical records, and expert staff, and attract large numbers of visitors and volunteers. Networks of botanical gardens spanning biomes and continents can expand the value of these resources. Over the past decade, research at botanical gardens has advanced our understanding of climate change impacts on plant phenology, physiology, anatomy, and conservation. For example, researchers have utilized botanical garden networks to assess anatomical and functional traits associated with phenological responses to climate change. New methods have enhanced the pace and impact of this research, including phylogenetic and comparative methods, and online databases of herbarium specimens and photographs that allow studies to expand geographically, temporally, and taxonomically in scope. Botanical gardens have grown their community and citizen science programs, informing the public about climate change and monitoring plants more intensively than is possible with garden staff alone. Despite these advances, botanical gardens are still underutilized in climate change research. To address this, we review recent progress and describe promising future directions for research and public engagement at botanical gardens.

Using target sequence capture to improve the phylogenetic resolution of a rapid radiation in New Zealand Veronica

PREMISE

Recent, rapid radiations present a challenge for phylogenetic reconstruction. Fast successive speciation events typically lead to low sequence divergence and poorly resolved relationships with standard phylogenetic markers. Target sequence capture of many independent nuclear loci has the potential to improve phylogenetic resolution for rapid radiations.

METHODS

Here we applied target sequence capture with 353 protein-coding genes (Angiosperms353 bait kit) to Veronica sect. Hebe (common name hebe) to determine its utility for improving the phylogenetic resolution of rapid radiations. Veronica section Hebe originated 5-10 million years ago in New Zealand, forming a monophyletic radiation of ca 130 extant species.

RESULTS

We obtained approximately 150 kbp of 353 protein-coding exons and an additional 200 kbp of flanking noncoding sequences for each of 77 hebe and two outgroup species. When comparing coding, noncoding, and combined data sets, we found that the latter provided the best overall phylogenetic resolution. While some deep nodes in the radiation remained unresolved, our phylogeny provided broad and often improved support for subclades identified by both morphology and standard markers in previous studies. Gene-tree discordance was nonetheless widespread, indicating that additional methods are needed to disentangle fully the history of the radiation.

CONCLUSIONS

Phylogenomic target capture data sets both increase phylogenetic signal and deliver new insights into the complex evolutionary history of rapid radiations as compared with traditional markers. Improving methods to resolve remaining discordance among loci from target sequence capture is now important to facilitate the further study of rapid radiations.

Physcraper: a Python package for continually updated phylogenetic trees using the Open Tree of Life

BACKGROUND

Phylogenies are a key part of research in many areas of biology. Tools that automate some parts of the process of phylogenetic reconstruction, mainly molecular character matrix assembly, have been developed for the advantage of both specialists in the field of phylogenetics and non-specialists. However, interpretation of results, comparison with previously available phylogenetic hypotheses, and selection of one phylogeny for downstream analyses and discussion still impose difficulties to one that is not a specialist either on phylogenetic methods or on a particular group of study.

RESULTS

Physcraper is a command-line Python program that automates the update of published phylogenies by adding public DNA sequences to underlying alignments of previously published phylogenies. It also provides a framework for straightforward comparison of published phylogenies with their updated versions, by leveraging upon tools from the Open Tree of Life project to link taxonomic information across databases. The program can be used by the nonspecialist, as a tool to generate phylogenetic hypotheses based on publicly available expert phylogenetic knowledge. Phylogeneticists and taxonomic group specialists will find it useful as a tool to facilitate molecular dataset gathering and comparison of alternative phylogenetic hypotheses (topologies).

CONCLUSION

The Physcraper workflow showcases the benefits of doing open science for phylogenetics, encouraging researchers to strive for better scientific sharing practices. Physcraper can be used with any OS and is released under an open-source license. Detailed instructions for installation and usage are available at https://physcraper.readthedocs.io.

Repeated evolution of a reproductive polyphenism in plants is strongly associated with bilateral flower symmetry

Polyphenisms are a special type of phenotypic plasticity in which the products of development are not continuous but instead are separate and distinct phenotypes produced in the same genetic background. One of the most widespread polyphenisms in the flowering plants is cleistogamy, in which the same individual plant produces both open, cross-pollinated flowers as well as highly reduced and closed, self-pollinated (cleistogamous) flowers.^1-5 Cleistogamy is not a rare evolutionary phenomenon. It has evolved independently at least 41 times.¹ But what favors the evolution of cleistogamy is still largely unknown.¹ Darwin⁶ proposed a hypothesis that has never been properly tested. He observed that cleistogamy is more common in taxa with bilaterally symmetric (zygomorphic) flowers than in those with radially symmetric (actinomorphic) flowers. Moreover, Darwin suggested that cleistogamous flowers help to ensure pollination, which he postulated is less certain in zygomorphic taxa that rely on more specialized groups of pollinators. Here, we combined the largest datasets on floral symmetry and cleistogamy and used phylogenetic approaches to show that cleistogamy is indeed disproportionately associated with zygomorphic flowers and that zygomorphic species are more likely to evolve cleistogamy than actinomorphic species. We also show that zygomorphic species are less capable of autonomous open-flower self-pollination (lower autofertility), suggesting that selection of cleistogamy via reproductive assurance in zygomorphic taxa could help account for Darwin's observation. Our results provide support for the hypothesis that polyphenisms are favored when organisms encounter contrasting environments.

Evolution of Disease Severity and Susceptibility in the Asteraceae to the Powdery Mildew Golovinomyces latisporus: Major Phylogenetic Structure Coupled With Highly Variable Disease Severity at Fine Scales

Pathogen host range and pathogen severity are dependent on interactions with their hosts and are hypothesized to have evolved as products of a coevolutionary arms race. An understanding of the factors that affect host range and pathogen severity is especially crucial in introduced pathogens that infect evolutionarily naïve hosts and cause substantial damage to ecosystems. Powdery mildews are detrimental pathogens found worldwide in managed and natural systems. Golovinomyces latisporus is a powdery mildew species that is especially damaging to plants within Asteraceae and to plants within the genus Helianthus in particular. In this study, we evaluated 126 species within Asteraceae to measure the role of host plant morphophysiological traits and evolutionary history on susceptibility to G. latisporus and disease severity. We observed phylogenetic signal in both susceptibility and severity within and among major clades of the Asteraceae. In general, there was a major phylogenetic structure of host severity to G. latisporus; however, there was some fine-scale phylogenetic variability. Phylogenetic statistical methods showed that chlorophyll content, biomass, stomatal index, and trichome density were not associated with disease severity, thus providing evidence that phylogenetic structure, rather than observed plant morphophysiological traits, is the most reliable predictor of pathogen severity. This work sheds light on the role that evolutionary history plays in plant susceptibility and severity to disease and underscores the relative unimportance of commonly assessed host plant traits in powdery mildew severity.

Do Alignment and Trimming Methods Matter for Phylogenomic (UCE) Analyses?

Alignment is a crucial issue in molecular phylogenetics because different alignment methods can potentially yield very different topologies for individual genes. But it is unclear if the choice of alignment methods remains important in phylogenomic analyses, which incorporate data from hundreds or thousands of genes. For example, problematic biases in alignment might be multiplied across many loci, whereas alignment errors in individual genes might become irrelevant. The issue of alignment trimming (i.e., removing poorly aligned regions or missing data from individual genes) is also poorly explored. Here, we test the impact of 12 different combinations of alignment and trimming methods on phylogenomic analyses. We compare these methods using published phylogenomic data from ultraconserved elements (UCEs) from squamate reptiles (lizards and snakes), birds, and tetrapods. We compare the properties of alignments generated by different alignment and trimming methods (e.g., length, informative sites, missing data). We also test whether these data sets can recover well-established clades when analyzed with concatenated (RAxML) and species-tree methods (ASTRAL-III), using the full data ($\sim $5000 loci) and subsampled data sets (10% and 1% of loci). We show that different alignment and trimming methods can significantly impact various aspects of phylogenomic data sets (e.g., length, informative sites). However, these different methods generally had little impact on the recovery and support values for well-established clades, even across very different numbers of loci. Nevertheless, our results suggest several "best practices" for alignment and trimming. Intriguingly, the choice of phylogenetic methods impacted the phylogenetic results most strongly, with concatenated analyses recovering significantly more well-established clades (with stronger support) than the species-tree analyses. [Alignment; concatenated analysis; phylogenomics; sequence length heterogeneity; species-tree analysis; trimming].

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

Building a Robust, Densely-Sampled Spider Tree of Life for Ecosystem Research

Phylogenetic relatedness is a key diversity measure for the analysis and understanding of how species and communities evolve across time and space. Understanding the nonrandom loss of species with respect to phylogeny is also essential for better-informed conservation decisions. However, several factors are known to influence phylogenetic reconstruction and, ultimately, phylogenetic diversity metrics. In this study, we empirically tested how some of these factors (topological constraint, taxon sampling, genetic markers and calibration) affect phylogenetic resolution and uncertainty. We built a densely sampled, species-level phylogenetic tree for spiders, combining Sanger sequencing of species from local communities of two biogeographical regions (Iberian Peninsula and Macaronesia) with a taxon-rich backbone matrix of Genbank sequences and a topological constraint derived from recent phylogenomic studies. The resulting tree constitutes the most complete spider phylogeny to date, both in terms of terminals and background information, and may serve as a standard reference for the analysis of phylogenetic diversity patterns at the community level. We then used this tree to investigate how partial data affect phylogenetic reconstruction, phylogenetic diversity estimates and their rankings, and, ultimately, the ecological processes inferred for each community. We found that the incorporation of a single slowly evolving marker (28S) to the DNA barcode sequences from local communities, had the highest impact on tree topology, closely followed by the use of a backbone matrix. The increase in missing data resulting from combining partial sequences from local communities only had a moderate impact on the resulting trees, similar to the difference observed when using topological constraints. Our study further revealed substantial differences in both the phylogenetic structure and diversity rankings of the analyzed communities estimated from the different phylogenetic treatments, especially when using non-ultrametric trees (phylograms) instead of time-stamped trees (chronograms). Finally, we provide some recommendations on reconstructing phylogenetic trees to infer phylogenetic diversity within ecological studies.