Maximum likelihood inference of protein phylogeny and the origin of chloroplasts

Mixed, not stirred: Genomic data confirm the first case of interspecific hybridization in planarian triclads (Platyhelminthes: Tricladida) and raise questions about a possibly novel form of hybrid speciation

Speciation is a complex process where many evolutionary forces interplay. The Mediterranean is acknowledged as one of the most relevant biodiverse areas in the Palearctic region and researchers have long studied the species inhabiting it to pursue the goals of evolutionary biology. Here, we study a complex of freshwater flatworm species of the genus Dugesia from Corsica and Sardinia using restriction site-associated DNA sequencing (specifically, 3RAD) data to unravel their evolutionary history and tackle the processes driving it. We assess the phylogenetic relationships and population structure within the group and evaluate new species boundaries using multispecies coalescent approaches. Furthermore, we offer insights into the environmental niche model of the group and use said model to guide our sampling efforts and collect and present molecular evidence for the first time of Dugesia leporii specimens, endemic from Sardinia last spotted in 1999. Our results indicate that paleoclimatic conditions rather than microplate tectonic dynamics were likely an important driver of diversification for the Corso-Sardinian group. Furthermore, our results warrant the taxonomic re-evaluation of the group as eight primary species candidates are established based on molecular data. Our study also reveals the first case of interspecific natural hybridization reported in Dugesiidae and, to our knowledge, in Tricladida. Finally, we discuss how this hybridization might constitute a new form of hybrid speciation.

The Early Evolution of Tudor Genes in Holozoa and How Their Distribution Was Influenced by Life History Traits in Metazoa

Early metazoan evolution was characterized by the expansion of multiple gene families, such as the Tudor family, involved in novel multicellularity-related functions. In eukaryotes, Tudor genes (i.e. genes including at least one Tudor domain) are numerous, heterogeneous, and mostly associated with gene expression regulation. However, they underwent an animal-specific expansion, with novel elements almost exclusively involved in retrotransposon regulation through Piwi-interacting RNAs, as spatiotemporal regulators of the key-element Piwi, another previously considered animal-specific gene. Here, we used online-available proteomes covering 25 major taxonomic groups to characterize the Tudor gene family at a holozoan-wide level, confirming the apomorphic metazoan expansion of Piwi-interacting RNA-related Tudor genes. However, we also annotated elements of the Piwi-interacting RNA pathway (Tudor and Piwi genes) in Ichthyosporea species, suggesting that elements of the Piwi-interacting RNA pathway were already present in the holozoan common ancestors. We observed an outstanding variability (34-fold) of Tudor gene number between and within metazoan phyla that could be associated with convergent genomic and phenotypic evolutions: expansions were usually sided by whole-genome duplications and/or life history traits such as parthenogenesis; reductions were mostly associated to overall phenotypic and genomic simplifications, like in almost all considered endoparasites. Lastly, we phylogenetically tested, and mostly (but not completely) confirmed, a previously proposed model for the evolution of the Tudor domain secondary structures.

Rediscovery and morpho-molecular characterization of three astome ciliates, with new insights into eco-evolutionary associations of astomes with their annelid hosts

Astome ciliates live in the digestive tract of a broad spectrum of marine, freshwater, and terricolous annelids. In aquatic lumbriculid and criodrilid oligochaetes collected in Central Europe, we rediscovered three insufficiently known astomes: Hoplitophrya secans, Mesnilella clavata, and Buchneriella criodrili. Their morphology was studied using in vivo observation, protargol, and dry silver nitrate impregnation. Multiple nuclear and mitochondrial molecular markers were used to determine their phylogenetic positions and reconstruct their evolutionary history. According to our phylogenetic analyses: (1) mouthless ciliates isolated from annelids form a robustly supported monophylum within the class Oligohymenophorea, (2) the progenitor of astomes invaded the digestive tract of marine polychaetes during the Paleozoic era, (3) lumbricid earthworms likely served as a source of astomes for criodrilid, almid, and megascolecid earthworms, (4) the ancestral host of the earthworm-dwelling astome clade led an endogeic lifestyle, and (5) there were multiple independent transfers of astomes from endogeic to epigeic and anecic earthworms. These findings support previous views of the annelid phylogeny, suggesting that astomes reside and evolve in tandem with annelids for several hundred million years.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00275-5.

Refining genetic classification of global porcine reproductive and respiratory syndrome virus type 1 (PRRSV-1) and investigating their geographic and temporal distributions

Porcine reproductive and respiratory syndrome virus type 1 (PRRSV-1) primarily circulates in Europe but is also detected in North America and Asia. Based on ORF5 sequences, previous studies classified PRRSV-1 into four subtypes. Subtype 1 was further classified into 12 clades (A-L) or into three lineages with lineage 1 including clades 1A-1G and lineage 3 including clades 3A-3G, but the systems are inconsistent and have not been adopted. In this study, we proposed a statistically supported PRRSV-1 genetic classification system based on 10,446 global PRRSV-1 ORF5 sequences spanning 1991-2023. We replaced the colloquial "subtype" designation with "lineage" to reflect evolutionary history and, subsequently, PRRSV-1 was classified into four lineages (L1-L4) with L1 including 18 sublineages (L1.1 to L1.18). The proposed classification system is flexible and may be amended if additional lineages, sublineages, or more granular classifications are needed to reflect contemporary PRRSV-1 detections and evolution. Geographic distributions of PRRSV-1 at lineage and sublineage levels were distinct, with L1 globally distributed and L2, L3 and L4 more restricted. Temporal dynamic changes in some countries were quantified. Classification and ORF5 nucleotide identity of six commercial PRRSV-1 vaccines to each lineage and sublineage and detection frequency of vaccine-like viruses were determined. The phylogenies based on whole-genome and ORF5 sequences demonstrated slightly different tree topologies. Recombination of PRRSV-1 was observed at within-sublineage and between-sublineage levels. A set of ORF5 reference sequences representing the refined classification is available for future diagnostic and epidemiological applications. This study provides a benchmark delineating the current genetic diversity of PRRSV-1 and introduces a refined classification system to support the global standardization and application of ORF5-based genetic classification for PRRSV-1.

Phylogenomic analyses unravel the tangled evolutionary history of Genisteae (Fabaceae)

Genisteae, a tribe in the subfamily Papilionoideae (Fabaceae), is characterized by the production of quinolizidine alkaloids that confer pest resistance in most of its members. Many relationships at the generic level remain unresolved due largely to a lack of modern attempts to reconstruct the phylogeny. Previous studies with limited taxon sampling and only a few molecular loci indicated the presence of three clades within the tribe: the Lupinus clade, the Cytisus-Genista complex and the Argyrolobium group. There are also two historical genera, Teline and Chamaecytisus, that have been reclassified over the years with some controversy. Species from Teline are currently classified in Genista, and Chamaecytisus species are placed inside of Cytisus. Sellocharis is another genus with vague placement inside of Genisteae near Anarthrophyllum, based mostly on morphology and cytology. Representative taxa from 24 of the 25 genera of Genisteae, along with species of historical genera Teline and Chamaecytisus, were sampled and utilized in a phylogenomic investigation using both plastid and nuclear data to resolve relationships at the generic level. Whole genomic DNA was sequenced and complete plastomes were assembled and annotated. Low-copy nuclear genes were retrieved from the genomic DNA sequences using a mapping-assembly-scaffold approach. Phylogenetic analyses using maximum likelihood, Bayesian and coalescence methods resulted in fully resolved and strongly supported trees for both nuclear and plastid data that show four major clades inside of Genisteae: Cytisus-Genista complex, Lupinus clade, Argyrolobium group and the novel Anarthrophyllum group. The resulting phylogenetic trees also supported the transfer of the Argyrolobium group from Crotalarieae to Genisteae, and the reclassification of Teline inside of Genista, both of which were previously suggested in literature. The phylogenetic trees also determined a placement for Sellocharis in the new Anarthrophyllum group. Although both nuclear and plastid trees were congruent with respect to the monophyly and relationships among the four major clades, incongruence was detected within some of the major clades and the potential causes are discussed.

Novel strategy for whole-genome sequencing of hepatitis A virus using NGS illumina technology and phylogenetic comparison with partial VP1/2A genomic region

Molecular epidemiology of hepatitis A virus (HAV) plays a critical role in identifying outbreak origin and conducting surveillance. Although it is mostly carried out using short partial VP1/2A genomic sequences, utilizing whole-genome sequences (WGS) provides more accurate and robust information. We developed an amplicon-based next-generation sequencing (NGS) strategy to obtain complete HAV genomes utilizing the COVIDSeq Test (Illumina). Twenty-five primer pairs were designed and used to amplify partial genomic fragments (400 bp) that comprise the entire HAV genome sequence from previously HAV positive serum and stool samples from Argentina. The DNA library was prepared using the Illumina COVIDSeq Test and sequenced in a MiSeq equipment. Phylogenetic analyses were performed with IQ-Tree using WGS and VP1/2A partial sequences of 1084pb and 422pb. Eleven samples were amplified and sequenced, with coverage between 79.3 and 100% (> 90% in 9 samples). Although phylogenetic analyses of partial sequences allowed genotype and subgenotype identification, WGS analyses yielded more accurate and reliable results for the phylogeny (phylogenetic definition). The amplicon-based NGS WGS tool developed by adapting the COVIDSeq test to HAV proved to be efficient. The study of partial VP1/2A regions (mainly the 1084 bp fragment) would constitute useful alternatives for outbreak investigation and surveillance when WGS could not be performed.

Diversity and phylogenetic position of the amphi-American lineages of the tapeworms of the genus Anthocephalum Linton, 1890 (Rhinebothriidea: Anthocephaliidae)

Among enduringly associated parasitic taxa, tapeworms (Platyhelminthes: Cestoda), particularly the family Anthocephaliidae Ruhnke, Caira & Cox, 2015, pose challenges to systematic history of the group. Within this family, the genus Anthocephalum Linton, 1890, remains insufficiently explored despite recent advancements. This study addresses the intricate taxonomy, phylogenetic relationships, diversity, and historical biogeography of tapeworms within the Anthocephaliidae, focusing on its type-genus Anthocephalum. To accomplish this objective, 15 specimens across various geographical regions including Alagoas and Pará in Brazil, Panama, and Senegal were selected for DNA extraction. For species of Anthocephalum, partial 18S and 28S rDNA sequences were amplified via PCR, purified, and sequenced using the ABI Big-Dye method. Three specimens of Alveobothrium grabatum Boudaya, Neifar & Euzet, 2018 were sequenced through Illumina technology. Phylogenetic analyses were conducted using IQTREE2 under two optimality criteria: Maximum Parsimony (MP) and Maximum Likelihood (ML). Biogeographic ancestral range estimations were performed using the R package BioGeoBEARS, incorporating multiple biogeographical models. Our phylogenetic analyses reaffirmed Anthocephaliidae's monophyly. However, both the relationships with or within Anthocephaliidae require further investigation. One example within the family is the positioning of Alveobothrium Boudaya, Neifar & Euzet, 2018, which challenged Anthocephalum's monophyly requiring taxonomic actions. Furthermore, when exploring new localities and/or new hosts for Anthocephalum, four independent lineages were identified, suggesting that Anthocephalum's diversity remains underestimated especially in unexplored regions and hosts. About the biogeographic ancestral range estimation, the analysis suggested a Central Indo-Pacific origin for early Anthocephalum lineages, with subsequent colonization events shaping the current diversity. The preliminary biogeographical framework presented here underscores the importance of refining phylogenetic hypotheses and enhancing taxonomic understanding. As taxonomic actions taken, the four new lineages were formally described and the genus Alveobothrium was synonymized with Anthocephalum, for which we proposed an amended diagnosis. This revision brings the total number of valid species of Anthocephalum to 30. Therefore, we suggest that future research initiatives should prioritize time-calibrated analyses, multiple genetic loci, and broader taxonomic representation for a detailed exploration of systematic history of anthocephaliid's.

nT4X and nT4M: Novel Time Non-reversible Mixture Amino Acid Substitution Models

One of the most important and difficult challenges in the research of molecular evolution is modeling the process of amino acid substitutions. Although single-matrix models, such as the LG model, are popular, their capability to properly capture the heterogeneity of the substitution process across sites is still questioned. Several mixture models with multiple matrices have been introduced and shown to offer advantages over single-matrix models. Current general mixture models assume the reversibility of the evolutionary process, implying that substitution rates between any two amino acids are equal in both forward and backward directions. This assumption is not based on biological properties but rather on computational simplicity. The well-known hypothesis is that more realistic models can yield more accurate evolutionary inferences; therefore, our aim is to estimate more biologically realistic models. To this end, we relax the assumption of reversibility and introduce two new general non-reversible 4-matrix mixture models, called nT4M and nT4X. Using alignments from HSSP and TreeBASE databases as data, our newly estimated models outperformed all single-matrix models and almost all reversible mixture models. Moreover, the new non-reversible mixture models enable us to infer rooted trees.

Two waves of photosymbiosis acquisition in extant planktonic foraminifera explained by ecological incumbency

Photosymbiosis, a mode of mixotrophy by algal endosymbiosis, provides key advantages to pelagic life in oligotrophic oceans. Despite its ecological importance, mechanisms underlying its emergence and association with the evolutionary success of photosymbiotic lineages remain unclear. We used planktonic foraminifera, a group of pelagic test-forming protists with an excellent fossil record, to reveal the history of symbiont acquisition among their three main extant clades. We used single-cell 18S rRNA gene amplicon sequencing to reveal symbiont identity and mapped the symbiosis on a phylogeny time-calibrated by fossil data. We show that the highly specific symbiotic interaction with dinoflagellates emerged in the wake of a major extinction of symbiont-bearing taxa at the end of the Eocene. In contrast, less specific and low-light-adapted symbioses with pelagophytes emerged 20 million years later, in multiple independent lineages in the Late Neogene, at a time when the vertical structure of pelagic ecosystems was transformed by global cooling. We infer that in foraminifera, photosymbiosis can evolve easily and that its establishment leads to diversification and ecological dominance to such an extent, that the proliferation of new symbioses is prevented by the incumbent lineages.

Convergent evolution revealed by paraphyly and polyphyly of many taxa of oribatid mites: A molecular approach

A reliable phylogeny is crucial for understanding the evolution and radiation of animal taxa. Phylogenies based on morphological data may be misleading due to frequent convergent evolution of traits-a problem from which molecular phylogenies suffer less. This may be particularly relevant in oribatid mites, an ancient soil animal taxon with more than 11,000 species, where the classification of species into high-ranking taxa such as superfamilies is equivocal. Here, we present a molecular phylogeny of 317 oribatid, 4 astigmatid and 17 endeostigmatid mite species/taxa based on 18S rDNA sequences. We aimed at testing the validity of the 41 superfamilies of oribatid mites recognized by Norton and Behan-Pelletier (in Krantz and Walter, A manual of Acarology, 3rd ed., Texas Technical University Press, Lubbock, 2009). The results indicate that 17 of the 41 oribatid mite superfamilies are monophyletic but that 18 superfamilies are paraphyletic or polyphyletic (5 superfamilies were only included with one species and Microzetoidea were not included). Our findings point to the importance of convergent evolution in polyphyletic oribatid mite taxa. Convergent evolution and the old age of mites likely resulted in a mosaic-like distribution of morphological characters impeding phylogenetic reconstructions based on morphology, calling for molecular approaches to improve oribatid mite systematics.

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.

Evidence of horizontal gene transfer and environmental selection impacting antibiotic resistance evolution in soil-dwelling Listeria

Soil is an important reservoir of antibiotic resistance genes (ARGs) and understanding how corresponding environmental changes influence their emergence, evolution, and spread is crucial. The soil-dwelling bacterial genus Listeria, including L. monocytogenes, the causative agent of listeriosis, serves as a key model for establishing this understanding. Here, we characterize ARGs in 594 genomes representing 19 Listeria species that we previously isolated from soils in natural environments across the United States. Among the five putatively functional ARGs identified, lin, which confers resistance to lincomycin, is the most prevalent, followed by mprF, sul, fosX, and norB. ARGs are predominantly found in Listeria sensu stricto species, with those more closely related to L. monocytogenes tending to harbor more ARGs. Notably, phylogenetic and recombination analyses provide evidence of recent horizontal gene transfer (HGT) in all five ARGs within and/or across species, likely mediated by transformation rather than conjugation and transduction. In addition, the richness and genetic divergence of ARGs are associated with environmental conditions, particularly soil properties (e.g., aluminum and magnesium) and surrounding land use patterns (e.g., forest coverage). Collectively, our data suggest that recent HGT and environmental selection play a vital role in the acquisition and diversification of bacterial ARGs in natural environments.

The Bayesian Phylogenetic Bootstrap and its Application to Short Trees and Branches

Felsenstein's bootstrap is the most commonly used method to measure branch support in phylogenetics. Current sequencing technologies can result in massive sampling of taxa (e.g. SARS-CoV-2). In this case, the sequences are very similar, the trees are short, and the branches correspond to a small number of mutations (possibly 0). Nevertheless, these trees contain a strong signal, with unresolved parts but a low rate of false branches. With such data, Felsenstein's bootstrap is not satisfactory. Due to the frequentist nature of bootstrap sampling, the expected support of a branch corresponding to a single mutation is ∼63%, even though it is highly likely to be correct. Here, we propose a Bayesian version of the phylogenetic bootstrap in which sites are assigned uninformative prior probabilities. The branch support can then be interpreted as a posterior probability. We do not view the alignment as a small subsample of a large sample of sites, but rather as containing all available information (e.g. as with complete viral genomes, which are becoming routine). We give formulas for expected supports under the assumption of perfect phylogeny, in both the frequentist and Bayesian frameworks, where a branch corresponding to a single mutation now has an expected support of ∼90%. Simulations show that these theoretical results are robust to realistic data. Analyses on low-homoplasy viral and nonviral datasets show that Bayesian bootstrap support is easier to interpret, with high supports for branches very likely to be correct. As homoplasy increases, the two supports become closer and strongly correlated.

Wētā Aotearoa-Polyphyly of the New Zealand Anostostomatidae (Insecta: Orthoptera)

The Anostostomatidae of Aotearoa New Zealand are well-characterized at the genus and species level, but the higher-level systematics of the family as a whole remain poorly resolved. We tested the hypothesis that the New Zealand anaostostomatid fauna consists of a single monophyletic group consistent with a single common ancestor. For phylogenetic analysis, we sampled the genera in Aotearoa New Zealand as well as representatives of the family from Australia and New Caledonia. Maximum likelihood analyses including topological comparison statistics with a DNA alignment of thirteen mitochondrial and four nuclear protein coding genes rejected the monophyly of lineages in New Zealand. We found phylogenetic support for four separate New Zealand lineages; three with their closest relatives in Australia and one in New Caledonia. The New Zealand genus Hemiandrus is paraphyletic and the establishment of a morphologically distinct genus is justified. We determined that six of the valid species previously placed in Hemiandrus form a distinct clade that we designated here as Anderus gen. nov. The putative Hemiandrus that we sampled from Australia was sister to neither of the New Zealand lineages.

Evolutionary origin and population diversity of a cryptic hybrid pathogen

Cryptic fungal pathogens pose disease management challenges due to their morphological resemblance to known pathogens. Here, we investigated the genomes and phenotypes of 53 globally distributed isolates of Aspergillus section Nidulantes fungi and found 30 clinical isolates-including four isolated from COVID-19 patients-were A. latus, a cryptic pathogen that originated via allodiploid hybridization. Notably, all A. latus isolates were misidentified. A. latus hybrids likely originated via a single hybridization event during the Miocene and harbor substantial genetic diversity. Transcriptome profiling of a clinical isolate revealed that both parental subgenomes are actively expressed and respond to environmental stimuli. Characterizing infection-relevant traits-such as drug resistance and growth under oxidative stress-revealed distinct phenotypic profiles among A. latus hybrids compared to parental and closely related species. Moreover, we identified four features that could aid A. latus taxonomic identification. Together, these findings deepen our understanding of the origin of cryptic pathogens.

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.

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.

Systematic revision of the South American "Nuncia" (Opiliones, Laniatores, Triaenonychidae)

The genus Nuncia has long been the most speciose within the Opiliones family Triaenonychidae, comprising 63 species and subspecies distributed across New Zealand and South America. Recent molecular studies utilizing Sanger sequencing and ultraconserved elements (UCEs) have indicated that this genus is not monophyletic, and true Nuncia are actually confined to New Zealand. Here, the morphology of all South American triaenonychids is re-examined and DNA sequence data compiled from three markers (18S rRNA, 28S rRNA and cytochrome c oxidase subunit I) for a large number of triaenonychid species, including specimens from all areas with species currently and formerly classified in Nuncia to reassess their phylogenetic position. Based on our findings we 1) revalidate the genus Chilenuncia (Muñoz-Cuevas, 1971) nom. rest.; 2) describe five new genera: Fresiax gen. nov., Mistralia gen. nov., Laftrachia gen. nov., Lautaria gen. nov., Nerudiella gen. nov.; 3) redescribe five species: Fresiaxspinulosa comb. nov., Mistraliaverrucosa comb. nov., Chilenunciachilensis comb. nov., Chilenunciarostrata comb. nov., Nerudiellaamericana comb. nov.; and 4) describe 22 new species of South American triaenonychids: Fresiaxconica sp. nov., Fresiaxfray sp. nov., Fresiaxmauryi sp. nov., Fresiaxpichicuy sp. nov., Mistraliaramirezi sp. nov., Laftrachiarobin sp. nov., Lautariaceachei sp. nov., Nerudiellacachai sp. nov., Nerudiellacaramavida sp. nov., Nerudiellacautin sp. nov., Nerudiellachoapa sp. nov., Nerudiellacuri sp. nov., Nerudiellagoroi sp. nov., Nerudiellajaimei sp. nov., Nerudiellamalleco sp. nov., Nerudiellapenco sp. nov., Nerudiellapichi sp. nov., Nerudiellaportai sp. nov., Nerudiellaquenes sp. nov., Nerudiellavilches sp. nov., Nerudiellawekufe sp. nov., and Nerudiellazapallar sp. nov. Furthermore, we provide detailed illustrations of all the South American species belonging to these lineages formerly classified in Nuncia.

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.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize diverse secondary metabolites, small organic bioactive molecules with key roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce a cadre of secondary metabolites, some of which are useful (e.g. the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and others harmful (e.g. the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. Fungal genomes often also encode resistance genes that confer protection against toxic secondary metabolites. Some Penicillium species, such as Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity. To investigate the evolutionary conservation of homologs of the gliotoxin BGC and of genes involved in gliotoxin resistance in Penicillium, we analyzed 35 Penicillium genomes from 23 species. Homologous, lesser fragmented gliotoxin BGCs were found in 12 genomes, mostly fragmented remnants of the gliotoxin BGC were found in 21 genomes, whereas the remaining 2 Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes was observed. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes, suggestive of constraint and potential functional significance or more recent decay. Gene tree-species tree reconciliation analyses suggested that the history of homologs in the gliotoxin BGC across the genus Penicillium likely involved multiple duplications, losses, and horizontal gene transfers. Our analyses suggest that genes encoded in BGCs can have complex evolutionary histories and be retained in genomes long after the loss of secondary metabolite biosynthesis.

Systematics and evolutionary dynamics of insect-fern interactions in the specialized fern-spore feeding Cuprininae (Lepidoptera, Stathmopodidae)

Fern-spore-feeding (FSF) is rare and found in only four families of Lepidoptera. Stathmopodidae is the most speciose family that contains FSF species, and its subfamily Cuprininae exclusively specializes on FSF. However, three species of Stathmopodinae also specialize on FSF. To better understand the evolutionary history of FSF and, more generally, the significance of specialization on a peculiar host, a phylogenetic and taxonomic revision for this group is necessary. We reconstructed the most comprehensive molecular phylogeny, including one mitochondrial and four nuclear genes, of Stathmopodidae to date, including 137 samples representing 62 species, with a particular focus on the FSF subfamily, Cuprininae, including 33 species (41% of named species) from 6 of the 7 Cuprininae genera. Species from two other subfamilies, Stathmopodinae and Atkinsoniinae, were also included. We found that FSF evolved only once in Stathmopodidae and that the previous hypothesis of multiple origins of FSF was misled by inadequate taxonomy. Moreover, we showed that (1) speciation/extinction rates do not differ significantly between FSF and non-FSF groups and that (2) oligophage is the ancestral character state in Cuprininae. We further revealed that a faster rate of accumulating specialists over time, and thus a higher number of specialists, was achieved by a higher transition rate from oligophagages to specialists compared to the transition rate in the opposite direction. We finish by describing three new genera, Trigonodagen. nov., Petalagen. nov., and Pediformisgen. nov., and revalidating five genera: Cuprina, Calicotis, Thylacosceles, Actinoscelis, Thylacosceloides in Cuprininae, and we provide an updated taxonomic key to genera and a revised global checklist of Cuprininae.

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.

Utilization of systematic error-assessment software to improve phylogenetic accuracy

Unlike classical systems based on the use of morphological data, modern phylogenetic analyses use genetic information to construct phylogenetic trees. Ongoing research in the field of phylogenetics is evaluating the accuracy of phylogenetic estimation results and the reliability of phylogenetic trees to explain evolutionary relationships. Recently, the probability of stochastic errors in large-scale phylogenetic datasets has decreased, while the probability of systematic errors has increased. Therefore, before constructing a phylogenetic tree, it is necessary to assess the causes of systematic bias to improve the accuracy of phylogenetic estimates. We performed analyses of three datasets (Terebelliformia, Daphniid, and Glires clades) using bioinformatics software to assess systematic error and improve phylogenetic tree accuracy. Then, we proposed a combination of systematic biases capable of discerning the most suitable gene markers within a series of taxa and generating conflicting phylogenetic topologies. Our findings will help improve the reliability of phylogenetic software to estimate phylogenies more accurately by exploiting systematic bias.

Cryptic diversity of cellulose-degrading gut bacteria in industrialized humans

Humans, like all mammals, depend on the gut microbiome for digestion of cellulose, the main component of plant fiber. However, evidence for cellulose fermentation in the human gut is scarce. We have identified ruminococcal species in the gut microbiota of human populations that assemble functional multienzymatic cellulosome structures capable of degrading plant cell wall polysaccharides. One of these species, which is strongly associated with humans, likely originated in the ruminant gut and was subsequently transferred to the human gut, potentially during domestication where it underwent diversification and diet-related adaptation through the acquisition of genes from other gut microbes. Collectively, these species are abundant and widespread among ancient humans, hunter-gatherers, and rural populations but are rare in populations from industrialized societies thus indicating potential disappearance in response to the westernized lifestyle.

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.

Comparative genomics reveals distinct diversification patterns among LysR-type transcriptional regulators in the ESKAPE pathogen Pseudomonas aeruginosa

Pseudomonas aeruginosa, a harmful nosocomial pathogen associated with cystic fibrosis and burn wounds, encodes for a large number of LysR-type transcriptional regulator proteins. To understand how and why LTTR proteins evolved with such frequency and to establish whether any relationships exist within the distribution we set out to identify the patterns underpinning LTTR distribution in P. aeruginosa and to uncover cluster-based relationships within the pangenome. Comparative genomic studies revealed that in the JGI IMG database alone ~86 000 LTTRs are present across the sequenced genomes (n=699). They are widely distributed across the species, with core LTTRs present in >93 % of the genomes and accessory LTTRs present in <7 %. Analysis showed that subsets of core LTTRs can be classified as either variable (typically specific to P. aeruginosa) or conserved (and found to be distributed in other Pseudomonas species). Extending the analysis to the more extensive Pseudomonas database, PA14 rooted analysis confirmed the diversification patterns and revealed PqsR, the receptor for the Pseudomonas quinolone signal (PQS) and 2-heptyl-4-quinolone (HHQ) quorum-sensing signals, to be amongst the most variable in the dataset. Successful complementation of the PAO1 pqsR - mutant using representative variant pqsR sequences suggests a degree of structural promiscuity within the most variable of LTTRs, several of which play a prominent role in signalling and communication. These findings provide a new insight into the diversification of LTTR proteins within the P. aeruginosa species and suggests a functional significance to the cluster, conservation and distribution patterns identified.

Many-core algorithms for high-dimensional gradients on phylogenetic trees

MOTIVATION

Advancements in high-throughput genomic sequencing are delivering genomic pathogen data at an unprecedented rate, positioning statistical phylogenetics as a critical tool to monitor infectious diseases globally. This rapid growth spurs the need for efficient inference techniques, such as Hamiltonian Monte Carlo (HMC) in a Bayesian framework, to estimate parameters of these phylogenetic models where the dimensions of the parameters increase with the number of sequences N. HMC requires repeated calculation of the gradient of the data log-likelihood with respect to (wrt) all branch-length-specific (BLS) parameters that traditionally takes O(N2) operations using the standard pruning algorithm. A recent study proposes an approach to calculate this gradient in O(N), enabling researchers to take advantage of gradient-based samplers such as HMC. The CPU implementation of this approach makes the calculation of the gradient computationally tractable for nucleotide-based models but falls short in performance for larger state-space size models, such as Markov-modulated and codon models. Here, we describe novel massively parallel algorithms to calculate the gradient of the log-likelihood wrt all BLS parameters that take advantage of graphics processing units (GPUs) and result in many fold higher speedups over previous CPU implementations.

RESULTS

We benchmark these GPU algorithms on three computing systems using three evolutionary inference examples exploring complete genomes from 997 dengue viruses, 62 carnivore mitochondria and 49 yeasts, and observe a >128-fold speedup over the CPU implementation for codon-based models and >8-fold speedup for nucleotide-based models. As a practical demonstration, we also estimate the timing of the first introduction of West Nile virus into the continental Unites States under a codon model with a relaxed molecular clock from 104 full viral genomes, an inference task previously intractable.

AVAILABILITY AND IMPLEMENTATION

We provide an implementation of our GPU algorithms in BEAGLE v4.0.0 (https://github.com/beagle-dev/beagle-lib), an open-source library for statistical phylogenetics that enables parallel calculations on multi-core CPUs and GPUs. We employ a BEAGLE-implementation using the Bayesian phylogenetics framework BEAST (https://github.com/beast-dev/beast-mcmc).

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.

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 comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Enigmatic Campyloxenus: Shedding light on the delayed origin of bioluminescence in ancient Gondwanan click beetles

Gondwanan elaterids, previously thought to be unrelated, include bioluminescent Campyloxenus earlier placed in bioluminescent Pyrophorinae. Genomic data suggest close relationships between Gondwanan groups. We maintain Morostomatinae and Hapatesinae and redefine Pityobiinae with Nearctic Pityobiini, Gondwanan Parablacini stat. nov., Campyloxenini stat. nov., and Tibionemini trib. nov. Their ancestors putatively underwent differentiation in Gondwana during the Cretaceous separation of southern continents. In contrast with their age, extant groups are species poor. Campyloxenus represents a recent origin of bioluminescence, no older than ∼53 my. Its large pronotal lanterns differ from Pyrophorini and resemble color patches of sympatric beetle co-mimics. This discovery highlights the fourth or fifth origin of bioluminescence in Elateroidea, alongside the lampyroid clade, click beetles Pyrophorini, Alampoides and Coctilelater in Anaissini (Pyrophorinae), and Balgus schnusei (Thylacosterninae). While our phylogenetic findings illuminate the phylogenetic aspects, the complete story awaits further field observations and in-depth genomic analyses of biochemical pathways used by bioluminescent elateroids.

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.

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.

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

The complete chloroplast genome of Diplodiscus trichospermus and phylogenetic position of Brownlowioideae within Malvaceae

BACKGROUND

Malvaceae is an economically important plant family of 4,225 species in nine subfamilies. Phylogenetic relationships among the nine subfamilies have always been controversial, especially for Brownlowioideae, whose phylogenetic position remains largely unknown due to the lack of samples in previous analysis datasets. To greatly clarify the phylogenetic relationship of Malvaceae, we newly sequenced and assembled the plastome of Diplodiscus trichospermus taxonomically located in Brownlowioideae, and downloaded the allied genomes from public database to build a dataset covering all subfamily members of Malvaceae.

RESULTS

The annotation results showed that the plastome of Diplodiscus trichospermus has a typical quadripartite structure, comprising 112 unique genes, namely 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The total length was 158,570 bp with 37.2% GC content. Based on the maximum likelihood method and Bayesian inference, a robust phylogenetic backbone of Malvaceae was reconstructed. The topology showed that Malvaceae was divided distinctly into two major branches which were previously recognized as Byttneriina and Malvadendrina. In the Malvadendrina clade, Malvoideae and Bombacoideae formed, as always, a close sister clade named as Malvatheca. Subfamily Helicteroideae occupied the most basal position and was followed by Sterculioideae which was sister to the alliance of Malvatheca, Brownlowioideae, Dombeyoideae, and Tilioideae. Brownlowioideae together with the clade comprising Dombeyoideae and Tilioideae formed a sister clade to Malvatheca. In addition, one specific conservation SSR and three specific palindrome sequences were observed in Brownlowioideae.

CONCLUSIONS

In this study, the phylogenetic framework of subfamilies in Malvaceae has been resolved clearly based on plastomes, which may contribute to a better understanding of the classification and plastome evolution for Malvaceae.

Likelihood-Based Tests of Species Tree Hypotheses

Likelihood-based tests of phylogenetic trees are a foundation of modern systematics. Over the past decade, an enormous wealth and diversity of model-based approaches have been developed for phylogenetic inference of both gene trees and species trees. However, while many techniques exist for conducting formal likelihood-based tests of gene trees, such frameworks are comparatively underdeveloped and underutilized for testing species tree hypotheses. To date, widely used tests of tree topology are designed to assess the fit of classical models of molecular sequence data and individual gene trees and thus are not readily applicable to the problem of species tree inference. To address this issue, we derive several analogous likelihood-based approaches for testing topologies using modern species tree models and heuristic algorithms that use gene tree topologies as input for maximum likelihood estimation under the multispecies coalescent. For the purpose of comparing support for species trees, these tests leverage the statistical procedures of their original gene tree-based counterparts that have an extended history for testing phylogenetic hypotheses at a single locus. We discuss and demonstrate a number of applications, limitations, and important considerations of these tests using simulated and empirical phylogenomic data sets that include both bifurcating topologies and reticulate network models of species relationships. Finally, we introduce the open-source R package SpeciesTopoTestR (SpeciesTopology Tests in R) that includes a suite of functions for conducting formal likelihood-based tests of species topologies given a set of input gene tree topologies.

Evolutionary origin, population diversity, and diagnostics for a cryptic hybrid pathogen

Cryptic fungal pathogens pose significant identification and disease management challenges due to their morphological resemblance to known pathogenic species while harboring genetic and (often) infectionrelevant trait differences. The cryptic fungal pathogen Aspergillus latus, an allodiploid hybrid originating from Aspergillus spinulosporus and an unknown close relative of Aspergillus quadrilineatus within section Nidulantes, remains poorly understood. The absence of accurate diagnostics for A. latus has led to misidentifications, hindering epidemiological studies and the design of effective treatment plans. We conducted an in-depth investigation of the genomes and phenotypes of 44 globally distributed isolates (41 clinical isolates and three type strains) from Aspergillus section Nidulantes. We found that 21 clinical isolates were A. latus; notably, standard methods of pathogen identification misidentified all A. latus isolates. The remaining isolates were identified as A. spinulosporus (8), A. quadrilineatus (1), or A. nidulans (11). Phylogenomic analyses shed light on the origin of A. latus, indicating one or two hybridization events gave rise to the species during the Miocene, approximately 15.4 to 8.8 million years ago. Characterizing the A. latus pangenome uncovered substantial genetic diversity within gene families and biosynthetic gene clusters. Transcriptomic analysis revealed that both parental genomes are actively expressed in nearly equal proportions and respond to environmental stimuli. Further investigation into infection-relevant chemical and physiological traits, including drug resistance profiles, growth under oxidative stress conditions, and secondary metabolite biosynthesis, highlight distinct phenotypic profiles of the hybrid A. latus compared to its parental and closely related species. Leveraging our comprehensive genomic and phenotypic analyses, we propose five genomic and phenotypic markers as diagnostics for A. latus species identification. These findings provide valuable insights into the evolutionary origin, genomic outcome, and phenotypic implications of hybridization in a cryptic fungal pathogen, thus enhancing our understanding of the underlying processes contributing to fungal pathogenesis. Furthermore, our study underscores the effectiveness of extensive genomic and phenotypic analyses as a promising approach for developing diagnostics applicable to future investigations of cryptic and emerging pathogens.

Exploring the Mitogenomes of Mantodea: New Insights from Structural Diversity and Higher-Level Phylogenomic Analyses

The recently reorganized classification of Mantodea has made significant progress in resolving past homoplasy problems, although some relationships among higher taxa remain uncertain. In the present study, we utilized newly sequenced mitogenomes and nuclear gene sequences of 23 mantid species, along with published data of 53 mantises, to perform familial-sampling structural comparisons of mantodean mitogenomes and phylogenomic studies. Our rstructural analysis revealed generally conserved mitogenome organizations, with a few cases of tRNA gene rearrangements, including the detection of trnL₂ duplication for the first time. In our phylogenetic analysis, we found a high degree of compositional heterogeneity and lineage-specific evolutionary rates among mantodean mitogenomes, which frequently corresponded to several unexpected groupings in the topologies under site-homogeneous models. In contrast, the topologies obtained using the site-heterogeneous mixture model fit the currently accepted phylogeny of Mantodea better. Topology tests and four-cluster likelihood mapping analyses further determined the preferred topologies. Our phylogenetic results confirm the monophyly of superfamilial groups Schizomantodea, Amerimantodea, Heteromantodea, Promantidea, and Mantidea and recover the early-branching relationships as (Mantoidoidea + (Amerimantodea + (Metallyticoidea + Cernomantodea))). Additionally, the results suggest that the long-unresolved phylogenetic position of Majangidae should be placed within Mantidea, close to Mantoidea, rather than within Epaphroditoidea. Our findings contribute to understanding the compositional and structural diversity in mantodean mitogenomes, underscore the importance of evolutionary model selection in phylogenomic studies, and provide new insights into the high-level phylogeny of Mantodea.

Statistical evaluation of character support reveals the instability of higher-level dinosaur phylogeny

The interrelationships of the three major dinosaur clades (Theropoda, Sauropodomorpha, and Ornithischia) have come under increased scrutiny following the recovery of conflicting phylogenies by a large new character matrix and its extensively modified revision. Here, we use tools derived from recent phylogenomic studies to investigate the strength and causes of this conflict. Using maximum likelihood as an overarching framework, we examine the global support for alternative hypotheses as well as the distribution of phylogenetic signal among individual characters in both the original and rescored dataset. We find the three possible ways of resolving the relationships among the main dinosaur lineages (Saurischia, Ornithischiformes, and Ornithoscelida) to be statistically indistinguishable and supported by nearly equal numbers of characters in both matrices. While the changes made to the revised matrix increased the mean phylogenetic signal of individual characters, this amplified rather than reduced their conflict, resulting in greater sensitivity to character removal or coding changes and little overall improvement in the ability to discriminate between alternative topologies. We conclude that early dinosaur relationships are unlikely to be resolved without fundamental changes to both the quality of available datasets and the techniques used to analyze them.

Divergent east-west lineages in an Australian fruit fly, (Bactrocera jarvisi), associated with the Carpentaria Basin divide

Bactrocera jarvisi is an endemic Australian fruit fly species (Diptera: Tephritidae). It occurs commonly across tropical and subtropical coastal Australia, from far-northern Western Australia, across the 'Top End' of the Northern Territory, and then down the Queensland east coast. Across this range, its distribution crosses several well documented biogeographic barriers. In order to better understand factors leading to the divergence of Australian fruit fly lineages, we carried out a population genetic study of B. jarvisi from across its range using genome-wide SNP analysis, utilising adult specimens gained from trapping and fruit rearing. Populations from the Northern Territory (NT) and Western Australia were genetically similar to each other, but divergent from the genetically uniform east-coast (= Queensland, QLD) population. Phylogenetic analysis demonstrated that the NT population derived from the QLD population. We infer a role for the Carpentaria Basin as a biogeographic barrier restricting east-west gene flow. The QLD populations were largely panmictic and recognised east-coast biogeographic barriers play no part in north-south population structuring. While the NT and QLD populations were genetically distinct, there was evidence for the historically recent translocation of flies from each region to the other. Flies reared from different host fruits collected in the same location showed no genetic divergence. While a role for the Carpentaria Basin as a barrier to gene flow for Australian fruit flies agrees with existing work on the related B. tryoni, the reason(s) for population panmixia for B. jarvisi (and B. tryoni) over the entire Queensland east coast, a linear north-south distance of >2000km, remains unknown.

Identification of a covert evolutionary pathway between two protein folds

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

Bacterial NLR-related proteins protect against phage

Bacteria use a wide range of immune pathways to counter phage infection. A subset of these genes shares homology with components of eukaryotic immune systems, suggesting that eukaryotes horizontally acquired certain innate immune genes from bacteria. Here, we show that proteins containing a NACHT module, the central feature of the animal nucleotide-binding domain and leucine-rich repeat containing gene family (NLRs), are found in bacteria and defend against phages. NACHT proteins are widespread in bacteria, provide immunity against both DNA and RNA phages, and display the characteristic C-terminal sensor, central NACHT, and N-terminal effector modules. Some bacterial NACHT proteins have domain architectures similar to the human NLRs that are critical components of inflammasomes. Human disease-associated NLR mutations that cause stimulus-independent activation of the inflammasome also activate bacterial NACHT proteins, supporting a shared signaling mechanism. This work establishes that NACHT module-containing proteins are ancient mediators of innate immunity across the tree of life.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

The Structure of Evolutionary Model Space for Proteins across the Tree of Life

The factors that determine the relative rates of amino acid substitution during protein evolution are complex and known to vary among taxa. We estimated relative exchangeabilities for pairs of amino acids from clades spread across the tree of life and assessed the historical signal in the distances among these clade-specific models. We separately trained these models on collections of arbitrarily selected protein alignments and on ribosomal protein alignments. In both cases, we found a clear separation between the models trained using multiple sequence alignments from bacterial clades and the models trained on archaeal and eukaryotic data. We assessed the predictive power of our novel clade-specific models of sequence evolution by asking whether fit to the models could be used to identify the source of multiple sequence alignments. Model fit was generally able to correctly classify protein alignments at the level of domain (bacterial versus archaeal), but the accuracy of classification at finer scales was much lower. The only exceptions to this were the relatively high classification accuracy for two archaeal lineages: Halobacteriaceae and Thermoprotei. Genomic GC content had a modest impact on relative exchangeabilities despite having a large impact on amino acid frequencies. Relative exchangeabilities involving aromatic residues exhibited the largest differences among models. There were a small number of exchangeabilities that exhibited large differences in comparisons among major clades and between generalized models and ribosomal protein models. Taken as a whole, these results reveal that a small number of relative exchangeabilities are responsible for much of the structure of the "model space" for protein sequence evolution. The clade-specific models we generated may be useful tools for protein phylogenetics, and the structure of evolutionary model space that they revealed has implications for phylogenomic inference across the tree of life.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize a diversity of secondary metabolites, small organic bioactive molecules that play diverse roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered sets of genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce diverse secondary metabolites that have been both useful (e.g., the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and harmful (e.g., the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. BGCs often also encode resistance genes that confer self-protection to the secondary metabolite-producing fungus. Some Penicillium species, such as Penicillium lilacinoechinulatum and Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity; however, an evolutionary characterization of the BGC responsible for gliotoxin biosynthesis among Penicillium species is lacking. Here, we examine the conservation of genes involved in gliotoxin biosynthesis and resistance in 35 Penicillium genomes from 23 species. We found homologous, less fragmented gliotoxin BGCs in 12 genomes, mostly fragmented remnants of the gliotoxin BGC in 21 genomes, whereas the remaining two Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, we observed broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes. Even though the gliotoxin BGC is fragmented to varying degrees in nearly all genomes examined, ancestral state reconstruction suggests that the ancestor of Penicillium species possessed the gliotoxin BGC. Our analyses suggest that genes that are part of BGCs can be retained in genomes long after the loss of secondary metabolite biosynthesis.

Exploring Evolutionary Relationships within Neodermata Using Putative Orthologous Groups of Proteins, with Emphasis on Peptidases

The phylogenetic relationships within Neodermata were examined based on putative orthologous groups of proteins (OGPs) from 11 species of Monogenea, Trematoda, and Cestoda. The dataset included OGPs from BUSCO and OMA. Additionally, peptidases were identified and evaluated as phylogenetic markers. Phylogenies were inferred using the maximum likelihood method. A network analysis and a hierarchical grouping analysis of the principal components (HCPC) of orthologous groups of peptidases were performed. The phylogenetic analyses showed the monopisthocotylean monogeneans as the sister-group of cestodes, and the polyopisthocotylean monogeneans as the sister-group of trematodes. However, the sister-group relationship between Monopisthocotylea and Cestoda was not statistically well supported. The network analysis and HCPC also showed a cluster formed by polyopisthocotyleans and trematodes. The present study supports the non-monophyly of Monogenea. An analysis of mutation rates indicated that secreted peptidases and inhibitors, and those with multiple copies, are under positive selection pressure, which could explain the expansion of some families such as C01, C19, I02, and S01. Whilst not definitive, our study presents another point of view in the discussion of the evolution of Neodermata, and we hope that our data drive further discussion and debate on this intriguing topic.

Unveiling the evolutionary relationships and the high cryptic diversity in Andean rainfrogs (Craugastoridae: Pristimantis myersi group)

Background

Pristimantis is the most diverse genus of terrestrial frogs. Historically, it has been divided into several phenetic groups in order to facilitate species identification. However, in light of phylogenetic analysis, many of these groups have been shown to be non-monophyletic, denoting a high degree of morphological convergence and limited number of diagnostic traits. In this study, we focus on the Pristimantis myersi group, an assemblage of small rainfrogs distributed throughout the Andes of Ecuador and Colombia, whose external morphology is highly conserved, and its species diversity and evolutionary relationships largely unknown.

Methods

We inferred a new phylogenetic hypothesis for the frog genus Pristimantis, including all available sequences of the mtDNA 16S rRNA, as well as new DNA sequences from 175 specimens. Our sampling included 19 of the 24 species currently recognized as part of the Pristimantis myersi group.

Results

Our new evolutionary hypothesis recovered the P. myersi group as non-monophyletic and composed of 16 species. Therefore, we exclude P. albujai, P. bicantus, P. sambalan, and P. nelsongalloi in order to preserve the monophyly of the group. We discovered at least eight candidate species, most of them hidden under the names of P. leoni, P. hectus, P. festae, P. gladiator, and P. ocreatus.

Discussion

Our results reveal the occurrence of a high level of cryptic diversity to the species level within the P. myersi group and highlight the need to redefine some of its species and reassess their conservation status. We suggest that the conservation status of six species within the group need to be re-evaluated because they exhibit smaller distributions than previously thought; these species are: P. festae, P. gladiator, P. hectus, P. leoni, P. ocreatus, and P. pyrrhomerus. Finally, given that the Pristimantis myersi group, as defined in this work, is monophyletic and morphologically diagnosable, and that Trachyphrynus is an available name for the clade containing P. myersi, we implement Trachyphrynus as a formal subgenus name for the Pristimantis myersi group.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.