Improving species delimitation for effective conservation: a case study in the endemic maple-leaf oak (Quercus acerifolia)

Introgression, Phylogeography, and Genomic Species Cohesion in the Eastern North American White Oak Syngameon

Hybridization and interspecific gene flow play a substantial role in the evolution of plant taxa. The eastern North American white oak syngameon, a group of approximately 15 ecologically, morphologically and genomically distinguishable species, has long been recognised as a model system for studying introgressive hybridization in temperate trees. However, the prevalence, genomic context and environmental correlates of introgression in this system remain largely unknown. To assess introgression in the eastern North American white oak syngameon and population structure within the widespread Quercus macrocarpa, we conducted a rangewide survey of Q. macrocarpa and four sympatric eastern North American white oak species. Using a Hyb-Seq approach, we assembled a dataset of 3412 thinned single-nucleotide polymorphisms (SNPs) in 445 enriched target loci including 62 genes putatively associated with various ecological functions, as well as associated intronic regions and some off-target intergenic regions (not associated with the exons). Admixture analysis and hybrid class inference demonstrated species coherence despite hybridization and introgressive gene flow (due to backcrossing of F1s to one or both parents). Additionally, we recovered a genetic structure within Q. macrocarpa associated with latitude. Generalised linear mixed models (GLMMs) indicate that proximity to range edge predicts interspecific admixture, but rates of genetic differentiation do not appear to vary between putative functional gene classes. Our study suggests that gene flow between eastern North American white oak species may not be as rampant as previously assumed and that hybridization is most strongly predicted by proximity to a species' range margin.

An integrative approach clarifies species delimitation and biogeographic history of Debregeasia (Urticaceae)

Integrative data from plastid and nuclear loci are increasingly utilized to resolve species boundaries and phylogenetic relationships within major angiosperm clades. Debregeasia (Urticaceae), an economically important genus, presents challenges in species delimitation due to its overlapping morphological traits and unstable taxonomic assignments. Here, we analyzed 14 morphological traits and generated 12 data matrices from the plastomes and nrDNA using genome skimming from the nine recognized morphospecies to clarify species boundaries and assess barcode performance in Debregeasia. We also used a universal set of 353 nuclear genes to explore reticulate evolution and biogeographic history of Debregeasia. Plastomes of Debregeasia exhibited the typical quadripartite structure with conserved gene content and marginal independent variations in the SC/IR boundary at inter- and intra-specific levels. Three Debregeasia species were non-monophyletic and could not be discerned by any barcode; however, ultra-barcodes identified the remaining six (67%), outperforming standard barcodes (56%). Our phylogenetic analyses placed Debregeasia wallichiana outside the genus and suggested six monophyletic clades in Debregeasia, although the placement between Debregeasia hekouensis and Debregeasia libera varied. There was extensive trait overlap in key morphologically diagnostic characters, with reticulation analysis showing potentially pervasive hybridization, likely influenced by speciation patterns and overlaps between species ranges. We inferred that Debregeasia crown diversification began at ca. 12.82 Ma (95% HPD: 11.54-14.63 Ma) in the mid-Miocene within Australia, followed by vicariance and later long-distance dispersal, mainly out of southern China. Our findings highlight the utility of genomic data with integrative lines of evidence to refine species delimitation and explore evolutionary relationships in complex plant lineages.

Deep learning to capture leaf shape in plant images: Validation by geometric morphometrics

Plant leaves play a pivotal role in automated species identification using deep learning (DL). However, achieving reproducible capture of leaf variation remains challenging due to the inherent "black box" problem of DL models. To evaluate the effectiveness of DL in capturing leaf shape, we used geometric morphometrics (GM), an emerging component of eXplainable Artificial Intelligence (XAI) toolkits. We photographed Ranunculus auricomus leaves directly in situ and after herbarization. From these corresponding leaf images, we automatically extracted DL features using a neural network and digitized leaf shapes using GM. The association between the extracted DL features and GM shapes was then evaluated using dimension reduction and covariation models. DL features facilitated the clustering of leaf images by source populations in both in situ and herbarized leaf image datasets, and certain DL features were significantly associated with biological leaf shape variation as inferred by GM. DL features also enabled leaf classification into morpho-phylogenomic groups within the intricate R. auricomus species complex. We demonstrated that simple in situ leaf imaging and DL reproducibly captured leaf shape variation at the population level, while combining this approach with GM provided key insights into the shape information extracted from images by computer vision, a necessary prerequisite for reliable automated plant phenotyping.

Divergent ecological selection maintains species boundaries despite gene flow in a rare endemic tree, Quercus acerifolia (maple-leaf oak)

Strong gene flow from outcrossing relatives tends to blur species boundaries, while divergent ecological selection can counteract gene flow. To better understand how these two forces affect the maintenance of species boundaries, we focused on a species complex including a rare species, maple-leaf oak (Quercus acerifolia), which is found in only four disjunct ridges in Arkansas. Its limited range and geographic proximity to co-occurring close relatives create the possibility for genetic swamping. In this study, we gathered genome-wide single nucleotide polymorphisms (SNPs) using restriction-site-associated DNA sequencing (RADseq) from 190 samples of Q. acerifolia and three of its close relatives, Q. shumardii, Q. buckleyi, and Q. rubra. We found that Q. shumardii and Q. acerifolia are reciprocally monophyletic with low support, suggesting incomplete lineage sorting, introgression between Q. shumardii and Q. acerifolia, or both. Analyses that model allele distributions demonstrate that admixture contributes strongly to this pattern. Populations of Q. acerifolia experience gene flow from Q. shumardii and Q. rubra, but we found evidence that divergent selection is likely maintaining species boundaries: 1) ex situ collections of Q. acerifolia have a higher proportion of hybrids compared to the mature trees of the wild populations, suggesting ecological selection against hybrids at the seed/seedling stage; 2) ecological traits co-vary with genomic composition; and 3) Q. acerifolia shows genetic differentiation at loci hypothesized to influence tolerance of radiation, drought, and high temperature. Our findings strongly suggest that in maple-leaf oak, selection results in higher divergence at regions of the genome despite gene flow from close relatives.

Geographic And Taxonomic Occurrence R-based Scrubbing (gatoRs): An R package and workflow for processing biodiversity data

Premise

Digitized biodiversity data offer extensive information; however, obtaining and processing biodiversity data can be daunting. Complexities arise during data cleaning, such as identifying and removing problematic records. To address these issues, we created the R package Geographic And Taxonomic Occurrence R-based Scrubbing (gatoRs).

Methods and Results

The gatoRs workflow includes functions that streamline downloading records from the Global Biodiversity Information Facility (GBIF) and Integrated Digitized Biocollections (iDigBio). We also created functions to clean downloaded specimen records. Unlike previous R packages, gatoRs accounts for differences in download structure between GBIF and iDigBio and allows for user control via interactive cleaning steps.

Conclusions

Our pipeline enables the scientific community to process biodiversity data efficiently and is accessible to the R coding novice. We anticipate that gatoRs will be useful for both established and beginning users. Furthermore, we expect our package will facilitate the introduction of biodiversity-related concepts into the classroom via the use of herbarium specimens.

Ex situ conservation of two rare oak species using microsatellite and SNP markers

Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.

Leaf morphometric analysis and potential distribution modelling contribute to taxonomic differentiation in the Quercus microphylla complex

Mexico is a major center of evolutionary radiation for the genus Quercus, with oak species occurring across different habitat types and showing a wide variation in morphology and growth form. Despite representing about 20% of Mexican species, scrub oaks have received little attention and even basic aspects of their taxonomy and geographic distribution remain unresolved. In this study, we analyzed the morphological and climatic niche differentiation of scrub oak populations forming a complex constituted by six named species, Quercus cordifolia, Quercus frutex, Quercus intricata, Quercus microphylla, Quercus repanda, Quercus striatula and a distinct morphotype of Q. striatula identified during field and herbarium work (hereafter named Q. striatula II). Samples were obtained from 35 sites covering the geographic distribution of the complex in northern and central Mexico. Morphological differentiation was analyzed through geometric morphometrics of leaf shape and quantification of trichome traits. Our results indicated the presence of two main morphological groups with geographic concordance. The first was formed by Q. frutex, Q. microphylla, Q. repanda and Q. striatula, distributed in the Trans-Mexican Volcanic Belt, the Sierra Madre Occidental and a little portion of the south of the Mexican Altiplano (MA). The second group consists of Q. cordifola, Q. intricata and Q. striatula II, found in the Sierra Madre Oriental and the MA. Therefore, our evidence supports the distinctness of the Q. striatula II morphotype, indicating the need for a taxonomic revision. Within the two groups, morphological differentiation among taxa varied from very clear to low or inexistent (i.e. Q. microphylla-Q. striatula and Q. cordifolia-Q. striatula II) but niche comparisons revealed significant niche differentiation in all pairwise comparisons, highlighting the relevance of integrative approaches for the taxonomic resolution of complicated groups such as the one studied here.

Achieving zero extinction for land plants

Despite the importance of plants for humans and the threats to their future, plant conservation receives far less support compared with vertebrate conservation. Plants are much cheaper and easier to conserve than are animals, but, although there are no technical reasons why any plant species should become extinct, inadequate funding and the shortage of skilled people has created barriers to their conservation. These barriers include the incomplete inventory, the low proportion of species with conservation status assessments, partial online data accessibility, varied data quality, and insufficient investment in both in and ex situ conservation. Machine learning, citizen science (CS), and new technologies could mitigate these problems, but we need to set national and global targets of zero plant extinction to attract greater support.