The Prevalence and Impact of Model Violations in Phylogenetic Analysis

Phylogenomics and biogeography of sawflies and woodwasps (Hymenoptera, Symphyta)

Phylogenomic approaches have recently helped elucidate various insect relationships, but large-scale comprehensive analyses on relationships within sawflies and woodwasps are still lacking. Here, we infer the relationships and long-term biogeographic history of these hymenopteran groups using a large dataset of 354 UCE loci collected from 385 species that represent all major lineages. Early Hymenoptera started diversifying during the Early Triassic ∼249 Ma and spread all over the ancient supercontinent Pangaea. We recovered Xyeloidea as a monophyletic sister group to other Hymenoptera and Pamphilioidea as sister to Unicalcarida. Within the diverse family Tenthredinidae, our taxonomically and geographically expanded taxon sampling highlights the non-monophyly of several traditionally defined subfamilies. In addition, the recent removal of Athalia and related genera from the Tenthredinidae into the separate family Athaliidae is supported. The deep historical biogeography of the group is characterised by independent dispersals and re-colonisations between the northern (Laurasia) and southern (Gondwana) palaeocontinents. The breakup of these landmasses led to ancient vicariance in several Gondwanan lineages, while interchange across the Northern Hemisphere has continued until the Recent. The little-studied African sawfly fauna is likewise a diverse mixture of groups with varying routes of colonization. Our results reveal interesting parallels in the evolution and biogeography of early hymenopterans and other ancient insect groups.

Phylogenomic analyses of Blattodea combining traditional methods, incremental tree-building, and quality-aware support

Despite the many advances of the genomic era, there is a persistent problem in assessing the uncertainty of phylogenomic hypotheses. We see this in the recent history of phylogenetics for cockroaches and termites (Blattodea), where huge advances have been made, but there are still major inconsistencies between studies. To address this, we present a phylogenetic analysis of Blattodea that emphasizes identification and quantification of uncertainty. We analyze 1183 gene domains using three methods (multi-species coalescent inference, concatenation, and a supermatrix-supertree hybrid approach) and assess support for controversial relationships while considering data quality. The hybrid approach-here dubbed "tiered phylogenetic inference"-incorporates information about data quality into an incremental tree building framework. Leveraging this method, we are able to identify cases of low or misleading support that would not be possible otherwise, and explore them more thoroughly with follow-up tests. In particular, quality annotations pointed towards nodes with high bootstrap support that later turned out to have large ambiguities, sometimes resulting from low-quality data. We also clarify issues related to some recalcitrant nodes: Anaplectidae's placement lacks unbiased signal, Ectobiidae s.s. and Anaplectoideini need greater taxon sampling, the deepest relationships among most Blaberidae lack signal. As a result, several previous phylogenetic uncertainties are now closer to being resolved (e.g., African and Malagasy "Rhabdoblatta" spp. are the sister to all other Blaberidae, and Oxyhaloinae is sister to the remaining Blaberidae). Overall, we argue for more approaches to quantifying support that take data quality into account to uncover the nature of recalcitrant nodes.

Whole-genome-based phylogenetic analyses provide new insights into the evolution of springtails (Hexapoda: Collembola)

Springtails (Collembola) stand as one of the most abundant, widespread, and ancient terrestrial arthropods on earth. However, their evolutionary history and deep phylogenetic relationships remain elusive. In this study, we employed phylogenomic approaches to elucidate the basal relationships among Collembola. We sampled whole-genome data representing all major collembolan lineages in proportion to their known diversity. To account for potential phylogenomic biases, we implemented various data extraction, locus sampling, and signal filtering strategies to generate matrices. Subsequently, we applied a diverse array of tree-searching and rate-modelling methods to reconstruct the phylogeny. Our analyses, utilizing different matrices and methods, converged on the same unrooted relationships among collembolan ingroups, supporting the current ordinal classification and challenging the monophyly of Arthropleona and Symphypleona s.l. However, discrepancies across analyses existed in the root of Collembola. Among various root positions, those based on more informative matrices and biologically realistic models, favoring a basal topology of Entomobryomorpha + (Symphypleona s.s. + (Neelipleona + Poduromorpha)), were supported by subsequent methodological assessment, topology tests, and rooting analyses. This optimal topology suggests multiple independent reduction of the pronotum in non-poduromorph orders and aligns with the plesiomorphic status of neuroendocrine organs and epicuticular structure of Entomobryomorpha. Fossil-calibrated dating analyses based on the optimal topology indicated late-Paleozoic to mid-Mesozoic origins of the crown Collembola and four orders. In addition, our results questioned the monophyly of Isotomidae and Neanuridae, underscoring the need for further attention to the systematics of these families. Overall, this study provides novel insights into the phylogenetic backbone of Collembola, which will inform future studies on the systematics, ecology, and evolution of this significant arthropod lineage.

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Genetic characterization of the first Deltacoronavirus from wild birds around Qinghai Lake

Deltacoronavirus, widely distributed among pigs and wild birds, pose a significant risk of cross-species transmission, including potential human epidemics. Metagenomic analysis of bird samples from Qinghai Lake, China in 2021 reported the presence of Deltacoronavirus. A specific gene fragment of Deltacoronavirus was detected in fecal samples from wild birds at a positive rate of 5.94% (6/101). Next-generation sequencing (NGS) identified a novel Deltacoronavirus strain, which was closely related to isolates from the United Arab Emirates (2018), China (2022), and Poland (2023). Subsequently the strain was named A/black-headed gull/Qinghai/2021(BHG-QH-2021) upon confirmation of the Cytochrome b gene of black-headed gull in the sample. All available genome sequences of avian Deltacoronavirus, including the newly identified BHG-QH-2021 and 5 representative strains of porcine Deltacoronavirus (PDCoV), were classified according to ICTV criteria. In contrast to Coronavirus HKU15, which infects both mammals and birds and shows the possibility of cross-species transmission from bird to mammal host, our analysis revealed that BHG-QH-2021 is classified as Putative species 4. Putative species 4 has been reported to infect 5 species of birds but not mammals, suggesting that cross-species transmission of Putative species 4 is more prevalent among birds. Recombination analysis traced BHG-QH-2021 origin to dut148cor1 and MW01_1o strains, with MW01_1o contributing the S gene. Surprisingly, SwissModle prediction showed that the optimal template for receptor-binding domain (RBD) of BHG-QH-2021 is derived from the human coronavirus 229E, a member of the Alphacoronavirus, rather than the anticipated RBD structure of PDCoV of Deltacoronavirus. Further molecular docking analysis revealed that substituting the loop 1-2 segments of HCoV-229E significantly enhanced the binding capability of BHG-QH-2021 with human Aminopeptidase N (hAPN), surpassing its native receptor-binding domain (RBD). Most importantly, this finding was further confirmed by co-immunoprecipitation experiment that loop 1-2 segments of HCoV-229E enable BHG-QH-2021 RBD binding to hAPN, indicating that the loop 1-2 segment of the RBD in Putative species 4 is a probable key determinant for the virus ability to spill over into humans. Our results summarize the phylogenetic relationships among known Deltacoronavirus, reveal an independent putative avian Deltacoronavirus species with inter-continental and inter-species transmission potential, and underscore the importance of continuous surveillance of wildlife Deltacoronavirus.

Synergistic innovations enabled the radiation of anglerfishes in the deep open ocean

Major ecological transitions are thought to fuel diversification, but whether they are contingent on the evolution of certain traits called key innovations1 is unclear. Key innovations are routinely invoked to explain how lineages rapidly exploit new ecological opportunities.1,2,3 However, investigations of key innovations often focus on single traits rather than considering trait combinations that collectively produce effects of interest.4 Here, we investigate the evolution of synergistic trait interactions in anglerfishes, which include one of the most species-rich vertebrate clades in the bathypelagic, or "midnight," zone of the deep sea: Ceratioidea.5 Ceratioids are the only vertebrates that possess sexual parasitism, wherein males temporarily attach or permanently fuse to females to mate.6,7 We show that the rapid transition of ancestrally benthic anglerfishes into pelagic habitats occurred during a period of major global warming 50-35 million years ago.8,9 This transition coincided with the origins of sexual parasitism, which is thought to increase the probability of successful reproduction once a mate is found in the midnight zone, Earth's largest habitat.5,6,7 Our reconstruction of the evolutionary history of anglerfishes and the loss of immune genes support that permanently fusing clades have convergently degenerated their adaptive immunity. We find that degenerate adaptive immune genes and sexual body size dimorphism, both variably present in anglerfishes outside the ceratioid radiation, likely promoted their transition into the bathypelagic zone. These results show how traits from separate physiological, morphological, and reproductive systems can interact synergistically to drive major transitions and subsequent diversification in novel environments.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata

Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.

A novel and diverse family of filamentous DNA viruses associated with parasitic wasps

Large dsDNA viruses from the Naldaviricetes class are currently composed of four viral families infecting insects and/or crustaceans. Since the 1970s, particles described as filamentous viruses (FVs) have been observed by electronic microscopy in several species of Hymenoptera parasitoids but until recently, no genomic data was available. This study provides the first comparative morphological and genomic analysis of these FVs. We analyzed the genomes of seven FVs, six of which were newly obtained, to gain a better understanding of their evolutionary history. We show that these FVs share all genomic features of the Naldaviricetes while encoding five specific core genes that distinguish them from their closest relatives, the Hytrosaviruses. By mining public databases, we show that FVs preferentially infect Hymenoptera with parasitoid lifestyle and that these viruses have been repeatedly integrated into the genome of many insects, particularly Hymenoptera parasitoids, overall suggesting a long-standing specialization of these viruses to parasitic wasps. Finally, we propose a taxonomical revision of the class Naldaviricetes in which FVs related to the Leptopilina boulardi FV constitute a fifth family. We propose to name this new family, Filamentoviridae.

Ordovician origin and subsequent diversification of the brown algae

Brown algae are the only group of heterokont protists exhibiting complex multicellularity. Since their origin, brown algae have adapted to various marine habitats, evolving diverse thallus morphologies and gamete types. However, the evolutionary processes behind these transitions remain unclear due to a lack of a robust phylogenetic framework and problems with time estimation. To address these issues, we employed plastid genome data from 138 species, including heterokont algae, red algae, and other red-derived algae. Based on a robust phylogeny and new interpretations of algal fossils, we estimated the geological times for brown algal origin and diversification. The results reveal that brown algae first evolved true multicellularity, with plasmodesmata and reproductive cell differentiation, during the late Ordovician Period (ca. 450 Ma), coinciding with a major diversification of marine fauna (the Great Ordovician Biodiversification Event) and a proliferation of multicellular green algae. Despite its early Paleozoic origin, the diversification of major orders within this brown algal clade accelerated only during the Mesozoic Era, coincident with both Pangea rifting and the diversification of other heterokont algae (e.g., diatoms), coccolithophores, and dinoflagellates, with their red algal-derived plastids. The transition from ancestral isogamy to oogamy was followed by three simultaneous reappearances of isogamy during the Cretaceous Period. These are concordant with a positive character correlation between parthenogenesis and isogamy. Our new brown algal timeline, combined with a knowledge of past environmental conditions, shed new light on brown algal diversification and the intertwined evolution of multicellularity and sexual reproduction.

Forty Years of Inferential Methods in the Journals of the Society for Molecular Biology and Evolution

We are launching a series to celebrate the 40th anniversary of the first issue of Molecular Biology and Evolution. In 2024, we will publish virtual issues containing selected papers published in the Society for Molecular Biology and Evolution journals, Molecular Biology and Evolution and Genome Biology and Evolution. Each virtual issue will be accompanied by a perspective that highlights the historic and contemporary contributions of our journals to a specific topic in molecular evolution. This perspective, the first in the series, presents an account of the broad array of methods that have been published in the Society for Molecular Biology and Evolution journals, including methods to infer phylogenies, to test hypotheses in a phylogenetic framework, and to infer population genetic processes. We also mention many of the software implementations that make methods tractable for empiricists. In short, the Society for Molecular Biology and Evolution community has much to celebrate after four decades of publishing high-quality science including numerous important inferential methods.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Adaptive RAxML-NG: Accelerating Phylogenetic Inference under Maximum Likelihood using Dataset Difficulty

Phylogenetic inferences under the maximum likelihood criterion deploy heuristic tree search strategies to explore the vast search space. Depending on the input dataset, searches from different starting trees might all converge to a single tree topology. Often, though, distinct searches infer multiple topologies with large log-likelihood score differences or yield topologically highly distinct, yet almost equally likely, trees. Recently, Haag et al. introduced an approach to quantify, and implemented machine learning methods to predict, the dataset difficulty with respect to phylogenetic inference. Easy multiple sequence alignments (MSAs) exhibit a single likelihood peak on their likelihood surface, associated with a single tree topology to which most, if not all, independent searches rapidly converge. As difficulty increases, multiple locally optimal likelihood peaks emerge, yet from highly distinct topologies. To make use of this information, we introduce and implement an adaptive tree search heuristic in RAxML-NG, which modifies the thoroughness of the tree search strategy as a function of the predicted difficulty. Our adaptive strategy is based upon three observations. First, on easy datasets, searches converge rapidly and can hence be terminated at an earlier stage. Second, overanalyzing difficult datasets is hopeless, and thus it suffices to quickly infer only one of the numerous almost equally likely topologies to reduce overall execution time. Third, more extensive searches are justified and required on datasets with intermediate difficulty. While the likelihood surface exhibits multiple locally optimal peaks in this case, a small proportion of them is significantly better. Our experimental results for the adaptive heuristic on 9,515 empirical and 5,000 simulated datasets with varying difficulty exhibit substantial speedups, especially on easy and difficult datasets (53% of total MSAs), where we observe average speedups of more than 10×. Further, approximately 94% of the inferred trees using the adaptive strategy are statistically indistinguishable from the trees inferred under the standard strategy (RAxML-NG).

Phylogeny of Rhynchium and Its Related Genera (Hymenoptera: Eumeninae) Based on Universal Single-Copy Orthologs and Ultraconserved Elements

The subfamily Eumeninae is a large group of fierce predatory insects that prey mainly on the larvae of Lepidoptera pests. Because of the highly similar morphologies of the genus Rhynchium and its related genera in the subfamily, including Rhynchium Spinola, Allorhynchium van der Vecht, Anterhynchium de Saussure, Pararrhynchium de Saussure, it is essential to delineate their relationships. A previous phylogenetic analysis based on mitochondrial genomes suggested the inconsistent relationships of these genera under traditional classification based on morphological characters. In this study, we first used single-copy orthologs [USCO] and ultraconserved elements [UCE] extracted from 10 newly sequenced low-coverage whole genomes to resolve the phylogenetic relationships of the above genera. The newly sequenced genomes are 152.99 Mb to 211.49 Mb in size with high completeness (BUSCO complete: 91.5-95.6%) and G + C content (36.31-38.76%). Based on extracted 5811 USCOs and 2312 UCEs, the phylogenetic relationships of Rhynchium and its related genera were: ((Allorhynchium + Lissodynerus) + (Pararrhynchium + (Anterhynchium + (Dirhynchium + Rhynchium)))), which was consistent with the mitochondrial genome results. The results supported the genus Rhynchium as monophyletic, whereas Anterhynchium was recovered as paraphyletic, with Anterhynchium (Dirhynchium) as a sister to Rhynchium and hence deserving generic status; In addition, in the genus Pararrhynchium, P. septemfasciatus feanus and P. venkataramani were separated, not clustered on a branch, just as P. septemfasciatus feanus was not together with P. striatum based on mitochondrial genomes. Since Lissodynerus septemfasciatus, the type species of the genus Lissodynerus, was transferred to Pararrhynchium, it is considered that the genus Lissodynerus should be restituted as a valid genus, not a synonym of Pararrhynchium.

A global catalog of whole-genome diversity from 233 primate species

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

Biogeography and evolution of social parasitism in Australian Myrmecia bulldog ants revealed by phylogenomics

Studying the historical biogeography and life history transitions from eusocial colony life to social parasitism contributes to our understanding of the evolutionary mechanisms generating biodiversity in eusocial insects. The ants in the genus Myrmecia are a well-suited system for testing evolutionary hypotheses about how their species diversity was assembled through time because the genus is endemic to Australia with the single exception of the species M. apicalis inhabiting the Pacific Island of New Caledonia, and because at least one social parasite species exists in the genus. However, the evolutionary mechanisms underlying the disjunct biogeographic distribution of M. apicalis and the life history transition(s) to social parasitism remain unexplored. To study the biogeographic origin of the isolated, oceanic species M. apicalis and to reveal the origin and evolution of social parasitism in the genus, we reconstructed a comprehensive phylogeny of the ant subfamily Myrmeciinae. We utilized Ultra Conserved Elements (UCEs) as molecular markers to generate a comprehensive molecular genetic dataset consisting of 2,287 loci per taxon on average for 66 out of the 93 known Myrmecia species as well as for the sister lineage Nothomyrmecia macrops and selected outgroups. Our time-calibrated phylogeny inferred that: (i) stem Myrmeciinae originated during the Paleocene ∼58 Ma ago; (ii) the current disjunct biogeographic distribution of M. apicalis was driven by long-distance dispersal from Australia to New Caledonia during the Miocene ∼14 Ma ago; (iii) the single social parasite species, M. inquilina, evolved directly from one of the two known host species, M. nigriceps, in sympatry via the intraspecific route of social parasite evolution; and (iv) 5 of the 9 previously established taxonomic species groups are non-monophyletic. We suggest minor changes to reconcile the molecular phylogenetic results with the taxonomic classification. Our study enhances our understanding of the evolution and biogeography of Australian bulldog ants, contributes to our knowledge about the evolution of social parasitism in ants, and provides a solid phylogenetic foundation for future inquiries into the biology, taxonomy, and classification of Myrmeciinae.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Chromosome-level genome assembly of bean flower thrips Megalurothrips usitatus (Thysanoptera: Thripidae)

Bean flower thrips Megalurothrips usitatus is a staple pest of cowpea and other legumes and causes dramatic economic losses. Its small size allows for easy concealment, and large reproductive capacity easily leads to infestations. Despite the importance of a genome in developing novel management strategies, genetic studies on M. usitatus remain limited. Thus, we generated a chromosome-level M. usitatus genome using a combination of PacBio long read and Hi-C technologies. The assembled genome was 238.14 Mb with a scaffold N50 of 13.85 Mb. The final genome was anchored into 16 pseudo-chromosomes containing 14,000 genes, of which 91.74% were functionally annotated. Comparative genomic analyses revealed that expanded gene families were enriched in fatty acid metabolism and detoxification metabolism (ABC transporters), and contracted gene families were strongly associated with chitin-based cuticle development and sensory perception of taste. In conclusion, this high-quality genome provides an invaluable resource for us to understand the thrips' ecology and genetics, contributing to pest management.

nRCFV: a new, dataset-size-independent metric to quantify compositional heterogeneity in nucleotide and amino acid datasets

MOTIVATION

Compositional heterogeneity-when the proportions of nucleotides and amino acids are not broadly similar across the dataset-is a cause of a great number of phylogenetic artefacts. Whilst a variety of methods can identify it post-hoc, few metrics exist to quantify compositional heterogeneity prior to the computationally intensive task of phylogenetic tree reconstruction. Here we assess the efficacy of one such existing, widely used, metric: Relative Composition Frequency Variability (RCFV), using both real and simulated data.

RESULTS

Our results show that RCFV can be biased by sequence length, the number of taxa, and the number of possible character states within the dataset. However, we also find that missing data does not appear to have an appreciable effect on RCFV. We discuss the theory behind this, the consequences of this for the future of the usage of the RCFV value and propose a new metric, nRCFV, which accounts for these biases. Alongside this, we present a new software that calculates both RCFV and nRCFV, called nRCFV_Reader.

AVAILABILITY AND IMPLEMENTATION

nRCFV has been implemented in RCFV_Reader, available at: https://github.com/JFFleming/RCFV_Reader . Both our simulation and real data are available at Datadryad: https://doi.org/10.5061/dryad.wpzgmsbpn .

Evolution of stridulatory mechanisms: vibroacoustic communication may be common in leaf-footed bugs and allies (Heteroptera: Coreoidea)

Intra- and interspecific communication is crucial to fitness via its role in facilitating mating, territoriality and defence. Yet, the evolution of animal communication systems is puzzling-how do they originate and change over time? Studying stridulatory morphology provides a tractable opportunity to deduce the origin and diversification of a communication mechanism. Stridulation occurs when two sclerotized structures rub together to produce vibratory and acoustic (vibroacoustic) signals, such as a cricket 'chirp'. We investigated the evolution of stridulatory mechanisms in the superfamily Coreoidea (Hemiptera: Heteroptera), a group of insects known for elaborate male fighting behaviours and enlarged hindlegs. We surveyed a large sampling of taxa and used a phylogenomic dataset to investigate the evolution of stridulatory mechanisms. We identified four mechanisms, with at least five evolutionary gains. One mechanism, occurring only in male Harmostini (Rhopalidae), is described for the first time. Some stridulatory mechanisms appear to be non-homoplastic apomorphies within Rhopalidae, while others are homoplastic or potentially homoplastic within Coreidae and Alydidae, respectively. We detected no losses of these mechanisms once evolved, suggesting they are adaptive. Our work sets the stage for further behavioural, evolutionary and ecological studies to better understand the context in which these traits evolve and change.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Roadmap to the study of gene and protein phylogeny and evolution-A practical guide

Developments in sequencing technologies and the sequencing of an ever-increasing number of genomes have revolutionised studies of biodiversity and organismal evolution. This accumulation of data has been paralleled by the creation of numerous public biological databases through which the scientific community can mine the sequences and annotations of genomes, transcriptomes, and proteomes of multiple species. However, to find the appropriate databases and bioinformatic tools for respective inquiries and aims can be challenging. Here, we present a compilation of DNA and protein databases, as well as bioinformatic tools for phylogenetic reconstruction and a wide range of studies on molecular evolution. We provide a protocol for information extraction from biological databases and simple phylogenetic reconstruction using probabilistic and distance methods, facilitating the study of biodiversity and evolution at the molecular level for the broad scientific community.

Fusarioid community diversity associated with conifer seedlings in forest nurseries across the contiguous USA

INTRODUCTION

Fusarioid fungi that cause damping-off and root diseases can result in significant losses to conifer crops produced in forest nurseries across the USA. These nurseries are vital to reforestation and forest restoration efforts. Understanding the diversity of Fusarioid fungi associated with damping-off and root diseases of conifer seedlings can provide an approach for targeted management techniques to limit seedling losses and pathogen spread to novel landscapes.

METHODS

This study identifies 26 Fusarium spp. (F. acuminatum, F. annulatum, F. avenaceum, F. brachygibbosum, F. clavus, F. commune, F. cugenangense, F. diversisporum, F. elaeagni, F. elaeidis, F. flocciferum, F. fredkrugeri, F. fujikuroi, F. grosmichelii, F. ipomoeae, F. lactis, F. languescens, F. luffae, F. odoratissimum, F. oxysporum, F. queenslandicum, F. redolens, F. torulosum, F. triseptatum, F. vanleeuwenii, & F. verticillioides), 15 potential species within Fusarium and Neocosmospora species complexes (two from F. fujikuroi species complex, nine from F. oxysporum species complex, three from F. tricinctum species complex, and one from Neocosmospora species complex), and four Neocosmospora spp. (N. falciforme, N. metavorans, N. pisi, & N. solani) and associated host information collected from conifer-producing nurseries across the contiguous USA.

RESULTS

Phylogenetic analyses identified Fusarioid fungi haplotypes that were associated with 1) host specificity, 2) localization to geographic regions, or 3) generalists found on multiple hosts across diverse geographic regions.

DISCUSSION

The haplotypes and novel species identified on conifer seedlings should be considered for further analysis to determine pathogenicity, pathogen spread, and assess management practices.

Pseudocercospora pancratii Causing Leaf Spots on Commercial Blackberry (Rubus sp.) in Florida

Blackberry (Rubus L. subgenus Rubus Watson) is a deciduous berry crop that is the fourth most economically important berry crop, and its production is expanding in the southeastern United States. However, since most commercially available cultivars were bred under temperate conditions, they are not always well adapted and could be threatened by new pathogen populations inhabiting subtropical areas. In 2017, plants showing purple or brown leaf spots and angular-to-irregular lesions on both leaf surfaces, with clusters of black conidiophores at the center, were observed in a field trial at the University of Florida's Gulf Coast Research and Education Center (UF/GCREC) in Wimauma, FL. A fungus resembling Cercospora/Pseudocercospora was isolated from the lesions. The ribosomal DNA internal transcribed spacers, the translation elongation factor 1-alpha, and the actin genes were amplified and sequenced. Based on the phylogenetic analysis, the closest related species was Pseudocercospora pancratii. Pathogenicity assays and subsequent reisolation confirmed that this species is the causal agent of the disease. Among eight cultivars screened, no complete resistance was found. However, 'Osage' was the least susceptible, and 'Kiowa' was the most susceptible. This study is the first report of P. pancratii causing leaf spots on blackberry worldwide, and it may help shape future research into disease epidemiology and management for a crop that is rapidly expanding but has very limited disease information currently available for Florida growers.

Testing the systematic status of Homalictus and Rostrohalictus with weakened cross-vein groups within Halictini (Hymenoptera: Halictidae) using low-coverage whole-genome sequencing

The halictid genus Lasioglossum, as one of the most species-rich bee groups with persistently contentious subgeneric boundaries, is one of the most challenging bee groups from a systematic standpoint. An enduring question is the relationship of Lasioglossum and Homalictus, whether all halictine bees with weakened distal wing venation comprise one or multiple genera. Here, we analyzed the phylogenetic relationships among the subgroups within Lasioglossum s.l. based on thousands of single-copy orthologs and ultraconserved elements, which were extracted from 23 newly sequenced low-coverage whole genomes alongside a published genome (22 ingroups plus 2 outgroups). Both marker sets provided consistent results across maximum likelihood and coalescent-based species tree approaches. The phylogenetic and topology test results show that the Lasioglossum and Hemihalictus series are reciprocally monophyletic and Homalictus and Rostrohalictus are valid subgenera of Lasioglossum. Consequently, we lower Homalictus to subgenus status within Lasioglossum again, and we also raise Rostrohalictus to subgenus status from its prior synonymy with subgenus Hemihalictus. Lasioglossum przewalskyi is also transferred to the subgenus Hemihalictus. Ultimately, we redefine Lasioglossum to include all halictine bees with weakened distal wing venation.

Radiation and hybridization underpin the spread of the fire ant social supergene

Supergenes are clusters of tightly linked genes that jointly produce complex phenotypes. Although widespread in nature, how such genomic elements are formed and how they spread are in most cases unclear. In the fire ant Solenopsis invicta and closely related species, a "social supergene controls whether a colony maintains one or multiple queens. Here, we show that the three inversions constituting the Social b (Sb) supergene emerged sequentially during the separation of the ancestral lineages of S. invicta and Solenopsis richteri. The two first inversions arose in the ancestral population of both species, while the third one arose in the S. richteri lineage. Once completely assembled in the S. richteri lineage, the supergene first introgressed into S. invicta, and from there into the other species of the socially polymorphic group of South American fire ant species. Surprisingly, the introgression of this large and important genomic element occurred despite recent hybridization being uncommon between several of the species. These results highlight how supergenes can readily move across species boundaries, possibly because of fitness benefits they provide and/or expression of selfish properties favoring their transmission.

Assessment of Absolute Substitution Model Fit Accommodating Time-Reversible and Non-Time-Reversible Evolutionary Processes

The loss of information accompanying assessment of absolute fit of substitution models to phylogenetic data negatively affects the discriminatory power of previous methods and can make them insensitive to lineage-specific changes in the substitution process. As an alternative, I propose evaluating absolute fit of substitution models based on a novel statistic which describes the observed data without information loss and which is unlikely to become zero-inflated with increasing numbers of taxa. This method can accommodate gaps and is sensitive to lineage-specific shifts in the substitution process. In simulation experiments, it exhibits greater discriminatory power than previous methods. The method can be implemented in both Bayesian and Maximum Likelihood phylogenetic analyses, and used to screen any set of models. Recently, it has been suggested that model selection may be an unnecessary step in phylogenetic inference. However, results presented here emphasize the importance of model fit assessment for reliable phylogenetic inference.

Recoding amino acids to a reduced alphabet may increase or decrease phylogenetic accuracy

Common molecular phylogenetic characteristics such as long branches and compositional heterogeneity can be problematic for phylogenetic reconstruction when using amino acid data. Recoding alignments to reduced alphabets before phylogenetic analysis has often been used both to explore and potentially decrease the effect of such problems. We tested the effectiveness of this strategy on topological accuracy using simulated data on four-taxon trees. We simulated alignments in phylogenetically challenging ways to test the phylogenetic accuracy of analyses using various recoding strategies together with commonly-used homogeneous models. We tested three recoding methods based on amino acid exchangeability, and another recoding method based on lowering the compositional heterogeneity among alignment sequences as measured by the Chi-squared statistic. Our simulation results show that on trees with long branches where sequences approach saturation, accuracy was not greatly affected by exchangeability-based recodings, but Chi-squared-based recoding decreased accuracy. We then simulated sequences with different kinds of compositional heterogeneity over the tree. Recoding often increased accuracy on such alignments. Exchangeability-based recoding was rarely worse than not recoding, and often considerably better. Recoding based on lowering the Chi-squared value improved accuracy in some cases but not in others, suggesting that low compositional heterogeneity by itself is not sufficient to increase accuracy in the analysis of these alignments. We also simulated alignments using site-specific amino acid profiles, making sequences that had compositional heterogeneity over alignment sites. Exchangeability-based recoding coupled with site-homogeneous models had poor accuracy for these datasets but Chi-squared-based recoding on these alignments increased accuracy. We then simulated datasets that were compositionally both site- and tree-heterogeneous, like many real datasets. The effect on accuracy of recoding such doubly problematic datasets varied widely, depending on the type of compositional tree-heterogeneity and on the recoding scheme. Interestingly, analysis of unrecoded compositionally heterogeneous alignments with the NDCH or CAT models was generally more accurate than homogeneous analysis, whether recoded or not. Overall, our results suggest that making trees for recoded amino acid datasets can be useful, but they need to be interpreted cautiously as part of a more comprehensive analysis. The use of better fitting models like NDCH and CAT, which directly account for the patterns in the data, may offer a more promising long-term solution for analysing empirical data.

Ecogenomics sheds light on diverse lifestyle strategies in freshwater CPR

BACKGROUND

The increased use of metagenomics and single-cell genomics led to the discovery of organisms from phyla with no cultivated representatives and proposed new microbial lineages such as the candidate phyla radiation (CPR or Patescibacteria). These bacteria have peculiar ribosomal structures, reduced metabolic capacities, small genome, and cell sizes, and a general host-associated lifestyle was proposed for the radiation. So far, most CPR genomes were obtained from groundwaters; however, their diversity, abundance, and role in surface freshwaters is largely unexplored. Here, we attempt to close these knowledge gaps by deep metagenomic sequencing of 119 samples of 17 different freshwater lakes located in Europe and Asia. Moreover, we applied Fluorescence in situ Hybridization followed by Catalyzed Reporter Deposition (CARD-FISH) for a first visualization of distinct CPR lineages in freshwater samples.

RESULTS

A total of 174 dereplicated metagenome-assembled genomes (MAGs) of diverse CPR lineages were recovered from the investigated lakes, with a higher prevalence from hypolimnion samples (162 MAGs). They have reduced genomes (median size 1 Mbp) and were generally found in low abundances (0.02-14.36 coverage/Gb) and with estimated slow replication rates. The analysis of genomic traits and CARD-FISH results showed that the radiation is an eclectic group in terms of metabolic capabilities and potential lifestyles, ranging from what appear to be free-living lineages to host- or particle-associated groups. Although some complexes of the electron transport chain were present in the CPR MAGs, together with ion-pumping rhodopsins and heliorhodopsins, we believe that they most probably adopt a fermentative metabolism. Terminal oxidases might function in O₂ scavenging, while heliorhodopsins could be involved in mitigation against oxidative stress.

CONCLUSIONS

A high diversity of CPR MAGs was recovered, and distinct CPR lineages did not seem to be limited to lakes with specific trophic states. Their reduced metabolic capacities resemble the ones described for genomes in groundwater and animal-associated samples, apart from Gracilibacteria that possesses more complete metabolic pathways. Even though this radiation is mostly host-associated, we also observed organisms from different clades (ABY1, Paceibacteria, Saccharimonadia) that appear to be unattached to any other organisms or were associated with 'lake snow' particles (ABY1, Gracilibacteria), suggesting a broad range of potential life-strategies in this phylum. Video Abstract.

Ant phylogenomics reveals a natural selection hotspot preceding the origin of complex eusociality

The evolution of eusociality has allowed ants to become one of the most conspicuous and ecologically dominant groups of organisms in the world. A large majority of the current ∼14,000 ant species belong to the formicoids,¹ a clade of nine subfamilies that exhibit the most extreme forms of reproductive division of labor, large colony size,² worker polymorphism,³ and extended queen longevity.⁴ The eight remaining non-formicoid subfamilies are less well studied, with few genomes having been sequenced so far and unclear phylogenetic relationships.⁵ By sequencing 65 genomes, we provide a robust phylogeny of the 17 ant subfamilies, retrieving high support to the controversial leptanillomorph clade (Leptanillinae and Martialinae) as the sister group to all other extant ants. Moreover, our genomic analyses revealed that the emergence of the formicoids was accompanied by an elevated number of positive selection events. Importantly, the top three gene functions under selection are linked to key features of complex eusociality, with histone acetylation being implicated in caste differentiation, gene silencing by RNA in worker sterility, and autophagy in longevity. These results show that the key pathways associated with eusociality have been under strong selection during the Cretaceous, suggesting that the molecular foundations of complex eusociality may have evolved rapidly in less than 20 Ma.

DEPP: Deep Learning Enables Extending Species Trees using Single Genes

Placing new sequences onto reference phylogenies is increasingly used for analyzing environmental samples, especially microbiomes. Existing placement methods assume that query sequences have evolved under specific models directly on the reference phylogeny. For example, they assume single-gene data (e.g., 16S rRNA amplicons) have evolved under the GTR model on a gene tree. Placement, however, often has a more ambitious goal: extending a (genome-wide) species tree given data from individual genes without knowing the evolutionary model. Addressing this challenging problem requires new directions. Here, we introduce Deep-learning Enabled Phylogenetic Placement (DEPP), an algorithm that learns to extend species trees using single genes without pre-specified models. In simulations and on real data, we show that DEPP can match the accuracy of model-based methods without any prior knowledge of the model. We also show that DEPP can update the multi-locus microbial tree-of-life with single genes with high accuracy. We further demonstrate that DEPP can combine 16S and metagenomic data onto a single tree, enabling community structure analyses that take advantage of both sources of data.

Chromosome-Level Genome Assembly of Anthidium xuezhongi Niu & Zhu, 2020 (Hymenoptera: Apoidea: Megachilidae: Anthidiini)

Anthidiini, a large bee tribe characterized by light-colored maculations, represents nearly 1,000 pollinator species, but no genomes are yet available for this tribe. Here, we report a chromosome-level genome assembly of Anthidium xuezhongi collected from the Tibetan Plateau. Using PacBio long reads and Hi-C data, we assembled a genome of 189.14 Mb with 99.94% of the assembly located in 16 chromosomes. Our assembly contains 23 scaffolds, with the scaffold N50 length of 12.53 Mb, and BUSCO completeness of 98.70% (n = 1,367). We masked 25.98 Mb (13.74%) of the assembly as repetitive elements, identified 385 noncoding RNAs, and predicted 10,820 protein-coding genes (99.20% BUSCO completeness). Gene family evolution analyses identified 9,251 gene families, of which 31 gene families experienced rapid evolution. Interspecific chromosomal variation among A. xuezhongi, Bombus terrestris, and Apis mellifera showed strong chromosomal syntenic relationships. This high-quality genome assembly is a valuable resource for evolutionary and comparative genomic analyses of bees.

Phylogenetic Relationships Within the Hyper-Diverse Genus Eugenia (Myrtaceae: Myrteae) Based on Target Enrichment Sequencing

Eugenia is one of the most taxonomically challenging lineages of flowering plants, in which morphological delimitation has changed over the last few years resulting from recent phylogenetic study based on molecular data. Efforts, until now, have been limited to Sanger sequencing of mostly plastid markers. These phylogenetic studies indicate 11 clades formalized as infrageneric groups. However, relationships among these clades are poorly supported at key nodes and inconsistent between studies, particularly along the backbone and within Eugenia sect. Umbellatae encompasses ca. 700 species. To resolve and better understand systematic discordance, 54 Eugenia taxa were subjected to phylogenomic Hyb-Seq using 353 low-copy nuclear genes. Twenty species trees based on coding and non-coding loci of nuclear and plastid datasets were recovered using coalescent and concatenated approaches. Concordant and conflicting topologies were assessed by comparing tree landscapes, topology tests, and gene and site concordance factors. The topologies are similar except between nuclear and plastid datasets. The coalescent trees better accommodate disparity in the intron dataset, which contains more parsimony informative sites, while concatenated trees recover more conservative topologies, as they have narrower distribution in the tree landscape. This suggests that highly supported phylogenetic relationships determined in previous studies do not necessarily indicate overwhelming concordant signal. Congruence must be interpreted carefully especially in concatenated datasets. Despite this, the congruence between the multi-species coalescent (MSC) approach and concatenated tree topologies found here is notable. Our analysis does not support Eugenia subg. Pseudeugenia or sect. Pilothecium, as currently circumscribed, suggesting necessary taxonomic reassessment. Five clades are further discussed within Eugenia sect. Umbellatae progress toward its division into workable clades. While targeted sequencing provides a massive quantity of data that improves phylogenetic resolution in Eugenia, uncertainty still remains in Eugenia sect. Umbellatae. The general pattern of higher site coefficient factor (CF) than gene CF in the backbone of Eugenia suggests stochastic error from limited signal. Tree landscapes in combination with concordance factor scores, as implemented here, provide a comprehensive approach that incorporates several phylogenetic hypotheses. We believe the protocols employed here will be of use for future investigations on the evolutionary history of Myrtaceae.

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.

Phylogenetic Analysis That Models Compositional Heterogeneity over the Tree

Molecular sequences in a phylogenetic analysis can differ in composition, and that shows that the process of evolution can change over time. However, models of evolution in common use are homogeneous over the tree, and if used in a phylogenetic analysis with compositionally tree-heterogeneous datasets these models can recover incorrect trees. The NDCH or Node-Discrete Compositional Heterogeneity model is able to model such data by accommodating differences in composition over the tree. Usage, problems, and limitations of this model are discussed, and a modification, the NDCH2 model, is described that can ameliorate some of these problems and limitations. Using these models can greatly increase the fit of the model to the data and can find better tree topologies. These models and various statistical tests are illustrated using a bacterial SSU rRNA dataset. These models are implemented in the software P4, and files for the analyses described here are made available.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

A unified evolutionary origin for the ubiquitous protein transporters SecY and YidC

BACKGROUND

Protein transporters translocate hydrophilic segments of polypeptide across hydrophobic cell membranes. Two protein transporters are ubiquitous and date back to the last universal common ancestor: SecY and YidC. SecY consists of two pseudosymmetric halves, which together form a membrane-spanning protein-conducting channel. YidC is an asymmetric molecule with a protein-conducting hydrophilic groove that partially spans the membrane. Although both transporters mediate insertion of membrane proteins with short translocated domains, only SecY transports secretory proteins and membrane proteins with long translocated domains. The evolutionary origins of these ancient and essential transporters are not known.

RESULTS

The features conserved by the two halves of SecY indicate that their common ancestor was an antiparallel homodimeric channel. Structural searches with SecY's halves detect exceptional similarity with YidC homologs. The SecY halves and YidC share a fold comprising a three-helix bundle interrupted by a helical hairpin. In YidC, this hairpin is cytoplasmic and facilitates substrate delivery, whereas in SecY, it is transmembrane and forms the substrate-binding lateral gate helices. In both transporters, the three-helix bundle forms a protein-conducting hydrophilic groove delimited by a conserved hydrophobic residue. Based on these similarities, we propose that SecY originated as a YidC homolog which formed a channel by juxtaposing two hydrophilic grooves in an antiparallel homodimer. We find that archaeal YidC and its eukaryotic descendants use this same dimerisation interface to heterodimerise with a conserved partner. YidC's sufficiency for the function of simple cells is suggested by the results of reductive evolution in mitochondria and plastids, which tend to retain SecY only if they require translocation of large hydrophilic domains.

CONCLUSIONS

SecY and YidC share previously unrecognised similarities in sequence, structure, mechanism, and function. Our delineation of a detailed correspondence between these two essential and ancient transporters enables a deeper mechanistic understanding of how each functions. Furthermore, key differences between them help explain how SecY performs its distinctive function in the recognition and translocation of secretory proteins. The unified theory presented here explains the evolution of these features, and thus reconstructs a key step in the origin of cells.

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

Relative model selection of evolutionary substitution models can be sensitive to multiple sequence alignment uncertainty

BACKGROUND

Multiple sequence alignments (MSAs) represent the fundamental unit of data inputted to most comparative sequence analyses. In phylogenetic analyses in particular, errors in MSA construction have the potential to induce further errors in downstream analyses such as phylogenetic reconstruction itself, ancestral state reconstruction, and divergence time estimation. In addition to providing phylogenetic methods with an MSA to analyze, researchers must also specify a suitable evolutionary model for the given analysis. Most commonly, researchers apply relative model selection to select a model from candidate set and then provide both the MSA and the selected model as input to subsequent analyses. While the influence of MSA errors has been explored for most stages of phylogenetics pipelines, the potential effects of MSA uncertainty on the relative model selection procedure itself have not been explored.

RESULTS

We assessed the consistency of relative model selection when presented with multiple perturbed versions of a given MSA. We find that while relative model selection is mostly robust to MSA uncertainty, in a substantial proportion of circumstances, relative model selection identifies distinct best-fitting models from different MSAs created from the same set of sequences. We find that this issue is more pervasive for nucleotide data compared to amino-acid data. However, we also find that it is challenging to predict whether relative model selection will be robust or sensitive to uncertainty in a given MSA.

CONCLUSIONS

We find that that MSA uncertainty can affect virtually all steps of phylogenetic analysis pipelines to a greater extent than has previously been recognized, including relative model selection.

Phylogenomics resolves ambiguous relationships within Aciculata (Errantia, Annelida)

Aciculata (Eunicida + Phyllodocida) is among the largest clades of annelids, comprising almost half of the known diversity of all marine annelids. Despite the group's large size and biological importance, most phylogenomic studies on Annelida to date have had a limited sampling of this clade. The phylogenetic placement of many clades within Phyllodocida in particular has remained poorly understood. To resolve the relationships within Aciculata we conducted a large-scale phylogenomic analysis based on 24 transcriptomes (13 new), chosen to represent many family-ranked taxa that have never been included in a broad phylogenomic study. Our sampling also includes several enigmatic taxa with challenging phylogenetic placement, such as Histriobdella, Struwela, Lacydonia, Pilargis and the holopelagic worms Lopadorrhynchus, Travisiopsis and Tomopteris. Our robust phylogeny allows us to name and place some of these problematic clades and has significant implications on the systematics of the group. Within Eunicida we reinstate the names Eunicoidea and Oenonoidea. Within Phyllodocida we delineate Phyllodociformia, Glyceriformia, Nereidiformia, Nephtyiformia and Aphroditiformia. Phyllodociformia now includes: Lacydonia, Typhloscolecidae, Lopadorrhynchidae and Phyllodocidae. Nephtyiformia includes Nephtyidae and Pilargidae. We also broaden the delineation of Glyceriformia to include Sphaerodoridae, Tomopteridae and Glyceroidea (Glyceridae + Goniadidae). Furthermore, our study demonstrates and explores how conflicting, yet highly supported topologies can result from confounding signals in gene trees.

Phylogeny of gracillariid leaf-mining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data

Gracillariidae is the most taxonomically diverse cosmopolitan leaf-mining moth family, consisting of nearly 2000 named species in 105 described genera, classified into eight extant subfamilies. The majority of gracillariid species are internal plant feeders as larvae, creating mines and galls in plant tissue. Despite their diversity and ecological adaptations, their phylogenetic relationships, especially among subfamilies, remain uncertain. Genomic data (83 taxa, 589 loci) were integrated with Sanger data (130 taxa, 22 loci), to reconstruct a phylogeny of Gracillariidae. Based on analyses of both datasets combined and analyzed separately, monophyly of Gracillariidae and all its subfamilies, monophyly of the clade "LAMPO" (subfamilies: Lithocolletinae, Acrocercopinae, Marmarinae, Phyllocnistinae, and Oecophyllembiinae) and relationships of its subclade "AMO" (subfamilies: Acrocercopinae, Marmarinae, and Oecophyllembiinae) were strongly supported. A sister-group relationship of Ornixolinae to the remainder of the family, and a monophyletic leaf roller lineage (Callicercops Vári + Parornichinae) + Gracillariinae, as sister to the "LAMPO" clade were supported by the most likely tree. Dating analyses indicate a mid-Cretaceous (105.3 Ma) origin of the family, followed by a rapid diversification into the nine subfamilies predating the Cretaceous-Palaeogene extinction. We hypothesize that advanced larval behaviours, such as making keeled or tentiform blotch mines, rolling leaves and galling, allowed gracillariids to better avoid larval parasitoids allowing them to further diversify. Finally, we stabilize the classification by formally re-establishing the subfamily ranks of Marmarinae stat.rev., Oecophyllembiinae stat.rev. and Parornichinae stat.rev., and erect a new subfamily, Callicercopinae Li, Ohshima and Kawahara to accommodate the enigmatic genus Callicercops.

Theoretical and practical considerations when using retroelement insertions to estimate species trees in the anomaly zone

A potential shortcoming of concatenation methods for species tree estimation is their failure to account for incomplete lineage sorting. Coalescent methods address this problem but make various assumptions that, if violated, can result in worse performance than concatenation. Given the challenges of analyzing DNA sequences with both concatenation and coalescent methods, retroelement insertions (RIs) have emerged as powerful phylogenomic markers for species tree estimation. Here, we show that two recently proposed quartet-based methods, SDPquartets and ASTRAL_BP, are statistically consistent estimators of the unrooted species tree topology under the coalescent when RIs follow a neutral infinite-sites model of mutation and the expected number of new RIs per generation is constant across the species tree. The accuracy of these (and other) methods for inferring species trees from RIs has yet to be assessed on simulated data sets, where the true species tree topology is known. Therefore, we evaluated eight methods given RIs simulated from four model species trees, all of which have short branches and at least three of which are in the anomaly zone. In our simulation study, ASTRAL_BP and SDPquartets always recovered the correct species tree topology when given a sufficiently large number of RIs, as predicted. A distance-based method (ASTRID_BP) and Dollo parsimony also performed well in recovering the species tree topology. In contrast, unordered, polymorphism, and Camin-Sokal parsimony typically fail to recover the correct species tree topology in anomaly zone situations with more than four ingroup taxa. Of the methods studied, only ASTRAL_BP automatically estimates internal branch lengths (in coalescent units) and support values (i.e. local posterior probabilities). We examined the accuracy of branch length estimation, finding that estimated lengths were accurate for short branches but upwardly biased otherwise. This led us to derive the maximum likelihood (branch length) estimate for when RIs are given as input instead of binary gene trees; this corrected formula produced accurate estimates of branch lengths in our simulation study, provided that a sufficiently large number of RIs were given as input. Lastly, we evaluated the impact of data quantity on species tree estimation by repeating the above experiments with input sizes varying from 100 to 100 000 parsimony-informative RIs. We found that, when given just 1 000 parsimony-informative RIs as input, ASTRAL_BP successfully reconstructed major clades (i.e clades separated by branches > 0.3 CUs) with high support and identified rapid radiations (i.e. shorter connected branches), although not their precise branching order. The local posterior probability was effective for controlling false positive branches in these scenarios.

Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae)

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].