Delineating species in the speciation continuum: A proposal

Genetic and morphological evidence support the specific status of the endemic Ericaandevalensis (Ericales, Ericaceae)

Assessing the taxonomic status of closely related taxa is crucial in plant systematics and can have important implications for conservation and human plant use. Ericaandevalensis Cabezudo & Rivera is a metallophyte endemic species from highly metal-polluted soils of SW Iberian Peninsula, an area with a mining history going back more than 5,000 years. Ericaandevalensis is closely related to Ericamackayana Bab., a northern Iberian species also present in western Ireland. The status of E.andevalensis as a species or subspecies subordinated to E.mackayana is subject to debate. Here, we assessed the genetic and phenotypic relationship between both species, including the population structure of E.andevalensis. We used high throughput sequencing to determine genome-wide Single Nucleotide Polymorphisms (SNPs), and morphometric analyses from 35 reproductive and vegetative traits. The morphological analysis showed at least eight characters that can discriminate the two species, from which ovary hairiness and the size of leaf glandular hairs were the most informative. Genetic analyses showed that each species formed a monophyletic cluster with full support, separated by an interspecific genetic distance >4-fold higher than intra-specific distance. Population genetic analyses of E.andevalensis shows that populations are highly structured, with the Portuguese one as the most isolated and less variable. These results support the recognition of E.andevalensis as a distinct species with a highly constrained ecological requirements and a narrow geographic distribution, but with a limited gene flow between populations. We discuss the implications of these outcomes in conservation policies and potential uses of E.andevalensis such as decontamination of polluted soils.

Biometric variability of inflorescence and flower traits among ex situ accessions of the neotropical oilseed palm Acrocomia Mart

The oilseed palm genus Acrocomia is suitable for sustainable oil production in South America. The high phenotypic diversity of wild populations poses a challenge for the delimitation of the genus. Comparing the inflorescence architecture, a first-order panicle, and staminate and pistillate flower traits could be a valuable tool in resolving the taxonomic disarray. Thus, this study aims to characterize the differences in the inflorescence architecture and floral structures of three common and economically significant Acrocomia species: A. aculeata, A. totai, and A. intumescens. Biometric traits of the inflorescence architecture and floral structures of various Acrocomia accessions in an ex situ germplasm collection in Brazil were assessed. The unweighted pair group method with arithmetic mean (UPGMA) cluster analysis based on the Gower distance was used to measure dissimilarities between the individual plants of the accessions. To our best knowledge, this study provides the first evidence of the presence of second-order rachillae in the genus Acrocomia. Evaluated traits showed a high level of variation within and between accessions, emphasizing the phenotypic diversity of the genus. The accessions of A. totai were distinguishable from those of the other two species by their inflorescence architecture and flower traits. The dissimilarities between A. aculeata and A. intumescens were not sufficient to differentiate both. In conclusion, the quantitative assessment of inflorescence and floral traits is a valuable tool for taxonomic resolution of the genus.

Phenotypic characterization of cryptic species in the fungal pathogen Histoplasma

Histoplasmosis is an endemic mycosis that often presents as a respiratory infection in immunocompromised patients. Hundreds of thousands of new infections are reported annually around the world. The etiological agent of the disease, Histoplasma, is a dimorphic fungus commonly found in the soil where it grows as mycelia. Humans can become infected by Histoplasma through inhalation of its spores (conidia) or mycelial particles. The fungi transition into the yeast phase in the lungs at 37°C. Once in the lungs, yeast cells reside and proliferate inside alveolar macrophages. Genomic work has revealed that Histoplasma is composed of at least five cryptic phylogenetic species that differ genetically. Three of those lineages have received new names. Here, we evaluated multiple phenotypic characteristics (colony morphology, secreted proteolytic activity, yeast size, and growth rate) of strains from five of the phylogenetic species of Histoplasma to identify phenotypic traits that differentiate between these species: Histoplasma capsulatum sensu stricto, Histoplasma ohiense, Histoplasma mississippiense, Histoplasma suramericanum, and an African lineage. We report diagnostic traits for three species. The other two species can be identified by a combination of traits. Our results suggest that (i) there are significant phenotypic differences among the cryptic species of Histoplasma and (ii) those differences can be used to positively distinguish those species in a clinical setting and for further study of the evolution of this fungal pathogen.IMPORTANCEIdentifying species boundaries is a critical component of evolutionary biology. Genome sequencing and the use of molecular markers have advanced our understanding of the evolutionary history of fungal pathogens, including Histoplasma, and have allowed for the identification of new species. This is especially important in organisms where morphological characteristics have not been detected. In this study, we revised the taxonomic status of the four named species of the genus Histoplasma, H. capsulatum sensu stricto (ss), H. ohiense, H. mississippiense, and H. suramericanum, and propose the use of species-specific phenotypic traits to aid their identification when genome sequencing is not available. These results have implications not only for evolutionary study of Histoplasma but also for clinicians, as the Histoplasma species could determine the outcome of disease and treatment needed.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae)

The Ficus erecta complex, characterized by its morphological diversity and frequent interspecific overlap, shares pollinating fig wasps among several species. This attribute, coupled with its intricate phylogenetic relationships, establishes it as an exemplary model for studying speciation and evolutionary patterns. Extensive researches involving RADseq (Restriction-site associated DNA sequencing), complete chloroplast genome data, and flow cytometry methods were conducted, focusing on phylogenomic analysis, genetic structure, and ploidy detection within the complex. Significantly, the findings exposed a pronounced nuclear-cytoplasmic conflict. This evidence, together with genetic structure analysis, confirmed that hybridization within the complex is a frequent occurrence. The ploidy detection revealed widespread polyploidy, with certain species exhibiting multiple ploidy levels, including 2×, 3×, and 4×. Of particular note, only five species (F. abelii, F. erecta, F. formosana, F. tannoensis and F. vaccinioides) in the complex were proved to be monophyletic. Species such as F. gasparriniana, F. pandurata, and F. stenophylla were found to encompass multiple phylogenetically distinct lineages. This discovery, along with morphological comparisons, suggests a significant underestimation of species diversity within the complex. This study also identified F. tannoensis as an allopolyploid species originating from F. vaccinioide and F. erecta. Considering the integration of morphological, molecular systematics, and cytological evidences, it is proposed that the scope of the F. erecta complex should be expanded to the entire subsect. Frutescentiae. This would redefine the complex as a continuously evolving group comprising at least 33 taxa, characterized by blurred species boundaries, frequent hybridization and polyploidization, and ambiguous genetic differentiation.

Integration of population genetics with oceanographic models reveals strong connectivity among coral reefs across Seychelles

Many countries with tropical reef systems face hard choices preserving coral reefs in the face of climate change on limited budgets. One approach to maximising regional reef resilience is targeting management efforts and resources at reefs that export large numbers of larvae to other reefs. However, this requires reef connectivity to be quantified. To map coral connectivity in the Seychelles reef system we carried out a population genomic study of the Porites lutea species complex using 241 sequenced colonies from multiple islands. To identify oceanographic drivers of this connectivity and quantify variability, we further used a 2 km resolution regional ocean simulation coupled with a larval dispersal model to predict the flow of coral larvae between reef sites. Patterns of admixture and gene flow are broadly supported by model predictions, but the realised connectivity is greater than that predicted from model simulations. Both methods detected a biogeographic dispersal barrier between the Inner and Outer Islands of Seychelles. However, this barrier is permeable and substantial larval transport is possible across Seychelles, particularly for one of two putative species found in our genomic study. The broad agreement between predicted connectivity and observed genetic patterns supports the use of such larval dispersal simulations in reef system management in Seychelles and the wider region.

Phenotypic characterization of cryptic species in the fungal pathogen Histoplasma

Histoplasmosis is an endemic mycosis that often presents as a respiratory infection in immunocompromised patients. Hundreds of thousands of new infections are reported annually around the world. The etiological agent of the disease, Histoplasma, is a dimorphic fungus commonly found in the soil where it grows as mycelia. Humans can become infected by Histoplasma through inhalation of its spores (conidia) or mycelial particles. The fungi transitions into the yeast phase in the lungs at 37°C. Once in the lungs, yeast cells reside and proliferate inside alveolar macrophages. We have previously described that Histoplasma is composed of at least five cryptic species that differ genetically, and assigned new names to the lineages. Here we evaluated multiple phenotypic characteristics of 12 strains from five phylogenetic species of Histoplasma to identify phenotypic traits that differentiate between these species: H. capsulatum sensu stricto, H. ohiense, H. mississippiense, H. suramericanum, and an African lineage. We report diagnostic traits for two species. The other three species can be identified by a combination of traits. Our results suggest that 1) there are significant phenotypic differences among the cryptic species of Histoplasma, and 2) that those differences can be used to positively distinguish those species in a clinical setting and for further study of the evolution of this fungal pathogen.

Genome-wide support for incipient Tula hantavirus species within a single rodent host lineage

Evolutionary divergence of viruses is most commonly driven by co-divergence with their hosts or through isolation of transmission after host shifts. It remains mostly unknown, however, whether divergent phylogenetic clades within named virus species represent functionally equivalent byproducts of high evolutionary rates or rather incipient virus species. Here, we test these alternatives with genomic data from two widespread phylogenetic clades in Tula orthohantavirus (TULV) within a single evolutionary lineage of their natural rodent host, the common vole Microtus arvalis. We examined voles from forty-two locations in the contact region between clades for TULV infection by reverse transcription (RT)-PCR. Sequencing yielded twenty-three TULV Central North and twenty-one TULV Central South genomes, which differed by 14.9-18.5 per cent at the nucleotide and 2.2-3.7 per cent at the amino acid (AA) level without evidence of recombination or reassortment between clades. Geographic cline analyses demonstrated an abrupt (<1 km wide) transition between the parapatric TULV clades in continuous landscape. This transition was located within the Central mitochondrial lineage of M. arvalis, and genomic single nucleotide polymorphisms showed gradual mixing of host populations across it. Genomic differentiation of hosts was much weaker across the TULV Central North to South transition than across the nearby hybrid zone between two evolutionary lineages in the host. We suggest that these parapatric TULV clades represent functionally distinct, incipient species, which are likely differently affected by genetic polymorphisms in the host. This highlights the potential of natural viral contact zones as systems for investigating the genetic and evolutionary factors enabling or restricting the transmission of RNA viruses.

DNA Barcoding in Species Delimitation: From Genetic Distances to Integrative Taxonomy

Over the past two decades, DNA barcoding has become the most popular exploration approach in molecular taxonomy, whether for identification, discovery, delimitation, or description of species. The present contribution focuses on the utility of DNA barcoding for taxonomic research activities related to species delimitation, emphasizing the following aspects:(1) To what extent DNA barcoding can be a valuable ally for fundamental taxonomic research, (2) its methodological and theoretical limitations, (3) the conceptual background and practical use of pairwise distances between DNA barcode sequences in taxonomy, and (4) the different ways in which DNA barcoding can be combined with complementary means of investigation within a broader integrative framework. In this chapter, we recall and discuss the key conceptual advances that have led to the so-called renaissance of taxonomy, elaborate a detailed glossary for the terms specific to this discipline (see Glossary in Chap. 35 ), and propose a newly designed step-by-step species delimitation protocol starting from DNA barcode data that includes steps from the preliminary elaboration of an optimal sampling strategy to the final decision-making process which potentially leads to nomenclatural changes.

A Glossary of DNA Barcoding Terms

This chapter provides a reference glossary for the protocols in this volume. We have chosen only the very basic terms in the DNA barcode lexicon to include, and provide clear and concise definitions of these terms. We hope the reader finds this glossary useful.

Reconstruction of the last bacterial common ancestor from 183 pangenomes reveals a versatile ancient core genome

BACKGROUND

Cumulative sequencing efforts have yielded enough genomes to construct pangenomes for dozens of bacterial species and elucidate intraspecies gene conservation. Given the diversity of organisms for which this is achievable, similar analyses for ancestral species are feasible through the integration of pangenomics and phylogenetics, promising deeper insights into the nature of ancient life.

RESULTS

We construct pangenomes for 183 bacterial species from 54,085 genomes and identify their core genomes using a novel statistical model to estimate genome-specific error rates and underlying gene frequencies. The core genomes are then integrated into a phylogenetic tree to reconstruct the core genome of the last bacterial common ancestor (LBCA), yielding three main results: First, the gene content of modern and ancestral core genomes are diverse at the level of individual genes but are similarly distributed by functional category and share several poorly characterized genes. Second, the LBCA core genome is distinct from any individual modern core genome but has many fundamental biological systems intact, especially those involving translation machinery and biosynthetic pathways to all major nucleotides and amino acids. Third, despite this metabolic versatility, the LBCA core genome likely requires additional non-core genes for viability, based on comparisons with the minimal organism, JCVI-Syn3A.

CONCLUSIONS

These results suggest that many cellular systems commonly conserved in modern bacteria were not just present in ancient bacteria but were nearly immutable with respect to short-term intraspecies variation. Extending this analysis to other domains of life will likely provide similar insights into more distant ancestral species.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Population Genomics for Insect Conservation

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Molecular Weevil Identification Project: A thoroughly curated barcode release of 1300 Western Palearctic weevil species (Coleoptera, Curculionoidea)

The Molecular Weevil Identification project (MWI) studies the systematics of Western Palearctic weevils (superfamily Curculionoidea) in an integrative taxonomic approach of DNA barcoding, morphology and ecology. This barcode release provides almost 3600 curated CO1 sequences linked to morphological vouchers in about 1300 weevil species. The dataset is presented in statistical distance tables and as a Neighbour-Joining tree. Bayesian Inference trees are computed for the subfamilies Cryptorhynchinae, Apioninae and Ceutorhynchinae. Altogether, 18 unresolved taxonomic issues are discussed. A new barcode primer set is presented. Finally, we establish group-specific genetic distances for many weevil genera to serve as a tool in species delineation. These values are statistically based on distances between "good species" and their congeners. With this morphologically calibrated approach, we could resolve most alpha-taxonomic questions within the MWI project.

Range‐wide population genomics of the spongy moth, Lymantria dispar (Erebidae): Implications for biosurveillance, subspecies classification and phylogeography of a destructive moth

The spongy moth, Lymantria dispar, is an irruptive forest pest native to Eurasia where its range extends from coast to coast and overspills into northern Africa. Accidentally introduced from Europe in Massachusetts in 1868-1869, it is now established in North America where it is considered a highly destructive invasive pest. A fine-scale characterization of its population genetic structure would facilitate identification of source populations for specimens intercepted during ship inspections in North America and would enable mapping of introduction pathways to help prevent future incursions into novel environments. In addition, detailed knowledge of L. dispar's global population structure would provide new insight into the adequacy of its current subspecies classification system and its phylogeographic history. To address these issues, we generated >2000 genotyping-by-sequencing-derived SNPs from 1445 contemporary specimens sampled at 65 locations in 25 countries/3 continents. Using multiple analytical approaches, we identified eight subpopulations that could be further partitioned into 28 groups, achieving unprecedented resolution for this species' population structure. Although reconciliation between these groupings and the three currently recognized subspecies proved to be challenging, our genetic data confirmed circumscription of the japonica subspecies to Japan. However, the genetic cline observed across continental Eurasia, from L. dispar asiatica in East Asia to L. d. dispar in Western Europe, points to the absence of a sharp geographical boundary (e.g., the Ural Mountains) between these two subspecies, as suggested earlier. Importantly, moths from North America and the Caucasus/Middle East displayed high enough genetic distances from other populations to warrant their consideration as separate subspecies of L. dispar. Finally, in contrast with earlier mtDNA-based investigations that identified the Caucasus as L. dispar's place of origin, our analyses suggest continental East Asia as its evolutionary cradle, from where it spread to Central Asia and Europe, and to Japan through Korea.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Australian Bryobia mites (Trombidiformes: Tetranychidae) form a complex of cryptic taxa with unique climatic niches and insecticide responses

BACKGROUND

Bryobia (Koch) mites belong to the economically important spider mite family, the Tetranychidae, with >130 species described worldwide. Due to taxonomic difficulties and most species being asexual, species identification relies heavily on genetic markers. Multiple putative Bryobia mite species have been identified attacking pastures and grain crops in Australia. In this study, we collected 79 field populations of Bryobia mites and combined these with 134 populations that were collected previously. We characterised taxonomic variation of mites using 28S rDNA amplicon-based DNA metabarcoding using next-generation sequencing approaches and direct Sanger sequencing. We then undertook species distribution modelling of the main genetic lineages and examined the chemical responses of multiple field populations.

RESULTS

We identified 47 unique haplotypes across all mites sampled that grouped into four distinct genetic lineages. These lineages have different distributions, with three of the four putative lineages showing different climatic envelopes, as inferred from species distribution modelling. Bryobia mite populations also showed different responses to a widely used insecticide (the organophosphate, omethoate), but not to another chemical (the pyrethroid, bifenthrin) when examined using laboratory bioassays.

CONCLUSION

Our findings indicate that cryptic diversity is likely to complicate the formulation of management strategies for Bryobia mites. Although focussed on Australia, this study demonstrates the challenges of studying Bryobia and highlights the importance of further research into this complex group of mites across the world. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Ten years of demographic modelling of divergence and speciation in the sea

Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the 'grey zone of speciation' for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.

The Diversity in the Genus Canis Challenges Conservation Biology: A Review of Available Data on Asian Wolves

Taxa belonging to the Genus Canis can challenge taxonomists because species boundaries and distribution ranges are often gradual. Species delineation within Canis is currently not based on consistent criteria, and is hampered by geographical bias and lack of taxonomic research. But a consistent taxonomy is critical, given its importance for assigning legal protection, conservation priorities, and financial resources. We carried out a qualitative review of the major wolf lineages so far identified from Asia from historical to contemporary time and considered relevant morphological, ecological, and genetic evidence. We present full mitochondrial phylogenies and genetic distances between these lineages. This review aims to summarize the available data on contemporary Asian wolf lineages within the context of the larger phylogenetic Canis group and to work toward a taxonomy that is consistent within the Canidae. We found support for the presence and taxon eligibility of Holarctic gray, Himalayan/Tibetan, Indian, and Arabian wolves in Asia and recommend their recognition at the taxonomic levels consistent within the group.

Morphologic, genetic, and biogeographic continua among subspecies hinder the conservation of threatened taxa: the case of Centaurea aspera ssp. scorpiurifolia (Asteraceae)

Subspecies are widely included as conservation units because of their potential to become new species. However, their practical recognition includes variable criteria, such as morphological, genetic, geographic and other differences. Centaurea aspera ssp. scorpiurifolia is a threatened taxon endemic to Andalusia (Spain), which coexists in most of its distribution area with similar taxa. Because of the difficulty to identify it using morphology alone, we aimed to sample all the populations cited as ssp. scorpiurifolia as exhaustively as possible, morphologically characterise them, and analyse their genetic structuring using microsatellites, to better understand difficulties when conserving subspecies. Three different Centaurea species were found which were easily identified. Within C. aspera, two genetic populations and some admixed individuals were observed, one including ssp. scorpiurifolia individuals and the other including individuals identified as subspecies aspera, stenophylla, and scorpiurifolia. A morphological continuum between these two genetic populations and a wide overlapping of their biogeographic distribution were also found. This continuum can affect the conservation of ssp. scorpiurifolia because of potential misidentifications and harmful effects of subspecific hybridization. Misidentifications could be partly overcome by using as many different traits as possible, and conservation priority should be given to populations representative of the ends of this continuum.

Molecular species delimitation refines the taxonomy of native and nonnative physinine snails in North America

Being able to associate an organism with a scientific name is fundamental to our understanding of its conservation status, ecology, and evolutionary history. Gastropods in the subfamily Physinae have been especially troublesome to identify because morphological variation can be unrelated to interspecific differences and there have been widespread introductions of an unknown number of species, which has led to a speculative taxonomy. To resolve uncertainty about species diversity in North America, we targeted an array of single-locus species delimitation methods at publically available specimens and new specimens collected from the Snake River basin, USA to generate species hypotheses, corroborated using nuclear analyses of the newly collected specimens. A total-evidence approach delineated 18 candidate species, revealing cryptic diversity within recognized taxa and a lack of support for other named taxa. Hypotheses regarding certain local endemics were confirmed, as were widespread introductions, including of an undescribed taxon likely belonging to a separate genus in southeastern Idaho for which the closest relatives are in southeast Asia. Overall, single-locus species delimitation was an effective first step toward understanding the diversity and distribution of species in Physinae and to guiding future investigation sampling and analyses of species hypotheses.

Sex chromosomes in meiotic, hemiclonal, clonal and polyploid hybrid vertebrates: along the 'extended speciation continuum'

We review knowledge about the roles of sex chromosomes in vertebrate hybridization and speciation, exploring a gradient of divergences with increasing reproductive isolation (speciation continuum). Under early divergence, well-differentiated sex chromosomes in meiotic hybrids may cause Haldane-effects and introgress less easily than autosomes. Undifferentiated sex chromosomes are more susceptible to introgression and form multiple (or new) sex chromosome systems with hardly predictable dominance hierarchies. Under increased divergence, most vertebrates reach complete intrinsic reproductive isolation. Slightly earlier, some hybrids (linked in 'the extended speciation continuum') exhibit aberrant gametogenesis, leading towards female clonality. This facilitates the evolution of various allodiploid and allopolyploid clonal ('asexual') hybrid vertebrates, where 'asexuality' might be a form of intrinsic reproductive isolation. A comprehensive list of 'asexual' hybrid vertebrates shows that they all evolved from parents with divergences that were greater than at the intraspecific level (K2P-distances of greater than 5-22% based on mtDNA). These 'asexual' taxa inherited genetic sex determination by mostly undifferentiated sex chromosomes. Among the few known sex-determining systems in hybrid 'asexuals', female heterogamety (ZW) occurred about twice as often as male heterogamety (XY). We hypothesize that pre-/meiotic aberrations in all-female ZW-hybrids present Haldane-effects promoting their evolution. Understanding the preconditions to produce various clonal or meiotic allopolyploids appears crucial for insights into the evolution of sex, 'asexuality' and polyploidy. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Species Delimitation and Conservation in Taxonomically Challenging Lineages: The Case of Two Clades of Capurodendron (Sapotaceae) in Madagascar

Capurodendron is the largest endemic genus of plants from Madagascar, with around 76% of its species threatened by deforestation and illegal logging. However, some species are not well circumscribed and many of them remain undescribed, impeding a confident evaluation of their conservation status. Here we focus on taxa delimitation and conservation of two species complexes within Capurodendron: the Arid and Western complexes, each containing undescribed morphologies as well as intermediate specimens alongside well-delimited taxa. To solve these taxonomic issues, we studied 381 specimens morphologically and selected 85 of them to obtain intergenic, intronic, and exonic protein-coding sequences of 794 nuclear genes and 227 microsatellite loci. These data were used to test species limits and putative hybrid patterns using different approaches such as phylogenies, PCA, structure analyses, heterozygosity level, FST, and ABBA-BABA tests. The potential distributions were furthermore estimated for each inferred species. The results show that the Capurodendron Western Complex contains three well-delimited species, C. oblongifolium, C. perrieri, and C. pervillei, the first two hybridizing sporadically with the last and producing morphologies similar to, but genetically distinct from C. pervillei. The Arid Complex shows a more intricate situation, as it contains three species morphologically well-delimited but genetically intermixed. Capurodendron mikeorum nom. prov. is shown to be an undescribed species with a restricted distribution, while C. androyense and C. mandrarense have wider and mostly sympatric distributions. Each of the latter two species contains two major genetic pools, one showing interspecific admixture in areas where both taxa coexist, and the other being less admixed and comprising allopatric populations having fewer contacts with the other species. Only two specimens out of 172 showed clear genetic and morphological signals of recent hybridization, while all the others were morphologically well-delimited, independent of their degree of genetic admixture. Hybridization between Capurodendron androyense and C. microphyllum, the sister species of the Arid Complex, was additionally detected in areas where both species coexist, producing intermediate morphologies. Among the two complexes, species are well-defined morphologically with the exception of seven specimens (1.8%) displaying intermediate patterns and genetic signals compatible with a F1 hybridization. A provisional conservation assessment for each species is provided.

Phylogenomic Assessment of Biodiversity Using a Reference-Based Taxonomy: An Example With Horned Lizards (Phrynosoma)

Phylogenomic investigations of biodiversity facilitate the detection of fine-scale population genetic structure and the demographic histories of species and populations. However, determining whether or not the genetic divergence measured among populations reflects species-level differentiation remains a central challenge in species delimitation. One potential solution is to compare genetic divergence between putative new species with other closely related species, sometimes referred to as a reference-based taxonomy. To be described as a new species, a population should be at least as divergent as other species. Here, we develop a reference-based taxonomy for Horned Lizards (Phrynosoma; 17 species) using phylogenomic data (ddRADseq data) to provide a framework for delimiting species in the Greater Short-horned Lizard species complex (P. hernandesi). Previous species delimitation studies of this species complex have produced conflicting results, with morphological data suggesting that P. hernandesi consists of five species, whereas mitochondrial DNA support anywhere from 1 to 10 + species. To help address this conflict, we first estimated a time-calibrated species tree for P. hernandesi and close relatives using SNP data. These results support the paraphyly of P. hernandesi; we recommend the recognition of two species to promote a taxonomy that is consistent with species monophyly. There is strong evidence for three populations within P. hernandesi, and demographic modeling and admixture analyses suggest that these populations are not reproductively isolated, which is consistent with previous morphological analyses that suggest hybridization could be common. Finally, we characterize the population-species boundary by quantifying levels of genetic divergence for all 18 Phrynosoma species. Genetic divergence measures for western and southern populations of P. hernandesi failed to exceed those of other Phrynosoma species, but the relatively small population size estimated for the northern population causes it to appear as a relatively divergent species. These comparisons underscore the difficulties associated with putting a reference-based approach to species delimitation into practice. Nevertheless, the reference-based approach offers a promising framework for the consistent assessment of biodiversity within clades of organisms with similar life histories and ecological traits.

Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow-up study required to verify their geographic or evolutionary limits. Here, we follow-up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species' range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome-wide divergence (F ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species' range. Within lineages, gene flow was high, with low differentiation (mean F ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non-contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.

Using near-infrared spectroscopy to discriminate closely related species: a case study of neotropical ferns

Identifying plant species requires considerable knowledge and can be difficult without complete specimens. Fourier-transform near-infrared spectroscopy (FT-NIR) is an effective technique for discriminating plant species, especially angiosperms. However, its efficacy has never been tested on ferns. Here we tested the accuracy of FT-NIR at discriminating species of the genus Microgramma. We obtained 16 spectral readings per individual from the adaxial and abaxial surfaces of 100 specimens belonging to 13 species. The analyses included all 1557 spectral variables. We tested different datasets (adaxial + abaxial, adaxial, and abaxial) to compare the correct identification of species through the construction of discriminant models (Linear discriminant analysis and partial least squares discriminant analysis) and cross-validation techniques (leave-one-out, K-fold). All analyses recovered an overall high percentage (> 90%) of correct predictions of specimen identifications for all datasets, regardless of the model or cross-validation used. On average, there was > 95% accuracy when using partial least squares discriminant analysis and both cross-validations. Our results show the high predictive power of FT-NIR at correctly discriminating fern species when using leaves of dried herbarium specimens. The technique is sensitive enough to reflect species delimitation problems and possible hybridization, and it has the potential of helping better delimit and identify fern species.

Flower evolution in the presence of heterospecific gene flow and its contribution to lineage divergence

The success of species depends on their ability to exploit ecological resources in order to optimize their reproduction. However, species are not usually found within single-species ecosystems but in complex communities. Because of their genetic relatedness, closely related lineages tend to cluster within the same ecosystem, rely on the same resources, and be phenotypically similar. In sympatry, they will therefore compete for the same resources and, in the case of flowering plants, exchange their genes through heterospecific pollen transfer. These interactions, nevertheless, pose significant challenges to species co-existence because they can lead to resource limitation and reproductive interference. In such cases, divergent selective pressures on floral traits will favour genotypes that isolate or desynchronize the reproduction of sympatric lineages. The resulting displacement of reproductive characters will, in turn, lead to pre-mating isolation and promote intraspecific divergence, thus initiating or reinforcing the speciation process. In this review, we discuss the current theoretical and empirical knowledge on the influence of heterospecific pollen transfer on flower evolution, highlighting its potential to uncover the ecological and genomic constraints shaping the speciation process.

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDI) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

What Is the Best Lens? Comparing the Resolution Power of Genome-Derived Markers and Standard Barcodes

Fungal species delimitation was traditionally carried out with multicopy ribosomal RNA (rRNA) genes, principally for their ease of amplification. Since the efficacy of these markers has been questioned, single-copy protein-encoding genes have been proposed alone or in combination for Multi-Locus Sequence Typing (MLST). In this context, the role of the many sequences obtained with Next-Generation Sequencing (NGS) techniques, in both genomics and metagenomics, further pushes toward an analysis of the efficacy of NGS-derived markers and of the metrics to evaluate the marker efficacy in discriminating fungal species. This paper aims at proposing MeTRe (Mean Taxonomic Resolution), a novel index that could be used both for measuring marker efficacy and for assessing the actual resolution (i.e., the level of separation) between species obtained with different markers or their combinations. In this paper, we described and then employed this index to compare the efficacy of two rRNAs and four single-copy markers obtained from public databases as both an amplicon-based approach and genome-derived sequences. Two different groups of species were used, one with a pathogenic species of Candida that was characterized by relatively well-separated taxa, whereas the other, comprising some relevant species of the sensu stricto group of the genus Saccharomyces, included close species and interspecific hybrids. The results showed the ability of MeTRe to evaluate marker efficacy in general and genome-derived markers specifically.

DILS: Demographic inferences with linked selection by using ABC

We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single-population or two-population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.

Neotype designation for Thymallus aeliani Valenciennes, 1848 from a museum topotype specimen and its affiliation with Adriatic grayling on the basis of mitochondrial DNA

In 1848, the grayling Thymallus aeliani (Valenciennes) was described from Lake Maggiore, Italy, in the north Adriatic basin. Genetic analyses of the mitochondrial control region showed a unique evolutionary history of grayling inhabiting the rivers of northern Adriatic basin, from the upper reaches of the Po River and its left tributaries in the west to the Soča River in the east, which resulted in the designation of this phylogenetic lineage as Adriatic grayling. Consequently, the name T. aeliani was connected to the Adriatic lineage, re-establishing the validity of this taxon. However, the mitochondrial haplotypes belonging to Adriatic grayling were never compared with the type specimens of T. aeliani, as their whereabouts were unknown. In this study, a neotype for T. aeliani was designated using topotypical specimens stored at the Natural History Museum in Vienna. The neotype (NMW 68027:2 labelled as "Lago Maggiore, Bellotti, 1880") was designated pursuant to the conditions stipulated in Article 75.3 of the International Code of Zoological Nomenclature. Furthermore, the mitochondrial control region of the neotype was compared to haplotypes of the Adriatic lineage and showed high genetic similarity, which therefore connects the species name T. aeliani to the Adriatic grayling. This crucial step in fixing nomenclatural status of this species is very important for its protection and management.

On the way to specificity - Microbiome reflects sponge genetic cluster primarily in highly structured populations

Most animals, including sponges (Porifera), have species-specific microbiomes. Which genetic or environmental factors play major roles structuring the microbial community at the intraspecific level in sponges is, however, largely unknown. In this study, we tested whether geographic location or genetic structure of conspecific sponges influences their microbial assembly. For that, we used three sponge species with different rates of gene flow, and collected samples along their entire distribution range (two from the Mediterranean and one from the Southern Ocean) yielding a total of 393 samples. These three sponge species have been previously analysed by microsatellites or single nucleotide polymorphisms, and here we investigate their microbiomes by amplicon sequencing of the microbial 16S rRNA gene. The sponge Petrosia ficiformis, with highly isolated populations (low gene flow), showed a stronger influence of the host genetic distance on the microbial composition than the spatial distance. Host-specificity was therefore detected at the genotypic level, with individuals belonging to the same host genetic cluster harbouring more similar microbiomes than distant ones. On the contrary, the microbiome of Ircinia fasciculata and Dendrilla antarctica - both with weak population structure (high gene flow) - seemed influenced by location rather than by host genetic distance. Our results suggest that in sponge species with high population structure, the host genetic cluster influence the microbial community more than the geographic location.

Shining a light on species delimitation in the tree genus Engelhardia Leschenault ex Blume (Juglandaceae)

Enhanced efficacy in species delimitation is critically important in biology given the pending biodiversity crisis under global warming and anthropogenic activity. In particular, delineation of traditional classifications in view of the complexity of species requires an integrative approach to effectively define species boundaries, and this is a major focus of systematic biology. Here, we explored species delimitation of Engelhardia in tropical and subtropical Asia. In total, 716 individuals in 71 populations were genotyped using five chloroplast regions, one nuclear DNA region (nrITS), and 11 nuclear simple sequence repeats (nSSR). Phylogenetic trees were constructed and relationships among species were assessed. Molecular analyses were then combined with 14 morphological characteristics of 720 specimens to further explore the species boundaries of Engelhardia. Integrating phylogenetic and morphological clusters provided well-resolved relationships to delineate seven species. The results suggested that: first, that E. fenzelii, E. roxburghiana, E. hainanensis, E. apoensis, and E. serrata are distinct species; second, E. spicata var. spicata, E. spicata var. aceriflora, E. spicata var. colebrookeana, and E. rigida should be combined under E. spicata and treated as a species complex; third, E. serrata var. cambodica should be raised to species level and named E. villosa. We illuminated that bias thresholds determining the cluster number for delimiting species boundaries were substantially reduced when morphological data were incorporated. Our results urge caution when using the concepts of subspecies and varieties in order to prevent confusion, particularly with respect to species delimitation for tropical and subtropical species. In some cases, re-ranking or combining subspecies and/or varieties may enable more accurate species delimitation.

Comparative genomics approach to evolutionary process connectivity

The influence of species life history traits and historical demography on contemporary connectivity is still poorly understood. However, these factors partly determine the evolutionary responses of species to anthropogenic landscape alterations. Genetic connectivity and its evolutionary outcomes depend on a variety of spatially dependent evolutionary processes, such as population structure, local adaptation, genetic admixture, and speciation. Over the last years, population genomic studies have been interrogating these processes with increasing resolution, revealing a large diversity of species responses to spatially structured landscapes. In parallel, multispecies meta-analyses usually based on low-genome coverage data have provided fundamental insights into the ecological determinants of genetic connectivity, such as the influence of key life history traits on population structure. However, comparative studies still lack a thorough integration of macro- and micro-evolutionary scales to fully realize their potential. Here, I present how a comparative genomics framework may provide a deeper understanding of evolutionary process connectivity. This framework relies on coupling the inference of long-term demographic and selective history with an assessment of the contemporary consequences of genetic connectivity. Standardizing this approach across several species occupying the same landscape should help understand how spatial environmental heterogeneity has shaped the diversity of historical and contemporary connectivity patterns in different taxa with contrasted life history traits. I will argue that a reasonable amount of genome sequence data can be sufficient to resolve and connect complex macro- and micro-evolutionary histories. Ultimately, implementing this framework in varied taxonomic groups is expected to improve scientific guidelines for conservation and management policies.

Oak symbolism in the light of genomics

Throughout the Northern Hemisphere, human societies, political systems, and religions have appropriated oaks in symbolic representations. In this review, we explore the possible associations between recent genetic and genomic findings and the symbolic representations of oaks. We first consider the ways in which evolutionary history during the Holocene has tightened links between humans and oaks in Europe, and how this may have led to symbolic representations. We then show how recent findings concerning the structure and evolution of the oak genome have provided additional knowledge about symbolic representations, such as longevity, cohesiveness, and robustness.

Twin introductions by independent invader mussel lineages are both associated with recent admixture with a native congener in Australia

Introduced species can impose profound impacts on the evolution of receiving communities with which they interact. If native and introduced taxa remain reproductively semi-isolated, human-mediated secondary contact may promote genetic exchange across newly created hybrid zones, potentially impacting native genetic diversity and invasive species spread. Here, we investigate the contributions of recent divergence histories and ongoing (post-introduction) gene flow between the invasive marine mussel, Mytilus galloprovincialis, and a morphologically indistinguishable and taxonomically contentious native Australian taxon, Mytilus planulatus. Using transcriptome-wide markers, we demonstrate that two contemporary M. galloprovincialis introductions into south-eastern Australia originate from genetically divergent lineages from its native range in the Mediterranean Sea and Atlantic Europe, where both introductions have led to repeated instances of admixture between introduced and endemic populations. Through increased genome-wide resolution of species relationships, combined with demographic modelling, we validate that mussels sampled in Tasmania are representative of the endemic Australian taxon (M. planulatus), but share strong genetic affinities to M. galloprovincialis. Demographic inferences indicate late-Pleistocene divergence times and historical gene flow between the Tasmanian endemic lineage and northern M. galloprovincialis, suggesting that native and introduced taxa have experienced a period of historical isolation of at least 100,000 years. Our results demonstrate that many genomic loci and sufficient sampling of closely related lineages in both sympatric (e.g. Australian populations) and allopatric (e.g. northern hemisphere Mytilus taxa) ranges are necessary to accurately (a) interpret patterns of intraspecific differentiation and to (b) distinguish contemporary invasive introgression from signatures left by recent divergence histories in high dispersal marine species. More broadly, our study fills a significant gap in systematic knowledge of native Australian biodiversity and sheds light on the intrinsic challenges for invasive species research when native and introduced species boundaries are not well defined.

The genome of the sea anemone Actinia equina (L.): Meiotic toolkit genes and the question of sexual reproduction

The beadlet anemone Actinia equina (L.) (Cnidaria: Anthozoa: Actiniaria: Actiniidae) is one of the most familiar organisms of the North European intertidal zone. Once considered a single, morphologically variable species across northern Europe, it is now recognised as one member of a variable species complex. Previous studies of distribution, aggression, allozymes and mitochondrial DNA suggest that the diversity in form and colour within A. equina may hide still unrecognised species diversity. To empower further study of A. equina population genetics and systematics, we sequenced (PacBio Sequel) the genome of a single A. equina individual to produce a high-quality genome assembly (contig N₅₀ = 492,607 bp, 1485 contigs, number of protein coding genes = 47,671, 97% BUSCO completeness). There is debate as to whether A. equina reproduces solely asexually, since no reliable, consistent evidence of sexual reproduction has been found. To gain further insight, we examined the genome for evidence of a 'meiotic toolkit' - genes believed to be found consistently in sexually reproducing organisms - and demonstrate that the A. equina genome appears not to have this full complement. Additionally, Smudgeplot analysis, coupled with high haplotype diversity, indicates this genome assembly to be of ambiguous ploidy, suggesting that A. equina may not be diploid. The suggested polyploid nature of this species coupled with the deficiency in meiotic toolkit genes, indicates that further field and laboratory studies of this species is warranted to understand how this species reproduces and what role ploidy may play in speciation within this speciose genus.