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.

Phylogenomics resolves a 100-year-old debate regarding the evolutionary history of caddisflies (Insecta: Trichoptera)

Trichoptera (caddisfly) phylogeny provides an interesting example of aquatic insect evolution, with rich ecological diversification, especially for underwater architecture. Trichoptera provide numerous critical ecosystem services and are also one of the most important groups of aquatic insects for assessing water quality. The phylogenetic relationships of Trichoptera have been debated for nearly a century. In particular, the phylogenetic position of the "cocoon-makers" within Trichoptera has long been contested. Here, we designed a universal single-copy orthologue and sets of ultraconserved element markers specific for Trichoptera for the first time. Simultaneously, we reconstructed the phylogenetic relationship of Trichoptera based on genome data from 111 species, representing 29 families and 71 genera. Our phylogenetic inference clarifies the probable phylogenetic relationships of "cocoon-makers" within Integripalpia. Hydroptilidae is considered as the basal lineage within Integripalpia, and the families Glossosomatidae, Hydrobiosidae, and Rhyacophilidae formed a monophyletic clade, sister to the integripalpian subterorder Phryganides. The resulting divergence time and ancestral state reconstruction suggest that the most recent common ancestor of Trichoptera appeared in the early Permian and that diversification was strongly correlated with habitat change.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Movement of the A-strain maize streak virus in and out of Madagascar

The maize streak virus belongs in the genus Mastrevirus, in the family Geminiviridae. The A-strain of the virus (MSV-A) is recognised as the principal causative agent of the most severe manifestation of maize streak disease (MSD). This disease continues to be a persistent limitation on maize output across sub-Saharan Africa and the nearby Indian Ocean islands. Irrespective of the causes behind the spread of MSV-A, we can determine the paths and speeds with which MSV-A spreads by analysing MSV genome sequence data along with information on when and where samples were taken. This information is valuable for identifying the geographical origins of viral strains that cause sporadic MSD epidemics in specific places and the geographical regions where viruses remain in reservoirs and contribute to prolonged epidemics during outbreaks. Our aim is to utilise these analyses to estimate the timing and origin of the MSV-A that arrived on the island of Madagascar in the Indian Ocean. Specifically, we employ model-based phylogeographic analyses on 524 complete MSV-A genome sequences, which consist of 56 newly obtained genomes from infected maize plants collected in Madagascar. These studies allow us to reconstruct the most likely paths of MSV-A to Madagascar. We found strong evidence for the existence of at least four separate movements of MSV-A variants from East and southern Africa to Madagascar. These movements took place between roughly 1979 (with a 95% highest probability density interval [HPD] ranging from 1976 to 1982) and 2003 (with a 95% HPD ranging from 2002 to 2003). While we inferred that MSV-A variants are spreading at an average rate of 38.9 km/year (with a 95% highest posterior density interval of 34.0-44.4) across their geographical range. Since their arrival in Madagascar, MSV-A variants have been migrating at an average rate of 47.6 km/year (with a 95% highest posterior density interval of 36.05-61.70). Human influences are likely significant contributors to both sporadic long-range movements of MSV-A between mainland Africa and Madagascar, as well as shorter to medium range movements within the island.

A molecular phylogeny for all 21 families within Chiroptera (bats)

Bats, members of the Chiroptera order, rank as the second most diverse group among mammals. Recent molecular systematic studies on bats have successfully classified 21 families within two suborders: Yinpterochiroptera and Yangochiroptera. Nevertheless, the phylogeny within these 21 families has remained a subject of controversy. In this study, we have employed a balanced approach to establish a robust family-level phylogenetic hypothesis for bats, utilizing a more comprehensive molecular dataset. This dataset includes representative species from all 21 bat families, resulting in a reduced level of missing genetic information. The resulting phylogenetic tree comprises 21 lineages that are strongly supported, each corresponding to one of the bat families. Our findings support to place the Emballonuroidea superfamily as the basal lineage of Yangochiroptera, and that Myzopodidae should be situated as a basal lineage of Emballonuroidea, forming a sister relationship with the clade consisting of Nycteridae and Emballonuridae. Finally, we have conducted dating analyses on this newly resolved phylogenetic tree, providing divergence times for each bat family. Collectively, our study has employed a relatively comprehensive molecular dataset to establish a more robust phylogeny encompassing all 21 bat families. This improved phylogenetic framework will significantly contribute to our understanding of evolutionary processes, ecological roles, disease dynamics, and biodiversity conservation in the realm of bats.

Proposal of Lactobacillus amylovorus subsp. animalis subsp. nov. and an emended description of Lactobacillus amylovorus

Seven novel lactic acid bacterial strains (BF125T, BF186, TKL145, YK3, YK6, YK10 and NSK) were isolated from the fresh faeces of Japanese black beef cattle and weanling piglets, spent mushroom substrates, or steeping water of a corn starch production plant. These strains are rod-shaped, Gram-stain-positive, non-motile, non-spore-forming, catalase-negative, cytochrome oxidase-negative, facultatively anaerobic, and homofermentative. Strain BF125T did not produce any gas from glucose; both d- and l-lactate were produced as end-products of glucose (D/L, 40 : 60). Growth occurred at 30-45 °C (optimum, 37 °C), pH 5.0-8.0 (optimum, pH 6.0), and with NaCl concentration of 1.0-3.0% (w/v). The G+C content of genomic DNA of strain BF125T was 37.8 mol% (whole-genome analysis). The major fatty acids were C16 : 0, C18 : 1 ω9c, C19 cyclopropane 9, 10, and summed feature 10. The 16S rRNA gene in strain BF125T showed high similarity to that of the type strain of Lactobacillus amylovorus (99.93%), and the other isolates were also identified as L. amylovorus based on these similarities. A phylogenetic tree based on the core genomes of L. amylovorus strains (n=54), including the seven isolates, showed that they could be divided into two clusters. Strains YK3, YK6, YK10, and NSK were in the first cluster, along with the type strain DSM 20531T, while the second cluster included isolates BF125T, BF186, TKL145, and other strains isolated from various animal origins. Phenotypic differences in fermentability were observed for lactose, salicin, and gentiobiose between these two groups. The intergroup digital DNA-DNA hybridization values (72.9-78.6%) and intergroup average nucleotide identity values (95.64-96.92%) were comparable to values calculated using datasets of other valid subspecies of the genus (ex-) Lactobacillus. In light of the physiological, genotypic, and phylogenetic evidence, we propose a novel subspecies of L. amylovorus, named Lactobacillus amylovorus subsp. animalis subsp. nov. (type strain BF125T=MAFF 212522T=DSM 115528T). Our findings also led to the automatic creation of Lactobacillus amylovorus subsp. amylovorus subsp. nov. and an emended description of the species L. amylovorus.

Antennal Transcriptome Screening and Identification of Chemosensory Proteins in the Double-Spine European Spruce Bark Beetle, Ips duplicatus (Coleoptera: Scolytinae)

The northern bark beetle, Ips duplicatus, is an emerging economic pest, reportedly infesting various species of spruce (Picea spp.), pine (Pinus spp.), and larch (Larix spp.) in Central Europe. Recent climate changes and inconsistent forest management practices have led to the rapid spread of this species, leaving the current monitoring strategies inefficient. As understanding the molecular components of pheromone detection is key to developing novel control strategies, we generated antennal transcriptomes from males and females of this species and annotated the chemosensory proteins. We identified putative candidates for 69 odorant receptors (ORs), 50 ionotropic receptors (IRs), 25 gustatory receptors (GRs), 27 odorant-binding proteins (OBPs), including a tetramer-OBP, 9 chemosensory proteins (CSPs), and 6 sensory neuron membrane proteins (SNMPs). However, no sex-specific chemosensory genes were detected. The phylogenetic analysis revealed conserved orthology in bark beetle chemosensory proteins, especially with a major forest pest and co-habitant, Ips typographus. Recent large-scale functional studies in I. typographus chemoreceptors add greater significance to the orthologous sequences reported here. Nevertheless, identifying chemosensory genes in I. duplicatus is valuable to understanding the chemosensory system and its evolution in bark beetles (Coleoptera) and, generally, insects.

Pinibarreniales, a new order of Sordariomycetes from pine barrens ecosystem

Pinibarrenia chlamydospora, sp. nov. isolated from the roots of highbush blueberry in the New Jersey Pine Barrens, is described and illustrated. Based on multigene phylogenetic analysis, as well as morphological and ecological characteristics, Pinibarreniales and Pinibarreniaceae are established to accommodate this novel lineage in Sordariomycetidae, Sordariomycetes. Pinibarreniales, Tracyllalales, and Vermiculariopsiellales are proposed to be included in the subclass Sordariomycetidae. Pinibarreniales likely have a wide distribution and forms association with Ericaceae plants that live in acidic and oligotrophic environments because its DNA barcode matches with environmental sequences from other independent ecological studies. The plant-fungal interaction experiment revealed negative impacts on Arabidopsis, indicating its pathogenicity. This uncovered new fungal lineage will contribute to a better understanding of the diversity and systematics of Sordariomycetes.

The complete chloroplast genome of Lindernia crustacea (L.) F. Muell 1882 (Linderniaceae) and its phylogenetic analysis

Lindernia crustacea (L.) F. Muell 1882, a species in the Linderniaceae family, holds traditional medicinal value in China. To investigate its genetic diversity, we assembled, annotated, and characterized the first complete chloroplast genome of L. crustacea using Illumina sequencing data and various bioinformatics tools. The genome is 153,647 bp in length, with a GC content of 37.6%. It exhibits a typical quadripartite structure, consisting of a large single-copy region (LSC) of 85,411 bp, a small single-copy region (SSC) of 18,724 bp, and two inverted repeat sequences (IRa and IRb) of 25,816 bp each. The genome was predicted to contain 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis indicated that L. crustacea is closely related to L. stricta. These findings provide a foundation for further research on the evolution and potential medicinal applications of the Linderniaceae family.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Vicariant speciation resulting from biogeographic barriers in the Australian tropics: The case of the red-cheeked dunnart (Sminthopsis virginiae)

Global biodiversity loss continues unabated, and in Australia, the rate of recent mammal extinctions is among the worst in the world. Meanwhile, the diversity among and within many endemic mammal species remains undescribed. This information is crucial to delineate species boundaries and thus inform decision-making for conservation. Sminthopsis virginiae (the red-cheeked dunnart) is a small, dasyurid marsupial found in four disjunct populations around the northern coast of Australia and New Guinea. There are three currently recognized subspecies, each occupying a distinct geographic location. Sminthopsis v. virginiae occurs in Queensland, S. v. rufigenis is distributed across New Guinea and the Aru Islands, and S. v. nitela has populations in the Top End of the Northern Territory and the Kimberley region of Western Australia. Previous molecular work has suggested the current subspecies definitions are not aligned with DNA sequence data, though the sampling was limited. We undertook a comprehensive genetic and morphological review of S. virginiae to clarify relationships within the species. This included mitochondrial (CR, 12S, and cytb) and nuclear (omega-globin, IRBP, and bfib7) loci, and morphometric analysis of skulls and whole wet-preserved specimens held in museums. Maximum Likelihood and Bayesian phylogenetic analyses resolved samples into two distinct clades, demarcated by the Gulf of Carpentaria in Australia's north. Sminthopsis. v. nitela was consistently separated from S. v. virginiae and S. v. rufigenis, based on the overall body and skull size and craniodental features, while S. v. virginiae and S. v. rufigenis were more difficult to distinguish from each other. Thus, we redescribed S. virginiae, recognizing two species, S. nitela (raised from subspecies) and S. virginiae (now comprising the subspecies S. v. virginiae and S. v. rufigenis). This study highlights the importance of recognizing cryptic mammal fauna to help address the gap in our knowledge about diagnosing diversity during a time of conservation crisis.

Functional Characterization of an Odorant Receptor Expressed in Newly Hatched Larvae of Fall Armyworm Spodoptera frugiperda

In the past decade, Spodoptera frugiperda has emerged as a significant invasive pest globally, posing a serious threat to agriculture due to its broad diet, migratory behavior, and ability to cause extensive plant damage. While extensive research has focused on the olfactory capabilities of adult S. frugiperda, understanding of the olfactory process in larvae remains limited, despite larvae playing a crucial role in crop damage. To address this gap, we identified an odorant receptor (OR), SfruOR40, expressed in the first-instar larvae through phylogenetic analysis. Using quantitative real-time PCR, we compared SfruOR40 expression levels in larvae and adults. We then characterized the function of SfruOR40 against 67 compounds using the Xenopus oocyte expression system and found that SfruOR40 responded to three plant volatiles. Further, behavioral experiments revealed a larval attraction to (-)-trans-Caryophyllene oxide. This study elucidates SfruOR40's role in the olfactory recognition of newly hatched S. frugiperda larvae, expanding our knowledge of such mechanisms in Noctuid moths. Furthermore, it highlights the potential of plant-derived natural products for biological pest control from a behavioral ecology perspective.

Genomic diversity in Paenibacillus polymyxa: unveiling distinct species groups and functional variability

BACKGROUND

Paenibacillus polymyxa is a bacterial species of high interest, as suggested by the increased number of publications on its functions in the past years. Accordingly, the number of described strains and sequenced genomes is also on the rise. While functional diversity of P. polymyxa has been suggested before, the available genomic data is now sufficient for robust comparative genomics analyses.

RESULTS

Using 157 genomes, we found significant disparities among strains currently affiliated to P. polymyxa. Multiple taxonomic groups were identified with conserved predicted functions putatively impacting their respective ecology. As strains of this species have been reported to exhibit considerable potential in agriculture, medicine, and bioremediation, it is preferable to clarify their taxonomic organization to facilitate reliable and durable approval as active ingredients.

CONCLUSIONS

Strains currently affiliated to P. polymyxa can be separated into two major species groups with differential potential in nitrogen fixation, plant interaction, secondary metabolism, and antimicrobial resistance, as inferred from genomic data.

Virus ecology and 7-year temporal dynamics across a permafrost thaw gradient

Soil microorganisms are pivotal in the global carbon cycle, but the viruses that affect them and their impact on ecosystems are less understood. In this study, we explored the diversity, dynamics, and ecology of soil viruses through 379 metagenomes collected annually from 2010 to 2017. These samples spanned the seasonally thawed active layer of a permafrost thaw gradient, which included palsa, bog, and fen habitats. We identified 5051 virus operational taxonomic units (vOTUs), doubling the known viruses for this site. These vOTUs were largely ephemeral within habitats, suggesting a turnover at the vOTU level from year to year. While the diversity varied by thaw stage and depth-related patterns were specific to each habitat, the virus communities did not significantly change over time. The abundance ratios of virus to host at the phylum level did not show consistent trends across the thaw gradient, depth, or time. To assess potential ecosystem impacts, we predicted hosts in silico and found viruses linked to microbial lineages involved in the carbon cycle, such as methanotrophy and methanogenesis. This included the identification of viruses of Candidatus Methanoflorens, a significant global methane contributor. We also detected a variety of potential auxiliary metabolic genes, including 24 carbon-degrading glycoside hydrolases, six of which are uniquely terrestrial. In conclusion, these long-term observations enhance our understanding of soil viruses in the context of climate-relevant processes and provide opportunities to explore their role in terrestrial carbon cycling.

A machine-learning-based alternative to phylogenetic bootstrap

MOTIVATION

Currently used methods for estimating branch support in phylogenetic analyses often rely on the classic Felsenstein's bootstrap, parametric tests, or their approximations. As these branch support scores are widely used in phylogenetic analyses, having accurate, fast, and interpretable scores is of high importance.

RESULTS

Here, we employed a data-driven approach to estimate branch support values with a probabilistic interpretation. To this end, we simulated thousands of realistic phylogenetic trees and the corresponding multiple sequence alignments. Each of the obtained alignments was used to infer the phylogeny using state-of-the-art phylogenetic inference software, which was then compared to the true tree. Using these extensive data, we trained machine-learning algorithms to estimate branch support values for each bipartition within the maximum-likelihood trees obtained by each software. Our results demonstrate that our model provides fast and more accurate probability-based branch support values than commonly used procedures. We demonstrate the applicability of our approach on empirical datasets.

AVAILABILITY AND IMPLEMENTATION

The data supporting this work are available in the Figshare repository at https://doi.org/10.6084/m9.figshare.25050554.v1, and the underlying code is accessible via GitHub at https://github.com/noaeker/bootstrap_repo.

Multisubstrate specificity shaped the complex evolution of the aminotransferase family across the tree of life

Aminotransferases (ATs) are an ancient enzyme family that play central roles in core nitrogen metabolism, essential to all organisms. However, many of the AT enzyme functions remain poorly defined, limiting our fundamental understanding of the nitrogen metabolic networks that exist in different organisms. Here, we traced the deep evolutionary history of the AT family by analyzing AT enzymes from 90 species spanning the tree of life (ToL). We found that each organism has maintained a relatively small and constant number of ATs. Mapping the distribution of ATs across the ToL uncovered that many essential AT reactions are carried out by taxon-specific AT enzymes due to wide-spread nonorthologous gene displacements. This complex evolutionary history explains the difficulty of homology-based AT functional prediction. Biochemical characterization of diverse aromatic ATs further revealed their broad substrate specificity, unlike other core metabolic enzymes that evolved to catalyze specific reactions today. Interestingly, however, we found that these AT enzymes that diverged over billion years share common signatures of multisubstrate specificity by employing different nonconserved active site residues. These findings illustrate that AT family enzymes had leveraged their inherent substrate promiscuity to maintain a small yet distinct set of multifunctional AT enzymes in different taxa. This evolutionary history of versatile ATs likely contributed to the establishment of robust and diverse nitrogen metabolic networks that exist throughout the ToL. The study provides a critical foundation to systematically determine diverse AT functions and underlying nitrogen metabolic networks across the ToL.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

Revision of the 'Acanthephyra purpurea' species complex (Crustacea: Decapoda), with an emphasis on species diversification in the Atlantic

We inventoried all nine species of the 'Acanthephyra purpurea' complex, one of the most abundant and cosmopolitan group of mesopelagic shrimps. We used 119 specimens at hand and genetic data for 124 specimens from GenBank and BOLD. Phylogenetic analysis of four genes (COI, 16S, NaK, and enolase) showed that the 'Acanthephyra purpurea' complex is polyphyletic and encompasses two species groups, 'A. purpurea' (mostly Atlantic) and 'A. smithi' (Indo-West Pacific). The 'A. purpurea' species group consists of two major molecular clades A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa. Molecular data suggest that hitherto accepted species A. acanthitelsonis, A. pelagica, and A. sica should be considered as synonyms. The Atlantic is inhabited by at least two cryptic genetic lineages of A. pelagica and A. quadrispinosa. Morphological analyses of qualitative and quantitative (900 measurements) characters resulted in a tabular key to species and in a finding of four evolutionary traits. Atlantic species showed various scenarios of diversification visible on mitochondrial gene level, nuclear gene level, and morphological level. We recorded and discussed similar phylogeographic trends in diversification and in distribution of genetic lineages within two different clades: A. pelagica and A. kingsleyi - A. purpurea - A. quadrispinosa.

Diversity of Cellular Slime Molds (Dictyostelids) in the Fanjing Mountain Nature Reserve and Geographical Distribution Comparisons with Other Representative Nature Reserves in Different Climate Zones of China

Protected areas are widely considered an essential strategy for biodiversity conservation. Dictyostelids are unique protists known to have important ecological functions in promoting soil and plant health through their top-down regulation of ecosystem processes, such as decomposition, that involve bacterial populations. But the relationship between dictyostelid diversity within protected areas remains poorly understood, especially on a large scale. Herein, we report data on the distribution of dictyostelids, identified with ITS + SSU rRNA molecular and morphology-based taxonomy, from soil samples collected in the Fanjing Mountain protected area of Guizhou Province, Southwest China. We compared the biodiversity data of dictyostelids in Fanjing Mountain with similar data from previously sampled sites in four other protected areas, including Changbai Mountain (CB), Gushan Mountain (GS), Baiyun Mountain (BY), and Qinghai-Tibet Plateau (QT) in China. We identified four species of dictyostelids belonging to three genera (Dictyostelium, Heterostelium, and Polysphondylium) and herein provide information on the taxonomy of these species. Two species (Heterostelium pallidum and Dictyostelium purpureum) are common and widely distributed throughout the world, but one species (Polysphondylium fuscans) was new to China. Our data indicate that there is no distinguishable significant correlation between the dictyostelid species studied and environmental factors. Overall, the similarity index between Baiyun Mountain in Henan Province and Fanjing Mountain in Guizhou Province, located at approximately the same longitude, is the highest, and the Jaccard similarity coefficients (Jaccard index) of family, genus, and species are 100%, 100%, and 12.5%, respectively. From a species perspective, species in the same climate zone are not closely related, but obvious geographical distributions are evident in different climate zones. This preliminary study provided evidence of the ecological adaptation of dictyostelids to different biological niches.

The complete plastome of Acalypha australis (Euphorbiaceae) and its phylogenetic analysis

Acalypha australis L. 1753 is a potherb popular among Asian populations and is a traditional herbal medicine. In the current study, the overall genetic diversity of A. australis still needs to be better. Here, we assembled and characterized the complete plastome of A. australis. The plastome is 168,885 bp in length with a large single-copy (LSC) of 94,576 bp, a small single-copy (SSC) of 19,715 bp, and two copies of inverted repeat region (IRs) of 27,297 bp each. The overall GC content is 34.9%. The plastome contains 127 genes, including 83 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis of the representative species of Euphorbiaceae showed that A. australis and A. hispida formed a monophyletic sister clade. The results of this study will support further research on the evolution and conservation of the Euphorbiaceae species; they will benefit pharmaceutical applications and ornamentation of the medicinal plant A. australis.

Population genomics of the invasive Northern Giant Hornet Vespa mandarinia in North America and across its native range

The northern giant hornet Vespa mandarinia (NGH) is a voracious predator of other insect species, including honey bees. NGH's native range spans subtropical and temperate regions across much of east and southeast Asia and, in 2019, exotic populations of the species were discovered in North America. Despite this broad range and invasive potential, investigation of the population genomic structure of NGH across its native and introduced ranges has thus far been limited to a small number of mitochondrial samples. Here, we present analyses of genomic data from NGH individuals collected across the species' native range and from exotic individuals collected in North America. We provide the first survey of whole-genome population variation for any hornet species, covering this species' native and invasive ranges, and in doing so confirm likely origins in Japan and South Korea for the two introductions. We additionally show that, while this introduced population exhibited strongly elevated levels of inbreeding, these signatures of inbreeding are also present in some long-standing native populations, which may indicate that inbreeding depression alone is insufficient to prevent the persistence of NGH populations. As well as highlighting the importance of ongoing monitoring and eradication efforts to limit the spread of this species outside of its natural range, our data will serve as a foundational database for future genomic studies into introduced hornet populations.

Genomic evidence for the widespread presence of GH45 cellulases among soil invertebrates

Lignocellulose is a major component of vascular plant biomass. Its decomposition is crucial for the terrestrial carbon cycle. Microorganisms are considered primary decomposers, but evidence increases that some invertebrates may also decompose lignocellulose. We investigated the taxonomic distribution and evolutionary origins of GH45 hydrolases, important enzymes for the decomposition of cellulose and hemicellulose, in a collection of soil invertebrate genomes. We found that these genes are common in springtails and oribatid mites. Phylogenetic analysis revealed that cellulase genes were acquired early in the evolutionary history of these groups. Domain architectures and predicted 3D enzyme structures indicate that these cellulases are functional. Patterns of presence and absence of these genes across different lineages prompt further investigation into their evolutionary and ecological benefits. The ubiquity of cellulase genes suggests that soil invertebrates may play a role in lignocellulose decomposition, independently or in synergy with microorganisms. Understanding the ecological and evolutionary implications might be crucial for understanding soil food webs and the carbon cycle.

Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

Genomic comparison of deep-sea hydrothermal genera related to Aeropyrum, Thermodiscus and Caldisphaera, and proposed emended description of the family Acidilobaceae

Deep-sea hydrothermal vents host archaeal and bacterial thermophilic communities, including taxonomically and functionally diverse Thermoproteota. Despite their prevalence in high-temperature submarine communities, Thermoproteota are chronically under-represented in genomic databases and issues have emerged regarding their nomenclature, particularly within the Aeropyrum-Thermodiscus-Caldisphaera. To resolve some of these problems, we identified 47 metagenome-assembled genomes (MAGs) within this clade, from 20 previously published deep-sea hydrothermal vent and submarine volcano metagenomes, and 24 MAGs from public databases. Using phylogenomic analysis, Genome Taxonomy Database Toolkit (GTDB-Tk) taxonomic assessment, 16S rRNA gene phylogeny, average amino acid identity (AAI) and functional gene patterns, we re-evaluated of the taxonomy of the Aeropyrum-Thermodiscus-Caldisphaera. At least nine genus-level clades were identified with two or more MAGs. In accordance with SeqCode requirements and recommendations, we propose names for three novel genera, viz. Tiamatella incendiivivens, Hestiella acidicharens and Calypsonella navitae. A fourth genus was also identified related to Thermodiscus maritimus, for which no available sequenced genome exists. We propose the novel species Thermodiscus eudorianus to describe our high-quality Thermodiscus MAG, which represents the type genome for the genus. All three novel genera and T. eudorianus are likely anaerobic heterotrophs, capable of fermenting protein-rich carbon sources, while some Tiamatella, Calypsonella and T. eudorianus may also reduce polysulfides, thiosulfate, sulfur and/or selenite, and the likely acidophile, Hestiella, may reduce nitrate and/or perchlorate. Based on phylogenomic evidence, we also propose the family Acidilobaceae be amended to include Caldisphaera, Aeropyrum, Thermodiscus and Stetteria and the novel genera described here.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

The efficacy of single mitochondrial genes at reconciling the complete mitogenome phylogeny-a case study on dwarf chameleons

Although genome-scale data generation is becoming more tractable for phylogenetics, there are large quantities of single gene fragment data in public repositories and such data are still being generated. We therefore investigated whether single mitochondrial genes are suitable proxies for phylogenetic reconstruction as compared to the application of full mitogenomes. With near complete taxon sampling for the southern African dwarf chameleons (Bradypodion), we estimated and compared phylogenies for the complete mitogenome with topologies generated from individual mitochondrial genes and various combinations of these genes. Our results show that the topologies produced by single genes (ND2, ND4, ND5, COI, and COIII) were analogous to the complete mitogenome, suggesting that these genes may be reliable markers for generating mitochondrial phylogenies in lieu of generating entire mitogenomes. In contrast, the short fragment of 16S commonly used in herpetological systematics, produced a topology quite dissimilar to the complete mitogenome and its concatenation with ND2 weakened the resolution of ND2. We therefore recommend the avoidance of this 16S fragment in future phylogenetic work.

Co-Occurrence of Borrelia burgdorferi Sensu Lato and Babesia spp. DNA in Ixodes ricinus Ticks Collected from Vegetation and Pets in the City of Poznań, Poland

Here, we described the prevalence of Borrelia burgdorferi s.l. and Babesia species found in mono- and double infections among Ixodes ricinus ticks occurring in urban areas of the city of Poznań, Poland. We tested 1029 host-seeking ticks and 1268 engorged ticks removed from pet animals. Borrelia afzelii and B. garinii prevailed both in ticks from vegetation (3.7% and 3.7%, respectively) and from pets (3.7% and 0.6%, respectively). Babesia canis and Ba. microti were the most prevalent in host-seeking (2.6% and 1.4%, respectively) and feeding ticks (2.8% and 2.2%, respectively). Babesia microti sequences proved to be identical to the human pathogenic Ba. microti genotype "Jena/Germany". Sequences of the rarest piroplasm Ba. venatorum (0.7%) were identical with those isolated from European patients. About 1.0% of tested ticks yielded dual infections; in host-seeking ticks, Ba. canis prevailed in co-infections with B. afzelii and B. garinii, whereas Ba. microti and B. afzelii dominated in double-infected feeding ticks. Dual infections, even with a low prevalence, pose a challenge for differential diagnosis in patients with acute febrile disease after a tick bite. The finding of Ba. canis in both tick groups suggests that I. ricinus could be involved in the circulation of this piroplasm.

Redefining the treponemal history through pre-Columbian genomes from Brazil

The origins of treponemal diseases have long remained unknown, especially considering the sudden onset of the first syphilis epidemic in the late 15th century in Europe and its hypothesized arrival from the Americas with Columbus' expeditions1,2. Recently, ancient DNA evidence has revealed various treponemal infections circulating in early modern Europe and colonial-era Mexico3-6. However, there has been to our knowledge no genomic evidence of treponematosis recovered from either the Americas or the Old World that can be reliably dated to the time before the first trans-Atlantic contacts. Here, we present treponemal genomes from nearly 2,000-year-old human remains from Brazil. We reconstruct four ancient genomes of a prehistoric treponemal pathogen, most closely related to the bejel-causing agent Treponema pallidum endemicum. Contradicting the modern day geographical niche of bejel in the arid regions of the world, the results call into question the previous palaeopathological characterization of treponeme subspecies and showcase their adaptive potential. A high-coverage genome is used to improve molecular clock date estimations, placing the divergence of modern T. pallidum subspecies firmly in pre-Columbian times. Overall, our study demonstrates the opportunities within archaeogenetics to uncover key events in pathogen evolution and emergence, paving the way to new hypotheses on the origin and spread of treponematoses.

Genetic Basis and Evolution of Structural Color Polymorphism in an Australian Songbird

Island organisms often evolve phenotypes divergent from their mainland counterparts, providing a useful system for studying adaptation under differential selection. In the white-winged fairywren (Malurus leucopterus), subspecies on two islands have a black nuptial plumage whereas the subspecies on the Australian mainland has a blue nuptial plumage. The black subspecies have a feather nanostructure that could in principle produce a blue structural color, suggesting a blue ancestor. An earlier study proposed independent evolution of melanism on the islands based on the history of subspecies divergence. However, the genetic basis of melanism and the origin of color differentiation in this group are still unknown. Here, we used whole-genome resequencing to investigate the genetic basis of melanism by comparing the blue and black M. leucopterus subspecies to identify highly divergent genomic regions. We identified a well-known pigmentation gene ASIP and four candidate genes that may contribute to feather nanostructure development. Contrary to the prediction of convergent evolution of island melanism, we detected signatures of a selective sweep in genomic regions containing ASIP and SCUBE2 not in the black subspecies but in the blue subspecies, which possesses many derived SNPs in these regions, suggesting that the mainland subspecies has re-evolved a blue plumage from a black ancestor. This proposed re-evolution was likely driven by a preexisting female preference. Our findings provide new insight into the evolution of plumage coloration in island versus continental populations, and, importantly, we identify candidate genes that likely play roles in the development and evolution of feather structural coloration.

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.

The Hypersaline Soils of the Odiel Saltmarshes Natural Area as a Source for Uncovering a New Taxon: Pseudidiomarina terrestris sp. nov

The hypersaline soils of the Odiel Saltmarshes Natural Area are an extreme environment with high levels of some heavy metals; however, it is a relevant source of prokaryotic diversity that we aim to explore. In this study, six strains related to the halophilic genus Pseudidiomarina were isolated from this habitat. The phylogenetic study based on the 16S rRNA gene sequence and the fingerprinting analysis suggested that they constituted a single new species within the genus Pseudidiomarina. Comparative genomic analysis based on the OGRIs indices and the phylogeny inferred from the core genome were performed considering all the members of the family Idiomarinaceae. Additionally, a completed phenotypic characterization, as well as the fatty acid profile, were also carried out. Due to the characteristics of the habitat, genomic functions related to salinity and high heavy metal concentrations were studied, along with the global metabolism of the six isolates. Last, the ecological distribution of the isolates was studied in different hypersaline environments by genome recruitment. To sum up, the six strains constitute a new species within the genus Pseudidiomarina, for which the name Pseudidiomarina terrestris sp. nov. is proposed. The low abundance in all the studied hypersaline habitats indicates that it belongs to the rare biosphere in these habitats. In silico genome functional analysis suggests the presence of heavy metal transporters and pathways for nitrate reduction and nitrogen assimilation in low availability, among other metabolic traits.

Out of the ancient lake: Multiple riverine colonizations and diversification of the freshwater snails in the genus Semisulcospira around Lake Biwa

Ancient lakes are a hotspot of biodiversity. Freshwater species often experience spectacular species radiation after colonizing lakes from riverine habitats. Therefore, the relationship between the fauna of the ancient lakes and the surrounding riverine system has a special significance in understanding their origin and evolutionary history. The study of ancient lake species often focused on the lake colonization of riverine species. In contrast, far less attention has been placed on the reverse direction: the riverine colonization of the lake species, despite its importance in disentangling their complex evolutionary history. The freshwater snails in the genus Semisulcospira involve endemic groups that radiated in the ancient Lake Biwa. Using genetics and fossil records, we inferred that the ancestors of these lake-endemic Semisulcospira snails historically colonized the riverine habitats at least three times during the Middle Pleistocene. Each colonization resulted in the formation of a new lineage that was genetically and morphologically distinct from other lineages. Further, one of these colonizations was followed by hybridization with a cosmopolitan riverine species, which potentially facilitated the population persistence of the colonizers in the new environment. Despite their complex histories, all these colonizers were currently grouped within a single species, Semisulcospira kurodai, suggesting cryptic diversity in this species. This study highlights the significance of the riverine colonizations of the lake species to fully understand the diversification history of freshwater fauna in and around the ancient lakes.

Sample design in biodiversity studies matters: a fine-scale study of Lawrence's velvet worm, Peripatopsis lawrencei (Onychophora: Peripatopsidae), reveals hidden diversity

A fine-scale phylogenetic and phylogeographic analysis of Peripatopsis lawrencei s.l. was conducted with both mitochondrial and nuclear DNA sequence data, using both external morphology and scanning electron microscopy of taxonomically important characters. A total of 119 sequences were used for the mitochondrial cytochrome c oxidase subunit I (COI ) whereas a single representative specimen from each locality was sequenced for the nuclear 18S rRNA locus. Phylogenetic analyses were conducted on the total COI data set and the combined COI + 18S rRNA data set using a Bayesian analysis and maximum likelihood analyses. For the combined DNA sequence data set, a divergence time estimation was further undertaken in BEAST and specimens placed in a phylogenetic framework including all the described Peripatopsis species from South Africa. In addition, a phylogeographic study was conducted exclusively on P. lawrencei s.s. (clade A) using an analysis of molecular variance and haplotype network. Phylogenetic results indicated that, at the Oubos sample locality, two highly distinct genetic lineages were present (clades A and B), whereas a divergence time estimation suggests a Miocene cladogenesis of the novel Oubos lineage. Marked phylogeographic structure was observed for P. lawrencei s.s. (restricted to clade A) across the distribution range with limited maternal dispersal. Morphologically, the two sympatric lineages at Oubos A and B differed in leg pair number, ventral colour and dorsal scale rank counts, as evident from scanning electron microscopy. Our results support the recognition of a distinct species that occurs in sympatry with P. lawrencei s.s. The new species, P. aereus sp. nov. (clade B) is described and the implication for fine-scale taxonomic studies on saproxylic taxa is discussed. ZooBank: urn:lsid:zoobank.org:pub:AB6E0BDA-7B5F-4FD3-A863-BA7C814E278C.

Cryogenian Origins of Multicellularity in Archaeplastida

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

Geastrum sanglinense, a new species from the Manghe Rhesus Monkey National Nature Reserve, China

A novel species of earthstar from China, Geastrum sanglinense is described. Phylogenetic analyses based on sequences of the nuclear ribosomal DNA internal transcribed spacer (ITS), large subunit nuclear ribosomal RNA (nrLSU), and subunit 6 of ATP synthase (atp6) regions showed that the species belongs to subsect. Epigaea in sect. Myceliostroma. The sequences of the new taxon formed a sister group to G. yanshanense and G. rubellum. This species was mainly characterized by scattered or clustered basidiomata (1.9-2.2 cm in width × 2.3-2.5 cm in height), small to medium-sized saccate exoperidium (1.9-4.3 cm diam. when expanded), smooth endoperidial bodies (1.2-2.7 cm diam.), and globose to subglobose basidiospores (3.7-4.1 μm diam.), surface with short columnar warts. The species can also be distinguished by ITS, nrLSU, and atp6 sequences. The new species was described in detail and can provide a reference for the investigation of macrofungi resources in Shanxi Province, China.

Mating Pheromone (Gamone 1) in Blepharisma: A Glycoprotein Responsible for Species Diversity in Unicellular Organisms (Alveolata, Ciliophora)

The genus Blepharisma (Alveolata, Ciliophora) is a unicellular organism distributed worldwide, even in extreme environments, and comprises numerous species. While usually proliferating through cell division, Blepharisma undergoes sexual reproduction (conjugation) when cells are moderately starved. Conjugation is initiated by mating pheromones (gamone 1 and gamone 2) secreted by complementary mating-type cells. Gamone 1, a glycoprotein, functions in a species-specific manner, while gamone 2, an amino acid derivative, is a common molecule across species. The specific function of gamone 1 suggests the possibility that mutations in gamone 1 might have led to reproductive isolation during the evolutionary process, triggering species diversification. In this study, by comparing the amino acid sequences of gamone 1 homologs from 15 strains (representing five species), we found that mutations resulting in distinct amino acid properties occur across species boundaries and are mainly concentrated at two specific regions within gamone 1. These mutations potentially alter the binding affinity of gamone 1 to its receptors, suggesting their effect in causing reproductive isolation. The interspecies artificial conjugation conducted previously and the molecular phylogenetic tree constructed using the gamone 1 homolog sequences in this study provide insights into the speciation process within the genus Blepharisma.

Phylogenomics of the psychoactive mushroom genus Psilocybe and evolution of the psilocybin biosynthetic gene cluster

Psychoactive mushrooms in the genus Psilocybe have immense cultural value and have been used for centuries in Mesoamerica. Despite the recent surge of interest in these mushrooms due to the psychotherapeutic potential of their natural alkaloid psilocybin, their phylogeny and taxonomy remain substantially incomplete. Moreover, the recent elucidation of the psilocybin biosynthetic gene cluster is known for only five of ~165 species of Psilocybe, four of which belong to only one of two major clades. We set out to improve the phylogeny of Psilocybe using shotgun sequencing of fungarium specimens, from which we obtained 71 metagenomes including from 23 types, and conducting phylogenomic analysis of 2,983 single-copy gene families to generate a fully supported phylogeny. Molecular clock analysis suggests the stem lineage of Psilocybe arose ~67 mya and diversified ~56 mya. We also show that psilocybin biosynthesis first arose in Psilocybe, with 4 to 5 possible horizontal transfers to other mushrooms between 40 and 9 mya. Moreover, predicted orthologs of the psilocybin biosynthetic genes revealed two distinct gene orders within the biosynthetic gene cluster that corresponds to a deep split within the genus, possibly a signature of two independent acquisitions of the cluster within Psilocybe.

Performance of tree-building methods using a morphological dataset and a well-supported Hexapoda phylogeny

Recently, many studies have addressed the performance of phylogenetic tree-building methods (maximum parsimony, maximum likelihood, and Bayesian inference), focusing primarily on simulated data. However, for discrete morphological data, there is no consensus yet on which methods recover the phylogeny with better performance. To address this lack of consensus, we investigate the performance of different methods using an empirical dataset for hexapods as a model. As an empirical test of performance, we applied normalized indices to effectively measure accuracy (normalized Robinson-Foulds metric, nRF) and precision, which are measured via resolution, one minus Colless' consensus fork index (1-CFI). Additionally, to further explore phylogenetic accuracy and support measures, we calculated other statistics, such as the true positive rate (statistical power) and the false positive rate (type I error), and constructed receiver operating characteristic plots to visualize the relationship between these statistics. We applied the normalized indices to the reconstructed trees from the reanalyses of an empirical discrete morphological dataset from extant Hexapoda using a well-supported phylogenomic tree as a reference. Maximum likelihood and Bayesian inference applying the k-state Markov (Mk) model (without or with a discrete gamma distribution) performed better, showing higher precision (resolution). Additionally, our results suggest that most available tree topology tests are reliable estimators of the performance measures applied in this study. Thus, we suggest that likelihood-based methods and tree topology tests should be used more often in phylogenetic tree studies based on discrete morphological characters. Our study provides a fair indication that morphological datasets have robust phylogenetic signal.

Comparative chloroplast genomics of 24 species shed light on the genome evolution and phylogeny of subtribe Coelogyninae (Orchidaceae)

BACKGROUND

The orchids of the subtribe Coelogyninae are among the most morphologically diverse and economically important groups within the subfamily Epidendroideae. Previous molecular studies have revealed that Coelogyninae is an unambiguously monophyletic group. However, intergeneric and infrageneric relationships within Coelogyninae are largely unresolved. There has been long controversy over the classification among the genera within the subtribe.

RESULTS

The complete chloroplast (cp.) genomes of 15 species in the subtribe Coelogyninae were newly sequenced and assembled. Together with nine available cp. genomes in GenBank from representative clades of the subtribe, we compared and elucidated the characteristics of 24 Coelogyninae cp. genomes. The results showed that all cp. genomes shared highly conserved structure and contained 135 genes arranged in the same order, including 89 protein-coding genes, 38 tRNAs, and eight rRNAs. Nevertheless, structural variations in relation to particular genes at the IR/SC boundary regions were identified. The diversification pattern of the cp. genomes showed high consistency with the phylogenetic placement of Coelogyninae. The number of different types of SSRs and long repeats exhibited significant differences in the 24 Coelogyninae cp. genomes, wherein mononucleotide repeats (A/T), and palindromic repeats were the most abundant. Four mutation hotspot regions (ycf1a, ndhF-rp132, psaC-ndhE, and rp132-trnL) were determined, which could serve as effective molecular markers. Selection pressure analysis revealed that three genes (ycf1a, rpoC2 and ycf2 genes) might have experienced apparent positive selection during the evolution. Using the alignments of whole cp. genomes and protein-coding sequences, this study presents a well-resolved phylogenetic framework of Coelogyninae.

CONCLUSION

The inclusion of 55 plastid genome data from a nearly complete generic-level sampling provide a comprehensive view of the phylogenetic relationships among genera and species in subtribe Coelogyninae and illustrate the diverse genetic variation patterns of plastid genomes in this species-rich plant group. The inferred relationships and informally recognized major clades within the subtribe are presented. The genetic markers identified here will facilitate future studies on the genetics and phylogeny of subtribe Coelogyninae.

The first released available genome of the common ice plant ( Mesembryanthemum crystallinum L.) extended the research region on salt tolerance, C 3-CAM photosynthetic conversion, and halophilism

Background

The common ice plant ( Mesembryanthemum crystallinum L.) is an annual herb belonging to the genus Mesembryanthemum of the family Aizoaceae, native to Southern Africa.

Methods

We performed shotgun genome paired-end sequencing using the Illumina platform to determine the genome sequence of the ice plants. We assembled the whole genome sequences using the genome assembler "ALGA" and "Redundans", then released them as available genomic information. Finally, we mainly estimated the potential genomic function by the homology search method.

Results

A draft genome was generated with a total length of 286 Mb corresponding to 79.2% of the estimated genome size (361 Mb), consisting of 49,782 contigs. It encompassed 93.49% of the genes of terrestrial higher plants, 99.5% of the ice plant transcriptome, and 100% of known DNA sequences. In addition, 110.9 Mb (38.8%) of repetitive sequences and untranslated regions, 971 tRNA, and 100 miRNA loci were identified, and their effects on stress tolerance and photosynthesis were investigated. Molecular phylogenetic analysis based on ribosomal DNA among 26 kinds of plant species revealed genetic similarity between the ice plant and poplar, which have salt tolerance. Overall, 35,702 protein-coding regions were identified in the genome, of which 56.05% to 82.59% were annotated and submitted to domain searches and gene ontology (GO) analyses, which found that eighteen GO terms stood out among five plant species. These terms were related to biological defense, growth, reproduction, transcription, post-transcription, and intermembrane transportation, regarded as one of the fundamental results of using the utilized ice plant genome.

Conclusions

The information that we characterized is useful for elucidation of the mechanism of growth promotion under salinity and reversible conversion of the photosynthetic type from C3 to Crassulacean Acid Metabolism (CAM).

Performance of Topology Tests under Extreme Selection Bias

Tree tests like the Kishino-Hasegawa (KH) test and chi-square test suffer a selection bias that tests like the Shimodaira-Hasegawa (SH) test and approximately unbiased test were intended to correct. We investigate tree-testing performance in the presence of severe selection bias. The SH test is found to be very conservative and, surprisingly, its uncorrected analog, the KH test has low Type I error even in the presence of extreme selection bias, leading to a recommendation that the SH test be abandoned. A chi-square test is found to usually behave well and but to require correction in extreme cases. We show how topology testing procedures can be used to get support values for splits and compare the likelihood-based support values to the approximate likelihood ratio test (aLRT) support values. We find that the aLRT support values are reasonable even in settings with severe selection bias that they were not designed for. We also show how they can be used to construct tests of topologies and, in doing so, point out a multiple comparisons issue that should be considered when looking at support values for splits.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus Complex)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

Refining PRRSV-2 genetic classification based on global ORF5 sequences and investigation of their geographic distributions and temporal changes

IMPORTANCE

In this study, comprehensive analysis of 82,237 global porcine reproductive and respiratory syndrome virus type 2 (PRRSV-2) open reading frame 5 sequences spanning from 1989 to 2021 refined PRRSV-2 genetic classification system, which defines 11 lineages and 21 sublineages and provides flexibility for growth if additional lineages, sublineages, or more granular classifications are needed in the future. Geographic distribution and temporal changes of PRRSV-2 were investigated in detail. This is a thorough study describing the molecular epidemiology of global PRRSV-2. In addition, the reference sequences based on the refined genetic classification system are made available to the public for future epidemiological and diagnostic applications worldwide. The data from this study will facilitate global standardization and application of PRRSV-2 genetic classification.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae)

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

Phylogenetic and taxonomic updates of Agaricales, with an emphasis on Tricholomopsis

The order Agaricales was divided into eight suborders. However, the phylogenetic relationships among some suborders are largely unresolved, and the phylogenetic positions and delimitations of some taxa, such as Sarcomyxaceae and Tricholomopsis, remain unsettled. In this study, sequence data of 38 genomes were generated through genome skimming on an Illumina sequencing system. To anchor the systematic position of Sarcomyxaceae and Tricholomopsis, a phylogenetic analysis based on 555 single-copy orthologous genes from the aforementioned genomes and 126 publicly accessible genomes was performed. The results fully supported the clustering of Tricholomopsis with Phyllotopsis and Pleurocybella within Phyllotopsidaceae, which formed a divergent monophyletic major lineage together with Pterulaceae, Radulomycetaceae, and Macrotyphula in Agaricales. The analysis also revealed that Sarcomyxaceae formed a unique major clade. Therefore, two new suborders, Phyllotopsidineae and Sarcomyxineae, are proposed for the two major lineages. Analyses of 450 single-copy orthologous genes and four loci suggested that Tricholomopsis consisted of at least four clades. Tricholomopsis is subsequently subdivided into four distinct sections. Seventeen Tricholomopsis species in China, including six new species, are reported. Conoloma is established to accommodate T. mucronata. The substrate preference of Tricholomopsis species and the transitions of the pileate ornamentations among the species within the genus are discussed.