Gymnosperms on the EDGE

Localization and activity of lipoxygenase in the ovule of Larix kaempferi (Lamb.) Carr. during female gametophyte maturation

KEY MESSAGE

Lipoxygenase activity and localization vary throughout the development of Larix kaempferi ovules, with the highest enzyme activity observed in ovules at the cellular stage and the most intense immunogold reaction noted at the mature archegonium stage of gametophyte development. Lipoxygenases are a family of oxidoreductases with a significant role in biological systems, widespread in living organisms e.g. mammals, fish, corals, plants, mosses, algae, fungi, yeasts, and bacteria. Lipoxygenase activity in plants leads to the formation of phytooxylipins, i.e. signaling molecules, which play a crucial role in many significant physiological processes such as male and female gametophyte maturation, germination and seedling growth, pathogen resistance, abiotic stress response, fruit ripening, and senescence. The activity and localization of lipoxygenase change during plant growth and development. The localization of lipoxygenase in a developing ovule of Larix kaempferi was analyzed using the immunogold labeling method, and the activity was determined spectrophotometrically with linolenic acid as a substrate. Among the investigated stages, the immunogold reaction was the most intense at the mature archegonium stage in the ovule. Lipoxygenase was found in all parts of the L. kaempferi ovule. The largest number of immunogold particles was detected in the integument cells of all the analyzed stages of ovule development. Only one isoform of lipoxygenase with an optimum at pH 8 was active in the ovules during female gametophyte maturation. The highest enzyme activity was determined at the cellular stage, whereas the mature archegonium stage was characterized by its lowest level, which means that LOX activity in developing ovules of the Japanese larch is not correlated with the number of antibody-labeled molecules of the enzyme.

Linking the spatial and genomic structure of adaptive potential for conservation management: a review

We unified the recent literature with the goal to contribute to the discussion on how genetic diversity might best be conserved. We argue that this decision will be guided by how genomic variation is distributed among manageable populations (i.e., its spatial structure), the degree to which adaptive potential is best predicted by variation across the entire genome or the subset of that variation that is identified as putatively adaptive (i.e., its genomic structure), and whether we are managing species as single entities or as collections of diversifying lineages. The distribution of genetic variation and our ultimate goal will have practical implications for on-the-ground management. If adaptive variation is largely polygenic or responsive to change, its spatial structure might be broadly governed by the forces determining genome-wide variation (linked selection, drift, and gene flow), making measurement and prioritization straightforward. If we are managing species as single entities, then population-level prioritization schemes are possible so as to maximize future pooled genetic variation. We outline one such scheme based on the popular Shapley value from cooperative game theory that considers the relative genetic contribution of a population to an unknown future collection of populations.

Advancing EDGE Zones to identify spatial conservation priorities of tetrapod evolutionary history

The biodiversity crisis is pruning the Tree of Life in a way that threatens billions of years of evolutionary history and there is a need to understand where the greatest losses are predicted to occur. We therefore present threatened evolutionary history mapped for all tetrapod groups and describe patterns of Evolutionarily Distinct and Globally Endangered (EDGE) species. Using a complementarity procedure with uncertainty incorporated for 33,628 species, we identify 25 priority tetrapod EDGE Zones, which are insufficiently protected and disproportionately exposed to high human pressure. Tetrapod EDGE Zones are spread over five continents, 33 countries, and 117 ecoregions. Together, they occupy 0.723% of the world's surface but harbour one-third of the world's threatened evolutionary history and EDGE tetrapod species, half of which is endemic. These EDGE Zones highlight areas of immediate concern for researchers, practitioners, policymakers, and communicators looking to safeguard the tetrapod Tree of Life.

Spatial decoupling of taxon richness, phylogenetic diversity and threat status in the megagenus Erica (Ericaceae)

Estimates of the number of vascular plant species currently under threat of extinction are shockingly high, with the highest extinction rates reported for narrow-range, woody plants, especially in biodiversity hotspots with Mediterranean and tropical climates. The large genus Erica is a prime example, as a large proportion of its 851 species, all shrubs or small trees, are endemic to the Cape Floristic Region (CFR) of South Africa. Almost two hundred are known to be threatened and a further hundred are 'Data Deficient'. We need to target conservation efforts and research to fill the most problematic knowledge gaps. This can be especially challenging in large genera, such as Erica, with numerous threatened species that are closely related. One approach involves combining knowledge of phylogenetic diversity with that of IUCN threat status to identify the most Evolutionarily Distinct and Globally Endangered (EDGE) species. We present an expanded and improved phylogenetic hypothesis for Erica (representing 65% of described species diversity) and combine this with available threat and distribution data to identify species and geographic areas that could be targeted for conservation effort to maximise preservation of phylogenetic diversity (PD). The resulting 39 EDGE taxa include 35 from the CFR. A further 32 high PD, data deficient taxa are mostly from outside the CFR, reflecting the low proportion of assessed taxa outside South Africa. The most taxon-rich areas are found in the south-western CFR. They are not the most phylogenetically diverse, but do include the most threatened PD. These results can be cross-referenced to existing living and seed-banked ex situ collections and used to target new and updated threat assessments and conservation action.

The robustness of phylogenetic diversity indices to extinctions

Phylogenetic diversity indices provide a formal way to apportion evolutionary history amongst living species. Understanding the properties of these measures is key to determining their applicability in conservation biology settings. In this work, we investigate some questions posed in a recent paper by Fischer et al. (Syst Biol 72(3):606-615, 2023). In that paper, it is shown that under certain extinction scenarios, the ranking of the surviving species by their Fair Proportion index scores may be the complete reverse of their ranking beforehand. Our main results here show that this behaviour extends to a large class of phylogenetic diversity indices, including the Equal-Splits index. We also provide a necessary condition for reversals of Fair Proportion rankings to occur on phylogenetic trees whose edge lengths obey the ultrametric constraint. Specific examples of rooted phylogenetic trees displaying these behaviours are given and the impact of our results on the use of phylogenetic diversity indices more generally is discussed.

An assessment of methods to combine evolutionary history and conservation: A case study in the Brazilian campo rupestre

Premise

Conservation policies typically focus on biodiversity hotspots. An alternative approach involves analyzing the evolutionary history of lineages in geographic areas along with their threat levels to guide conservation efforts. Mountains exhibit high levels of plant species richness and micro-endemism, and biogeographic studies commonly point to recent and rapid evolutionary radiations in these areas. Using a nearly endemic clade of legumes, our study evaluates conservation prioritization approaches in the campo rupestre, a Neotropical ecosystem associated with mountaintops that is located between two biodiversity hotspots.

Methods

We compared the EDGE and EDGE2 metrics, which combine the evolutionary distinctiveness and the extinction risk of a species in a single value. These metrics are compared with traditional metrics used to assess conservation priority, such as phylogenetic diversity.

Results

The EDGE values reported are lower than those of other studies using this metric, mostly due to the prevalence of threatened species with short phylogenetic branch lengths (low values of evolutionary distinctiveness). Certain areas of campo rupestre with relatively high phylogenetic diversity and EDGE values do not correspond to areas with high species richness, agreeing with previous studies on biodiversity hotspots.

Discussion

Our study highlights the necessity of conservation of the campo rupestres as well as advantages and disadvantages of using EDGE, EDGE2, and phylogenetic diversity for appropriate selection of conservation areas with rapid evolutionary radiations. The selection of the metrics will depend primarily on the life history of the focus group and the data availability, as well as the conservation approach.

Extinction risk predictions for the world's flowering plants to support their conservation

More than 70% of all vascular plants lack conservation status assessments. We aimed to address this shortfall in knowledge of species extinction risk by using the World Checklist of Vascular Plants to generate the first comprehensive set of predictions for a large clade: angiosperms (flowering plants, c. 330 000 species). We used Bayesian Additive Regression Trees (BART) to predict the extinction risk of all angiosperms using predictors relating to range size, human footprint, climate, and evolutionary history and applied a novel approach to estimate uncertainty of individual species-level predictions. From our model predictions, we estimate 45.1% of angiosperm species are potentially threatened with a lower bound of 44.5% and upper bound of 45.7%. Our species-level predictions, with associated uncertainty estimates, do not replace full global, or regional Red List assessments, but can be used to prioritise predicted threatened species for full Red List assessment and fast-track predicted non-threatened species for Least Concern assessments. Our predictions and uncertainty estimates can also guide fieldwork, inform systematic conservation planning and support global plant conservation efforts and targets.

Global conservation status of the jawed vertebrate Tree of Life

Human-driven extinction threatens entire lineages across the Tree of Life. Here we assess the conservation status of jawed vertebrate evolutionary history, using three policy-relevant approaches. First, we calculate an index of threat to overall evolutionary history, showing that we expect to lose 86-150 billion years (11-19%) of jawed vertebrate evolutionary history over the next 50-500 years. Second, we rank jawed vertebrate species by their EDGE scores to identify the highest priorities for species-focused conservation of evolutionary history, finding that chondrichthyans, ray-finned fish and testudines rank highest of all jawed vertebrates. Third, we assess the conservation status of jawed vertebrate families. We found that species within monotypic families are more likely to be threatened and more likely to be in decline than other species. We provide a baseline for the status of families at risk of extinction to catalyse conservation action. This work continues a trend of highlighting neglected groups-such as testudines, crocodylians, amphibians and chondrichthyans-as conservation priorities from a phylogenetic perspective.

Conserving the evolutionary history of birds

In the midst of the sixth mass extinction, limited resources are forcing conservationists to prioritize which species and places will receive conservation action. Evolutionary distinctiveness measures the isolation of a species on its phylogenetic tree. Combining a species' evolutionary distinctiveness with its globally endangered status creates an EDGE score. We use EDGE scores to prioritize the places and species that should be managed to conserve bird evolutionary history. We analyzed all birds in all countries and important bird areas. We examined parrots, raptors, and seabirds in depth because these groups are especially threatened and relatively speciose. The three focal groups had greater median threatened evolutionary history than other taxa, making them important for conserving bird evolutionary history. Australia, Brazil, Indonesia, Madagascar, New Zealand, and the Philippines were especially critical countries for bird conservation because they had the most threatened evolutionary history for endemic birds and are important for parrots, raptors, and seabirds. Increased enforcement of international agreements for the conservation of parrots, raptors, and seabirds is needed because these agreements protect hundreds of millions of years of threatened bird evolutionary history. Decisive action is required to conserve the evolutionary history of birds into the Anthropocene.

Indicators to monitor the status of the tree of life

Following the failure to fully achieve any of the 20 Aichi biodiversity targets, the future of biodiversity rests in the balance. The Convention on Biological Diversity's Kunming-Montreal Global Biodiversity Framework (GBF) presents the opportunity to preserve nature's contributions to people (NCPs) for current and future generations by conserving biodiversity and averting extinctions. There is a need to safeguard the tree of life-the unique and shared evolutionary history of life on Earth-to maintain the benefits it bestows into the future. Two indicators have been adopted within the GBF to monitor progress toward safeguarding the tree of life: the phylogenetic diversity (PD) indicator and the evolutionarily distinct and globally endangered (EDGE) index. We applied both to the world's mammals, birds, and cycads to show their utility at the global and national scale. The PD indicator can be used to monitor the overall conservation status of large parts of the evolutionary tree of life, a measure of biodiversity's capacity to maintain NCPs for future generations. The EDGE index is used to monitor the performance of efforts to conserve the most distinctive species. The risk to PD of birds, cycads, and mammals increased, and mammals exhibited the greatest relative increase in threatened PD over time. These trends appeared robust to the choice of extinction risk weighting. EDGE species had predominantly worsening extinction risk. A greater proportion of EDGE mammals (12%) had increased extinction risk compared with threatened mammals in general (7%). By strengthening commitments to safeguarding the tree of life, biodiversity loss can be reduced and thus nature's capacity to provide benefits to humanity now and in the future can be preserved.

Re-evaluating the importance of threatened species in maintaining global phytoregions

Anthropogenic introductions are known to be changing the structure of global phytogeographical regions (phytoregions), but previous studies have been limited by incomplete or biased data sets that are likely to underestimate the importance of threatened species. In this work, we analyse a comprehensive data set of all known species and their occurrences (at botanical country resolution) to quantify the impact of potential future extinction scenarios. We used Infomap, a network-based community detection algorithm, to generate phytoregional delineations for six species-distribution scenarios (native, introduced and extinctions of species that are either documented as threatened or likely to be threatened, as well as combinations thereof). We compared the numbers and sizes of phytoregions to characterise the amount and spatial distribution of changes in global phytoregions under each scenario. Extinctions of species that are predicted to be threatened had a greater homogenising effect on phytoregions than introductions, and there was some evidence that introductions may even mitigate the homogenisation caused by extinctions, though this interaction is complex. This research provides the first evidence that the loss of threatened species would have significant ramifications for global phytoregions and demonstrates the need to consider extinction processes in studies of anthropogenic effects on biodiversity patterns.

rWCVP: a companion R package for the World Checklist of Vascular Plants

The World Checklist of Vascular Plants (WCVP) is an extremely valuable resource that is being used to address many fundamental and applied questions in plant science, conservation, ecology and evolution. However, databases of this size require data manipulation skills that pose a barrier to many potential users. Here, we present rWCVP, an open-source R package that aims to facilitate the use of the WCVP by providing clear, intuitive functions to execute many common tasks. These functions include taxonomic name reconciliation, geospatial integration, mapping and generation of multiple different summaries of the WCVP in both data and report format. We have included extensive documentation and tutorials, providing step-by-step guides that are accessible even to users with minimal programming experience. rWCVP is available on cran and GitHub.

Conserving avian evolutionary history can effectively safeguard future benefits for people

Phylogenetic diversity (PD)-the evolutionary history of a set of species-is conceptually linked to the maintenance of yet-to-be-discovered benefits from biodiversity or "option value." We used global phylogenetic and utilization data for birds to test the PD option value link, under the assumption that the performance of sets of PD-maximizing species at capturing known benefits is analogous to selecting the same species at a point in human history before these benefits were realized. PD performed better than random at capturing utilized bird species across 60% of tests, with performance linked to the phylogenetic dispersion and prevalence of each utilization category. Prioritizing threatened species for conservation by the PD they encapsulate performs comparably to prioritizing by their functional distinctiveness. However, species selected by each metric show low overlap, indicating that we should conserve both components of biodiversity to effectively conserve a variety of uses. Our findings provide empirical support for the link between evolutionary history and benefits for future generations.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

Genetic diversity, population structure, and genome-wide association analysis of ginkgo cultivars

Ginkgo biloba is an economically valuable tree worldwide. The species has nearly become extinct during the Quaternary, which has likely resulted in reduction of its genetic variability. The genetic variability is now conserved in few natural populations in China and a number of cultivars that are, however, derived from a few ancient trees, helping the species survive in China through medieval times. Despite the recent interest in ginkgo, however, detailed knowledge of its genetic diversity, conserved in cultivated trees and cultivars, has remained poor. This limits efficient conservation of its diversity as well as efficient use of the existing germplasm resources. Here we performed genotyping-by-sequencing (GBS) on 102 cultivated germplasms of ginkgo collected to explore their genetic structure, kinship, and inbreeding prediction. For the first time in ginkgo, a genome-wide association analysis study (GWAS) was used to attempt gene mapping of seed traits. The results showed that most of the germplasms did not show any obvious genetic relationship. The size of the ginkgo germplasm population expanded significantly around 1500 years ago during the Sui and Tang dynasties. Classification of seed cultivars based on a phylogenetic perspective does not support the current classification criteria based on phenotype. Twenty-four candidate genes were localized after performing GWAS on the seed traits. Overall, this study reveals the genetic basis of ginkgo seed traits and provides insights into its cultivation history. These findings will facilitate the conservation and utilization of the domesticated germplasms of this living fossil plant.

Are Mediterranean Island Mountains Hotspots of Taxonomic and Phylogenetic Biodiversity? The Case of the Endemic Flora of the Balearic Islands

The Mediterranean islands are exceptionally rich in endemism, most of which is narrowly distributed. Conservation measures, such as protected areas, have been prioritised, mainly on the basis of species richness and endemism, but phylogenetic information should also be taken into account. In this study, we calculated several taxonomic and phylogenetic metrics at a high resolution for the endemic flora of the Balearic Islands (154 taxa), in order to identify (i) the spatial patterns and environmental factors that explain this endemism, (ii) hotspots of species and phylogenetic endemism, and (iii) gaps in the protected areas. The taxonomic and phylogenetic metrics showed different distribution patterns, but the mountainous areas of Mallorca, and some coastal areas of the Balearic Islands, have the highest values. These values were positively related to elevation, precipitation, temperature, and slope, and negatively related to the distance from the coast, aspect, and the temperature of the wettest quarter. We identified top grid hotspots where all the metrics had the highest values, and we also identified nano-hotspots within these hotspots, in some of the highest peaks of Mallorca, where most of these metrics' maximum values coincided. This approach allowed us to identify some gaps in the conservation priority areas, and to highlight the need to review their boundaries and definition.

The diversity and abundance of bacterial and fungal communities in the rhizosphere of Cathaya argyrophylla are affected by soil physicochemical properties

Cathaya argyrophylla is an ancient Pinaceae species endemic to China that is listed on the IUCN Red List. Although C. argyrophylla is an ectomycorrhizal plant, the relationship between its rhizospheric soil microbial community and soil properties related to the natural habitat remains unknown. High-throughput sequencing of bacterial 16S rRNA genes and fungal ITS region sequences was used to survey the C. argyrophylla soil community at four natural spatially distributed points in Hunan Province, China, and functional profiles were predicted using PICRUSt2 and FUNGuild. The dominant bacterial phyla included Proteobacteria, Acidobacteria, Actinobacteria, and Chloroflexi, and the dominant genus was Acidothermus. The dominant fungal phyla were Basidiomycota and Ascomycota, while Russula was the dominant genus. Soil properties were the main factors leading to changes in rhizosphere soil bacterial and fungal communities, with nitrogen being the main driver of changes in soil microbial communities. The metabolic capacities of the microbial communities were predicted to identify differences in their functional profiles, including amino acid transport and metabolism, energy production and conversion, and the presence of fungi, including saprotrophs and symbiotrophs. These findings illuminate the soil microbial ecology of C. argyrophylla, and provide a scientific basis for screening rhizosphere microorganisms that are suitable for vegetation restoration and reconstruction for this important threatened species.

Cretaceous pollen cone with three-dimensional preservation sheds light on the morphological evolution of cycads in deep time

The Cycadales are an ancient and charismatic group of seed plants. However, their morphological evolution in deep time is poorly understood. While molecular divergence time analyses estimate a Cretaceous origin for most major living cycad clades, much of the extant diversity is inferred to be a result of Neogene diversifications. This leads to long branches throughout the cycadalean phylogeny that, with few exceptions, have yet to be rectified by unequivocal fossil cycads. We report a permineralized pollen cone from the Campanian Holz Shale located in Silverado Canyon, CA, USA (c. 80 million yr ago). This fossil was studied via serial sectioning, SEM, 3D reconstruction and phylogenetic analyses. Microsporophyll and pollen morphology indicate this cone is assignable to Skyttegaardia, a recently described genus based on disarticulated lignitized microsporophylls from the Early Cretaceous of Denmark. Data from this new species, including a simple cone architecture, anatomical details and vasculature organization, indicate cycadalean affinities for Skyttegaardia. Phylogenetic analyses support this assignment and recover Skyttegaardia as crown-group Cycadales, nested within Zamiaceae. Our findings support a Cretaceous diversification for crown-group Zamiaceae, which included the evolution of morphological divergent extinct taxa with unique traits that have yet to be widely identified in the fossil record.

Global conservation prioritization for the Orchidaceae

Quantitative assessments of endemism, evolutionary distinctiveness and extinction threat underpin global conservation prioritization for well-studied taxa, such as birds, mammals, and amphibians. However, such information is unavailable for most of the world's taxa. This is the case for the Orchidaceae, a hyperdiverse and cosmopolitan family with incomplete phylogenetic and threat information. To define conservation priorities, we present a framework based on phylogenetic and taxonomic measures of distinctiveness and rarity based on the number of regions and the area of occupancy. For 25,434 orchid species with distribution data (89.3% of the Orchidaceae), we identify the Neotropics as hotspots for richness, New Guinea as a hotspot for evolutionary distinctiveness, and several islands that contain many rare and distinct species. Orchids have a similar proportion of monotypic genera as other Angiosperms, however, more taxonomically distinct orchid species are found in a single region. We identify 278 species in need of immediate conservation actions and find that more than 70% of these do not currently have an IUCN conservation assessment and are not protected in ex-situ collections at Botanical Gardens. Our study highlights locations and orchid species in urgent need of conservation and demonstrates a framework that can be applied to other data-deficient taxa.

The EDGE2 protocol: Advancing the prioritisation of Evolutionarily Distinct and Globally Endangered species for practical conservation action

The conservation of evolutionary history has been linked to increased benefits for humanity and can be captured by phylogenetic diversity (PD). The Evolutionarily Distinct and Globally Endangered (EDGE) metric has, since 2007, been used to prioritise threatened species for practical conservation that embody large amounts of evolutionary history. While there have been important research advances since 2007, they have not been adopted in practice because of a lack of consensus in the conservation community. Here, building from an interdisciplinary workshop to update the existing EDGE approach, we present an "EDGE2" protocol that draws on a decade of research and innovation to develop an improved, consistent methodology for prioritising species conservation efforts. Key advances include methods for dealing with uncertainty and accounting for the extinction risk of closely related species. We describe EDGE2 in terms of distinct components to facilitate future revisions to its constituent parts without needing to reconsider the whole. We illustrate EDGE2 by applying it to the world's mammals. As we approach a crossroads for global biodiversity policy, this Consensus View shows how collaboration between academic and applied conservation biologists can guide effective and practical priority-setting to conserve biodiversity.

Integrating deep-time palaeontology in conservation prioritisation

Halting biodiversity loss under growing anthropogenic pressure is arguably the greatest environmental challenge we face. Given that not all species are equally threatened and that resources are always limited, establishing robust prioritisation schemes is critical for implementing effective conservation actions. To this end, the International Union for Conservation of Nature (IUCN) Red List of Threatened Species has become a widely used source of information on species’ extinction risk. Various metrics have been proposed that combine IUCN status with different aspects of biodiversity to identify conservation priorities. However, current strategies do not take full advantage of palaeontological data, with conservation palaeobiology often focussing on the near-time fossil record (the last 2 million years). Here, we make a case for the value of the deep-time (over 2 million years ago), as it can offer tangible parallels with today’s biodiversity crisis and inform on the intrinsic traits that make species prone to extinction. As such, palaeontological data holds great predictive power, which could be harnessed to flag species likely to be threatened but that are currently too poorly known to be identified as such. Finally, we identify key IUCN-based prioritisation metrics and outline opportunities for integrating palaeontological data to validate their implementation. Although the human signal of the current extinction crisis makes direct comparisons with the geological past challenging, the deep-time fossil record has more to offer to conservation than is currently recognised.

Correlations between leaf economics, mechanical resistance and drought tolerance across 41 cycad species

BACKGROUND AND AIMS

We conducted a comprehensive analysis of the functional traits of leaves (leaflets) of cycads. The aim of this study was to clarify the functional divergence between the earlier origin Cycadaceae and the later differentiated Zamiaceae, and the differences in trait associations between cycads and angiosperms.

METHODS

We selected 20 Cycadaceae species and 21 Zamiaceae species from the same cycad garden in South China, and measured their leaf structure, economic traits, mechanical resistance (Fp) and leaf water potential at the turgor loss point (πtlp). In addition, we compiled a dataset of geographical distribution along with climatic variables for these cycad species, and some leaf traits of tropical-sub-tropical angiosperm woody species from the literature for comparison.

KEY RESULTS

The results showed significantly contrasting leaf trait syndromes between the two families, with Zamiaceae species exhibiting thicker leaves, higher carbon investments and greater Fp than Cycadaceae species. Leaf thickness (LT) and πtlp were correlated with mean climatic variables in their native distribution ranges, indicating their evolutionary adaptation to environmental conditions. Compared with the leaves of angiosperms, the cycad leaves were thicker and tougher, and more tolerant to desiccation. Greater Fp was associated with a higher structural investment in both angiosperms and cycads; however, cycads showed lower Fp at a given leaf mass per area or LT than angiosperms. Enhancement of Fp led to more negative πtlp in angiosperms, but the opposite trend was observed in cycads.

CONCLUSIONS

Our results reveal that variations in leaf traits of cycads are mainly influenced by taxonomy and the environment of their native range. We also demonstrate similar leaf functional associations in terms of economics, but different relationships with regard to mechanics and drought tolerance between cycads and angiosperms. This study expands our understanding of the ecological strategies and likely responses of cycads to future climate change.

Contrasting Responses of Multispatial Soil Fungal Communities of Thuja sutchuenensis Franch., an Extremely Endangered Conifer in Southwestern China

Thuja sutchuenensis Franch. is an endangered species in southwest China, distributed sporadically in mountainous areas. Soil property and soil fungal community play a crucial role in plant growth and survival. Nevertheless, understanding soil properties and the soil fungal community in the areas where T. sutchuenensis is distributed is extremely limited. Hence, this study collected a total of 180 soil samples from five altitudinal distribution areas (altitudinal gradients) and three vertical depths throughout four horizontal distances from the base of each tree. The results found that altitudinal gradients and vertical depths altered soil properties, including pH, organic matter content, water content, total nitrogen, phosphorus, and potassium, and available nitrogen, phosphorus, and potassium. The fungal alpha diversity indexes (Chao1 and Shannon) and beta diversity were dramatically decreased with elevation. In addition, high altitudes (2,119 m) harbored the highest relative abundance of ectomycorrhizal fungi (27.57%) and the lowest relative abundance of plant-pathogenic fungi (1.81%). Meanwhile, we identified a series of fungal communities, such as Tomentella, Piloderma, Cortinarius, Sebacina, and Boletaceae, that play an essential role in the survival of T. sutchuenensis. The correlation analysis and random forest model identified that water content and total phosphorus showed strong relationships with fungal characteristics and were the primary variables for Zygomycota and Rozellomycota. Collectively, the findings of this integrated analysis provide profound insights into understanding the contrasting responses of T. sutchuenensis soil fungal communities and provide a theoretical basis for T. sutchuenensis habitat restoration and species conservation from multispatial perspectives. IMPORTANCE The present study highlights the importance of fungal communities in an endangered plant, T. sutchuenensis. Comparative analysis of soil samples in nearly all extant T. sutchuenensis populations identified that soil properties, especially soil nutrients, might play critical roles in the survival of T. sutchuenensis. Our findings prove that a series of fungal communities (e.g., Tomentella, Piloderma, and Cortinarius) could be key indicators for T. sutchuenensis survival. In addition, this is the first time that large-scale soil property and fungal community investigations have been carried out in southwest China, offering important values for exploring the distribution pattern of regional soil microorganisms. Collectively, our findings display a holistic picture of soil microbiome and environmental factors associated with T. sutchuenensis.

Mummified Seed Cones of Pinus prehwangshanensis sp. nov. (Subgenus Pinus, Pinaceae) From the Upper Pleistocene of Guangdong, South China: Taxonomical Significance and Implication for Phytogeography and Ecology

Anatomical characters of Cenozoic pine seed cones are known mainly from North American fossils, while data on cone anatomy of Cenozoic species from Asia remain scarce. To date, only one seed cone of Pinus from the Miocene of eastern China has been studied using micro-computed tomography (micro-CT). A new fossil-species, Pinus prehwangshanensis sp. nov., of mummified seed cones from the upper Pleistocene of South China is described using a combination of scanning electron microscopy (SEM) and micro-CT. The new fossil-species combines a mosaic of seed cone morphological and anatomical characters observed in the group of closely related East Asian extant species of subgenus Pinus, section Pinus, subsection Pinus, comprising Pinus taiwanensis, Pinus hwangshanensis, Pinus luchuensis, Pinus thunbergii, and Pinus densiflora. The data obtained indicate that the characteristic anatomical features of this group were formed no later than the end of the Pleistocene. Based on the external seed cone morphology, the East Asian pine fossils confirm the existence of floristic exchange between continental Asia and the Japan archipelago prior to the formation of the Sea of Japan and later, in the middle Miocene to the late Pliocene, when the connection between the Japanese islands and Eurasian continent became re-established. Pollen grains associated with the new fossil-species are similar to those of some extant pine species related to P. thunbergii. A taxonomic and ecological analysis of the Pleistocene plant taxa from the Maoming Basin suggests that the regional climate was a humid subtropical monsoon with hot wet summers and cool dry winters, similar to the present-day climate of northeastern Vietnam.

Fleshy or dry: transcriptome analyses reveal the genetic mechanisms underlying bract development in Ephedra

BACKGROUND

Gnetales have a key phylogenetic position in the evolution of seed plants. Among the Gnetales, there is an extraordinary morphological diversity of seeds, the genus Ephedra, in particular, exhibits fleshy, coriaceous or winged (dry) seeds. Despite this striking diversity, its underlying genetic mechanisms remain poorly understood due to the limited studies in gymnosperms. Expanding the genomic and developmental data from gymnosperms contributes to a better understanding of seed evolution and development.

RESULTS

We performed transcriptome analyses on different plant tissues of two Ephedra species with different seed morphologies. Anatomical observations in early developing ovules, show that differences in the seed morphologies are established early in their development. The transcriptomic analyses in dry-seeded Ephedra californica and fleshy-seeded Ephedra antisyphilitica, allowed us to identify the major differences between the differentially expressed genes in these species. We detected several genes known to be involved in fruit ripening as upregulated in the fleshy seed of Ephedra antisyphilitica.

CONCLUSIONS

This study allowed us to determine the differentially expressed genes involved in seed development of two Ephedra species. Furthermore, the results of this study of seeds with the enigmatic morphology in Ephedra californica and Ephedra antisyphilitica, allowed us to corroborate the hypothesis which suggest that the extra envelopes covering the seeds of Gnetales are not genetically similar to integument. Our results highlight the importance of carrying out studies on less explored species such as gymnosperms, to gain a better understanding of the evolutionary history of plants.

The Common Bean V Gene Encodes Flavonoid 3'5' Hydroxylase: A Major Mutational Target for Flavonoid Diversity in Angiosperms

The classic V (violet, purple) gene of common bean (Phaseolus vulgaris) functions in a complex genetic network that controls seed coat and flower color and flavonoid content. V was cloned to understand its role in the network and the evolution of its orthologs in the Viridiplantae. V mapped genetically to a narrow interval on chromosome Pv06. A candidate gene was selected based on flavonoid analysis and confirmed by recombinational mapping. Protein and domain modeling determined V encodes flavonoid 3'5' hydroxylase (F3'5'H), a P450 enzyme required for the expression of dihydromyricetin-derived flavonoids in the flavonoid pathway. Eight recessive haplotypes, defined by mutations of key functional domains required for P450 activities, evolved independently in the two bean gene pools from a common ancestral gene. V homologs were identified in Viridiplantae orders by functional domain searches. A phylogenetic analysis determined F3'5'H first appeared in the Streptophyta and is present in only 41% of Angiosperm reference genomes. The evolutionarily related flavonoid pathway gene flavonoid 3' hydroxylase (F3'H) is found nearly universally in all Angiosperms. F3'H may be conserved because of its role in abiotic stress, while F3'5'H evolved as a major target gene for the evolution of flower and seed coat color in plants.

The Angiosperm Terrestrial Revolution and the origins of modern biodiversity

Biodiversity today has the unusual property that 85% of plant and animal species live on land rather than in the sea, and half of these live in tropical rainforests. An explosive boost to terrestrial diversity occurred from c. 100-50 million years ago, the Late Cretaceous and early Palaeogene. During this interval, the Earth-life system on land was reset, and the biosphere expanded to a new level of productivity, enhancing the capacity and species diversity of terrestrial environments. This boost in terrestrial biodiversity coincided with innovations in flowering plant biology and evolutionary ecology, including their flowers and efficiencies in reproduction; coevolution with animals, especially pollinators and herbivores; photosynthetic capacities; adaptability; and ability to modify habitats. The rise of angiosperms triggered a macroecological revolution on land and drove modern biodiversity in a secular, prolonged shift to new, high levels, a series of processes we name here the Angiosperm Terrestrial Revolution.

Conserving evolutionarily distinct species is critical to safeguard human well-being

Although there is growing interest in safeguarding the Tree of Life to preserve the human benefits that are directly provided by biodiversity, their evolutionary distribution remains unknown, which has hampered our understanding of the potential of phylodiversity indicators to evince them. Here, I drew on a global review of plant benefits and comprehensive phylogenetic information to breakdown their evolutionary distribution and thereby show why the commonly used Phylogenetic Diversity and Evolutionary Distinctiveness indicators can unequivocally help to preserve these natural services. Beneficial species clumped within phylogenetically overdispersed genera and closely related species often contributed very few and redundant benefits, suggesting that multiple plant lineages are required to maintain a wide variety of services. Yet, a reduced number of species stood out as multi-beneficial and evolutionarily distinct plants relative to both the entire phylogeny and the subset of beneficial species, and they collectively contributed a higher-than-expected number of records for most types of benefits. In addition to providing a clear mechanistic understanding for the recently proved success of Phylogenetic Diversity in capturing plant benefits, these findings stress the decisive role that conservation programmes aimed at protecting evolutionarily distinct taxa will play in safeguarding the beneficial potential of biodiversity for the future.

Taxus wallichiana var. chinensis (Pilg.) Florin Aqueous Extract Suppresses the Proliferation and Metastasis in Lung Carcinoma via JAK/STAT3 Signaling Pathway

As one of the most common neoplasms globally, lung cancer (LC) is the leading cause of cancer-related mortality. Recurrence and metastasis negatively influencing therapeutic efficacy and overall survival demand new strategies in LC treatment. The advantages of TCM are increasingly highlighted. In this study, we obtained the major chemical components and their ratios in the aqueous extract of Taxus wallichiana var. chinensis (Pilg.) Florin (AETW) by UPLC-Q/TOF-MS/MS detection. The CCK-8 assay revealed that AETW could selectively inhibit the growth of A549 and HCC827 cells in a dose-dependent manner with little effect on normal human lung cells. Moreover, both in vitro and in vivo experiments showed that AETW was able to suppress the capacities of cell migration and invasion and downregulate the EMT and the JAK/STAT3 signaling pathway. To further probe into the molecular mechanism, the overexpression of STAT3 was performed into LC cells with AETW treatment, which counteracted the inhibitory effect on malignant behaviors of A549 and HCC827 cells with the decline in the expressions of p-JAK and p-STAT3. Taken together, we propose that AETW may inhibit the proliferation and metastasis by inactivating the JAK/STAT3 axis.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

Deciphering the evolution of the ovule genetic network through expression analyses in Gnetum gnemon

BACKGROUND AND AIMS

The ovule is a synapomorphy of all seed plants (gymnosperms and angiosperms); however, there are some striking differences in ovules among the major seed plant lineages, such as the number of integuments or the orientation of the ovule. The genetics involved in ovule development have been well studied in the model species Arabidopsis thaliana, which has two integuments and anatropous orientation. This study is approached from what is known in arabidopsis, focusing on the expression patterns of homologues of four genes known to be key for the proper development of the integuments in arabidopsis: AINTEGUMENTA (ANT), BELL1, (BEL1), KANADIs (KANs) and UNICORN (UCN).

METHODS

We used histology to describe the morphoanatomical development from ovules to seeds in Gnetum gnemon. We carried out spatiotemporal expression analyses in G. gnemon, a gymnosperm, which has a unique ovule morphology with an integument covering the nucellus, two additional envelopes where the outermost becomes fleshy as the seed matures, and an orthotropous orientation.

KEY RESULTS

Our anatomical and developmental descriptions provide a framework for expression analyses in the ovule of G. gnemon. Our expression results show that although ANT, KAN and UCN homologues are expressed in the inner integument, their spatiotemporal patterns differ from those found in angiosperms. Furthermore, all homologues studied here are expressed in the nucellus, revealing major differences in seed plants. Finally, no expression of the studied homologues was detected in the outer envelopes.

CONCLUSIONS

Altogether, these analyses provide significant comparative data that allows us to better understand the functional evolution of these gene lineages, providing a compelling framework for evolutionary and developmental studies of seeds. Our findings suggest that these genes were most likely recruited from the sporangium development network and became restricted to the integuments of angiosperm ovules.

Arabinogalactan-proteins from non-coniferous gymnosperms have unusual structural features

Arabinogalactan-proteins (AGPs), important signalling molecules of the plant cell wall, are structurally extensively investigated in angiosperms, but information on AGPs in gymnosperms is still limited. We characterized AGPs from the gymnosperms Ginkgo biloba, Ephedra distachya, Encephalartos longifolius and Cycas revoluta. The protein contents are comparable to that of angiosperm AGPs. Hydroxyproline is the site of linking the carbohydrate part and was detected in all AGPs with highest concentration in Cycas AGP (1.1 % of the AGP). Interestingly, with the exception of Cycas, all AGPs contained the monosaccharide 3-O-methylrhamnose not present in angiosperm polysaccharides. The carbohydrate moieties of Cycas and Ephredra showed the main components 1,3,6-linked galactose and terminal arabinose typical of angiosperm AGPs, whereas that of Ginkgo AGP was unique with 1,4-linked galactose as dominant structural element. Bioinformatic search for glycosyltransferases in Ginkgo genome also revealed a lower number of galactosyltransferases responsible for biosynthesis of the 1,3-Gal/1,6-Gal AGP backbone.

Evolutionary analysis of chloroplast tRNA of Gymnosperm revealed the novel structural variation and evolutionary aspect

Gymnosperms such as ginkgo, conifers, cycads, and gnetophytes are vital components of land ecosystems, and they have significant economic and ecologic value, as well as important roles as forest vegetation. In this study, we investigated the structural variation and evolution of chloroplast transfer RNAs (tRNAs) in gymnosperms. Chloroplasts are important organelles in photosynthetic plants. tRNAs are key participants in translation where they act as adapter molecules between the information level of nucleic acids and functional level of proteins. The basic structures of gymnosperm chloroplast tRNAs were found to have family-specific conserved sequences. The tRNAΨ -loop was observed to contain a conforming sequence, i.e., U-U-C-N-A-N2. In gymnosperms, tRNAIle was found to encode a "CAU" anticodon, which is usually encoded by tRNAMet. Phylogenetic analysis suggested that plastid tRNAs have a common polyphyletic evolutionary pattern, i.e., rooted in abundant common ancestors. Analyses of duplication and loss events in chloroplast tRNAs showed that gymnosperm tRNAs have experienced little more gene loss than gene duplication. Transition and transversion analysis showed that the tRNAs are iso-acceptor specific and they have experienced unequal evolutionary rates. These results provide new insights into the structural variation and evolution of gymnosperm chloroplast tRNAs, which may improve our comprehensive understanding of the biological characteristics of the tRNA family.

Energy-water and seasonal variations in climate underlie the spatial distribution patterns of gymnosperm species richness in China

Studying the pattern of species richness is crucial in understanding the diversity and distribution of organisms in the earth. Climate and human influences are the major driving factors that directly influence the large-scale distributions of plant species, including gymnosperms. Understanding how gymnosperms respond to climate, topography, and human-induced changes is useful in predicting the impacts of global change. Here, we attempt to evaluate how climatic and human-induced processes could affect the spatial richness patterns of gymnosperms in China. Initially, we divided a map of the country into grid cells of 50 × 50 km2 spatial resolution and plotted the geographical coordinate distribution occurrence of 236 native gymnosperm taxa. The gymnosperm taxa were separated into three response variables: (a) all species, (b) endemic species, and (c) nonendemic species, based on their distribution. The species richness patterns of these response variables to four predictor sets were also evaluated: (a) energy-water, (b) climatic seasonality, (c) habitat heterogeneity, and (d) human influences. We performed generalized linear models (GLMs) and variation partitioning analyses to determine the effect of predictors on spatial richness patterns. The results showed that the distribution pattern of species richness was highest in the southwestern mountainous area and Taiwan in China. We found a significant relationship between the predictor variable set and species richness pattern. Further, our findings provide evidence that climatic seasonality is the most important factor in explaining distinct fractions of variations in the species richness patterns of all studied response variables. Moreover, it was found that energy-water was the best predictor set to determine the richness pattern of all species and endemic species, while habitat heterogeneity has a better influence on nonendemic species. Therefore, we conclude that with the current climate fluctuations as a result of climate change and increasing human activities, gymnosperms might face a high risk of extinction.

Formal Links between Feature Diversity and Phylogenetic Diversity

The extent to which phylogenetic diversity (PD) captures feature diversity (FD) is a topical and controversial question in biodiversity conservation. In this short paper, we formalize this question and establish a precise mathematical condition for FD (based on discrete characters) to coincide with PD. In this way, we make explicit the two main reasons why the two diversity measures might disagree for given data; namely, the presence of certain patterns of feature evolution and loss, and using temporal branch lengths for PD in settings that may not be appropriate (e.g., due to rapid evolution of certain features over short periods of time). Our article also explores the relationship between the "Fair Proportion" index of PD and a simple index of FD (both of which correspond to Shapley values in cooperative game theory). In a second mathematical result, we show that the two indices can take identical values for any phylogenetic tree, provided the branch lengths in the tree are chosen appropriately. [Evolutionary distinctiveness; feature diversity; phylogenetic diversity; shapley value.].

Steep topography buffers threatened gymnosperm species against anthropogenic pressures in China

China is one of the most species-rich countries in the world, harboring many rare gymnosperms. Following recent human-led loss of forests, China is now experiencing increases in forest cover resulting from efforts of reforestation schemes. As anthropogenic activities have previously been found to interact with topography in shaping forest cover in China and considering the large human population and the ongoing population increase of the country, it is important to understand the role of anthropogenic pressures relative to environmental drivers for shaping species distributions here. Based on the well-established relationship between human population density and topography, we propose a hypothesis for explaining species distributions in a country dominated by human activities, predicting that species are more likely to occur in areas of steep topography under medium human population densities compared to low and high human population densities. Using species occurrence data from the Chinese Vascular Plant Distribution Database along with a common SDM method (maximum entropy modeling), we tested this hypothesis. Our results show that steep topography has the highest importance for predicting Chinese gymnosperm species occurrences in general, and threatened species specifically, in areas of medium human population densities. Consequently, these species are more often found in areas of steep terrain, supporting the proposed hypothesis. Results from this study highlight the need to include topographically heterogeneous habitats when planning new protected areas for species conservation.

Anti-Obesity Effect of Ginkgo Vinegar, a Fermented Product of Ginkgo Seed Coat, in Mice Fed a High-Fat Diet and 3T3-L1 Preadipocyte Cells

Ginkgo seed coat is rarely used and is typically discarded, due to its offensive odor and its toxicity. Ginkgo vinegar is a fermented product of ginkgo seed coat, and fermentation removes the bad smell and most of the toxicity. Thus, ginkgo vinegar contains very low concentrations of toxic components. The present study examined the anti-obesity effect of ginkgo vinegar in mice fed a high-fat diet and its inhibition of adipogenesis in 3T3-L1 cells. Ginkgo vinegar suppressed high-fat diet-induced body weight gain and reduced the size of fat cells in mice. Ginkgo vinegar suppressed the expression of C/EBPδ and PPARγ, key proteins in adipogenesis, and inhibited lipid accumulation in 3T3-L1 cells that were induced to become adipocytes. These results suggested that ginkgo vinegar inhibited adipocyte differentiation. On the other hand, a corresponding concentration of acetic acid had significantly less effect on lipid accumulation and virtually no effect on adipogenic gene expression. These results suggested that, similar to Ginkgo biloba extract, ginkgo vinegar might prevent and improve adiposity. Therefore, ginkgo seed coat could be a useful material for medicinal ingredients.

Reconstructing Krassilovia mongolica supports recognition of a new and unusual group of Mesozoic conifers

Previously unrecognized anatomical features of the cone scales of the enigmatic Early Cretaceous conifer Krassilovia mongolica include the presence of transversely oriented paracytic stomata, which is unusual for all other extinct and extant conifers. Identical stomata are present on co-occurring broad, linear, multiveined leaves assigned to Podozamites harrisii, providing evidence that K. mongolica and P. harrisii are the seed cones and leaves of the same extinct plant. Phylogenetic analyses of the relationships of the reconstructed Krassilovia plant place it in an informal clade that we name the Krassilovia Clade, which also includes Swedenborgia cryptomerioides-Podozamites schenkii, and Cycadocarpidium erdmanni-Podozamites schenkii. All three of these plants have linear leaves that are relatively broad compared to most living conifers, and that are also multiveined with transversely oriented paracytic stomata. We propose that these may be general features of the Krassilovia Clade. Paracytic stomata, and other features of this new group, recall features of extant and fossil Gnetales, raising questions about the phylogenetic homogeneity of the conifer clade similar to those raised by phylogenetic analyses of molecular data.

Cold spot microrefugia hold the key to survival for Brazil's Critically Endangered Araucaria tree

Brazil's Araucaria tree (Araucaria angustifolia) is an iconic living fossil and a defining element of the Atlantic Forest global biodiversity hotspot. But despite more than two millennia as a cultural icon in southern Brazil, Araucaria is on the brink of extinction, having lost 97% of its extent to 20th-century logging. Although logging is now illegal, 21st-century climate change constitutes a new-but so far unevaluated-threat to Araucaria's future survival. We use a robust ensemble modelling approach, using recently developed climate data, high-resolution topography and fine-scale vegetation maps, to predict the species' response to climate change and its implications for conservation on meso- and microclimate scales. We show that climate-only models predict the total disappearance of Araucaria's most suitable habitat by 2070, but incorporating topographic effects allows potential highland microrefugia to be identified. The legacy of 20th-century destruction is evident-more than a third of these likely holdouts have already lost their natural vegetation-and 21st-century climate change will leave just 3.5% of remnant forest and 28.4% of highland grasslands suitable for Araucaria. Existing protected areas cover only 2.5% of the surviving microrefugia for this culturally important species, and none occur in any designated indigenous territory. Our results suggest that anthropogenic climate change is likely to commit Araucaria to a second consecutive century of significant losses, but targeted interventions could help ensure its survival in the wild.

A Universal Probe Set for Targeted Sequencing of 353 Nuclear Genes from Any Flowering Plant Designed Using k-Medoids Clustering

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5-15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself.