Incongruence in the phylogenomics era

Phylogenomic analysis reveals the evolutionary history of Paleartic needle-leaved junipers

Needle-leaved junipers (Juniperus sect. Juniperus, Cupressaceae) are coniferous trees and shrubs with red or blue fleshy cones. They are distributed across Asia, Macaronesia and the Mediterranean Basin, with one species (J. communis) having a circumboreal distribution. Here we aim to resolve the phylogeny of this clade to infer its intricate evolutionary history. To do so, we built a comprehensive, time-calibrated phylogeny using genotyping-by-sequencing (GBS) and combine it with species occurrence using phylogeographic tools. Our results provide solid phylogenetic resolution to propose a new taxonomic classification and a biogeographical history of the section. Specifically, we confirm the monophyly of two groups within J. sect. Juniperus: the Asian (blue-cone) species including the circumboreal J. communis, and the Mediterranean-Macaronesian (red-cone) species. In addition, we provide strong phylogenetic evidence for three distinct species (J. badia, J. conferta, J. lutchuensis) previously considered subspecies or varieties, as well as for the differentiation between the eastern and western Mediterranean lineages of J. macrocarpa. Our findings suggest that the Mediterranean basin was the primary center of diversification for Juniperus sect. Juniperus, followed by an East Asian-Tethyan disjunction resulting from uplifts of the Qinghai-Tibetan Plateau and climatic shifts. The colonization history of Macaronesia by red-cone junipers from the western Mediterranean appears to have taken place independently in two different geological periods: the Miocene (Azores) and the Pliocene (Madeira-Canary Islands). Overall, genomic data and phylogenetic analysis are key to consider a new taxonomic proposal and reconstruct the biogeographical history of the iconic needle-leaved junipers across the Paleartic.

Increased microgastropoda sampling give new insights into the phylogenetic relationships of Littorinoidea (Littorinimorpha)

Littorinoidea is one of the most diverse radiations and the most successful group that evolutionary transitions from marine to terrestrial within Littorinimorpha. With such an unmatched diversity, few phylogenetic investigations have attempted to understand their evolutionary relationships, and existing research has primarily focused on typical intertidal species. To address this gap, we conducted the first phylogenomic analysis of the Littorinoidea, leveraging 35 transcriptomes to investigate their internal relationships. Our analyses revealed significant revisions necessary within the Littorinoidea: 1) Pomatias appears distantly related to Littorinidae, suggesting a potential ancestral origin outside of Littorinoidea, challenging traditional classification. The homology of penial innervation within Littorinoidea warrants reevaluation. 2) Lacuna's placement indicates a close relationship with Naticidae, prompting consideration for its removal from Littorinidae. 3) Based on the current phylogenetic research, Peasiella may belong to a distinct family separate from Littorinidae. 4) Our findings support revising the placement of Pteropods within the Littorinimorpha, which is situated phylogenetically between the families Littorinoidea and Naticoidea. Additionally, we highlight the impact of site heterogeneity and evolutionary rate variation on phylogenetic inference. Our study provides a robust phylogenomic framework for the Littorinoidea, emphasizing the importance of including microgastropoda taxa in molecular phylogenetic reconstructions of gastropod subgroups.

Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

The genus Cortinarius should not (yet) be split

The genus Cortinarius (Agaricales, Basidiomycota) is one of the most species-rich fungal genera, with thousands of species reported. Cortinarius species are important ectomycorrhizal fungi and form associations with many vascular plants globally. Until recently Cortinarius was the single genus of the family Cortinariaceae, despite several attempts to provide a workable, lower-rank hierarchical structure based on subgenera and sections. The first phylogenomic study for this group elevated the old genus Cortinarius to family level and the family was split into ten genera, of which seven were described as new. Here, by careful re-examination of the recently published phylogenomic dataset, we detected extensive gene-tree/species-tree conflicts using both concatenation and multispecies coalescent approaches. Our analyses demonstrate that the Cortinarius phylogeny remains unresolved and the resulting phylogenomic hypotheses suffer from very short and unsupported branches in the backbone. We can confirm monophyly of only four out of ten suggested new genera, leaving uncertain the relationships between each other and the general branching order. Thorough exploration of the tree space demonstrated that the topology on which Cortinarius revised classification relies on does not represent the best phylogenetic hypothesis and should not be used as constrained topology to include additional species. For this reason, we argue that based on available evidence the genus Cortinarius should not (yet) be split. Moreover, considering that phylogenetic uncertainty translates to taxonomic uncertainty, we advise for careful evaluation of phylogenomic datasets before proposing radical taxonomic and nomenclatural changes.

Polyphyly of Boehmeria (Urticaceae) congruent with plastome structural variation

Boehmeria is a taxonomically challenging group within the nettle family (Urticaceae). The polyphyly of the genus has been proposed by previous studies with respect to five genera (Debregeasia, Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus). Extensive homoplasy of morphological characters has made generic delimitation problematic. Previous studies in other plant groups suggest that plastome structural variations have the potential to provide characters useful in reconstructing evolutionary relationships. We aimed to test this across Boehmeria and its allied genera by mapping plastome structural variations onto a resolved strongly supported phylogeny. In doing so, we expanded the sampling of the plastome to include Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus for the first time. The results of our phylogenomic analyses provide strong support for Sarcochlamys as being more closely related to Leucosyke puya than to Boehmeria and for the clustering of Boehmeria s.l. into four subclades. The sizes of the plastomes in Boehmeria s.l. ranged from 142,627 bp to 170,958 bp. The plastomes recovered a typical quadripartite structure comprising 127~146 genes. We observe several obvious structural variations across the taxa such as gene loss and multiple gene duplication, inverted repeat (IR) contraction and wide expansions, and inversions. Moreover, we recover a trend for these variations that the early clades were relatively conserved in evolution, whereas the later diverging clades were variable. We propose that the structural variations documented may be linked to the adaptation of Boehmeria s.l. to a wide range of habitats, from moist broadleaf forests in Asia to xeric shrublands and deserts in Africa. This study confirms that variation in plastome gene loss/duplication, IR contraction/expansion, and inversions can provide evidence useful for the reconstruction of evolutionary relationships.

Phylogenomic perspectives on speciation and reproductive isolation in a North American biodiversity hotspot: an example using California sages (Salvia subgenus Audibertia: Lamiaceae)

BACKGROUND AND AIMS

The California Floristic Province (CA-FP) is the most species-rich region of North America north of Mexico. One of several proposed hypotheses explaining the exceptional diversity of the region is that the CA-FP harbours myriad recently diverged lineages with nascent reproductive barriers. Salvia subgenus Audibertia is a conspicuous element of the CA-FP, with multiple sympatric and compatible species.

METHODS

Using 305 nuclear loci and both organellar genomes, we reconstruct species trees, examine genomic discordance, conduct divergence-time estimation, and analyse contemporaneous patterns of gene flow and mechanical reproductive isolation.

KEY RESULTS

Despite strong genomic discordance, an underlying bifurcating tree is supported. Organellar genomes capture additional introgression events not detected in the nuclear genome. Most interfertility is found within clades, indicating that reproductive barriers arise with increasing genetic divergence. Species are generally not mechanically isolated, suggesting that it is unlikely to be the primary factor leading to reproductive isolation.

CONCLUSIONS

Rapid, recent speciation with some interspecific gene flow in conjunction with the onset of a Mediterranean-like climate is the underlying cause of extant diversity in Salvia subgenus Audibertia. Speciation has largely not been facilitated by gene flow. Its signal in the nuclear genome seems to mostly be erased by backcrossing, but organellar genomes each capture different instances of historical gene flow, probably characteristic of many CA-FP lineages. Mechanical reproductive isolation appears to be only part of a mosaic of factors limiting gene flow.

Discovering Fragile Clades and Causal Sequences in Phylogenomics by Evolutionary Sparse Learning

Phylogenomic analyses of long sequences, consisting of many genes and genomic segments, reconstruct organismal relationships with high statistical confidence. But, inferred relationships can be sensitive to excluding just a few sequences. Currently, there is no direct way to identify fragile relationships and the associated individual gene sequences in species. Here, we introduce novel metrics for gene-species sequence concordance and clade probability derived from evolutionary sparse learning models. We validated these metrics using fungi, plant, and animal phylogenomic datasets, highlighting the ability of the new metrics to pinpoint fragile clades and the sequences responsible. The new approach does not necessitate the investigation of alternative phylogenetic hypotheses, substitution models, or repeated data subset analyses. Our methodology offers a streamlined approach to evaluating major inferred clades and identifying sequences that may distort reconstructed phylogenies using large datasets.

Candida glabrata maintains two Hap1 homologs, Zcf27 and Zcf4, for distinct roles in ergosterol gene regulation to mediate sterol homeostasis under azole and hypoxic conditions

Candida glabrata exhibits innate resistance to azole antifungal drugs but also has the propensity to rapidly develop clinical drug resistance. Azole drugs, which target Erg11, is one of the three major classes of antifungals used to treat Candida infections. Despite their widespread use, the mechanism controlling azole-induced ERG gene expression and drug resistance in C. glabrata has primarily revolved around Upc2 and/or Pdr1. In this study, we determined the function of two zinc cluster transcription factors, Zcf27 and Zcf4, as direct but distinct regulators of ERG genes. Our phylogenetic analysis revealed C. glabrata Zcf27 and Zcf4 as the closest homologs to Saccharomyces cerevisiae Hap1. Hap1 is a known zinc cluster transcription factor in S. cerevisiae in controlling ERG gene expression under aerobic and hypoxic conditions. Interestingly, when we deleted HAP1 or ZCF27 in either S. cerevisiae or C. glabrata, respectively, both deletion strains showed altered susceptibility to azole drugs, whereas the strain deleted for ZCF4 did not exhibit azole susceptibility. We also determined that the increased azole susceptibility in a zcf27Δ strain is attributed to decreased azole-induced expression of ERG genes, resulting in decreased levels of total ergosterol. Surprisingly, Zcf4 protein expression is barely detected under aerobic conditions but is specifically induced under hypoxic conditions. However, under hypoxic conditions, Zcf4 but not Zcf27 was directly required for the repression of ERG genes. This study provides the first demonstration that Zcf27 and Zcf4 have evolved to serve distinct roles allowing C. glabrata to adapt to specific host and environmental conditions.

Genomic incongruence accompanies the evolution of flower symmetry in Eudicots: a case study in the poppy family (Papaveraceae, Ranunculales)

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

Complete Organelle Genome of the Desiccation-Tolerant (DT) Moss Tortula atrovirens and Comparative Analysis of the Pottiaceae Family

Tortula atrovirens (Sm.) Lindb. is an important component of biological soil crusts and possesses an extraordinary tolerance against desiccation in dryland habitats. However, knowledge of the organelle genome of this desiccation-tolerant (DT) moss is still lacking. Here, we assembled the first reported Tortula organelle genome and conducted a comprehensive analysis within the Pottiaceae family. T. atrovirens exhibited the second largest chloroplast genome (129,646 bp) within the Pottiaceae, whereas its mitogenome (105,877 bp) and those of other mosses were smaller in size compared to other land plants. The chloroplast and mitochondrial genomes of T. atrovirens were characterized by the expansion of IR boundaries and the absence of homologous recombination-mediated by large repeats. A total of 57 RNA editing sites were detected through mapping RNA-seq data. Moreover, the gene content and order were highly conserved among the Pottiaceae organelle genomes. Phylogenetic analysis showed that bryophytes are paraphyletic, with their three lineages (hornworts, mosses, and liverworts) and vascular plants forming successive sister clades. Timmiella anomala is clearly separated from the monophyletic Pottiaceae, and T. atrovirens is closely related to Syntrichia filaris within the Pottioideae. In addition, we detected four hypervariable regions for candidate-molecular markers. Our findings provide valuable insights into the organelle genomes of T. atrovirens and the evolutionary relationships within the Pottiaceae family, facilitating future discovery of DT genetic resources from bryophytes.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae)

The evolution of the latitudinal diversity gradient (LDG), characterized by a peak in diversity toward the tropics, has captured significant attention in evolutionary biology and ecology. However, the inverse LDG (i-LDG) mechanism, wherein species richness increases toward the poles, remains inadequately explored. Cycads are among one of the oldest lineages of extant seed plants and have undergone extensive diversification in the tropics. Intriguingly, the extant cycad abundance exhibits an i-LDG pattern, and the underlying causes for this phenomenon remain largely elusive. Here, using 1,843 nuclear genes from a nearly complete sampling, we conducted comprehensive phylogenomic analyses to establish a robust species-level phylogeny for Cycas, the largest genus within cycads. We then reconstructed the spatial-temporal dynamics and integrated global environmental data to evaluate the roles of species ages, diversification rates, contemporary environment, and conservatism to ancestral niches in shaping the i-LDG pattern. We found Cycas experienced decreased diversification rates, coupled with the cooling temperature since its origin in the Eocene from continental Asia. Different regions have distinctively contributed to the formation of i-LDG for Cycas, with the northern hemisphere acting as evolutionary museums and the southern hemisphere serving as cradles. Moreover, water-related climate variables, specifically precipitation seasonality and potential evapotranspiration, were identified as paramount factors constraining Cycas species richness in the rainforest biome near the equator. Notably, the adherence to ancestral monsoonal climates emerges as a critical factor in sustaining the diversity pattern. This study underscores the imperative of integrating both evolutionary and ecological approaches to comprehensively unravel the mechanisms underpinning global biodiversity patterns.

Common Methods for Phylogenetic Tree Construction and Their Implementation in R

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

The complete plastome and phylogenetic analysis of Commelina benghalensis L.1753 (Commelinaceae)

Commelina benghalensis L. 1753, a member of the Commelinaceae family, holds significant medicinal and culinary value. This study represents the first documentation of the sequencing and assembly of the entire plastome of C. benghalensis. The genome spans a total length of 160,663 bp, exhibiting a conventional quadripartite architecture that comprises a large single-copy (LSC) region (87,750 bp), a small single-copy (SSC) region (18,417 bp), and two inverted repeats (IR) regions (both 27,248 bp). In its entirety, the C. benghalensis plastome encompasses 129 genes (with 108 being unique), incorporating 77 individual protein-coding genes, 37 unique tRNA genes, and four unique rRNA genes. Phylogenetic analysis revealed a close resemblance between C. benghalensis and C. communis. The sequencing of this plastome stands to expedite the development of molecular markers and significantly contribute to genetic assays involving this distinctive plant.

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.

Putative genome contamination has minimal impact on the GTDB taxonomy

The Genome Taxonomy Database (GTDB) provides a species to domain classification of publicly available genomes based on average nucleotide identity (ANI) (for species) and a concatenated gene phylogeny normalized by evolutionary rates (for genus to phylum), which has been widely adopted by the scientific community. Here, we use the Genome UNClutterer (GUNC) software to identify putatively contaminated genomes in GTDB release 07-RS207. We found that GUNC reported 35,723 genomes as putatively contaminated, comprising 11.25 % of the 317,542 genomes in GTDB release 07-RS207. To assess the impact of this high level of inferred contamination on the delineation of taxa, we created 'clean' versions of the 34,846 putatively contaminated bacterial genomes by removing the most contaminated half. For each clean half, we re-calculated the ANI and concatenated gene phylogeny and found that only 77 (0.22 %) of the genomes were not consistent with their original classification. We conclude that the delineation of taxa in GTDB is robust to the putative contamination detected by GUNC.

The contours of evolution: In defence of Darwin's tree of life paradigm

Both the concept of a Darwinian tree of life (TOL) and the possibility of its accurate reconstruction have been much criticized. Criticisms mostly revolve around the extensive occurrence of lateral gene transfer (LGT), instances of uptake of complete organisms to become organelles (with the associated subsequent gene transfer to the nucleus), as well as the implications of more subtle aspects of the biological species concept. Here we argue that none of these criticisms are sufficient to abandon the valuable TOL concept and the biological realities it captures. Especially important is the need to conceptually distinguish between organismal trees and gene trees, which necessitates incorporating insights of widely occurring LGT into modern evolutionary theory. We demonstrate that all criticisms, while based on important new findings, do not invalidate the TOL. After considering the implications of these new insights, we find that the contours of evolution are best represented by a TOL.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus

Modern taxonomic classification is often based on phylogenetic analyses of a few molecular markers, although single-gene studies are still common. Here, we leverage genome-scale molecular phylogenetics (phylogenomics) of species and populations to reconstruct evolutionary relationships in a dense data set of 710 fungal genomes from the biomedically and technologically important genus Aspergillus. To do so, we generated a novel set of 1,362 high-quality molecular markers specific for Aspergillus and provided profile Hidden Markov Models for each, facilitating their use by others. Examining the resulting phylogeny helped resolve ongoing taxonomic controversies, identified new ones, and revealed extensive strain misidentification (7.59% of strains were previously misidentified), underscoring the importance of population-level sampling in species classification. These findings were corroborated using the current standard, taxonomically informative loci. These findings suggest that phylogenomics of species and populations can facilitate accurate taxonomic classifications and reconstructions of the Tree of Life.IMPORTANCEIdentification of fungal species relies on the use of molecular markers. Advances in genomic technologies have made it possible to sequence the genome of any fungal strain, making it possible to use genomic data for the accurate assignment of strains to fungal species (and for the discovery of new ones). We examined the usefulness and current limitations of genomic data using a large data set of 710 publicly available genomes from multiple strains and species of the biomedically, agriculturally, and industrially important genus Aspergillus. Our evolutionary genomic analyses revealed that nearly 8% of publicly available Aspergillus genomes are misidentified. Our work highlights the usefulness of genomic data for fungal systematic biology and suggests that systematic genome sequencing of multiple strains, including reference strains (e.g., type strains), of fungal species will be required to reduce misidentification errors in public databases.

Microbial Diversity and Open Questions about the Deep Tree of Life

In this perspective, we explore the transformative impact and inherent limitations of metagenomics and single-cell genomics on our understanding of microbial diversity and their integration into the Tree of Life. We delve into the key challenges associated with incorporating new microbial lineages into the Tree of Life through advanced phylogenomic approaches. Additionally, we shed light on enduring debates surrounding various aspects of the microbial Tree of Life, focusing on recent advances in some of its deepest nodes, such as the roots of bacteria, archaea, and eukaryotes. We also bring forth current limitations in genome recovery and phylogenomic methodology, as well as new avenues of research to uncover additional key microbial lineages and resolve the shape of the Tree of Life.

Persistent homology reveals strong phylogenetic signal in 3D protein structures

Changes that occur in proteins over time provide a phylogenetic signal that can be used to decipher their evolutionary history and the relationships between organisms. Sequence comparison is the most common way to access this phylogenetic signal, while those based on 3D structure comparisons are still in their infancy. In this study, we propose an effective approach based on Persistent Homology Theory (PH) to extract the phylogenetic information contained in protein structures. PH provides efficient and robust algorithms for extracting and comparing geometric features from noisy datasets at different spatial resolutions. PH has a growing number of applications in the life sciences, including the study of proteins (e.g. classification, folding). However, it has never been used to study the phylogenetic signal they may contain. Here, using 518 protein families, representing 22,940 protein sequences and structures, from 10 major taxonomic groups, we show that distances calculated with PH from protein structures correlate strongly with phylogenetic distances calculated from protein sequences, at both small and large evolutionary scales. We test several methods for calculating PH distances and propose some refinements to improve their relevance for addressing evolutionary questions. This work opens up new perspectives in evolutionary biology by proposing an efficient way to access the phylogenetic signal contained in protein structures, as well as future developments of topological analysis in the life sciences.

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Genomic, transcriptomic, and ecological diversity of Penicillium species in cheese rind microbiomes

Although Penicillium molds can have significant impacts on agricultural, industrial, and biomedical systems, the ecological roles of Penicillium species in many microbiomes are not well characterized. Here we utilized a collection of 35 Penicillium strains isolated from cheese rinds to broadly investigate the genomic potential for secondary metabolism in cheese-associated Penicillium species, the impact of Penicillium on bacterial community assembly, and mechanisms of Penicillium-bacteria interactions. Using antiSMASH, we identified 1558 biosynthetic gene clusters, 406 of which were mapped to known pathways, including several mycotoxins and antimicrobial compounds. By measuring bacterial abundance and fungal mRNA expression when culturing representative Penicillium strains with a cheese rind bacterial community, we observed divergent impacts of different Penicillium strains, from strong inhibitors of bacterial growth to those with no impact on bacterial growth or community composition. Through differential mRNA expression analyses, Penicillium strains demonstrated limited differential gene expression in response to the bacterial community. We identified a few shared responses between the eight tested Penicillium strains, primarily upregulation of nutrient metabolic pathways, but we did not identify a conserved fungal response to growth in a multispecies community. These results in tandem suggest high variation among cheese-associated Penicillium species in their ability to shape bacterial community development and highlight important ecological diversity within this iconic genus.

Genetics and Genomics of Infectious Diseases in Key Aquaculture Species

Diseases pose a significant and pressing concern for the sustainable development of the aquaculture sector, particularly as their impact continues to grow due to climatic shifts such as rising water temperatures. While various approaches, ranging from biosecurity measures to vaccines, have been devised to combat infectious diseases, their efficacy is disease and species specific and contingent upon a multitude of factors. The fields of genetics and genomics offer effective tools to control and prevent disease outbreaks in aquatic animal species. In this study, we present the key findings from our recent research, focusing on the genetic resistance to three specific diseases: White Spot Syndrome Virus (WSSV) in white shrimp, Bacterial Necrotic Pancreatitis (BNP) in striped catfish, and skin fluke (a parasitic ailment) in yellowtail kingfish. Our investigations reveal that all three species possess substantial heritable genetic components for disease-resistant traits, indicating their potential responsiveness to artificial selection in genetic improvement programs tailored to combat these diseases. Also, we observed a high genetic association between disease traits and survival rates. Through selective breeding aimed at enhancing resistance to these pathogens, we achieved substantial genetic gains, averaging 10% per generation. These selection programs also contributed positively to the overall production performance and productivity of these species. Although the effects of selection on immunological traits or immune responses were not significant in white shrimp, they yielded favorable results in striped catfish. Furthermore, our genomic analyses, including shallow genome sequencing of pedigreed populations, enriched our understanding of the genomic architecture underlying disease resistance traits. These traits are primarily governed by a polygenic nature, with numerous genes or genetic variants, each with small effects. Leveraging a range of advanced statistical methods, from mixed models to machine and deep learning, we developed prediction models that demonstrated moderate-to-high levels of accuracy in forecasting these disease-related traits. In addition to genomics, our RNA-seq experiments identified several genes that undergo upregulation in response to infection or viral loads within the populations. Preliminary microbiome data, while offering limited predictive accuracy for disease traits in one of our studied species, underscore the potential for combining such data with genome sequence information to enhance predictive power for disease traits in our populations. Lastly, this paper briefly discusses the roles of precision agriculture systems and AI algorithms and outlines the path for future research to expedite the development of disease-resistant genetic lines tailored to our target species. In conclusion, our study underscores the critical role of genetics and genomics in fortifying the aquaculture sector against the threats posed by diseases, paving the way for more sustainable and resilient aquaculture development.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

Patchwork: Alignment-Based Retrieval and Concatenation of Phylogenetic Markers from Genomic Data

Low-coverage whole-genome sequencing (also known as "genome skimming") is becoming an increasingly affordable approach to large-scale phylogenetic analyses. While already routinely used to recover organellar genomes, genome skimming is rather rarely utilized for recovering single-copy nuclear markers. One reason might be that only few tools exist to work with this data type within a phylogenomic context, especially to deal with fragmented genome assemblies. We here present a new software tool called Patchwork for mining phylogenetic markers from highly fragmented short-read assemblies as well as directly from sequence reads. Patchwork is an alignment-based tool that utilizes the sequence aligner DIAMOND and is written in the programming language Julia. Homologous regions are obtained via a sequence similarity search, followed by a "hit stitching" phase, in which adjacent or overlapping regions are merged into a single unit. The novel sliding window algorithm trims away any noncoding regions from the resulting sequence. We demonstrate the utility of Patchwork by recovering near-universal single-copy orthologs within a benchmarking study, and we additionally assess the performance of Patchwork in comparison with other programs. We find that Patchwork allows for accurate retrieval of (putatively) single-copy genes from genome skimming data sets at different sequencing depths with high computational speed, outperforming existing software targeting similar tasks. Patchwork is released under the GNU General Public License version 3. Installation instructions, additional documentation, and the source code itself are all available via GitHub at https://github.com/fethalen/Patchwork.

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.