The impact of taxon sampling on phylogenetic inference: a review of two decades of controversy

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Impact of environmental factors on spider silk properties

Spider orb webs have evolved to stop flying prey, fast and slow alike. One of the main web elements dissipating impact energy is the radial fibers, or major ampullate silks, which possess a toughness surpassing most man-made materials. Orb webs are extended phenotypes, and as such their architectural elements, including major ampullate silks, have been selected to optimize prey capture under the respective environmental conditions. In this study, we investigated the correlation of three landscape scales and three microhabitat characteristics with intrinsic silk properties (elastic modulus, yield stress, tensile strength, extensibility, and toughness) to understand underlying ecological patterns. For this purpose, we collected and mechanically tested major ampullate silks from 50 spider species inhabiting large altitudinal and climatic gradients in Colombia. Using regression analysis and model selection, we investigated the environmental drivers of inter- and intra-specific patterns of major ampullate silk properties, taking into account phylogenetic relatedness based on newly sequenced mitochondrial genomes. We found that the total amount of energy absorbed, i.e., toughness and tensile strength, is higher for fibers from species inhabiting regions where heavy rainfall is common. Interestingly, we observe the same general trend between individuals of the same species, stressing the importance of this environmental driver. We also observe a phylogenetic conservation in the relation of environmental variables with silk tensile strength and yield stress. In conclusion, the increase in major ampullate silk tensile strength and toughness may reflect an adaptation to prevent frequent rain damage to orb webs and the associated energetic loss.

A supermatrix phylogeny of the world's bees (Hymenoptera: Anthophila)

The increasing availability of large molecular phylogenies has provided new opportunities to study the evolution of species traits, their origins and diversification, and biogeography; yet there are limited attempts to synthesise existing phylogenetic information for major insect groups. Bees (Hymenoptera: Anthophila) are a large group of insect pollinators that have a worldwide distribution, and a wide variation in ecology, morphology, and life-history traits, including sociality. For these reasons, as well as their major economic importance as pollinators, numerous molecular phylogenetic studies of family and genus-level relationships have been published, providing an opportunity to assemble a bee 'tree-of-life'. We used publicly available genetic sequence data, including phylogenomic data, reconciled to a taxonomic database, to produce a concatenated supermatrix phylogeny for the Anthophila comprising 4,586 bee species, representing 23% of species and 82% of genera. At family, subfamily, and tribe levels, support for expected relationships was robust, but between and within some genera relationships remain uncertain. Within families, sampling of genera ranged from 67 to 100% but species coverage was lower (17-41%). Our phylogeny mostly reproduces the relationships found in recent phylogenomic studies with a few exceptions. We provide a summary of these differences and the current state of molecular data available and its gaps. We discuss the advantages and limitations of this bee supermatrix phylogeny (available online at beetreeoflife.org), which may enable new insights into long standing questions about evolutionary drivers in bees, and potentially insects more generally.

Systemic Caryospora-like coccidiosis in a clutch of hatchling red-eared slider turtles (Trachemys scripta elegans)

Caryospora-like organisms (CLOs) form a clade of at least 11 genotypes of related coccidia that can cause epizootic mortality in marine turtles. The biology, transmission, host species range, and host cell tropism of these organisms are still largely unknown. The goal of this study was to characterize the host cell tropism, pathologic and ultrastructural features, and phylogeny associated with the first report of a mortality event due to CLO in the freshwater red-eared slider turtle (Trachemys scripta elegans). Sudden mortalities within a clutch of captive-raised red-eared slider hatchlings (n = 8) were recorded, and deceased animals had severe segmental to diffuse, transmural, fibrinonecrotic enterocolitis and multifocal to coalescing hepatic necrosis, among other lesions associated with numerous intracytoplasmic developing stages of intralesional coccidia. Among the different developmental stages, merozoites were ultrastructurally characterized by an apical complex. A pan-apicomplexan polymerase chain reaction (PCR) yielded a 347 bp-amplicon matching the Schellackia/Caryospora-like clade with 99.1% identity to the US3 strain from green sea turtles (Chelonia mydas) and 99.1% identity to Schellackia sp. Isolate OC116. Surviving hatchlings were treated with toltrazuril sulfone (ponazuril) but were subsequently euthanized due to the risk of spreading the parasite to other chelonids in the collection. The ponazuril-treated hatchlings (n = 4) had mild proliferative anterior enteritis, with few intraepithelial coccidia in one hatchling confirmed as CLO by PCR. This is the first report of Caryospora-like coccidiosis in non-cheloniid turtles, highlighting the relevance of this disease as an emerging highly pathogenic intestinal and extra-intestinal form of coccidiosis of turtles with potential cross-species infectivity.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

Phylogeny and diversification of planthoppers (Hemiptera Fulgoromorpha) based on a comprehensive molecular dataset and large taxon sampling

Our understanding of the evolution of Fulgoromorpha (Insects, Hemiptera) has relied on molecular studies that have only considered either a limited number of taxa where all the families were not represented simultaneously, or a reduced number of genes.The absence of a global analysis comparing all the available data has thus led to significant biases in the analyzes, as evidenced by the incongruence of the results reported for planthopper phylogeny. Here we provide a phylogenetic and dating analysis of the Fulgoromorpha with a large sampling of 531 ingroup taxa, representing about 80% of the currently described suprageneric taxonomic diversity in this group. This study is based on most of the molecular sequences available to date and duly verified, for a set of nuclear and mitochondrial genes from a taxonomic sampling as complete as possible. The most significant results of our study are: (1) the unexpected paraphyly of Delphacidae whose Protodelphacida seem more related to Cixiidae than to other Delphacidae;(2) the group Meenoplidae-Kinnaridae recovered sister to the remaining Fulgoroidea families; (3) the early branching node of Tettigometridae sister of all the other families;(4) the Achilidae-Derbidae clade with Achilidae Plectoderini including Achilixiidae recovered as monophyletic as well as theFulgoridae-Dictyopharidae clade; and (5) the Tropiduchidae placed sister to the other so called 'higher' families (sec. Shcherbakov, 2006).Our divergence times analysis, calibrated with a set of duly verified fossils, suggests that the first diversification of planthoppers occurred in the Early Triassic around 240 Mya and those of the superfamilies Delphacoidea and Fulgoroidea in the Middle-Late Triassic around 210 Mya and 230 Mya, respectively. By the end of the Jurassic, all major planthopper lineages were originated, and all families, around 125 Mya, might havebeen driven in their distribution and evolution (in their first subfamilial divisions) by the geographical constraints of the Gondwanan break-up.Rapid evolutionary radiations occurred particularly in Fulgoridae around 125-130 Mya. Our results stress the importance of the good quality of the sequences used in the molecular analyzes and the primordial importance of a large sampling when analyzing the phylogeny of the group.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

Bizarre Morphology Obscures Real Affiliation: An Integrative Study of Enigmatic Cephalaspid Philine denticulata from Arctic Waters Reveals Its Unique Phylogenetic Position

The biodiversity of Cephalaspidea (Gastropoda: Heterobranchia) is poorly studied, and novel findings often lead to revisions at different taxonomic levels. The family Philinidae has a distinct set of defining characters in the shell, copulatory apparatus, and gizzard morphology, but several species, considered part of the family, deviate from it significantly. Philine denticulata (J. Adams, 1800) was considered to be a Philinidae despite the species morphology not fitting well with the family diagnosis. This species has an oval cylindrical external shell, jaws, and a gizzard lined with a chitinous layer with three thickened ridges. We studied Philine denticulata morphology on samples from the White and Barents Seas using the light and scanning electron microscopy methods in addition to computer microtomography. We also reconstructed its phylogenetic position using COI, 16S, 28S, and H3 genetic markers. Our integrative analysis revealed close relationships of this species to the family Aglajidae. Thus, we describe a new genus Philinissima gen. nov., including a single species Philinissima denticulata (J. Adams, 1800) comb. nov. which is the first aglajid registered from the Arctic waters. Our findings highlight possible problems in the Aglajidae taxonomical composition and an overall need for a dedicated integrative revision of the Cephalaspidea.

Plastid Phylogenomic Insights into the Inter-Tribal Relationships of Plantaginaceae

Plantaginaceae, consisting of 12 tribes, is a diverse, cosmopolitan family. To date, the inter-tribal relationships of this family have been unresolved, and the plastome structure and composition within Plantaginaceae have seldom been comprehensively investigated. In this study, we compared the plastomes from 41 Plantaginaceae species (including 6 newly sequenced samples and 35 publicly representative species) representing 11 tribes. To clarify the inter-tribal relationships of Plantaginaceae, we inferred phylogenic relationships based on the concatenated and coalescent analyses of 68 plastid protein-coding genes. PhyParts analysis was performed to assess the level of concordance and conflict among gene trees across the species tree. The results indicate that most plastomes of Plantaginaceae are largely conserved in terms of genome structure and gene content. In contrast to most previous studies, a robust phylogeny was recovered using plastome data, providing new insights for better understanding the inter-tribal relationships of Plantaginaceae. Both concatenated and coalescent phylogenies favored the sister relationship between Plantagineae and Digitalideae, as well as between Veroniceae and Hemiphragmeae. Sibthorpieae diverged into a separate branch which was sister to a clade comprising the four tribes mentioned above. Furthermore, the sister relationship between Russelieae and Cheloneae is strongly supported. The results of PhyParts showed gene tree congruence and conflict to varying degrees, but most plastid genes were uninformative for phylogenetic nodes, revealing the defects of previous studies using single or multiple plastid DNA sequences to infer the phylogeny of Plantaginaceae.

Not the Last Piece of the Puzzle: Niphargus Phylogeny in Hungary

The Palaearctic genus Niphargus is a promising model system to understand subterranean fauna genesis in Europe. The Pannonian Plain (mainly covered by Hungary) in Central Europe, once being the area of the Paratethys, is a key area for Niphargus diversification. However, our knowledge on Hungarian species of Niphargus is primarily based on sporadic taxonomical works from the pre-molecular era. Here, we studied 14 localities, covering the eight valid Hungarian species of Niphargus and including nine previously unstudied populations. Based on sequences of three gene fragments, we reconstructed their phylogeny using maximum likelihood and Bayesian approaches. We found that not all Hungarian species of Niphargus are closely related, and even species sampled at the same localities can belong to different clades. Some Hungarian species form monophyletic clades, while others are nested in various non-Hungarian lineages. The new populations are all genetically distinct from the known species. Our results suggest that the Hungarian Niphargus fauna has originated from seven unrelated clades and its diversity is underestimated due to unknown populations and cryptic species. The detection of genetically distinct species of Niphargus from non-carbonate regions calls for further research efforts. The high diversity and the number of putative new species in the N. tatrensis clade warrants further, high-resolution phylogenetic studies.

Best practice for wildlife gut microbiome research: A comprehensive review of methodology for 16S rRNA gene investigations

Extensive research in well-studied animal models underscores the importance of commensal gastrointestinal (gut) microbes to animal physiology. Gut microbes have been shown to impact dietary digestion, mediate infection, and even modify behavior and cognition. Given the large physiological and pathophysiological contribution microbes provide their host, it is reasonable to assume that the vertebrate gut microbiome may also impact the fitness, health and ecology of wildlife. In accordance with this expectation, an increasing number of investigations have considered the role of the gut microbiome in wildlife ecology, health, and conservation. To help promote the development of this nascent field, we need to dissolve the technical barriers prohibitive to performing wildlife microbiome research. The present review discusses the 16S rRNA gene microbiome research landscape, clarifying best practices in microbiome data generation and analysis, with particular emphasis on unique situations that arise during wildlife investigations. Special consideration is given to topics relevant for microbiome wildlife research from sample collection to molecular techniques for data generation, to data analysis strategies. Our hope is that this article not only calls for greater integration of microbiome analyses into wildlife ecology and health studies but provides researchers with the technical framework needed to successfully conduct such investigations.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

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

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

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

Phylogenetic resolution of the fly superfamily Ephydroidea-Molecular systematics of the enigmatic and diverse relatives of Drosophilidae

The schizophoran superfamily Ephydroidea (Diptera: Cyclorrhapha) includes eight families, ranging from the well-known vinegar flies (Drosophilidae) and shore flies (Ephydridae), to several small, relatively unusual groups, the phylogenetic placement of which has been particularly challenging for systematists. An extraordinary diversity in life histories, feeding habits and morphology are a hallmark of fly biology, and the Ephydroidea are no exception. Extreme specialization can lead to "orphaned" taxa with no clear evidence for their phylogenetic position. To resolve relationships among a diverse sample of Ephydroidea, including the highly modified flies in the families Braulidae and Mormotomyiidae, we conducted phylogenomic sampling. Using exon capture from Anchored Hybrid Enrichment and transcriptomics to obtain 320 orthologous nuclear genes sampled for 32 species of Ephydroidea and 11 outgroups, we evaluate a new phylogenetic hypothesis for representatives of the superfamily. These data strongly support monophyly of Ephydroidea with Ephydridae as an early branching radiation and the placement of Mormotomyiidae as a family-level lineage sister to all remaining families. We confirm placement of Cryptochetidae as sister taxon to a large clade containing both Drosophilidae and Braulidae-the latter a family of honeybee ectoparasites. Our results reaffirm that sampling of both taxa and characters is critical in hyperdiverse clades and that these factors have a major influence on phylogenomic reconstruction of the history of the schizophoran fly radiation.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Phylogenies in yeast species descriptions: in defense of neighbor-joining

The neighbor-joining (NJ) method of tree inference is examined, with special attention to its use in yeast species descriptions. How the often-vilified method works is often misunderstood. More importantly, given the right kind of data, its output is a phylogram that illustrates a hypothetical phylogeny that is just as credible as that obtained by any other method. And as with any other method, the result is greatly affected by sampling intensity, particularly the number of aligned positions used for analysis. I address various allegations, including the claim that the method is phenetic, and therefore not phylogenetic. I argue that NJ is the most suitable tree inference method to use in yeast species descriptions, primarily because it is best at visually preserving the extent of sequence divergence between close relatives, which continues to be the primary criterion for yeast species delineation. The relevance of bootstraps in application of the phylogenetic species concept is discussed. This article is protected by copyright. All rights reserved.

A tale of two bellies: systematics of the oval frogs (Anura: Microhylidae: Elachistocleis Parker, 1927)

Oval frogs (Elachistocleis) have a broad geographic distribution covering nearly all of South America and parts of Central America. They also have a large inter- and intraspecific variation of the few morphological characters commonly used as diagnostic traits among species of the genus. Based on molecular data, we provide the most complete phylogeny of Elachistocleis to date, and explore its genetic diversity using distance-based and tree-based methods for putative species delimitation. Our results show that at least two of the most relevant traditional characters used in the taxonomy of this group (belly pattern and dorsal median white line) carry less phylogenetic information than previously thought. Based on our results, we propose some synonymizations and some candidate new species. This study is a first major step in disentangling the current systematics of Elachistocleis. Yet, a comprehensive review of morphological data is needed before any new species descriptions can be properly made.

Phylogenomic Analysis Reconstructed the Order Matoniales from Paleopolyploidy Veil

Phylogenetic conflicts limit our understanding of the evolution of terrestrial life under multiple whole genome duplication events, and the phylogeny of early terrestrial plants remains full of controversy. Although much incongruence has been solved with so-called robust topology based on single or lower copy genes, the evolutionary mechanisms behind phylogenetic conflicts such as polyploidization remain poorly understood. Here, through decreasing the effects of polyploidization and increasing the samples of species, which represent all four orders and eight families that comprise early leptosporangiate ferns, we have reconstructed a robust phylogenetic tree and network with 1125 1-to-1 orthologs based on both coalescent and concatenation methods. Our data consistently suggest that Matoniales, as a monophyletic lineage including Matoniaceae and Dipteridaceae, should be redefined as an ordinal rank. Furthermore, we have identified and located at least 11 whole-genome duplication events within the evolutionary history of four leptosporangiates lineages, and associated polyploidization with higher speciation rates and mass extinction events. We hypothesize that paleopolyploidization may have enabled leptosporangiate ferns to survive during mass extinction events at the end Permian period and then flourish throughout the Mesozoic era, which is supported by extensive fossil records. Our results highlight how ancient polyploidy can result in rapid species radiation, thus causing phylogenetic conflicts yet allowing plants to survive and thrive during mass extinction events.

Mitochondrial Genomes Provide New Phylogenetic and Evolutionary Insights into Psilidae (Diptera: Brachycera)

Psilidae (Diptera: Brachycera) is a moderate-sized family currently placed in the superfamily Diopsoidea and contains some destructive agricultural and forestry pests. The systematic position and intrafamilial classification of rust flies are in need of further study, and the available molecular data of Psilidae are still limited. In this study, we present the mitochondrial genomes of 6 Psilidae species (Chamaepsilatestudinaria Wang and Yang, Chyliza bambusae Wang and Yang, Chy. chikuni Wang, Loxocera lunata Wang and Yang, L. planivena Wang and Yang and L. sinica Wang and Yang). Comparative analyses show a conserved genome structure, in terms of gene composition and arrangement, and a highly Adenine plus Thymine biased nucleotide composition of the 6 psilid mitogenomes. Mitochondrial evolutionary rates vary among the 6 species, with species of Chylizinae exhibiting a slower average rate than species of Psilinae. The length, the nucleotide composition, and the copy number of repeat units of the control region are variable among the 6 species, which may offer useful information for phylogenetic and evolutionary studies of Psilidae. Phylogenetic analyses based on 4 mitogenomic datasets (AA, PCG, PCG12RNA, and PCGRNA) support the monophyly of Psilidae, and the sister relationship between Chylizinae and Psilinae, while Diopsoidea is suggested to be non-monophyletic. Our study enlightens the future application of mitogenomic data in the phylogenetic and evolutionary studies of Psilidae, based on denser taxon sampling.

A Linear Time Solution to the Labeled Robinson-Foulds Distance Problem

MOTIVATION

A large variety of pairwise measures of similarity or dissimilarity have been developed for comparing phylogenetic trees, e.g. species trees or gene trees. Due to its intuitive definition in terms of tree clades and bipartitions and its computational efficiency, the Robinson-Foulds (RF) distance is the most widely used for trees with unweighted edges and labels restricted to leaves (representing the genetic elements being compared). However, in the case of gene trees, an important information revealing the nature of the homologous relation between gene pairs (orthologs, paralogs, xenologs) is the type of event associated to each internal node of the tree, typically speciations or duplications, but other types of events may also be considered, such as horizontal gene transfers. This labeling of internal nodes is usually inferred from a gene tree/species tree reconciliation method. Here, we address the problem of comparing such event-labeled trees. The problem differs from the classical problem of comparing uniformly labeled trees (all labels belonging to the same alphabet) that may be done using the Tree Edit Distance (TED) mainly due to the fact that, in our case, two different alphabets are considered for the leaves and internal nodes of the tree, and leaves are not affected by edit operations.

RESULTS

We propose an extension of the RF distance to event-labeled trees, based on edit operations comparable to those considered for TED: node insertion, node deletion and label substitution. We show that this new Labeled Robinson Foulds (LRF) distance can be computed in linear time, in addition of maintaining other desirable properties: being a metric, reducing to RF for trees with no labels on internal nodes and maintaining an intuitive interpretation. The algorithm for computing the LRF distance enables novel analyses on event-label trees such as reconciled gene trees. Here, we use it to study the impact of taxon sampling on labeled gene tree inference, and conclude that denser taxon sampling yields trees with better topology but worse labeling.

Dactylogyridae 2022: a meta-analysis of phylogenetic studies and generic diagnoses of parasitic flatworms using published genetic and morphological data

Dactylogyridae is one of the most studied families of parasitic flatworms with more than 1000 species and 166 genera described to date including ecto- and endoparasites. Dactylogyrid monogeneans were suggested as model organisms for host-parasite macroevolutionary and biogeographical studies due to the scientific and economic importance of some of their host lineages. Consequently, an array of phylogenetic research into different dactylogyrid lineages has been produced over the past years but the last family-wide study was published 16 years ago. Here, we provide a meta-analysis of the phylogenetic relationships of Dactylogyridae including representatives of all genera with available molecular data (n=67). First, we investigate the systematic informativeness of morphological characters widely used to diagnose dactylogyrid genera through a parsimony analysis of the characters, character mapping, and phylogenetic comparative methods. Second, we provide an overview of the current state of the systematics of the family and its subfamilies, and summarise potentially poly- and paraphyletic genera. Third, we elaborate on the implications of taxonomic, citation, and confirmation bias in past studies. Fourth, we discuss host range, biogeographical, and freshwater-marine patterns. We found two well-supported macroclades which we assigned to the subfamilies Dactylogyrinae and Ancyrocephalinae. These subfamilies further include 16 well-supported clades with only a few synapomorphies that could be deduced from generic diagnoses in the literature. Furthermore, few morphological characters considered systematically informative at the genus level display a strong phylogenetic signal. However, the parsimony analysis suggests that these characters provide little information on the relationships between genera. We conclude that a strong taxonomic bias and low coverage of DNA sequences and regions limit knowledge on morphological and biogeographical evolutionary patterns that can be inferred from these results. We propose addressing potential citation and confirmation biases through a 'level playing field' multiple sequence alignment as provided by this study.

Plastid phylogenomics of the Sansevieria Clade of Dracaena (Asparagaceae) resolves a recent radiation

Best known as low maintenance houseplants, sansevierias are a diverse group of flowering plants native to Africa, Madagascar, the Arabian Peninsula, and the Indian subcontinent. Traditionally recognised as a distinct genus, Sansevieria was recently merged with the larger genus Dracaena based on molecular phylogenetic data. Within the Sansevieria Clade of Dracaena, taxonomic uncertainties remain despite attempts to unravel the relationships between the species. To investigate the evolutionary relationships, morphological evolution and biogeographical history in the group, we aim to reconstruct a robust dated phylogenetic hypothesis. Using genome skimming, a chloroplast genome (cpDNA) dataset and a nuclear ribosomal (nrDNA) dataset were generated. The sampling included representatives of all sections and informal groups previously described in Sansevieria based on morphology. Analysis of the cpDNA dataset using a maximum likelihood approach resulted in a well-supported phylogeny. The time-calibrated phylogeny indicated a recent radiation with five main clades emerging in the Pliocene. Two strongly supported clades align with previously defined groups, i.e., Sansevieria section Dracomima, characterised by the Dracomima-type inflorescence, and the Zeylanica informal group, native to the Indian subcontinent. Other previously defined groups were shown to be polyphyletic; a result of convergent evolution of the identifying characters. Switches between flat and cylindrical leaves occurred multiple times in the evolution of the Sansevieria Clade. Similarly, the Cephalantha-type inflorescence has originated multiple times from an ancestor with a Sansevieria-type inflorescence. Analysis of the nrDNA dataset resulted in a phylogenetic hypothesis with low resolution, yet it supported the same two groups confirmed by the cpDNA dataset. This study furthers our understanding of the evolution of the Sansevieria Clade, which will benefit taxonomic and applied research, and aid conservation efforts.

Re-evaluating Deep Neural Networks for Phylogeny Estimation: The Issue of Taxon Sampling

Deep neural networks (DNNs) have been recently proposed for quartet tree phylogeny estimation. Here, we present a study evaluating recently trained DNNs in comparison to a collection of standard phylogeny estimation methods on a heterogeneous collection of datasets simulated under the same models that were used to train the DNNs, and also under similar conditions but with higher rates of evolution. Our study shows that using DNNs with quartet amalgamation is less accurate than several standard phylogeny estimation methods we explore (e.g., maximum likelihood and maximum parsimony). We further find that simple standard phylogeny estimation methods match or improve on DNNs for quartet accuracy, especially, but not exclusively, when used in a global manner (i.e., the tree on the full dataset is computed and then the induced quartet trees are extracted from the full tree). Thus, our study provides evidence that a major challenge impacting the utility of current DNNs for phylogeny estimation is their restriction to estimating quartet trees that must subsequently be combined into a tree on the full dataset. In contrast, global methods (i.e., those that estimate trees from the full set of sequences) are able to benefit from taxon sampling, and hence have higher accuracy on large datasets.

The Azurin Coding Gene: Origin and Phylogenetic Distribution

Azurin is a bacterial-derived cupredoxin, which is mainly involved in electron transport reactions. Interest in azurin protein has risen in recent years due to its anticancer activity and its possible applications in anticancer therapies. Nevertheless, the attention of the scientific community only focused on the azurin protein found in Pseudomonas aeruginosa (Proteobacteria, Gammaproteobacteria). In this work, we performed the first comprehensive screening of all the bacterial genomes available in online repositories to assess azurin distribution in the three domains of life. The Azurin coding gene was not detected in the domains Archaea and Eucarya, whereas it was detected in phyla other than Proteobacteria, such as Bacteroidetes, Verrucomicrobia and Chloroflexi, and a phylogenetic analysis of the retrieved sequences was performed. Observed patchy distribution and phylogenetic data suggest that once it appeared in the bacterial domain, the azurin coding gene was lost in several bacterial phyla and/or anciently horizontally transferred between different phyla, even though a vertical inheritance appeared to be the major force driving the transmission of this gene. Interestingly, a shared conserved domain has been found among azurin members of all the investigated phyla. This domain is already known in P. aeruginosa as p28 domain and its importance for azurin anticancer activity has been widely explored. These findings may open a new and intriguing perspective in deciphering the azurin anticancer mechanisms and to develop new tools for treating cancer diseases.

Power and Weakness of Repetition - Evaluating the Phylogenetic Signal From Repeatomes in the Family Rosaceae With Two Case Studies From Genera Prone to Polyploidy and Hybridization (Rosa and Fragaria)

Plant genomes consist, to a considerable extent, of non-coding repetitive DNA. Several studies showed that phylogenetic signals can be extracted from such repeatome data by using among-species dissimilarities from the RepeatExplorer2 pipeline as distance measures. Here, we advanced this approach by adjusting the read input for comparative clustering indirectly proportional to genome size and by summarizing all clusters into a main distance matrix subjected to Neighbor Joining algorithms and Principal Coordinate Analyses. Thus, our multivariate statistical method works as a "repeatomic fingerprint," and we proved its power and limitations by exemplarily applying it to the family Rosaceae at intrafamilial and, in the genera Fragaria and Rosa, at the intrageneric level. Since both taxa are prone to hybridization events, we wanted to show whether repeatome data are suitable to unravel the origin of natural and synthetic hybrids. In addition, we compared the results based on complete repeatomes with those from ribosomal DNA clusters only, because they represent one of the most widely used barcoding markers. Our results demonstrated that repeatome data contained a clear phylogenetic signal supporting the current subfamilial classification within Rosaceae. Accordingly, the well-accepted major evolutionary lineages within Fragaria were distinguished, and hybrids showed intermediate positions between parental species in data sets retrieved from both complete repeatomes and rDNA clusters. Within the taxonomically more complicated and particularly frequently hybridizing genus Rosa, we detected rather weak phylogenetic signals but surprisingly found a geographic pattern at a population scale. In sum, our method revealed promising results at larger taxonomic scales as well as within taxa with manageable levels of reticulation, but success remained rather taxon specific. Since repeatomes can be technically easy and comparably inexpensively retrieved even from samples of rather poor DNA quality, our phylogenomic method serves as a valuable alternative when high-quality genomes are unavailable, for example, in the case of old museum specimens.

Phylogenetic analysis of the family Megalopodidae (Coleoptera: Chrysomeloidea): better taxon-sampling facilitates detection of new relationships and new taxa

The systematics of Megalopodidae is not adequately known, in spite of it being a relatively small group of phytophagous beetles. The first phylogenetic analysis of Megalopodidae with a comprehensive generic representation (25 genera of 30 described, 10 subgenera and 77 species) is undertaken. A parsimony analysis under equal and implied weights was carried out based on 147 adult and larval morphological characters. Subfamilies Palophaginae and Zeugophorinae were recovered as monophyletic, by contrast with Megalopodinae, which proved to be paraphyletic. Atelederinae are proposed as a new subfamily. Also, three tribes and three subtribes within Megalopodinae are proposed: Leucasteini trib.n., Sphondyliini trib.n. and Megalopodini, the latter including Macrolophina subtrib.n., Temnaspidina subtrib.n. and Megalopodina. The genera Macrolopha, Kuilua, Poecilomorpha, Temnaspis, Antonaria, Agathomerus, Megalopus and Bothromegalopus were recovered as non-monophyletic. New delimitations of the polyphyletic genera Poecilomorpha and Macrolopha are proposed, Clythraxeloma is resurrected, and the subgenera of Agathomerus are suppressed. The following new combinations are proposed: Kuilua apicata (Fairmaire), K. nyassae (Jacoby), Poecilomorpha cribricollis (Pic), P. minuta (Pic), Clythraxeloma assamensis (Jacoby), C. bipartita (Lacordaeri), C. discolineata (Pic), C. downesii (Baly), C. gerstaeckeri (Westwood), C. laosensis (Pic), C. maculata (Pic), C. mouhoti (Baly), C. nigrocyanea (Motschulsky), C. pretiosa (Reineck), Temnaspis tricoloripes (Pic) and Barticaria faciatus (Dalman). Clythraxeloma cyanipennis Kraatz is a restored combination. Distribution patterns of Megalopodidae largely conform to the breakup of Gondwanaland, with its main clades having particular distributions: Andean-Australian (Palophaginae), Ethiopian (Leucasteini, Sphondyliini, and Macrolophina), Neotropical (Ateledrinae and Megalopodina) and Ethiopian-Oriental-Palaearctic (Temnaspidina the result of a secondary expansion. Zeugophorinae present a worldwide distribution, except for the Neotropical and Andean regions, which may be the result of geodispersal. The findings of the present study also shed light on groups with taxonomic issues, where phylogenetic analyses are strongly needed.

Global diversification and evolutionary history of onchidiid slugs (Gastropoda, Pulmonata)

Many marine species are specialized to specific parts of a habitat. In a mangrove forest, for instance, species may be restricted to the mud surface, the roots and trunks of mangrove trees, or rotting logs, which can be regarded as distinct microhabitats. Shifts to new microhabitats may be an important driver of sympatric speciation. However, the evolutionary history of these shifts is still poorly understood in most groups of marine organisms, because it requires a well-supported phylogeny with relatively complete taxon sampling. Onchidiid slugs are an ideal case study for the evolutionary history of habitat and microhabitat shifts because onchidiid species are specialized to different tidal zones and microhabitats in mangrove forests and rocky shores, and the taxonomy of the family in the Indo-West Pacific has been recently revised in a series of monographs. Here, DNA sequences for onchidiid species from the North and East Pacific, the Caribbean, and the Atlantic are used to reconstruct phylogenetic relationships among Onchidella species, and are combined with new data for Indo-West Pacific species to reconstruct a global phylogeny of the family. The phylogenetic relationships of onchidiid slugs are reconstructed based on three mitochondrial markers (COI, 12S, 16S) and three nuclear markers (28S, ITS2, H3) and nearly complete taxon sampling (all 13 genera and 62 of the 67 species). The highly-supported phylogeny presented here suggests that ancestral onchidiids most likely lived in the rocky intertidal, and that a lineage restricted to the tropical Indo-West Pacific colonized new habitats, including mudflats, mangrove forests, and high-elevation rainforests. Many onchidiid species in the Indo-West Pacific diverged during the Miocene, around the same time that a high diversity of mangrove plants appears in the fossil record, while divergence among Onchidella species occurred earlier, likely beginning in the Eocene. It is demonstrated that ecological specialization to microhabitats underlies the divergence between onchidiid genera, as well as the diversification through sympatric speciation in the genera Wallaconchis and Platevindex. The geographic distributions of onchidiid species also indicate that allopatric speciation played a key role in the diversification of several genera, especially Onchidella and Peronia. The evolutionary history of several morphological traits (penial gland, rectal gland, dorsal eyes, intestinal loops) is examined in relation to habitat and microhabitat evolutionary transitions and that the rectal gland of onchidiids is an adaptation to high intertidal and terrestrial habitats.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L

BACKGROUNDS

Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae.

RESULTS

In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists.

CONCLUSIONS

The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

Transmission dynamics of an antimicrobial resistant Campylobacter jejuni lineage in New Zealand's commercial poultry network

Understanding the relative contribution of different between-farm transmission pathways is essential in guiding recommendations for mitigating disease spread. This study investigated the association between contact pathways linking poultry farms in New Zealand and the genetic relatedness of antimicrobial resistant Campylobacter jejuni Sequence Type 6964 (ST-6964), with the aim of identifying the most likely contact pathways that contributed to its rapid spread across the industry. Whole-genome sequencing was performed on 167C. jejuni ST-6964 isolates sampled from across 30 New Zealand commercial poultry enterprises. The genetic relatedness between isolates was determined using whole genome multilocus sequence typing (wgMLST). Permutational multivariate analysis of variance and distance-based linear models were used to explore the strength of the relationship between pairwise genetic associations among the C. jejuni isolates and each of several pairwise distance matrices, indicating either the geographical distance between farms or the network distance of transportation vehicles. Overall, a significant association was found between the pairwise genetic relatedness of the C. jejuni isolates and the parent company, the road distance and the network distance of transporting feed vehicles. This result suggests that the transportation of feed within the commercial poultry industry as well as other local contacts between flocks, such as the movements of personnel, may have played a significant role in the spread of C. jejuni. However, further information on the historical contact patterns between farms is needed to fully characterise the risk of these pathways and to understand how they could be targeted to reduce the spread of C. jejuni.

Phylogenetic analysis and development of molecular markers for five medicinal Alpinia species based on complete plastome sequences

BACKGROUND

Alpinia species are widely used as medicinal herbs. To understand the taxonomic classification and plastome evolution of the medicinal Alpinia species and correctly identify medicinal products derived from Alpinia species, we systematically analyzed the plastome sequences from five Alpinia species. Four of the Alpinia species: Alpinia galanga (L.) Willd., Alpinia hainanensis K.Schum., Alpinia officinarum Hance, and Alpinia oxyphylla Miq., are listed in the Chinese pharmacopeia. The other one, Alpinia nigra (Gaertn.) Burtt, is well known for its medicinal values.

RESULTS

The four Alpinia species: A. galanga, A. nigra, A. officinarum, and A. oxyphylla, were sequenced using the Next-generation sequencing technology. The plastomes were assembled using Novoplasty and annotated using CPGAVAS2. The sizes of the four plastomes range from 160,590 bp for A. galanga to 164,294 bp for A. nigra, and display a conserved quadripartite structure. Each of the plastomes encodes a total of 111 unique genes, including 79 protein-coding, 28 tRNA, and four rRNA genes. In addition, 293-296 SSRs were detected in the four plastomes, of which the majority are mononucleotides Adenine/Thymine and are found in the noncoding regions. The long repeat analysis shows all types of repeats are contained in the plastomes, of which palindromic repeats occur most frequently. The comparative genomic analyses revealed that the pair of the inverted repeats were less divergent than the single-copy region. Analysis of sequence divergence on protein-coding genes showed that two genes (accD and ycf1) had undergone positive selection. Phylogenetic analysis based on coding sequence of 77 shared plastome genes resolves the molecular phylogeny of 20 species from Zingiberaceae. In particular, molecular phylogeny of four sequenced Alpinia species (A. galanga, A. nigra, A. officinarum, and A. oxyphylla) based on the plastome and nuclear sequences showed congruency. Furthermore, a comparison of the four newly sequenced Alpinia plastomes and one previously reported Alpinia plastomes (accession number: NC_048461) reveals 59 highly divergent intergenic spacer regions. We developed and validated two molecular markers Alpp and Alpr, based on two regions: petN-psbM and psaJ-rpl33, respectively. The discrimination success rate was 100 % in validation experiments.

CONCLUSIONS

The results from this study will be invaluable for ensuring the effective and safe uses of Alpinia medicinal products and for the exploration of novel Alpinia species to improve human health.

A Repertory of Rearrangements and the Loss of an Inverted Repeat Region in Passiflora Chloroplast Genomes

Chloroplast genomes (cpDNA) in angiosperms are usually highly conserved. Although rearrangements have been observed in some lineages, such as Passiflora, the mechanisms that lead to rearrangements are still poorly elucidated. In the present study, we obtained 20 new chloroplast genomes (18 species from the genus Passiflora, and Dilkea retusa and Mitostemma brevifilis from the family Passifloraceae) in order to investigate cpDNA evolutionary history in this group. Passiflora cpDNAs vary in size considerably, with ∼50 kb between shortest and longest. Large inverted repeat (IR) expansions were identified, and at the extreme opposite, the loss of an IR was detected for the first time in Passiflora, a rare event in angiosperms. The loss of an IR region was detected in Passiflora capsularis and Passiflora costaricensis, a species in which occasional biparental chloroplast inheritance has previously been reported. A repertory of rearrangements such as inversions and gene losses were detected, making Passiflora one of the few groups with complex chloroplast genome evolution. We also performed a phylogenomic study based on all the available cp genomes and our analysis implies that there is a need to reconsider the taxonomic classifications of some species in the group.

Remotely sensed between-individual functional trait variation in a temperate forest

Trait-based ecology holds the promise to explain how plant communities work, for example, how functional diversity may support community productivity. However, so far it has been difficult to combine field-based approaches assessing traits at the level of plant individuals with limited spatial coverage and approaches using remote sensing (RS) with complete spatial coverage but assessing traits at the level of vegetation pixels rather than individuals. By delineating all individual-tree crowns within a temperate forest site and then assigning RS-derived trait measures to these trees, we combine the two approaches, allowing us to use general linear models to estimate the influence of taxonomic or environmental variation on between- and within-species variation across contiguous space.We used airborne imaging spectroscopy and laser scanning to collect individual-tree RS data from a mixed conifer-angiosperm forest on a mountain slope extending over 5.5 ha and covering large environmental gradients in elevation as well as light and soil conditions. We derived three biochemical (leaf chlorophyll, carotenoids, and water content) and three architectural traits (plant area index, foliage-height diversity, and canopy height), which had previously been used to characterize plant function, from the RS data. We then quantified the contributions of taxonomic and environmental variation and their interaction to trait variation and partitioned the remaining within-species trait variation into smaller-scale spatial and residual variation. We also investigated the correlation between functional trait and phylogenetic distances at the between-species level. The forest consisted of 13 tree species of which eight occurred in sufficient abundance for quantitative analysis.On average, taxonomic variation between species accounted for more than 15% of trait variation in biochemical traits but only around 5% (still highly significant) in architectural traits. Biochemical trait distances among species also showed a stronger correlation with phylogenetic distances than did architectural trait distances. Light and soil conditions together with elevation explained slightly more variation than taxonomy across all traits, but in particular increased plant area index (light) and reduced canopy height (elevation). Except for foliage-height diversity, all traits were affected by significant interactions between taxonomic and environmental variation, the different responses of the eight species to the within-site environmental gradients potentially contributing to the coexistence of the eight abundant species.We conclude that with high-resolution RS data it is possible to delineate individual-tree crowns within a forest and thus assess functional traits derived from RS data at individual level. With this precondition fulfilled, it is then possible to apply tools commonly used in field-based trait ecology to partition trait variation among individuals into taxonomic and potentially even genetic variation, environmental variation, and interactions between the two. The method proposed here presents a promising way of assessing individual-based trait information with complete spatial coverage and thus allowing analysis of functional diversity at different scales. This information can help to better understand processes shaping community structure, productivity, and stability of forests.

Phylogenetic uncertainty and the inference of patterns in community ecology and comparative studies

Progress in phylogenetic community ecology is often limited by the availability of phylogenetic information and the lack of appropriate methods and solutions to deal with this problem. We estimate the effect of the lack of phylogenetic information on the relations among taxa measured by commonly used phylogenetic metrics in comparative studies and community ecology, namely: Blomberg's K phylogenetic signal, Faith's Phylogenetic Diversity (PD), Mean Phylogenetic Distance (MPD) and Mean Nearest Taxon Distance (MNTD). To overcome this problem, we tested two possible solutions: Polytomic trees and Operational trees. Our results show that the effects on K values strongly depended on the level of phylogenetic signal. In the case of the community metrics, the effects were insensitive to the patterns of species distribution in the communities. Community metrics tended to be overestimated with both Polytomic and Operational trees, but the overestimation was higher with Polytomic trees. PD and MPD metrics were less biased than MNTD metric. We show that the lack of phylogenetic resolution is not necessarily problematic for all analyses and that its effect will depend on the chosen metric and on the solutions used to deal with the problem. Based on our results, we suggest that ecologists should prefer the Operational tree solution to remove polytomies in the phylogenetic tree and take careful consideration while designing experiments, and analyzing and interpreting the results of phylogenetic metrics.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

Historical biogeography, systematics, and integrative taxonomy of the non-Ethiopian speckled pelage brush-furred rats (Lophuromys flavopunctatus group)

BACKGROUND

The speckled-pelage brush-furred rats (Lophuromys flavopunctatus group) have been difficult to define given conflicting genetic, morphological, and distributional records that combine to obscure meaningful accounts of its taxonomic diversity and evolution. In this study, we inferred the systematics, phylogeography, and evolutionary history of the L. flavopunctatus group using maximum likelihood and Bayesian phylogenetic inference, divergence times, historical biogeographic reconstruction, and morphometric discriminant tests. We compiled comprehensive datasets of three loci (two mitochondrial [mtDNA] and one nuclear) and two morphometric datasets (linear and geometric) from across the known range of the genus Lophuromys.

RESULTS

The mtDNA phylogeny supported the division of the genus Lophuromys into three primary groups with nearly equidistant pairwise differentiation: one group corresponding to the subgenus Kivumys (Kivumys group) and two groups corresponding to the subgenus Lophuromys (L. sikapusi group and L. flavopunctatus group). The L. flavopunctatus group comprised the speckled-pelage brush-furred Lophuromys endemic to Ethiopia (Ethiopian L. flavopunctatus members [ETHFLAVO]) and the non-Ethiopian ones (non-Ethiopian L. flavopunctatus members [NONETHFLAVO]) in deeply nested relationships. There were distinctly geographically structured mtDNA clades among the NONETHFLAVO, which were incongruous with the nuclear tree where several clades were unresolved. The morphometric datasets did not systematically assign samples to meaningful taxonomic units or agree with the mtDNA clades. The divergence dating and ancestral range reconstructions showed the NONETHFLAVO colonized the current ranges over two independent dispersal events out of Ethiopia in the early Pleistocene.

CONCLUSION

The phylogenetic associations and divergence times of the L. flavopunctatus group support the hypothesis that paleoclimatic impacts and ecosystem refugia during the Pleistocene impacted the evolutionary radiation of these rodents. The overlap in craniodental variation between distinct mtDNA clades among the NONETHFLAVO suggests unraveling underlying ecomorphological drivers is key to reconciling taxonomically informative morphological characters. The genus Lophuromys requires a taxonomic reassessment based on extensive genomic evidence to elucidate the patterns and impacts of genetic isolation at clade contact zones.

Morphological Characters Can Strongly Influence Early Animal Relationships Inferred from Phylogenomic Data Sets

There are considerable phylogenetic incongruencies between morphological and phylogenomic data for the deep evolution of animals. This has contributed to a heated debate over the earliest-branching lineage of the animal kingdom: the sister to all other Metazoa (SOM). Here, we use published phylogenomic data sets ($\sim $45,000-400,000 characters in size with $\sim $15-100 taxa) that focus on early metazoan phylogeny to evaluate the impact of incorporating morphological data sets ($\sim $15-275 characters). We additionally use small exemplar data sets to quantify how increased taxon sampling can help stabilize phylogenetic inferences. We apply a plethora of common methods, that is, likelihood models and their "equivalent" under parsimony: character weighting schemes. Our results are at odds with the typical view of phylogenomics, that is, that genomic-scale data sets will swamp out inferences from morphological data. Instead, weighting morphological data 2-10$\times $ in both likelihood and parsimony can in some cases "flip" which phylum is inferred to be the SOM. This typically results in the molecular hypothesis of Ctenophora as the SOM flipping to Porifera (or occasionally Placozoa). However, greater taxon sampling improves phylogenetic stability, with some of the larger molecular data sets ($>$200,000 characters and up to $\sim $100 taxa) showing node stability even with $\geqq100\times $ upweighting of morphological data. Accordingly, our analyses have three strong messages. 1) The assumption that genomic data will automatically "swamp out" morphological data is not always true for the SOM question. Morphological data have a strong influence in our analyses of combined data sets, even when outnumbered thousands of times by molecular data. Morphology therefore should not be counted out a priori. 2) We here quantify for the first time how the stability of the SOM node improves for several genomic data sets when the taxon sampling is increased. 3) The patterns of "flipping points" (i.e., the weighting of morphological data it takes to change the inferred SOM) carry information about the phylogenetic stability of matrices. The weighting space is an innovative way to assess comparability of data sets that could be developed into a new sensitivity analysis tool. [Metazoa; Morphology; Phylogenomics; Weighting.].

Global Diversification Dynamics Since the Jurassic: Low Dispersal and Habitat-Dependent Evolution Explain Hotspots of Diversity and Shell Disparity in River Snails (Viviparidae)

The Viviparidae, commonly known as River Snails, is a dominant group of freshwater snails with a nearly worldwide distribution that reaches its highest taxonomic and morphological diversity in Southeast Asia. The rich fossil record is indicative of a probable Middle Jurassic origin on the Laurasian supercontinent where the group started to diversify during the Cretaceous. However, it remains uncertain when and how the biodiversity hotspot in Southeast Asia was formed. Here, we used a comprehensive genetic data set containing both mitochondrial and nuclear markers and comprising species representing 24 out of 28 genera from throughout the range of the family. To reconstruct the spatiotemporal evolution of viviparids on a global scale, we reconstructed a fossil-calibrated phylogeny. We further assessed the roles of cladogenetic and anagenetic events in range evolution. Finally, we reconstructed the evolution of shell features by estimating ancestral character states to assess whether the appearance of sculptured shell morphologies was driven by major habitat shifts. The molecular phylogeny supports the monophyly of the three subfamilies, the Bellamyinae, Lioplacinae, and Viviparinae, but challenges the currently accepted genus-level classification in several cases. The almost global distribution of River Snails has been influenced both by comparatively ancient vicariance and more recent founder events. In Southeast Asia, Miocene dispersal was a main factor in shaping the modern species distributions. A recurrent theme across different viviparid taxa is that many species living in lentic waters exhibit sculptured shells, whereas only one strongly sculptured species is known from lotic environments. We show that such shell sculpture is habitat-dependent and indeed evolved several times independently in lentic River Snails. Considerably high transition rates between shell types in lentic habitats probably caused the co-occurrence of morphologically distinct shell types in several lakes. In contrast, directional evolution toward smooth shells in lotic habitats, as identified in the present analyses, explains why sculptured shells are rarely found in these habitats. However, the specific factors that promoted changes in shell morphology require further work. [biogeographical analyses; fossil-calibrated phylogeny; fossil-constrained analyses; Southeast Asia; stochastic character mapping.].

An ocean yet to be discovered: increasing systematic knowledge of Indo-Pacific Okenia Menke, 1830 (Nudibranchia:Goniodorididae)

Numerous faunistic and ecological studies have been conducted throughout the Indo-Pacific Ocean to assess its biodiversity. Despite the abundance of research, studies on the species that inhabit the Indo-Pacific are still necessary due to its extent and high species richness. The major species richness of the genus Okenia Menke, 1830 (Nudibranchia, Goniodorididae) is found in the Indo-Pacific Ocean, including 38 of 60 valid species. Nevertheless, this number does not represent the real biodiversity, since at least 20 more species are already reported in field-guides as undescribed species belonging to this genus. The systematics of the genus Okenia are still unclear since it has been the subject of only a few and incomplete studies. In the present paper, we describe five new Okenia species from the coastlines of Japan, Mozambique and Australia: Okenia aurorapapillata sp. nov., Okenia elisae sp. nov., Okenia nakanoae sp. nov., Okenia siderata sp. nov. and Okenia tenuifibrata sp. nov. Moreover, anatomical details not previously described of Okenia atkinsonorum, Okenia barnardi, Okenia cf. echinata, Okenia hallucigenia, Okenia hiroi, Okenia japonica, Okenia pellucida, Okenia pilosa and Okenia rhinorma are provided. New partial sequences of standard markers (COI, 16S rRNA and H3) were obtained and a phylogenetic analysis that included all species with available data was performed. ZOOBANK urn:lsid:zoobank.org:pub:28AE2536-A264–4194–8AE3-C430620572E7

Diverse phylogenomic datasets uncover a concordant scenario of laurasiatherian interordinal relationships

Resolving the interordinal relationships in the mammalian superorder Laurasiatheria has been among the most intractable problems in higher-level mammalian systematics, with many conflicting hypotheses having been proposed. The present study collected three different sources of genome-scale data with comprehensive taxon sampling of laurasiatherian species, including two protein-coding datasets (4,186 protein-coding genes for an amino acid dataset comprising 2,761,247 amino acid residues and a nucleotide dataset comprising 5,516,340 nucleotides from 1st and 2nd codon positions), an intronic dataset (1,210 introns comprising 1,162,723 nucleotides) and an ultraconserved elements (UCEs) dataset (1,246 UCEs comprising 1,946,472 nucleotides) from 40 species representing all six laurasiatherian orders and 7 non-laurasiatherian outgroups. Remarkably, phylogenetic trees reconstructed with the four datasets using different tree-building methods (RAxML, FastTree, ASTRAL and MP-EST) all supported the relationship (Eulipotyphla, (Chiroptera, ((Carnivora, Pholidota), (Cetartiodactyla, Perissodactyla)))). We find a resolution of interordinal relationships of Laurasiatheria among all types of markers used in the present study, and the likelihood ratio tests for tree comparisons confirmed that the present tree topology is the optimal hypothesis compared to other examined hypotheses. Jackknifing subsampling analyses demonstrate that the results of laurasiatherian tree reconstruction varied with the number of loci and ordinal representatives used, which are likely the two main contributors to phylogenetic disagreements of Laurasiatheria seen in previous studies. Our study provides significant insight into laurasiatherian evolution, and moreover, an important methodological strategy and reference for resolving phylogenies of adaptive radiation, which have been a long-standing challenge in the field of phylogenetics.

Enabling evolutionary studies at multiple scales in Apocynaceae through Hyb-Seq

PREMISE

Apocynaceae is the 10th largest flowering plant family and a focus for study of plant-insect interactions, especially as mediated by secondary metabolites. However, it has few genomic resources relative to its size. Target capture sequencing is a powerful approach for genome reduction that facilitates studies requiring data from the nuclear genome in non-model taxa, such as Apocynaceae.

METHODS

Transcriptomes were used to design probes for targeted sequencing of putatively single-copy nuclear genes across Apocynaceae. The sequences obtained were used to assess the success of the probe design, the intrageneric and intraspecific variation in the targeted genes, and the utility of the genes for inferring phylogeny.

RESULTS

From 853 candidate nuclear genes, 835 were consistently recovered in single copy and were variable enough for phylogenomics. The inferred gene trees were useful for coalescent-based species tree analysis, which showed all subfamilies of Apocynaceae as monophyletic, while also resolving relationships among species within the genus Apocynum. Intraspecific comparison of Elytropus chilensis individuals revealed numerous single-nucleotide polymorphisms with potential for use in population-level studies.

DISCUSSION

Community use of this Hyb-Seq probe set will facilitate and promote progress in the study of Apocynaceae across scales from population genomics to phylogenomics.

Competition for access to mates predicts female-specific ornamentation and male investment in relative testis size

Sexually selected ornaments are highly variable and the factors that drive variation in ornament expression are not always clear. Rare instances of female-specific ornament evolution (such as in some dance fly species) are particularly puzzling. While some evidence suggests that such rare instances represent straightforward reversals of sexual selection intensity, the distinct nature of trade-offs between ornaments and offspring pose special constraints in females. To examine whether competition for access to mates generally favors heightened ornament expression, we built a phylogeny and conducted a comparative analysis of Empidinae dance fly taxa that display female-specific ornaments. We show that species with more female-biased operational sex ratios in lek-like mating swarms have greater female ornamentation, and in taxa with more ornate females, male relative testis investment is increased. These findings support the hypothesis that ornament diversity in dance flies depends on female receptivity to mates, which is associated with contests for nutritious nuptial gifts provided by males. Moreover, our results suggest that increases in female receptivity lead to higher levels of sperm competition among males. The incidence of both heightened premating sexual selection on females and postmating selection on males contradicts assertions that sex roles are straightforwardly reversed in dance flies.

Molecular phylogenetic analysis of Punctoidea (Gastropoda, Stylommatophora)

A phylogenetic analysis using a combination of mitochondrial (COI, 16S) and nuclear markers (ITS2, 28S) indicated that Punctoidea, as previously interpreted, is polyphyletic. It comprises two main groups, containing northern hemisphere (Laurasian) and predominantly southern hemisphere (Gondwanan) taxa respectively, treated here as separate superfamilies. Within Punctoidea sensu stricto, Punctidae, Cystopeltidae and Endodontidae form separate monophyletic clades, but Charopidae, as currently interpreted, is paraphyletic. Most of the charopid taxa that we sequenced, including Charopa coma (Gray, 1843) and other Charopinae, grouped in a clade with Punctidae but some charopid taxa from Australia and South America grouped with Cystopeltidae. Cystopeltidae previously contained a single Australia-endemic genus, Cystopelta Tate, 1881, but our analysis suggests that it is considerably more diverse taxonomically and has a much wider distribution. For taxonomic stability, we suggest that Charopidae be retained as a family-level group for now, pending further study of the systematic relationships of its constituent taxa. A new superfamily, Discoidea, is erected here for two Northern Hemisphere families, Discidae and Oreohelicidae, which were previously assigned to Punctoidea. The North American species Radiodomus abietum, previously in Charopidae, is also here assigned to Discoidea. The phylogenetic relationships of Helicodiscidae, previously assigned to Punctoidea, were not fully resolved in our analysis, but the family is apparently closely related to Arionoidea Gray, 1840 and infraorder Limacoidei.

The Mitochondrial Genome of the Phytopathogenic Fungus Bipolaris sorokiniana and the Utility of Mitochondrial Genome to Infer Phylogeny of Dothideomycetes

A number of species in Bipolaris are important plant pathogens. Due to a limited number of synapomorphic characters, it is difficult to perform species identification and to estimate phylogeny of Bipolaris based solely on morphology. In this study, we sequenced the complete mitochondrial genome of Bipolaris sorokiniana, and presented the detailed annotation of the genome. The B. sorokiniana mitochondrial genome is 137,775 bp long, and contains two ribosomal RNA genes, 12 core protein-coding genes, 38 tRNA genes. In addition, two ribosomal protein genes (rps3 gene and rps5 gene) and the fungal mitochondrial RNase P gene (rnpB) are identified. The large genome size is mostly determined by the presence of numerous intronic and intergenic regions. A total of 28 introns are inserted in eight core protein-coding genes. Together with the published mitochondrial genome sequences, we conducted a preliminary phylogenetic inference of Dothideomycetes under various datasets and substitution models. The monophyly of Capnodiales, Botryosphaeriales and Pleosporales are consistently supported in all analyses. The Venturiaceae forms an independent lineage, with a distant phylogenetic relationship to Pleosporales. At the family level, the Mycosphaerellaceae, Botryosphaeriaceae. Phaeosphaeriaceae, and Pleosporaceae are recognized in the majority of trees.