Phylogenomic species delimitation dramatically reduces species diversity in an Antarctic adaptive radiation

Integrative taxonomy clarifies the evolution of a cryptic primate clade

Global biodiversity is under accelerating threats, and species are succumbing to extinction before being described. Madagascar's biota represents an extreme example of this scenario, with the added complication that much of its endemic biodiversity is cryptic. Here we illustrate best practices for clarifying cryptic diversification processes by presenting an integrative framework that leverages multiple lines of evidence and taxon-informed cut-offs for species delimitation, while placing special emphasis on identifying patterns of isolation by distance. We systematically apply this framework to an entire taxonomically controversial primate clade, the mouse lemurs (genus Microcebus, family Cheirogaleidae). We demonstrate that species diversity has been overestimated primarily due to the interpretation of geographic variation as speciation, potentially biasing inference of the underlying processes of evolutionary diversification. Following a revised classification, we find that crypsis within the genus is best explained by a model of morphological stasis imposed by stabilizing selection and a neutral process of niche diversification. Finally, by clarifying species limits and defining evolutionarily significant units, we provide new conservation priorities, bridging fundamental and applied objectives in a generalizable framework.

A Diverging Species within the Stewartia gemmata (Theaceae) Complex Revealed by RAD-Seq Data

Informed species delimitation is crucial in diverse biological fields; however, it can be problematic for species complexes. Showing a peripatric distribution pattern, Stewartia gemmata and S. acutisepala (the S. gemmata complex) provide us with an opportunity to study species boundaries among taxa undergoing nascent speciation. Here, we generated genomic data from representative individuals across the natural distribution ranges of the S. gemmata complex using restriction site-associated DNA sequencing (RAD-seq). Based on the DNA sequence of assembled loci containing 41,436 single-nucleotide polymorphisms (SNPs) and invariant sites, the phylogenetic analysis suggested strong monophyly of both the S. gemmata complex and S. acutisepala, and the latter was nested within the former. Among S. gemmata individuals, the one sampled from Mt. Tianmu (Zhejiang) showed the closest evolutionary affinity with S. acutisepala (which is endemic to southern Zhejiang). Estimated from 2996 high-quality SNPs, the genetic divergence between S. gemmata and S. acutisepala was relatively low (an Fst of 0.073 on a per-site basis). Nevertheless, we observed a proportion of genomic regions showing relatively high genetic differentiation on a windowed basis. Up to 1037 genomic bins showed an Fst value greater than 0.25, accounting for 8.31% of the total. After SNPs subject to linkage disequilibrium were pruned, the principal component analysis (PCA) showed that S. acutisepala diverged from S. gemmata along the first and the second PCs to some extent. By applying phylogenomic analysis, the present study determines that S. acutisepala is a variety of S. gemmata and is diverging from S. gemmata, providing empirical insights into the nascent speciation within a species complex.

Chromosome-level genome assembly and annotation of the Antarctica whitefin plunderfish Pogonophryne albipinna

The Antarctic whitefin plunderfish Pogonophryne albipinna belongs to the family Artedidraconidae, a key component of Antarctic benthic ecosystems within the order Perciformes and the suborder Notothenioidei. While genome research on P. albipinna using short-read sequencing is available, high-quality genome assembly and annotation employing long-read sequencing have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for P. albipinna, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The resulting genome assembly spans approximately 1.07 Gb, with a longest scaffold measuring 59.39 Mb and an N50 length of 41.76 Mb. Of the 1,111 Hi-C scaffolds, 23 exceeded 10 Mb and were thus classified as chromosome-level. BUSCO completeness was assessed at 95.6%. The assembled genome comprises 50.68% repeat sequences, and a total of 31,128 protein-coding genes were predicted. This study will enhance our understanding of the genomic characteristics of cryonotothenioids and facilitate comparative analyses of their adaptation and evolution in extreme environments.

Unsupervised machine learning for species delimitation, integrative taxonomy, and biodiversity conservation

Integrative taxonomy, combining data from multiple axes of biologically relevant variation, is a major goal of systematics. Ideally, such taxonomies will derive from similarly integrative species-delimitation analyses. Yet, most current methods rely solely or primarily on molecular data, with other layers often incorporated only in a post hoc qualitative or comparative manner. A major limitation is the difficulty of devising quantitative parametric models linking different datasets in a unified ecological and evolutionary framework. Machine Learning (ML) methods offer flexibility in this arena by easily learning high-dimensional associations between observations (e.g., individual specimens) across a wide array of input features (e.g., genetics, geography, environment, and phenotype) to delimit statistically meaningful clusters. Here, I implement an unsupervised method using Self-Organizing (or "Kohonen") Maps (SOMs) for such purposes. Recent extensions called "SuperSOMs" can integrate multiple layers, each of which exerts independent influence on a two-dimensional output grid via empirically estimated weights. The grid cells are then delimited into K distinct units that can be interpreted as species or other entities. I show empirical examples in salamanders (Desmognathus) and snakes (Storeria) with layers representing alleles, space, climate, and traits. Simulations reveal that the SuperSOM approach can detect K = 1, tends not to over-split, reflects contributions from all layers, and limits large layers (e.g., genetic matrices) from overwhelming other datasets, desirable properties addressing major concerns from previous studies. Finally, I suggest that these and similar methods could integrate conservation-relevant layers such as population trends and human encroachment to delimit management units from an explicitly quantitative framework grounded in the ecology and evolution of species limits and boundaries.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.

Evolutionary time and species diversity in aquatic ecosystems worldwide

The latitudinal diversity gradient (LDG) is frequently described as the most dramatic biodiversity pattern on Earth, yet ecologists and biogeographers have failed to reach consensus on its primary cause. A key problem in explaining the LDG involves collinearity between multiple factors that are predicted to affect species richness in the same direction. In terrestrial systems, energy input, geographic area, and evolutionary time for species accumulation tend to covary positively with species richness at the largest spatial scales, such that their individual contributions to the LDG are confounded in global analyses. I review three diversity patterns from marine and freshwater systems that break this collinearity and which may thus provide stronger tests of the influence of time on global richness gradients. Specifically, I contrast biodiversity patterns along oceanic depth gradients, in geologically young versus ancient lakes, and in the north versus south polar marine biomes. I focus primarily on fishes due to greater data availability but synthesize patterns for invertebrates where possible. I find that regional-to-global species richness generally declines with depth in the oceans, despite the great age and stability of the deep-sea biome. Geologically ancient lakes generally do not contain more species than young lakes, and the Antarctic marine biome is not appreciably more species rich than the much younger Arctic marine biome. However, endemism is consistently higher in older systems. Patterns for invertebrate groups are less clear than for fishes and reflect a critical need for primary biodiversity data. In summary, the available data suggest that species richness is either decoupled from or only weakly related to the amount of time for diversification. These results suggest that energy, productivity, or geographic area are the primary drivers of large-scale diversity gradients. To the extent that marine and terrestrial diversity gradients result from similar processes, these examples provide evidence against a primary role for evolutionary time as the cause of the LDG.

Taxonomy of Aspergillus series Versicolores: species reduction and lessons learned about intraspecific variability

Aspergillus series Versicolores members occur in a wide range of environments and substrates such as indoor environments, food, clinical materials, soil, caves, marine or hypersaline ecosystems. The taxonomy of the series has undergone numerous re-arrangements including a drastic reduction in the number of species and subsequent recovery to 17 species in the last decade. The identification to species level is however problematic or impossible in some isolates even using DNA sequencing or MALDI-TOF mass spectrometry indicating a problem in the definition of species boundaries. To revise the species limits, we assembled a large dataset of 518 strains. From these, a total of 213 strains were selected for the final analysis according to their calmodulin (CaM) genotype, substrate and geography. This set was used for phylogenetic analysis based on five loci (benA, CaM, RPB2, Mcm7, Tsr1). Apart from the classical phylogenetic methods, we used multispecies coalescence (MSC) model-based methods, including one multilocus method (STACEY) and five single-locus methods (GMYC, bGMYC, PTP, bPTP, ABGD). Almost all species delimitation methods suggested a broad species concept with only four species consistently supported. We also demonstrated that the currently applied concept of species is not sustainable as there are incongruences between single-gene phylogenies resulting in different species identifications when using different gene regions. Morphological and physiological data showed overall lack of good, taxonomically informative characters, which could be used for identification of such a large number of existing species. The characters expressed either low variability across species or significant intraspecific variability exceeding interspecific variability. Based on the above-mentioned results, we reduce series Versicolores to four species, namely A. versicolor, A. creber, A. sydowii and A. subversicolor, and the remaining species are synonymized with either A. versicolor or A. creber. The revised descriptions of the four accepted species are provided. They can all be identified by any of the five genes used in this study. Despite the large reduction in species number, identification based on phenotypic characters remains challenging, because the variation in phenotypic characters is high and overlapping among species, especially between A. versicolor and A. creber. Similar to the 17 narrowly defined species, the four broadly defined species do not have a specific ecology and are distributed worldwide. We expect that the application of comparable methodology with extensive sampling could lead to a similar reduction in the number of cryptic species in other extensively studied Aspergillus species complexes and other fungal genera. Citation: Sklenář F, Glässnerová K, Jurjević Ž, Houbraken J, Samson RA, Visagie CM, Yilmaz N, Gené J, Cano J, Chen AJ, Nováková A, Yaguchi T, Kolařík M, Hubka V (2022). Taxonomy of Aspergillus series Versicolores: species reduction and lessons learned about intraspecific variability. Studies in Mycology 102 : 53-93. doi: 10.3114/sim.2022.102.02.

Reducing the number of accepted species in Aspergillus series Nigri

The Aspergillus series Nigri contains biotechnologically and medically important species. They can produce hazardous mycotoxins, which is relevant due to the frequent occurrence of these species on foodstuffs and in the indoor environment. The taxonomy of the series has undergone numerous rearrangements, and currently, there are 14 species accepted in the series, most of which are considered cryptic. Species-level identifications are, however, problematic or impossible for many isolates even when using DNA sequencing or MALDI-TOF mass spectrometry, indicating a possible problem in the definition of species limits or the presence of undescribed species diversity. To re-examine the species boundaries, we collected DNA sequences from three phylogenetic markers (benA, CaM and RPB2) for 276 strains from series Nigri and generated 18 new whole-genome sequences. With the three-gene dataset, we employed phylogenetic methods based on the multispecies coalescence model, including four single-locus methods (GMYC, bGMYC, PTP and bPTP) and one multilocus method (STACEY). From a total of 15 methods and their various settings, 11 supported the recognition of only three species corresponding to the three main phylogenetic lineages: A. niger, A. tubingensis and A. brasiliensis. Similarly, recognition of these three species was supported by the GCPSR approach (Genealogical Concordance Phylogenetic Species Recognition) and analysis in DELINEATE software. We also showed that the phylogeny based on benA, CaM and RPB2 is suboptimal and displays significant differences from a phylogeny constructed using 5 752 single-copy orthologous proteins; therefore, the results of the delimitation methods may be subject to a higher than usual level of uncertainty. To overcome this, we randomly selected 200 genes from these genomes and performed ten independent STACEY analyses, each with 20 genes. All analyses supported the recognition of only one species in the A. niger and A. brasiliensis lineages, while one to four species were inconsistently delimited in the A. tubingensis lineage. After considering all of these results and their practical implications, we propose that the revised series Nigri includes six species: A. brasiliensis, A. eucalypticola, A. luchuensis (syn. A. piperis), A. niger (syn. A. vinaceus and A. welwitschiae), A. tubingensis (syn. A. chiangmaiensis, A. costaricensis, A. neoniger and A. pseudopiperis) and A. vadensis. We also showed that the intraspecific genetic variability in the redefined A. niger and A. tubingensis does not deviate from that commonly found in other aspergilli. We supplemented the study with a list of accepted species, synonyms and unresolved names, some of which may threaten the stability of the current taxonomy. Citation: Bian C, Kusuya Y, Sklenář F, D'hooge E, Yaguchi T, Ban S, Visagie CM, Houbraken J, Takahashi H, Hubka V (2022). Reducing the number of accepted species in Aspergillus series Nigri. Studies in Mycology 102: 95-132. doi: 10.3114/sim.2022.102.03.

The mitochondrial genome of the red icefish (Channichthys rugosus) casts doubt on its species status

Antarctic notothenioid fishes are recognised as one of the rare examples of adaptive radiation in the marine system. Withstanding the freezing temperatures of Antarctic waters, these fishes have diversified into over 100 species within no more than 10–20 million years. However, the exact species richness of the radiation remains contested. In the genus Channichthys, between one and nine species are recognised by different authors. To resolve the number of Channichthys species, genetic information would be highly valuable; however, so far, only sequences of a single species, C. rhinoceratus, are available. Here, we present the nearly complete sequence of the mitochondrial genome of C. rugosus, obtained from a formalin-fixed museum specimen sampled in 1974. This sequence differs from the mitochondrial genome of C. rhinoceratus in no more than 27 positions, suggesting that the two species may be synonymous.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

Taxonomic delimitation and molecular identification of clusters within the species Zanthoxylumnitidum (Rutaceae) in China

Zanthoxylumnitidum, known as Liang-Mian-Zhen in China, is a traditional Chinese medicinal plant used to treat traumatic injury, rheumatism, paralysis, toothache, stomach ache, and venomous snake bites. Two varieties of the species have been described and three morphological types have been reported within the original variety. However, taxonomic delimitation and molecular markers for distinguishing these varieties and types within this species remain unknown. Since different populations exhibit varying chemical compositions, easy identification of intraspecific taxa is crucial. We collected 420 individuals from 38 natural populations, 3 samples of standard medicinal material, and 17 folk-medicine samples to perform classification and identification within Zanthoxylumnitidum. Four distinct genetic clusters (A, B, C, and D) were highly supported by the nuclear barcode. Two distinct chloroplast clusters (A1 and A2) were further detected within A, and three others had one-to-one correspondence with the remaining nuclear clusters. Molecular identification showed that the 17 folk samples comprised A1, A2, B, and D, while the 3 standard samples belonged to A2. The internal transcribed spacer (ITS) region and rbcL gene are proposed as barcodes for rapid and accurate identification of the different Liang-Mian-Zhen lineages in China. This study highlights the importance of accurate taxonomic delimitation in combination with rapid and accurate molecular identification of medicinal plants.