phylotaR: An Automated Pipeline for Retrieving Orthologous DNA Sequences from GenBank in R

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.

Evolutionary signatures of a trade-off in direct and indirect defenses across the wild grape genus, Vitis

Evolutionary correlations between chemical defense and protection by mutualist bodyguards have been long predicted, but tests of these patterns remain rare. We use a phylogenetic framework to test for evolutionary correlations indicative of trade-offs or synergisms between direct defense in the form of plant secondary metabolism and indirect defense in the form of leaf domatia, across 33 species in the wild grape genus, Vitis. We also performed a bioassay with a generalist herbivore to associate our chemical phenotypes with herbivore palatability. Finally, we tested whether defensive traits correlated with the average abiotic characteristics of each species' contemporary range and whether these correlations were consistent with plant defense theory. We found a negative evolutionary correlation between domatia size and the diversity of secondary metabolites in Vitis leaf tissue across the genus, and also that leaves with a higher diversity and richness of secondary metabolites were less palatable to a generalist herbivore, consistent with a trade-off in chemical and mutualistic defense investment. Predictions from plant defense theory were not supported by associations between investment in defense phenotypes and abiotic variables. Our work demonstrates an evolutionary pattern indicative of a trade-off between indirect and direct defense strategies across the Vitis genus.

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.

Evolutionary modularity, integration and disparity in an accretionary skeleton: analysis of venerid Bivalvia

Modular evolution, the relatively independent evolution of body parts, may promote high morphological disparity in a clade. Conversely, integrated evolution via stronger covariation of parts may limit disparity. However, integration can also promote high disparity by channelling morphological evolution along lines of least resistance-a process that may be particularly important in the accumulation of disparity in the many invertebrate systems having accretionary growth. We use a time-calibrated phylogenetic hypothesis and high-density, three-dimensional semilandmarking to analyse the relationship between modularity, integration and disparity in the most diverse extant bivalve family: the Veneridae. In general, venerids have a simple, two-module parcellation of their body that is divided into features of the calcium carbonate shell and features of the internal soft anatomy. This division falls more along developmental than functional lines when placed in the context of bivalve anatomy and biomechanics. The venerid body is tightly integrated in absolute terms, but disparity appears to increase with modularity strength among subclades and ecologies. Thus, shifts towards more mosaic evolution beget higher morphological variance in this speciose family.

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.

DNA barcoding reveals different cestode helminth species in northern European marine and freshwater ringed seals

Three subspecies of the ringed seal (Pusa hispida) are found in northeastern Europe: P. h. botnica in the Baltic Sea, P. h saimensis in Lake Saimaa in Finland, and P. h. ladogensis in Lake Ladoga in Russia. We investigated the poorly-known cestode helminth communities of these closely related but ecologically divergent subspecies using COI barcode data. Our results show that, while cestodes from the Baltic Sea represent Schistocephalus solidus, all worms from the two lakes are identified as Ligula intestinalis, a species that has previously not been reported from seals. The observed shift in cestode communities appears to be driven by differential availability of intermediate fish host species in marine vs. freshwater environments. Both observed cestode species normally infect fish-eating birds, so further work is required to elucidate the health and conservation implications of cestode infections in European ringed seals, whether L. intestinalis occurs also in marine ringed seals, and whether the species is able to reproduce in seal hosts. In addition, a deep barcode divergence found within S. solidus suggests the presence of cryptic diversity under this species name.

•

COI barcoding reveals different cestodes in marine and freshwater ringed seals.

•

Ligula intestinalis is reported for the first time from seals.

•

A deep barcode divergence is found within Schistocephalus solidus in the Baltic Sea.

Genome sequencing reveals extraordinary cephalic horns in the Madagascan dung beetle genus Helictopleurus (Coleoptera, Scarabaeinae): insight from a revision of fungicola species group

In this study, we test and corroborate the phylogenetic position of Heterosyphus within Helictopleurus using mitogenomes and nuclear loci. Our recent samplings revealed that males of the former Heterosyphus sicardi Paulian, 1975 (today under Helictopleurus d'Orbigny, 1915) have extraordinary bilateral clypeal horns which are exclusive within the genus. We provide a taxonomic review of the fungicola species group of Helictopleurus and discuss the systematic position of H. sicardi within the group. The male phenotype of H. sicardi is described and photographs of the body and genitalia of the members of the fungicola group are given, as well as a diagnostic key to species of the group. Helictopleurus fungicola peyrierasi is considered to be a distinct species within the genus (H. peyrierasi stat. rest.). Helictopleurus pluristriatus d'Orbigny, 1915 syn. nov. is established as a junior synonym of H. fungicola (Fairmaire, 1899).

Integrative taxonomic revision of the woodlouse-hunter spider genus Dysdera (Araneae: Dysderidae) in the Madeira archipelago with notes on its conservation status

Dysdera is a highly speciose genus of mid-sized, nocturnal hunting spiders, mostly circumscribed to the Mediterranean. The genus managed to colonize all Macaronesian archipelagos, and underwent major diversification in the Canary Islands. Here, we report on an independent diversification event on the Madeira archipelago. Based on the integration of morphological and molecular evidence, we describe 8 new species to science, Dysdera dissimilis sp. nov., Dysdera exigua sp. nov., Dysdera isambertoi sp. nov., Dysdera precaria sp. nov., Dysdera recondita sp. nov., Dysdera sandrae sp. nov., Dysdera teixeirai sp. nov., Dysdera titanica sp. nov. and redescribe Dysdera coiffaiti, Dysdera diversa and Dysdera portisancti. We synonymize Dysdera longibulbis and Dysdera vandeli under D. coiffaiti and D. diversa, respectively. Additionally, we use a multilocus target gene phylogeny to support a single colonization event of the archipelago followed by in situ diversification. We further discuss the discovered diversity patterns and their drivers. We conclude to that many of the species inhabit disturbed or fragile habitats and should be considered of high conservation concern.

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].

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.

Phylogenetic revision of Gymnotidae (Teleostei: Gymnotiformes), with descriptions of six subgenera

The diversity of gymnotid electric fishes has been intensely studied over the past 25 years, with 35 species named since 1994, compared to 11 species in the previous 236 years. Substantial effort has also been applied in recent years to documenting gymnotid interrelationships, with seven systematic studies published using morphological and molecular datasets. Nevertheless, until now, all gymnotids have been assigned to one of just two supraspecific taxa, the subfamily Electrophorinae with one genus Electrophorus and three valid species and the subfamily Gymnotine also with one genus Gymnotus and 43 valid species. This simple classification has obscured the substantial phenotypic and lineage diversity within the subfamily Gymnotine and hampered ecological and evolutionary studies of gymnotid biology. Here we present the most well-resolved and taxon-complete phylogeny of the Gymnotidae to date, including materials from all but one of the valid species. This phylogeny was constructed using a five-gene molecular dataset and a 115-character morphological dataset, enabling the inclusion of several species for which molecular data are still lacking. This phylogeny was time-calibrated using biogeographical priors in the absence of a fossil record. The tree topology is similar to those of previous studies, recovering all the major clades previously recognized with informal names. We propose a new gymnotid classification including two subfamilies (Electrophorinae and Gymnotinae) and six subgenera within the genus Gymnotus. Each subgenus exhibits a distinctive biogeographic distribution, within which most species have allopatric distributions and the subgenera are diverged from one another by an estimated 5–35 million years. We further provide robust taxonomic diagnoses, descriptions and identification keys to all gymnotid subgenera and all but four species. This new taxonomy more equitably partitions species diversity among supra-specific taxa, employing the previously vacant subgenus and subfamily ranks. This new taxonomy renders known gymnotid diversity more accessible to study by highlighting the deep divergences (chronological, geographical, genetic and morphological) among its several clades.

Bony Patchwork: Mosaic Patterns of Evolution in the Skull of Electric Fishes (Apteronotidae: Gymnotiformes)

Mosaic evolution refers to the pattern whereby different organismal traits exhibit differential rates of evolution typically due to reduced levels of trait covariation through deep time (i.e., modularity). These differences in rates can be attributed to variation in responses to selective pressures between individual traits. Differential responses to selective pressures also have the potential to facilitate functional specialization, allowing certain traits to track environmental stimuli more closely than others. The teleost skull is a multifunctional structure comprising a complex network of bones and thus an excellent system for which to study mosaic evolution. Here we construct an ultrametric phylogeny for a clade of Neotropical electric fishes (Apteronotidae: Gymnotiformes) and use three-dimensional geometric morphometrics to investigate patterns of mosaic evolution in the skull and jaws. We find strong support for a developmental, three-module hypothesis that consists of the face, braincase, and mandible, and we find that the mandible has evolved four times faster than its neighboring modules. We hypothesize that the functional specialization of the mandible in this group of fishes has allowed it to outpace the face and braincase and evolve in a more decoupled manner. We also hypothesize that this pattern of mosaicism may be widespread across other clades of teleost fishes.