New Methods to Calculate Concordance Factors for Phylogenomic Datasets

Morphology, genetic characterization and molecular phylogeny of Xystrognathus priscus (Nematoda: Rhigonematomorpha)

The infraorder Rhigonematomorpha is a unique group of zooparasitic nematodes commonly occurring in the hindgut of millipedes. Our present knowledge of the composition of the Rhigonematomorpha nematode fauna in China is far from complete. The current genetic data base for Rhigonematomorpha remains very limited, and the molecular phylogeny of Rhigonematomorpha is still in its beginning phase. In the present study, Xystrognathus priscus was redescribed using light microscopy and, for the first time, scanning electron microscopy based on the freshly collected specimens from Litostrophus scaber (Spirobolida: Pachybolidae) in China. SEM observations showed the detailed morphology of some structures for the first time, including cephalic structures bearing 1 pair of peg-like submedian cephalic papillae and 1 pair of amphidial apertures in ventral jaw, cuticular spines (microtrichs) varied from saw-tooth blade-like at anterior of body to tooth-like at posterior of body, excretory pore without cuticular spines surrounded, caudal papillae nipple-shaped, vulva without anterior vulval flap, egg with smooth surface and oval phasmid. The present study is only the second record of X. priscus, and also represents the first record of the species of the genus Xystrognathus and family Ichthyocephalidae in China. The partial 18S, 28S, ITS, cox1 and cox2 sequences of X. priscus were provided for the first time to enrich the genetic data of Rhigonematomorpha, and enabled the molecular identification of Xystrognathus nematodes. Additionally, phylogenetic results based on the partial 18S + 28S sequence data indicated that the families Ichthyocephalidae, Carnoyidae, Rhigonematidae and Xustromatidae in Rhigonematomorpha are all non-monophyletic, but strongly supported the monophyly of the genus Xystrognathus, and suggested a close affinity between Xystrognathus and the genus Obainia belonging to the Xustrostomatidae.

Multi-omics insights into the response of Aspergillus parasiticus to long-chain alkanes in relation to polyethylene modification

Plastic pollution presents a global challenge, with polyethylene (PE) being among the most persistent plastics due to its durability and environmental resilience. Long-chain alkane (lcAlk) degrading microbes are a potential source of PE-degrading enzymes, as both lcAlk and PE are large hydrophobic compounds that consist exclusively of C-C and C-H bonds. In this work, we employed a multi-omics approach to study the ability of Aspergillus parasiticus MM36, an isolate derived from Tenebrio molitor intestines, to metabolize lcAlk and secrete enzymes that are potentially capable of modifying PE. The fungus was grown with hexadecane (C16) or a mixture of lcAlk (C24 to C36) as carbon sources and culture supernatants were tested daily for their ability to modify PE. Proteomic analysis identified induced oxidases hypothetically involved in lcAlk and PE functionalization. Key enzymes include multicopper oxidases, peroxidases, an unspecific peroxygenase and FAD-dependent monooxygenases. Surfactant proteins facilitating enzymatic and cellular interaction with hydrophobic substrates, such as one hydrophobin, three hydrophobic surface-binding proteins (HsbA) and one cerato platanin, were present in all secretomes. Transcriptomic analysis comparing lcAlk to C16 cultures highlighted the enrichment of oxidoreductase activities and carboxylic acid metabolism in both lcAlk incubation days, with transmembrane transporters and transferases predominating on day 2 and biosynthetic processes on day 3. In C16 cultures, hydrolytic enzymes, including esterases, were upregulated alongside Baeyer-Villiger monooxygenases, suggesting a shift toward sub-terminal hydroxylation. Integrating transcriptomic and secretomic data, we propose a mechanism for lcAlk assimilation by A. parasiticus MM36, involving extracellular oxyfunctionalization, hydrocarbon uptake via surface-modifying proteins and channeling through membrane transporters for energy consumption and biosynthetic processes. This study provides insights into fungal mechanisms for alkane metabolism and highlights their potential relevance to plastic biotransformation.

Frequent Hybridisation Between Parapatric Lekking Bird-of-Paradise Species

Hybridisation is known to occur between a wide range of taxa, including species for which strong sexual selection has led to markedly different sexual phenotypes and lek-mating behaviours. To what extent occasional hybridisation can overcome the reproductive barriers in such systems and, for example, lead to the establishment of hybrid zones is poorly known. In this study, we address this question by focusing on one of the most well-known avian radiations in which sexual selection has resulted in an extraordinary assemblage of phenotypic diversity and lek-mating behaviours: the birds-of-paradise (Paradisaeidae). We quantify the genome-wide distribution of introgression and find multiple signals of recent and historical gene flow between and within two genera of birds-of-paradise, Astrapia and Paradigalla. In addition, we present the first empirical genomic indication of a putative hybrid zone between two lekking bird-of-paradise species that differ substantially in their sexually selected traits and behaviours. Our findings are consistent with the idea that behavioural and phenotypic traits may constitute weaker pre- and post-zygotic barriers to gene flow than generally thought in lek-mating species.

Genotyping-by-sequencing elucidates the species limits of Ukrainian taxa within Veronica subg. Pseudolysimachium

The subgenus Veronica subg. Pseudolysimachium contains mostly species occurring in natural and seminatural grasslands throughout Eurasia and is characterized by one of the highest diversification rates within Veronica, making it difficult to outline species boundaries, especially in taxa from its diversity hotspots. To resolve species limits among its Ukrainian members, the European diversity hotspot of the group, we combine evidence from GBS, niche modeling, palynology and morphometrics and compare the climatic niches of three closely related species. We confirm that Veronica maeotica and V.steppacea are distinct species. Both taxa, even though previously largely ignored, represent independent lineages in the phylogeny and are characterized each by a unique combination of morphological traits. Veronica steppacea is adapted to a drier, more continental climate than its closest relative, V.barrelieri, highlighting the importance of the zonality of the Eurasian steppe belt in shaping the diversity of the subgenus. Veronicaincana is split into two geographically distinct clades, providing evidence that plants with a cobwebby indumentum arose twice within the subgenus. Two other Ukrainian taxa, V.gryniana and V.incana subsp.hololeuca, may be of hybrid origin.

Multigene-based establishment of a novel Pleosporales family (Leandriaceae) and neotypification of Leandria momordicae Rangel

Leandria momordicae Rangel is the etiological agent of net spot, which is a very important disease of cucumber (Cucumis sativus), chayote (Sicyos edule), and in other Cucurbitaceae species. Although this pathogen was described more than a century ago infecting Momordica charantia, its taxonomic positioning has not yet been clarified. Moreover, a standard isolate of this fungus is not available in any reference mycological collection. In this context, the aims of the present work were (i) to clarify the taxonomic position of the fungus via morphological and phylogenetic analyses; (ii) to deposit a neotype of L. momordicae in a public mycological collection; (iii) to provide the first DNA sequences of the pathogen in GenBank. The morphometrical characteristics of distinct L. momordicae isolates were similar to those originally described. The phylogenetic analyses were carried out using information of five genomic regions. Using the BLASTn algorithm, it was observed that the L. momordicae-derived sequences displayed higher identity levels to fungal sequences of the order Pleosporales, which were subsequently included in the analyses. All sequences from the nine isolates of L. momordicae grouped within a single phylogenetic clade distinct from the other taxa in Pleosporales. Therefore, it is suggested the establishment of a new family of Pleosporales to include L. momordicae to which the name Leandriaceae is proposed. In addition, an ex-neotype and a neotype were illustrated and deposited in the Herbarium of the Federal University of Viçosa (VIC4754) and in the culture collection "Octávio Almeida Drummond" (COAD3403). The present study reinforces the view that many "hidden" fungal taxa are yet to be discovery under Neotropical conditions.

Genomic insights into novel extremotolerant bacteria isolated from the NASA Phoenix mission spacecraft assembly cleanrooms

BACKGROUND

Human-designed oligotrophic environments, such as cleanrooms, harbor unique microbial communities shaped by selective pressures like temperature, humidity, nutrient availability, cleaning reagents, and radiation. Maintaining the biological cleanliness of NASA's mission-associated cleanrooms, where spacecraft are assembled and tested, is critical for planetary protection. Even with stringent controls such as regulated airflow, temperature management, and rigorous cleaning, resilient microorganisms can persist in these environments, posing potential risks for space missions.

RESULTS

During the Phoenix spacecraft mission, genomes of 215 bacterial isolates were sequenced and based on overall genome-related indices, 53 strains belonging to 26 novel species were recognized. Metagenome mapping indicated less than 0.1% of the reads associated with novel species, suggesting their rarity. Genes responsible for biofilm formation, such as BolA (COG0271) and CvpA (COG1286), were predominantly found in proteobacterial members but were absent in other non-spore-forming and spore-forming species. YqgA (COG1811) was detected in most spore-forming members but was absent in Paenibacillus and non-spore-forming species. Cell fate regulators, COG1774 (YaaT), COG3679 (YlbF, YheA/YmcA), and COG4550 (YmcA, YheA/YmcA), controlling sporulation, competence, and biofilm development processes, were observed in all spore-formers but were missing in non-spore-forming species. COG analyses further revealed resistance-conferring proteins in all spore-formers (n = 13 species) and eight actinobacterial species, responsible for enhanced membrane transport and signaling under radiation (COG3253), transcription regulation under radiation stress (COG1108), and DNA repair and stress responses (COG2318). Additional functional analysis revealed that Agrococcus phoenicis, Microbacterium canaveralium, and Microbacterium jpeli contained biosynthetic gene clusters (BGCs) for ε-poly-L-lysine, beneficial in food preservation and biomedical applications. Two novel Sphingomonas species exhibited for zeaxanthin, an antioxidant beneficial for eye health. Paenibacillus canaveralius harbored genes for bacillibactin, crucial for iron acquisition. Georgenia phoenicis had BGCs for alkylresorcinols, compounds with antimicrobial and anticancer properties used in food preservation and pharmaceuticals.

CONCLUSION

Despite stringent decontamination and controlled environmental conditions, cleanrooms harbor unique bacterial species that form biofilms, resist various stressors, and produce valuable biotechnological compounds. The reduced microbial competition in these environments enhances the discovery of novel microbial diversity, contributing to the mitigation of microbial contamination and fostering biotechnological innovation. Video Abstract.

Phylogenomics establishes an Early Miocene reconstruction of reef vertebrate diversity

Oceans blanket more than two-thirds of Earth's surface, yet marine biodiversity is disproportionately concentrated in coral reefs. Investigating the origins of this exceptional diversity is crucial for predicting how reefs will respond to anthropogenic disturbances. Here, we use a genome-scale dataset to reconstruct the evolutionary history of the wrasses and parrotfishes (Labridae), which rank among the most species-rich and ecologically diverse lineages of reef fishes. We show that major labrid clades experienced pulses of evolutionary innovation and accelerated diversification during the Miocene approximately 20 to 15 million years ago that the origin of no single phenotypic trait can explain. These results draw parallels to the evolutionary histories of many clades after mass extinctions and corroborate recent fossil evidence for an Early Miocene extinction event in oceanic vertebrates and changes in coral reef faunal composition. Our data provide genomic evidence for a major Early Miocene reassembly of reef faunas.

Complete mitochondrial genome of the Korean endemic earthworm Amynthas bubonis (Clitellata: Megascolecidae): mitogenome characterization and phylogenetic positioning

Amynthas bubonis Hong & James, 2001 is an earthworm species endemic to Korea. This species is typically found in mountainous forests at low altitudes. Specimens were collected from Mt. Deogyu in South Korea, and the complete A. bubonis mitogenome was sequenced, assembled, and annotated. The A. bubonic mitogenome is a 15,095 bp circular DNA molecule with 64.85% A + T content. It contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 non-coding region (control region). Phylogenetic analysis revealed that A. bubonis clustered with A. jiriensis, A. yunoshimensis, and M. hilgendorfi in the well-supported Megascolecidae family.

Phylogenomic conflict analyses of plastid and mitochondrial genomes of Impatiens (Balsaminaceae) reveal its complex evolutionary history

Impatiens is among the most diverse angiosperm genera, comprising more than 1000 species. The phylogenetic relationship among sections of Impatiens remains unclear. Plastomes and mitogenomes are useful for resolving problematic relationships in plant phylogenetics; however, conflicts between the organellar genomes have been reported. Here, we reconstructed the phylogeny of Impatiens using concatenated and multispecies coalescent (MSC) methods based on coding and noncoding regions of the plastome and coding regions of the mitogenome from 139 species, representing all major clades in Impatiens. Conflict analyses were conducted to test and visualize the incongruences between the organellar genomes and within the plastome. The analyses supported the monophyly of all subgenera and sections and identified a new clade (clade Longlinensis), but the relationships among these sections are inconsistent. There was incongruence between the organellar genome trees regarding the relationships among sect. Semeiocardium, sect. Racemosae and clade Longlinensis. There was also incongruence within the plastome regarding the relationships among sections Fasciculatae, Impatiens, Tuberosae, Scorpioidae, and Uniflorae. Our results show that incongruence between organellar genomes likely results from the complex evolutionary history of the genus, involving mixed inheritance of organellar genomes and hybridization. The incongruence within the plastome may result from the limited phylogenetic signal in plastome data, which could be due to the rapid radiation between 15.07 and 12.93 Ma. Specific genes and regions that led to such incongruence have been identified. By confirming the monophyly of Impatiens sections and detecting phylogenetic conflicts, this study provides a unique plastid and mitogenomic perspective on Impatiens phylogeny.

Disentangling a genome-wide mosaic of conflicting phylogenetic signals in Western Rattlesnakes

Species tree inference is often assumed to be more accurate as datasets increase in size, with whole genomes representing the best-case-scenario for estimating a single, most-likely speciation history with high confidence. However, genomes may harbor a complex mixture of evolutionary histories among loci, which amplifies the opportunity for model misspecification and impacts phylogenetic inference. Accordingly, multiple distinct and well-supported phylogenetic trees are often recovered from genome-scale data, and approaches for biologically interpreting these distinct signatures are a major challenge for evolutionary biology in the age of genomics. Here, we analyze 32 whole genomes of nine taxa and two outgroups from the Western Rattlesnake species complex. Using concordance factors, topology weighting, and concatenated and species tree analyses with a chromosome-level reference genome, we characterize the distribution of phylogenetic signal across the genomic landscape. We find that concatenated and species tree analyses of autosomes, the Z (sex) chromosome, and mitochondrial genome yield distinct, yet strongly supported phylogenies. Analyses of site-specific likelihoods show additional patterns consistent with rampant model misspecification, a likely consequence of several evolutionary processes. Together, our results suggest that a combination of historic and recent introgression, along with natural selection, recombination rate variation, and cytonuclear co-evolution of nuclear-encoded mitochondrial genes, underlie genome-wide variation in phylogenetic signal. Our results highlight both the power and complexity of interpreting whole genomes in a phylogenetic context and illustrate how patterns of phylogenetic discordance can reveal the impacts of different evolutionary processes that contribute to genome-wide variation in phylogenetic signal.

Unraveling the extensive phylogenetic discordance and evolutionary history of spurless taxa within the Aquilegia ecalcarata complex

Parallel evolution of the same, or at least very similar, phenotype(s) in different lineages is often interpreted as evidence for the action of natural selection. However, caution is required when inferring parallel evolution based on uncertain or potentially incorrect phylogenetic frameworks. Here, by conducting extensive phylogenomic and population genetic analyses, we aim to clarify the evolutionary history of spurless taxa within the Aquilegia ecalcarata complex. We observed substantial discordance in the phylogenetic patterns across the entire genome, primarily attributed to ancient introgression and incomplete lineage sorting. Additionally, we identified several spurless lineages whose phylogenetic positions were distorted by admixture events. Using a backbone tree and demographic modeling, we determined that these spurless taxa independently originated twice within this group. Intriguingly, our investigation revealed that the spurless taxa experienced population expansion during global cooling, while their spurred sister groups underwent population contraction. The parallel losses of petal spurs, therefore, may be linked to adaptations for low-temperature conditions. These findings emphasize the importance of comprehensive population-level analyses in phylogenetic inference and provide valuable insights into the dynamics of trait loss and its implications for the adaptive strategies.

Incomplete lineage sorting and introgression among genera and species of Liliaceae tribe Tulipeae: insights from phylogenomics

BACKGROUND

Phylogenetic research in Tulipa (Liliaceae), a genus of significant economic and horticultural value, has relied on limited nuclear (mostly nuclear ribosomal internal transcribed spacer, nrITS) and plastid DNA sequences, resulting in low-resolution phylogenetic trees and uncertain intrageneric classifications. The genus, noted for its large genome, presents discordant relationships among Amana, Erythronium, and Tulipa, likely due to incomplete lineage sorting (ILS) and/or reticulate evolution. Thus, phylogenomic approaches are needed to clarify these relationships and the conflicting signals within the tribe Tulipeae.

RESULTS

We newly sequenced 50 transcriptomes of 46 species of tribe Tulipeae (including multiple accessions of all four genera) and one outgroup species of the sister tribe Lilieae (Notholirion campanulatum), and downloaded 15 previously published transcriptomes of tribe Tulipeae to supplement the sampling. One plastid dataset (74 plastid protein-coding genes, PCGs) and one nuclear dataset (2594 nuclear orthologous genes, OGs) were constructed, with the latter used for species tree inference based on maximum likelihood (ML) and multi-species coalescent (MSC) methods. To investigate causes of gene tree discordance, "site con/discordance factors" (sCF and sDF1/sDF2) were calculated first, after which phylogenetic nodes displaying high or imbalanced sDF1/2 were selected for phylogenetic network analyses and polytomy tests to determine whether ILS or reticulate evolution best explain incongruence. Key relationships not resolved by this technique, especially those among Amana, Erythronium, and Tulipa, were further investigated by applying D-statistics and QuIBL.

CONCLUSIONS

We failed to reconstruct a reliable and unambiguous evolutionary history among Amana, Erythronium, and Tulipa due to especially pervasive ILS and reticulate evolution, likely caused either by obscured minority phylogenetic signal or differing signals among genomic compartments. However, within Tulipa we confirmed the monophyly of most subgenera, with the exception of two species in the small subgenus Orithyia, of which Tulipa heterophylla was recovered as sister to the remainder of the genus, whereas T. sinkiangensis clustered within subgenus Tulipa. In contrast, most traditional sections of Tulipa were found to be non-monophyletic.

Reconstructing the phylogeny and evolutionary history of freshwater fishes (Nemacheilidae) across Eurasia since early Eocene

Eurasia has undergone substantial tectonic, geological, and climatic changes throughout the Cenozoic, primarily associated with tectonic plate collisions and a global cooling trend. The evolution of present-day biodiversity unfolded in this dynamic environment, characterised by intricate interactions of abiotic factors. However, comprehensive, large-scale reconstructions illustrating the extent of these influences are lacking. We reconstructed the evolutionary history of the freshwater fish family Nemacheilidae across Eurasia and spanning most of the Cenozoic on the base of 471 specimens representing 279 species and 37 genera plus outgroup samples. Molecular phylogeny using six genes uncovered six major clades within the family, along with numerous unresolved taxonomic issues. Dating of cladogenetic events and ancestral range estimation traced the origin of Nemacheilidae to Indochina around 48 mya. Subsequently, one branch of Nemacheilidae colonised eastern, central, and northern Asia, as well as Europe, while another branch expanded into the Burmese region, the Indian subcontinent, the Near East, and northeast Africa. These expansions were facilitated by tectonic connections, favourable climatic conditions, and orogenic processes. Conversely, aridification emerged as the primary cause of extinction events. Our study marks the first comprehensive reconstruction of the evolution of Eurasian freshwater biodiversity on a continental scale and across deep geological time.

Acanthocephalans from freshwater fishes in northeast Thailand with the description of a new species of the subgenus Acanthosentis Verma & Dutta, 1929 (Acanthocephala, Quadrigyridae)

During an ichthyoparasitological survey in northeast Thailand in 2015, four species of acanthocephalans were found in four species of freshwater fishes. Adult stages of Pallisentisrexus Wongkham & Whitfield, 1999 and Raosentis sp. (Eoacanthocephala, Quadrigyridae) were collected from Channastriata (Anabantiformes, Channidae) and Mystusmysticetus (Siluriformes, Bagridae), respectively, and cystacanths of Arhythmorhynchus sp. (Palaeacanthocephala, Polymorphidae) were found on the visceral organs of M.albolineatus (Siluriformes, Bagridae). Adult acanthocephalans of the subgenus Acanthosentis Verma & Dutta, 1929 isolated from Barbonymusschwanenfeldii (Cypriniformes, Cyprinidae) were morphologically distinct from all described species of the subgenus in the arrangement of rings of tegumental spines in two fields with a more or less pronounced distance between them, and by the presence of a dome-shaped process with a ring of small spines at the base at the posterior end in females. Molecular data were generated for three species and phylogenetic analysis based on the 18S rDNA clearly supported the generic identification of Acanthogyrus (Acanthosentis) barbonymisp. nov. and P.rexus. While the phylogenetic position of the former species within the genus could not be clarified, the latter species formed a well-supported sister lineage in a clade with isolates of four congeneric species. Acanthogyrus (Acanthosentis) barbonymisp. nov. is formally described, the first genetic data for P.rexus are generated, a species of the genus Raosentis Datta, 1947 is reported for the first time outside the Indian subcontinent, and a key to the species of the latter genus is presented.

Road to Extinction? Past and Present Population Structure and Genomic Diversity in the Koala

Koalas are arboreal herbivorous marsupials, endemic to Australia. During the late 1800s and early 1900s, the number of koalas declined dramatically due to hunting for their furs. In addition, anthropogenic activities have further decimated their available habitat, and decreased population numbers. Here, we utilize 37 historic and 25 modern genomes sampled from across their historic and present geographic range, to gain insights into how their population structure and genetic diversity have changed across time; assess the genetic consequences of the period of intense hunting, and the current genetic status of this iconic Australian species. Our analyses reveal how genome-wide heterozygosity has decreased through time and unveil previously uncharacterized mitochondrial haplotypes and nuclear genotypes in the historic dataset, which are absent from today's koala populations.

Origins and Distribution of Panicum Mosaic Virus and Sugarcane Mosaic Virus on Stenotaphrum secundatum in Australia

Stenotaphrum secundatum is a premium turf grass in warm temperate and subtropical regions of the world and is the most important turf species in Australia based on the value of its production. A new disease called buffalo grass yellows (BGY) has become a problem on turf farms in Australia. We surveyed turf farms in New South Wales, Queensland, and Western Australia to determine whether panicum mosaic virus (PMV) and sugarcane mosaic virus (SCMV) were associated with BGY. PMV was only found on three farms, two located in the Hawkesbury Valley near Sydney and a third at Echuca, about 800 km to the southwest of the former location. SCMV was more prevalent, present in all major cultivars and states surveyed. We analyzed phylogenetic relationships for SCMV and found that isolates infecting S. secundatum in Australia belonged to three clades. The first included Australian isolates typical of the population of viruses circulating in Digitaria didactyla. The second included a single New South Wales isolate from S. secundatum 'SS100' that grouped with otherwise American isolates of SCMV recorded in S. secundatum and Saccharum officinale from Florida and Zea mays from Ohio. Finally, an isolate of SCMV from S. secundatum originating from a turf farm in South East Queensland grouped with viruses mostly infecting maize; this record is potentially the first maize-adapted strain of SCMV in Australia. Our study sheds light on the etiology of the BGY disease syndrome and invasion history of PMV and SCMV in Australia.

Biogeographic History of Pigeons and Doves Drives the Origin and Diversification of Their Parasitic Body Lice

Despite their extensive diversity and ecological importance, the history of diversification for most groups of parasitic organisms remains relatively understudied. Elucidating broad macroevolutionary patterns of parasites is challenging, often limited by the availability of samples, genetic resources, and knowledge about ecological relationships with their hosts. In this study, we explore the macroevolutionary history of parasites by focusing on parasitic body lice from doves. Building on extensive knowledge of ecological relationships and previous phylogenomic studies of their avian hosts, we tested specific questions about the evolutionary origins of the body lice of doves, leveraging whole genome data sets for phylogenomics. Specifically, we sequenced whole genomes from 68 samples of dove body lice, including representatives of all body louse genera from 51 host taxa. From these data, we assembled > 2300 nuclear genes to estimate dated phylogenetic relationships among body lice and several outgroup taxa. The resulting phylogeny of body lice was well supported, although some branches had conflicting signals across the genome. We then reconstructed ancestral biogeographic ranges of body lice and compared the body louse phylogeny to the phylogeny of doves, and also to a previously published phylogeny of the wing lice of doves. Divergence estimates placed the origin of body lice in the late Oligocene. Body lice likely originated in Australasia and dispersed with their hosts during the early Miocene, with subsequent codivergence and host switching throughout the world. Notably, this evolutionary history is very similar to that of dove wing lice, despite the stronger dispersal capabilities of wing lice compared to body lice. Our results highlight the central role of the biogeographic history of host organisms in driving the evolutionary history of their parasites across time and geographic space.

A haplotype-resolved reference genome of Quercus alba sheds light on the evolutionary history of oaks

White oak (Quercus alba) is an abundant forest tree species across eastern North America that is ecologically, culturally, and economically important. We report the first haplotype-resolved chromosome-scale genome assembly of Q. alba and conduct comparative analyses of genome structure and gene content against other published Fagaceae genomes. We investigate the genetic diversity of this widespread species and the phylogenetic relationships among oaks using whole genome data. Despite strongly conserved chromosome synteny and genome size across Quercus, certain gene families have undergone rapid changes in size, including defense genes. Unbiased annotation of resistance (R) genes across oaks revealed that the overall number of R genes is similar across species - as are the chromosomal locations of R gene clusters - but, gene number within clusters is more labile. We found that Q. alba has high genetic diversity, much of which predates its divergence from other oaks and likely impacts divergence time estimations. Our phylogenetic results highlight widespread phylogenetic discordance across the genus. The white oak genome represents a major new resource for studying genome diversity and evolution in Quercus. Additionally, we show that unbiased gene annotation is key to accurately assessing R gene evolution in Quercus.

Phylogenomic insights into the taxonomy, ecology, and mating systems of the lorchel family Discinaceae (Pezizales, Ascomycota)

Lorchels, also known as false morels (Gyromitra sensu lato), are iconic due to their brain-shaped mushrooms and production of gyromitrin, a deadly mycotoxin. Molecular phylogenetic studies have hitherto failed to resolve deep-branching relationships in the lorchel family, Discinaceae, hampering our ability to settle longstanding taxonomic debates and to reconstruct the evolution of toxin production. We generated 75 draft genomes from cultures and ascomata (some collected as early as 1960), conducted phylogenomic analyses using 1542 single-copy orthologs to infer the early evolutionary history of lorchels, and identified genomic signatures of trophic mode and mating-type loci to better understand lorchel ecology and reproductive biology. Our phylogenomic tree was supported by high gene tree concordance, facilitating taxonomic revisions in Discinaceae. We recognized 10 genera across two tribes: tribe Discineae (Discina, Maublancomyces, Neogyromitra, Piscidiscina, and Pseudodiscina) and tribe Gyromitreae (Gyromitra, Hydnotrya, Paragyromitra, Pseudorhizina, and Pseudoverpa); Piscidiscina was newly erected and 26 new combinations were formalized. Paradiscina melaleuca and Marcelleina donadinii formed their own family-level clade sister to Morchellaceae, which merits further taxonomic study. Genome size and CAZyme content were consistent with a mycorrhizal lifestyle for the truffle species (Hydnotrya spp.), whereas the other Discinaceae genera possessed genomic properties of a saprotrophic habit. Lorchels were found to be predominantly heterothallic-either MAT1-1 or MAT1-2-but a single occurrence of colocalized mating-type idiomorphs indicative of homothallism was observed in Gyromitra esculenta strain CBS101906 and requires additional confirmation and follow-up study. Lastly, we confirmed that gyromitrin has a phylogenetically discontinuous distribution, having been detected exclusively in two distantly related genera (Gyromitra and Piscidiscina) belonging to separate tribes. Our genomic dataset will facilitate further investigations into the gyromitrin biosynthesis genes and their evolutionary history. With additional sampling of Geomoriaceae and Helvellaceae-two closely related families with no publicly available genomes-these data will enable comprehensive studies on the independent evolution of truffles and ecological diversification in an economically important group of pezizalean fungi.

Exon Capture Museomics Deciphers the Nine-Banded Armadillo Species Complex and Identifies a New Species Endemic to the Guiana Shield

The nine-banded armadillo (Dasypus novemcinctus) is the most widespread xenarthran species across the Americas. Recent studies have suggested it is composed of 4 morphologically and genetically distinct lineages of uncertain taxonomic status. To address this issue, we used a museomic approach to sequence 80 complete mitogenomes and capture 997 nuclear loci for 71 Dasypus individuals sampled across the entire distribution. We carefully cleaned up potential genotyping errors and cross-contaminations that could blur species boundaries by mimicking gene flow. Our results unambiguously support 4 distinct lineages within the D. novemcinctus complex. We found cases of mito-nuclear phylogenetic discordance but only limited contemporary gene flow confined to the margins of the lineage distributions. All available evidence including the restricted gene flow, phylogenetic reconstructions based on both mitogenomes and nuclear loci, and phylogenetic delimitation methods consistently supported the 4 lineages within D. novemcinctus as 4 distinct species. Comparable genetic differentiation values to other recognized Dasypus species further reinforced their status as valid species. Considering congruent morphological results from previous studies, we provide an integrative taxonomic view to recognize 4 species within the D. novemcinctus complex: D. novemcinctus, D. fenestratus, D. mexicanus, and D. guianensis sp. nov., a new species endemic of the Guiana Shield that we describe here. The 2 available individuals of D. mazzai and D. sabanicola were consistently nested within D. novemcinctus lineage and their status remains to be assessed. The present work offers a case study illustrating the power of museomics to reveal cryptic species diversity within a widely distributed and emblematic species of mammals.

Description of a new species of Stauropathes (Anthozoa, Antipatharia, Schizopathidae) from Puerto Rico

A new species of black coral, Stauropathesmonopinnata sp. nov., represented by two specimens collected 738 m and 1604 m deep off Puerto Rico and Hawaii, respectively, is recognized in the family Schizopathidae. The new species is characterized by a monopodial, unbranched corallum; simple, suboppositely arranged pinnules in two anterolateral rows along the stem with nearly 90° distal angles, spaced 12-17 mm apart in a row, and with smooth and triangular spines 0.05-0.08 mm tall; and polyps 4-9 mm in transverse diameter. A phylogeny composed of 90 taxa representing species in Schizopathidae and Cladopathidae (rooted in Cladopathidae) was reconstructed from 794 nuclear loci to show their systematic relationships. Herein, we provide morphological and molecular evidence to show that this new species is distinct from other species in the genus Stauropathes.

Annotated Bioinformatic Pipelines for Phylogenomic Placement of Mitochondrial Genomes

The limited standards for the rigorous and objective use of mitochondrial genomes (mitogenomes) can lead to uncertainties regarding the phylogenetic relationships of taxa under varying evolutionary constraints. The mitogenome exhibits heterogeneity in base composition, and evolutionary rates may vary across different regions, which can cause empirical data to violate assumptions of the applied evolutionary models. Consequently, the unique evolutionary signatures of the dataset must be carefully evaluated before selecting an appropriate approach for phylogenomic inference. Here, we present the bioinformatic pipeline and code used to expand the mitogenome phylogeny of the order Carcharhiniformes (groundsharks), with a focus on houndsharks (Chondrichthyes: Triakidae). We present a rigorous approach for addressing difficult-to-resolve phylogenies, incorporating multi-species coalescent modelling (MSCM) to address gene/species tree discordance. The protocol describes carefully designed approaches for preparing alignments, partitioning datasets, assigning models of evolution, inferring phylogenies based on traditional site-homogenous concatenation approaches as well as under multispecies coalescent and site heterogenous models, and generating statistical data for comparison of different topological outcomes. The datasets required to run our analyses are available on GitHub and Dryad repositories. Key features • An extensive statistical framework to conduct model selection and data partitioning and tackle difficult-to-resolve phylogenies. • Instructions for generating statistical data for comparison of different topological outcomes. • Tips for selecting mitochondrial phylogenomic (mitophylogenomic) approaches to suit unique datasets. • Access to the scripts, data files, and pipelines used to enable replication of all analyses.

Morphology, complete mitochondrial genome and molecular phylogeny of Heterakis pucrasia sp. n. (Nematoda: Ascaridida) from the koklass pheasant Pucrasia macrolopha (Lesson) (Galliformes: Phasianidae) in Pakistan

Species of Heterakis (Ascaridida: Heterakoidea) are commonest nematode parasites occurring in the alimentary canal of wild and domestic birds, which are of major socio-economic importance, due to some Heterakis species causing Heterakidosis in wild birds and poultry. In the present study, a new species of Heterakis, H. pucrasia sp. n., was described using integrated methods based on specimens collected from the koklass pheasant Pucrasia macrolopha (Lesson) (Galliformes: Phasianidae) in Pakistan. The complete mitochondrial genome of H. pucrasia sp. n. was sequenced and annotated for the first time to enrich the mitogenomic data, and reveal the pattern of mitogenomic evolution of the family Heterakidae. Moreover, phylogenetic analyses of the orders Ascaridida, Spirurida, Oxyurida and Rhigonematida based on the amino acid sequences of 12 protein coding genes (PCGs) of mitochondrial genomes, revealed that the order Ascaridida is not monophyletic, and the superfamily Heterakoidea has a closer affinity with Rhigonematida + Oxyurida + Spirurida, than the superfamily Ascaridoidea in Ascaridida. The present findings enriched the global species composition of heterakid nematodes and their mitogenomic data, and also provided novel insight on the phylogenetic relationships between Heterakoidea and its related groups.

Phylogenomic analyses re-evaluate the backbone of Corylus and unravel extensive signals of reticulate evolution

Phylogenomic analyses have shown that reticulate evolution greatly affects the accuracy of phylogenetic inferences, and thus may challenge the authority of bifurcating phylogenetic trees. In this study, we re-evaluated the phylogenetic backbone of the genus Corylus based on complete taxon sampling and genomic data. We assembled 581 single-copy nuclear genes and whole plastomes from 64 genome resequencing datasets to elucidate the reticulate relationships within Corylus. Nuclear coalescent and concatenation phylogenies revealed identical and fully supported backbone, clarifying the sisterhood between sect. Acanthochlamys and sect. Siphonochlamys as well as the phylogenetic position of C. fargesii and C. wangii, which have yet been addressed in previous phylogenetic studies. However, the monophyly of C. jacquemontii and C. kwechowensis and the distinction between C. ferox and C. ferox var. thibetica were not supported. Gene trees-species tree conflicts and cytonuclear discordance were identified, with multiple evidences supporting that hybridization/introgression, coupled with incomplete lineage sorting, have led to substantial phylogenetic incongruence in Corylus. Moreover, typical geographical clustering rather than strict monophyletic pattern in plastome phylogeny implies chloroplast capture within Corylus and offers evidence of cytoplasmic introgression. Overall, this study provides a robust phylogenomic backbone for Corylus and unravels that reticulate evolution can greatly shape taxonomic revision.

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.

Rapid Evolution of Host Repertoire and Geographic Range in a Young and Diverse Genus of Montane Butterflies

Evolutionary changes in geographic distribution and larval host plants may promote the rapid diversification of montane insects, but this scenario has been rarely investigated. We studied the rapid radiation of the butterfly genus Colias, which has diversified in mountain ecosystems in Eurasia, Africa, and the Americas. Based on a data set of 150 nuclear protein-coding genetic loci and mitochondrial genomes, we constructed a time-calibrated phylogenetic tree of Colias species with broad taxon sampling. We then inferred their ancestral geographic ranges, historical diversification rates, and the evolution of host use. We found that the most recent common ancestor of Colias was likely geographically widespread and originated ~3.5 Ma. The group subsequently diversified in different regions across the world, often in tandem with geographic expansion events. No aspect of elevation was found to have a direct effect on diversification. The genus underwent a burst of diversification soon after the divergence of the Neotropical lineage, followed by an exponential decline in diversification rate toward the present. The ancestral host repertoire included the legume genera Astragalus and Trifolium but later expanded to include a wide range of Fabaceae genera and plants in more distantly related families, punctuated with periods of host range expansion and contraction. We suggest that the widespread distribution of the ancestor of all extant Colias lineages set the stage for diversification by isolation of populations that locally adapted to the various different environments they encountered, including different host plants. In this scenario, elevation is not the main driver but might have accelerated diversification by isolating populations.

Phylogenetic Inferences and Historical Biogeography of Onocleaceae

The family Onocleaceae represents a small family of terrestrial ferns, with four genera and around five species. It has a circumboreal to north temperate distribution, and exhibits a disjunct distribution between Eurasia and North America, including Mexico. Historically, the taxonomy and classification of this family has been subject to debate and contention among scholars, leading to contradictory classifications and disagreements on the number of genera and species within the family. Furthermore, due to this disjunct intercontinental distribution and the lack of detailed study across its wide range, this family merits further study to clarify its distributional pattern. Maximum likelihood and Bayesian phylogenetic reconstructions were based on a concatenated sequence dataset for 17 plastid loci and one nuclear locus, which were generated from 106 ingroup and six outgroup taxa from three families. Phylogenetic analyses support that Onocleaceae is composed of four main clades, and Pentarhizidium was recovered as the first branching lineages in Onocleaceae. Molecular dating and ancestral area reconstruction analyses suggest that the stem group of Onocleaceae originated in Late Cretaceous, with subsequent diversification and establishment of the genera Matteuccia, Onoclea, Onocleopsis, and Pentarhizidium during the Paleogene and Neogene. The ancestors of Matteuccia, Onoclea, and Onocleopsis could have migrated to North America via the Beringian land bridge or North Atlantic land bridge which suggests that the diversification of Matteuccia + Onoclea + Onocleopsis closely aligns with the Paleocene-Eocene Thermal Maximum (PETM). In addition, these results suggest that Onocleaceae species diversity peaks during the late Neogene to Quaternary. Studies such as this enhance our understanding of the mechanisms and climatic conditions shaping disjunct distribution in ferns and lycophytes of eastern Asia, North America, and Mexico and contribute to a growing body of evidence from other taxa, to advance our understanding of the origins and migration of plants across continents.

The Genomic Landscape, Causes, and Consequences of Extensive Phylogenomic Discordance in Murine Rodents

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents-a large and ecologically diverse group that gave rise to the laboratory mouse and rat model systems. Combining recently published linked-read genome assemblies for seven murine species with other available rodent genomes, we first used ultraconserved elements (UCEs) to infer a robust time-calibrated species tree. We then used whole genomes to examine finer-scale patterns of discordance across ∼12 million years of divergence. We found that proximate chromosomal regions tended to have more similar phylogenetic histories. There was no clear relationship between local tree similarity and recombination rates in house mice, but we did observe a correlation between recombination rates and average similarity to the species tree. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to appreciably less discordance. Finally, we show that assuming a single species tree can result in substantial deviation from the results with gene trees when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Reassessing the evolutionary relationships of tropical wandering spiders using phylogenomics: A UCE-based phylogeny of Ctenidae (Araneae) with the discovery of a new lycosoid family

Tropical wandering spiders (Ctenidae) are a diverse family of cursorial predators whose species richness peaks in the tropics. The phylogeny of Ctenidae has been examined using morphology and Sanger-based sequencing data, but these studies have been limited by taxon sampling and have often recovered low branch support for many intrafamilial phylogenetic relationships. Herein, we present the most extensive phylogenetic sampling of this family using genome-scale data, leveraging museum collections of all ctenid subfamilies from across the world. We obtained a well-resolved phylogeny of Ctenidae, with the majority of nodes showing maximal nodal support and topological congruence across different phylogenetic analyses. For the first time, we show with high support that Ancylometes is not within Ctenidae but is the sister lineage to all the remaining lycosoid families. Therefore, we propose Ancylometidae as a new family. We assess the phylogenetic position of Ctenidae within Lycosoidea using a variety of phylogenetic methods and tests, demonstrating that the previously proposed position of Ctenidae as the sister clade of Psechridae, based on phylotranscriptomic analyses, lacks phylogenetic support. As a new finding, this study shows that the subfamily Acantheinae, as currently delimited, is polyphyletic. Therefore, we erect the new ctenid subfamily Enoplocteninae to accommodate the Neotropical genera Enoploctenus, Chococtenus, and Phymatoctenus. Our phylogenomic results using UCE data resolve the position of several problematic genera (e.g., Califorctenus and Acantheis) and add support to other parts of the tree that received low support in the most recent Sanger-based phylogeny. We discuss some of the putative morphological synapomorphies of the main ctenid lineages within the phylogenetic framework provided by the molecular phylogenetic results of this study.

Clade I mpox virus genomic diversity in the Democratic Republic of the Congo, 2018-2024: Predominance of zoonotic transmission

Recent reports raise concerns on the changing epidemiology of mpox in the Democratic Republic of the Congo (DRC). High-quality genomes were generated for 337 patients from 14/26 provinces to document whether the increase in number of cases is due to zoonotic spillover events or viral evolution, with enrichment of APOBEC3 mutations linked to human adaptation. Our study highlights two patterns of transmission contributing to the source of human cases. All new sequences from the eastern South Kivu province (n = 17; 4.8%) corresponded to the recently described clade Ib, associated with sexual contact and sustained human-to-human transmission. By contrast, all other genomes are clade Ia, which exhibits high genetic diversity with low numbers of APOBEC3 mutations compared with clade Ib, suggesting multiple zoonotic introductions. The presence of multiple clade I variants in urban areas highlights the need for coordinated international response efforts and more studies on the transmission and the reservoir of mpox.

Whole genome data confirm pervasive gene discordance in the evolutionary history of Coenonympha (Nymphalidae) butterflies

Phylogenetic inference is challenged by genealogical heterogeneity amongst molecular markers. Such discordance is driven predominantly by incomplete lineage sorting (ILS) and interspecific gene flow, and bears attendant consequences for the accurate resolution of species relationships. Understanding the distribution of gene conflict in organismal genomes is, therefore, a key aspect of phylogenetic analysis. In this study, three large phylogenomic datasets (i.e., whole mitogenomes, conserved nuclear protein-coding loci, and genomic windows) are used to probe the extent to which discordance pervades the unresolved phylogeny of Coenonympha (Nymphalidae) butterflies. Gene tree discordance is found to be elevated at multiple historically recalcitrant phylogenetic positions. In particular, species relationships near the crown of Coenonympha and within a rapidly diversifying subclade (the hero group) remain difficult to resolve, suggesting that ILS and gene flow have obscured the evolution of this genus. These findings have implications for the taxonomy of this butterfly group and the study of its diversification history. In addition, this work lends support to a growing body of evidence that gene conflict driven by biological processes stands to confound phylogeny, even when extensive data are used.

Population Genomics of Adaptive Radiation

Adaptive radiations are rich laboratories for exploring, testing, and understanding key theories in evolution and ecology because they offer spectacular displays of speciation and ecological adaptation. Particular challenges to the study of adaptive radiation include high levels of species richness, rapid speciation, and gene flow between species. Over the last decade, high-throughput sequencing technologies and access to population genomic data have lessened these challenges by enabling the analysis of samples from many individual organisms at whole-genome scales. Here we review how population genomic data have facilitated our knowledge of adaptive radiation in five key areas: (1) phylogenetics, (2) hybridization, (3) timing and rates of diversification, (4) the genomic basis of trait evolution, and (5) the role of genome structure in divergence. We review current knowledge in each area, highlight outstanding questions, and focus on methods that facilitate detection of complex patterns in the divergence and demography of populations through time. It is clear that population genomic data are revolutionising the ability to reconstruct evolutionary history in rapidly diversifying clades. Additionally, studies are increasingly emphasising the central role of gene flow, re-use of standing genetic variation during adaptation, and structural genomic elements as facilitators of the speciation process in adaptive radiations. We highlight hybridization-and the hypothesized processes by which it shapes diversification-and questions seeking to bridge the divide between microevolutionary and macroevolutionary processes as rich areas for future study. Overall, access to population genomic data has facilitated an exciting era in adaptive radiation research, with implications for deeper understanding of fundamental evolutionary processes across the tree of life.

Phylogenomic resolution of lampreys reveals the recent evolution of an ancient vertebrate lineage

Jawless vertebrates once dominated Palaeozoic waters, but just two lineages have persisted to the present day: lampreys and hagfishes. Living lampreys are a relatively small clade, with just over 50 species described, but knowledge of their evolutionary relationships has always been based on either a few mitochondrial genes or a small number of taxa. Biogeographers have noted the disjunct antitropical distribution of living lamprey families. Here, we present a comprehensive phylogenomic analysis of living and fossil lampreys, sampling 36 species with phylogenomic data and 46 in total with genetic data. We present new divergence time estimates based on comprehensive nuclear data and analysis of their diversification dynamics. Our analysis indicates a central role for extreme global warming during the Late Cretaceous Cenomanian-Turonian Boundary Event as a likely cause for the antitropical distribution of living lampreys, and a notable increase in lineage diversification in Northern Hemisphere lampreys during the Miocene corresponding with a period of global cooling.

Insect-specific RNA viruses detection in Field-Caught Aedes aegypti mosquitoes from Argentina using NGS technology

Mosquitoes are the primary vectors of arthropod-borne pathogens. Aedes aegypti is one of the most widespread mosquito species worldwide, responsible for transmitting diseases such as Dengue, Zika, and Chikungunya, among other medically significant viruses. Characterizing the array of viruses circulating in mosquitoes, particularly in Aedes aegypti, is a crucial tool for detecting and developing novel strategies to prevent arbovirus outbreaks. In this study, we address the implementation of a sequencing and analysis pipeline based on the Oxford Nanopore Technologies MinION Mk1b system, for arboviral detection in field-caught mosquitoes from Argentina. Full genome of Humaita Tubiacanga Virus (HTV), Phasi Charoen-like Phasivirus (PCLV), Aedes aegypti totivirus (AaeTV) has been sequenced in three distinct regions of Argentina comprising Buenos Aires province, Santa Fe province and the northern province of Salta. Viral sequences enriched by SISPA and coupled with Nanopore sequencing can be a useful tool for viral surveillance, not only for detecting viruses that have a high impact on human and animal health, but also for detecting insect-specific viruses that could promote the transmission of arboviruses.

Leaf transcriptomes from C3, C3-C4 intermediate, and C4Neurachne species give insights into C4 photosynthesis evolution

The C4 photosynthetic pathway is hypothesized to have evolved from the ancestral C3 pathway through progressive changes in leaf anatomy and biochemistry with extant C3-C4 photosynthetic intermediate species representing phenotypes between species demonstrating full C3 and full C4 states. The Australian endemic genus Neurachne is the only known grass group that contains distinct, closely related species that carry out C3, C3-C4 intermediate, or C4 photosynthesis. To explore and understand the molecular mechanisms underlying C4 photosynthesis evolution in this genus, leaf transcriptomes were generated from two C3, three photosynthetic intermediate (proto-Kranz, C2-like, and C2), and two C4Neurachne species. The data were used to reconstruct phylogenetic relationships in Neurachne, which confirmed two independent C4 origins in the genus. Relative transcript abundances substantiated the photosynthetic phenotypes of individual species and highlighted transcriptional investment differences between species, including between the two C4 species. The data also revealed proteins potentially involved in C4 cycle intermediate transport and identified molecular mechanisms responsible for the evolution of C4-associated proteins in the genus.

The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms

The rise of angiosperms to ecological dominance and the breakup of Gondwana during the Mesozoic marked major transitions in the evolutionary history of insect-plant interactions. To elucidate how contemporary trophic interactions were influenced by host plant shifts and palaeogeographical events, we integrated molecular data with information from the fossil record to construct a time tree for ancient phytophagous weevils of the beetle family Belidae. Our analyses indicate that crown-group Belidae originated approximately 138 Ma ago in Gondwana, associated with Pinopsida (conifer) host plants, with larvae likely developing in dead/decaying branches. Belids tracked their host plants as major plate movements occurred during Gondwana's breakup, surviving on distant, disjunct landmasses. Some belids shifted to Angiospermae and Cycadopsida when and where conifers declined, evolving new trophic interactions, including brood-pollination mutualisms with cycads and associations with achlorophyllous parasitic angiosperms. Extant radiations of belids in the genera Rhinotia (Australian region) and Proterhinus (Hawaiian Islands) have relatively recent origins.

Revisiting evolutionary relationships of Antrodiaetus (Araneae, Mygalomorphae, Antrodiaetidae) using phylogenomics; implications for species diversity and biogeography of a persistent Holarctic lineage

Antrodiaetus is a lineage of mygalomorph spider (Mygalomorphae: Antrodiaetidae) that has persisted since the late Cretaceous and has a disjunct Holarctic distribution and strong morphological conservatism. These folding-door spiders possess a life history (i.e., limited dispersal, conserved environmental niche) that closely ties their evolution to geology. This study produces a robust, well-supported phylogenomic inference of all currently recognized Antrodiaetus species using UCEs (Ultraconserved Elements), corroborates previous biogeographical hypotheses, and proposes new hypotheses about diversification patterns. We also confirm that previously suspected cryptic diversity within A. pacificus is underestimated, as this nominal species comprises multiple divergent and cryptic lineages. Our phylogeny now serves as a foundation for understanding Antrodiaetus species relationships, biogeography, and speciation.

A global blueberry phylogeny: Evolution, diversification, and biogeography of Vaccinieae (Ericaceae)

Vaccinieae is a morphologically diverse and species-rich (∼1430 species) tribe in Ericaceae. Although the majority of diversity is tropical, Vaccinieae are best known for temperate crops (i.e., blueberries, cranberries, and lingonberries) in Vaccinium. Vaccinium itself (∼500 species) has been previously suggested as highly polyphyletic and taxonomic boundaries among many of the other genera in the tribe remain uncertain. We assessed the evolutionary history of Vaccinieae with phylogenomic analyses based on a target-enrichment dataset containing 256 low-copy nuclear loci and 210 species representing 30 of the 35 genera in the tribe and 25 of the 29 sections of Vaccinium. We conducted time-calibrated biogeographic analyses and diversification analyses to explore the area of origin and global dispersal history of the tribe. The analysis recovered a temperate North American origin for Vaccinieae approximately 30 million years ago. Tropical diversity of Vaccinieae was inferred to result from multiple, independent movements into the tropics from north-temperate ancestors. Diversification rate increases corresponded to radiation into the Andes and SE Asia. The pseudo-10-locular ovary evolved once in the tribe from the five-locular state, coinciding with the diversification of a major clade that includes most Asian Vaccinium and the group from which commercial blueberries are derived (V. sect. Cyanococcus). A reconstruction from available chromosome counts suggests that a major polyploid event predated the evolution of nearly half the diversity of Vaccinieae. The extent of polyphyly in Vaccinium documented here supports the need for a generic reclassification of the tribe.

Phyllosticta paracitricarpa is synonymous with the EU quarantine fungus P. citricarpa based on phylogenomic analyses

Phyllosticta citricarpa is an important citrus-pathogen and a quarantine organism in the European Union. Its recently described relative, P. paracitricarpa, is very closely related and not listed as a quarantine organism. P. paracitricarpa is very difficult to distinguish from P. citricarpa, since its morphological features overlap and the barcoding gene sequences that were originally used to delimit them as distinct species have a low number of species-specific polymorphisms that have subsequently been shown to overlap between the two clades. Therefore, we performed extensive genomic analyses to determine whether the genetic variation between P. citricarpa and P. paracitricarpa strains should be considered to represent infraspecific variation within P. citricarpa, or whether it is indicative of distinct species. Using a phylogenomic analysis with 3,000 single copy ortholog genes and whole-genome comparisons, we determined that the variation between P. citricarpa and P. paracitricarpa can be considered as infraspecies variation within P. citricarpa. We also determined the level of variation in mitochondrial assemblies of several Phyllosticta species and concluded there are only minimal differences between the assemblies of P. citricarpa and P. paracitricarpa. Thus, using several orthogonal approaches, we here demonstrate that variation within the nuclear and mitochondrial genomes of other Phyllosticta species is larger than variation between genomes obtained from P. citricarpa and P. paracitricarpa strains. Thus, P. citricarpa and P. paracitricarpa should be considered as conspecific.

Practical Guidance and Workflows for Identifying Fast Evolving Non-Coding Genomic Elements Using PhyloAcc

Comparative genomics provides ample ways to study genome evolution and its relationship to phenotypic traits. By developing and testing alternate models of evolution throughout a phylogeny, one can estimate rates of molecular evolution along different lineages in a phylogeny and link these rates with observations in extant species, such as convergent phenotypes. Pipelines for such work can help identify when and where genomic changes may be associated with, or possibly influence, phenotypic traits. We recently developed a set of models called PhyloAcc, using a Bayesian framework to estimate rates of nucleotide substitution on different branches of a phylogenetic tree and evaluate their association with pre-defined or estimated phenotypic traits. PhyloAcc-ST and PhyloAcc-GT both allow users to define a priori a set of target lineages and then compare different models to identify loci accelerating in one or more target lineages. Whereas ST considers only one species tree across all input loci, GT considers alternate topologies for every locus. PhyloAcc-C simultaneously models molecular rates and rates of continuous trait evolution, allowing the user to ask whether the two are associated. Here, we describe these models and provide tips and workflows on how to prepare the input data and run PhyloAcc.

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

The Meaning and Measure of Concordance Factors in Phylogenomics

As phylogenomic datasets have grown in size, researchers have developed new ways to measure biological variation and to assess statistical support for specific branches. Larger datasets have more sites and loci and therefore less sampling variance. While we can more accurately measure the mean signal in these datasets, lower sampling variance is often reflected in uniformly high measures of branch support-such as the bootstrap and posterior probability-limiting their utility. Larger datasets have also revealed substantial biological variation in the topologies found across individual loci, such that the single species tree inferred by most phylogenetic methods represents a limited summary of the data for many purposes. In contrast to measures of statistical support, the degree of underlying topological variation among loci should be approximately constant regardless of the size of the dataset. "Concordance factors" (CFs) and similar statistics have therefore become increasingly important tools in phylogenetics. In this review, we explain why CFs should be thought of as descriptors of topological variation rather than as measures of statistical support, and argue that they provide important information about the predictive power of the species tree not contained in measures of support. We review a growing suite of statistics for measuring concordance, compare them in a common framework that reveals their interrelationships, and demonstrate how to calculate them using an example from birds. We also discuss how measures of topological variation might change in the future as we move beyond estimating a single "tree of life" toward estimating the myriad evolutionary histories underlying genomic variation.

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

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

Phylogenetic analyses show the Select Agent Coniothyrium glycines represents a single species that has significant morphological and genetic variation

Soybean red leaf blotch (RLB), caused by the fungus Coniothyrium glycines, represents a foliar disease of soybean that is thus far restricted to Africa. The fungus is listed as a Select Agent by the Federal Select Agent Program because it could pose a severe threat to plant health were it to establish in the United States. Previous work uncovered tremendous molecular diversity at the internal transcribed spacer region, suggesting that there may be multiple species causing RLB. To determine whether multiple species cause RLB, we reconstructed the phylogeny of C. glycines and taxonomic allies using sequence data from four genes. We included 33 C. glycines isolates collected from six African countries and determined that all isolates form a well-supported, monophyletic lineage. Within this lineage there are at least six well-supported clades that largely correspond to geography, with one clade exclusively composed of isolates from Ethiopia, another exclusively composed of isolates from Uganda, and four composed of isolates from southern Africa. However, we did not detect any concordance for these clades between the four genes, indicating that all isolates included in this analysis are representative of a single species. Isolates in the Ethiopia clade are morphologically distinct from isolates in the other clades, as they produce larger sclerotia and smaller pycnida and more sclerotia in planta. Additionally, ancestral range estimations suggest that the C. glycines lineage emerged in southern Africa. These results show that there is significantly more genetic and morphological diversity than was initially suspected with this high-consequence fungal plant pathogen.

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

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

Phylogenetics and biogeography of the olive family (Oleaceae)

BACKGROUND AND AIMS

Progress in the systematic studies of the olive family (Oleaceae) during the last two decades provides the opportunity to update its backbone phylogeny and to investigate its historical biogeography. We also aimed to understand the factors underlying the disjunct distribution pattern between East Asia and both West Asia and Europe that is found more commonly in this family than in any other woody plant family.

METHODS

Using a sampling of 298 species out of ~750, the largest in a phylogenetic study of Oleaceae thus far, with a set of 36 plastid and nuclear markers, we reconstructed and dated a new phylogenetic tree based on maximum likelihood and Bayesian methods and checked for any reticulation events. We also assessed the relative support of four competing hypotheses [Qinghai-Tibet Plateau uplift (QTP-only hypothesis); climatic fluctuations (climate-only hypothesis); combined effects of QTP uplift and climate (QTP-climate hypothesis); and no effects (null hypothesis)] in explaining these disjunct distributions.

KEY RESULTS

We recovered all tribes and subtribes within Oleaceae as monophyletic, but uncertainty in the position of tribe Forsythieae remains. Based on this dataset, no reticulation event was detected. Our biogeographical analyses support the QTP-climate hypothesis as the likely main explanation for the East-West Eurasian disjunctions in Oleaceae. Our results also show an earlier origin of Oleaceae at ~86 Mya and the role of Tropical Asia as a main source of species dispersals.

CONCLUSION

Our new family-wide and extensive phylogenetic tree highlights both the stable relationships within Oleaceae, including the polyphyly of the genus Chionanthus, and the need for further systematic studies within the largest and most undersampled genera of the family (Chionanthus and Jasminum). Increased sampling will also help to fine-tune biogeographical analyses across spatial scales and geological times.

Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting (ILS) and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and single nucleotide polymorphism-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that ILS, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera)

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

A consequential one-night stand: Episodic historical hybridization leads to mitochondrial takeover in sympatric desert ant-eating spiders

Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.