Unravelling phylogenetic relationships of the tribe Cereeae using target enrichment sequencing.

BACKGROUND AND AIMS: Cactaceae are succulent plants, quasi-endemic to the American continent, and one of the most endangered plant groups in the world. Molecular phylogenies have been key to unravelling phylogenetic relationships among major cactus groups, previously hampered by high levels of morphological convergence. Phylogenetic studies using plastid markers have not provided adequate resolution for determining generic relationships within cactus groups. This is the case for the tribe Cereeae s.l., a highly diverse group from tropical America. Here we aimed to reconstruct a well-resolved phylogenetic tree of tribe Cereeae and update the circumscription of suprageneric and generic groups in this tribe. METHODS: We integrated sequence data from public gene and genomic databases with new target sequences (generated using the customized Cactaceae591 probe set) across representatives of this tribe, with a denser taxon sampling of the subtribe Cereinae. We inferred concatenated and coalescent phylogenetic trees and compared the performance of both approaches. KEY RESULTS: Six well-supported suprageneric clades were identified using different datasets. However, only genomic datasets, especially the Cactaceae591, were able to resolve the contentious relationships within the subtribe Cereinae. CONCLUSIONS: We propose a new taxonomic classification within Cereeae based on well-resolved clades, including new subtribes (Aylosterinae subtr. nov., Uebelmanniinae subtr. nov. and Gymnocalyciinae subtr. nov.) and revised subtribes (Trichocereinae, Rebutiinae and Cereinae). We emphasize the importance of using genomic datasets allied with coalescent inference to investigate evolutionary patterns within the tribe Cereeae.

The iconic cactus of the Caatinga dry forest, Cereus jamacaru (Cactaceae) has high sphingophily specialization and pollinator dependence.

Cereus jamacaru is a cactus distributed in Northeastern Brazil, with high symbolic value to this region. However, the interaction, behavior and the role of pollinators remains poorly understood. Here, we investigate the reproductive biology, addressing the ecological significance of floral attributes, including details about floral signaling. The study was carried at three areas of the Caatinga, in 2015, 2017 and 2021. We analyzed the floral morphometry, volume and concentration of the nectar, and characterized the colour and scent of flowers. Additionally, we described the pollinator behavior and performed controlled pollination experiments. The 'Mandacaru' is self-incompatible, has nocturnal anthesis and the nectar is accumulated as droplets in a long hypanthial tube. The flowers have a reflective pattern with a dark outer surface and a white inner surface. (E)-nerolidol is the major component (87.4%) of its floral perfume. We registered the sphingid moth Cocytius antaeus visiting the flowers. The floral attributes, attractants and rewards drives to a sphingophily, and the pollination treatments showed the dependence to fruit set by C. antaeus, the pollinator registered. In this case, if the apparent lack of pollinator diversity encompasses its entire range, the loss of the hawkmoth could severely impact the reproductive success of the cactus.

Revisiting phylogeny, systematics, and biogeography of a Pleistocene radiation.

PREMISE: Pilosocereus (Cactaceae) is an important dry forest element in all subregions and transitional zones of the neotropics, with the highest diversity in eastern Brazil. The genus is subdivided into informal taxonomic groups; however, most of these are not supported by recent molecular phylogenetic inferences. This lack of confidence is probably due to the use of an insufficient number of loci and the complexity of cactus diversification. Here, we explored the species relationships in Pilosocereus in more detail, integrating multilocus phylogenetic approaches with the assessment of the ancestral range and the effect of geography on diversification shifts. METHODS: We used 28 nuclear, plastid, and mitochondrial loci from 54 plant samples of 31 Pilosocereus species for phylogenetic analyses. We used concatenated and coalescent phylogenetic trees and Bayesian models to estimate the most likely ancestral range and diversification shifts. RESULTS: All Pilosocereus species were clustered in the same branch, except P. bohlei. The phylogenetic relationships were more associated with the geographic distribution than taxonomic affinities among taxa. The genus began diversifying during the Plio-Pleistocene transition in the Caatinga domain and experienced an increased diversification rate during the Calabrian age. CONCLUSIONS: We recovered a well-supported multispecies coalescent phylogeny. Our results refine the pattern of rapid diversification of Pilosocereus species across neotropical drylands during the Pleistocene and highlight the need for taxonomic rearrangements in the genus. We recovered a pulse of diversification during the Pleistocene that was likely driven by multiple dispersal and vicariance events within and among the Caatinga, Cerrado, and Atlantic Forest domains.

Coalescent-based species delimitation meets deep learning: Insights from a highly fragmented cactus system.

Delimiting species boundaries is a major goal in evolutionary biology. An increasing volume of literature has focused on the challenges of investigating cryptic diversity within complex evolutionary scenarios of speciation, including gene flow and demographic fluctuations. New methods based on model selection, such as approximate Bayesian computation, approximate likelihoods, and machine learning are promising tools arising in this field. Here, we introduce a framework for species delimitation using the multispecies coalescent model coupled with a deep learning algorithm based on convolutional neural networks (CNNs). We compared this strategy with a similar ABC approach. We applied both methods to test species boundary hypotheses based on current and previous taxonomic delimitations as well as genetic data (sequences from 41 loci) in Pilosocereus aurisetus, a cactus species complex with a sky-island distribution and taxonomic uncertainty. To validate our method, we also applied the same strategy on data from widely accepted species from the genus Drosophila. The results show that our CNN approach has a high capacity to distinguish among the simulated species delimitation scenarios, with higher accuracy than ABC. For the cactus data set, a splitter hypothesis without gene flow showed the highest probability in both CNN and ABC approaches, a result agreeing with previous taxonomic classifications and in line with the sky-island distribution and low dispersal of P. aurisetus. Our results highlight the cryptic diversity within the P. aurisetus complex and show that CNNs are a promising approach for distinguishing complex evolutionary histories, even outperforming the accuracy of other model-based approaches such as ABC.

Cross-genera SSR transferability in cacti revealed by a case study using Cereus (Cereeae, Cactaceae).

The study of transferability of simple sequence repeats (SSR) among closely related species is a well-known strategy in population genetics, however transferability among distinct genera is less common. We tested cross-genera SSR amplification in the family Cactaceae using a total of 20 heterologous primers previously developed for the genera Ariocarpus, Echinocactus, Polaskia and Pilosocereus, in four taxa of the genus Cereus: C. fernambucensis subsp. fernambucensis, C. fernambucensis subsp. sericifer, C. jamacaru and C. insularis. Nine microsatellite loci were amplified in Cereus resulting in 35.2% of success in transferability, which is higher than the average rate of 10% reported in the literature for cross-genera transferability in eudicots. The genetic variation in the transferred markers was sufficient to perform standard clustering analysis, indicating each population as a cohesive genetic cluster. Overall, the amount of genetic variation found indicates that the transferred SSR markers might be useful in large-scale population studies within the genus Cereus.

Lineage-specific evolutionary rate in plants: Contributions of a screening for Cereus (Cactaceae).

PREMISE OF THE STUDY: Predictable chloroplast DNA (cpDNA) sequences have been listed for the shallowest taxonomic studies in plants. We investigated whether plastid regions that vary between closely allied species could be applied for intraspecific studies and compared the variation of these plastid segments with two nuclear regions. METHODS: We screened 16 plastid and two nuclear intronic regions for species of the genus Cereus (Cactaceae) at three hierarchical levels (species from different clades, species of the same clade, and allopatric populations). RESULTS: Ten plastid regions presented interspecific variation, and six of them showed variation at the intraspecific level. The two nuclear regions showed both inter- and intraspecific variation, and in general they showed higher levels of variability in almost all hierarchical levels than the plastid segments. DISCUSSION: Our data suggest no correspondence between variation of plastid regions at the interspecific and intraspecific level, probably due to lineage-specific variation in cpDNA, which appears to have less effect in nuclear data. Despite the heterogeneity in evolutionary rates of cpDNA, we highlight three plastid segments that may be considered in initial screenings in plant phylogeographic studies.

High proportion of cactus species threatened with extinction.

A high proportion of plant species is predicted to be threatened with extinction in the near future. However, the threat status of only a small number has been evaluated compared with key animal groups, rendering the magnitude and nature of the risks plants face unclear. Here we report the results of a global species assessment for the largest plant taxon evaluated to date under the International Union for Conservation of Nature (IUCN) Red List Categories and Criteria, the iconic Cactaceae (cacti). We show that cacti are among the most threatened taxonomic groups assessed to date, with 31% of the 1,478 evaluated species threatened, demonstrating the high anthropogenic pressures on biodiversity in arid lands. The distribution of threatened species and the predominant threatening processes and drivers are different to those described for other taxa. The most significant threat processes comprise land conversion to agriculture and aquaculture, collection as biological resources, and residential and commercial development. The dominant drivers of extinction risk are the unscrupulous collection of live plants and seeds for horticultural trade and private ornamental collections, smallholder livestock ranching and smallholder annual agriculture. Our findings demonstrate that global species assessments are readily achievable for major groups of plants with relatively moderate resources, and highlight different conservation priorities and actions to those derived from species assessments of key animal groups.

Interglacial microrefugia and diversification of a cactus species complex: phylogeography and palaeodistributional reconstructions for Pilosocereus aurisetus and allies.

The role of Pleistocene climate changes in promoting evolutionary diversification in global biota is well documented, but the great majority of data regarding this subject come from North America and Europe, which were greatly affected by glaciation. The effects of Pleistocene changes on cold- and/or dry-adapted species in tropical areas where glaciers were not present remain sparsely investigated. Many such species are restricted to small areas surrounded by unfavourable habitats, which may represent potential interglacial microrefugia. Here, we analysed the phylogeographic structure and diversification history of seven cactus species in the Pilosocereus aurisetus complex that are restricted to rocky areas with high diversity and endemism within the Neotropical savannas of eastern South America. We combined palaeodistributional estimates with standard phylogeographic approaches based on two chloroplast DNA regions (trnT-trnL and trnS-trnG), exon 1 of the nuclear gene PhyC and 10 nuclear microsatellite loci. Our analyses revealed a phylogeographic history marked by multiple levels of distributional fragmentation, isolation leading to allopatric differentiation and secondary contact among divergent lineages within the complex. Diversification and demographic events appear to have been affected by the Quaternary climatic cycles as a result of isolation in multiple patches of xerophytic vegetation. These small patches presently harbouring P. aurisetus populations seem to operate as microrefugia, both at present and during Pleistocene interglacial periods; the role of such microrefugia should be explored and analysed in greater detail.

Cross-species amplification of microsatellites reveals incongruence in the molecular variation and taxonomic limits of the Pilosocereus aurisetus group (Cactaceae).

The Pilosocereus aurisetus group contains eight cactus species restricted to xeric habitats in eastern and central Brazil that have an archipelago-like distribution. In this study, 5-11 microsatellite markers previously designed for Pilosocereus machrisii were evaluated for cross-amplification and polymorphisms in ten populations from six species of the P. aurisetus group. The genotypic information was subsequently used to investigate the genetic relationships between the individuals, populations, and species analyzed. Only the Pmac101 locus failed to amplify in all of the six analyzed species, resulting in an 88 % success rate. The number of alleles per polymorphic locus ranged from 2 to 12, and the most successfully amplified loci showed at least one population with a larger number of alleles than were reported in the source species. The population relationships revealed clear genetic clustering in a neighbor-joining tree that was partially incongruent with the taxonomic limits between the P. aurisetus and P. machrisii species, a fact which parallels the problematic taxonomy of the P. aurisetus group. A Bayesian clustering analysis of the individual genotypes confirmed the observed taxonomic incongruence. These microsatellite markers provide a valuable resource for facilitating large-scale genetic studies on population structures, systematics and evolutionary history in this group.

What does it take to resolve relationships and to identify species with molecular markers? An example from the epiphytic Rhipsalideae (Cactaceae).

PREMISE OF THE STUDY: The Cactaceae are a major New World plant family and popular in horticulture. Still, taxonomic units and species limits have been difficult to define, and molecular phylogenetic studies so far have yielded largely unresolved trees, so relationships within Cactaceae remain insufficiently understood. This study focuses on the predominantly epiphytic tribe Rhipsalideae and evaluates the utility of a spectrum of plastid genomic regions. • METHODS: We present a phylogenetic study including 52 of the 53 Rhipsalideae species and all the infraspecific taxa. Seven regions (trnK intron, matK, rbcL, rps3-rpl16, rpl16 intron, psbA-trnH, trnQ-rps16), ca. 5600 nucleotides (nt) were sequenced per sample. The regions used were evaluated for their phylogenetic performance and performance in DNA-based species recognition based on operational taxonomic units (OTUs) defined beforehand. • KEY RESULTS: The Rhipsalideae are monophyletic and contain five clades that correspond to the genera Rhipsalis, Lepismium, Schlumbergera, Hatiora, and Rhipsalidopsis. The species-level tree was well resolved and supported; the rpl16 and trnK introns yielded the best phylogenetic signal. Although the psbA-trnH and trnQ-rps16 spacers were the most successful individual regions for OTU identification, their success rate did not significantly exceed 70%. The highest OTU identification rate of 97% was found using the combination of psbA-trnH, rps3-rpl16, trnK intron, and trnQ-rps16 as a minimum possible marker length (ca. 1660 nt). • CONCLUSIONS: The phylogenetic performance of a marker is not determined by the level of sequence variability, and species discrimination power does not necessarily correlate with phylogenetic utility.

Isolation, characterization, and cross-species amplification of polymorphic microsatellite markers for Pilosocereus machrisii (Cactaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed in Pilosocereus machrisii, a columnar cactus with a patchy distribution in eastern tropical South America, to assess its genetic diversity and population structure. METHODS AND RESULTS: Eleven microsatellite markers were developed, of which one was monomorphic among 51 individuals from two populations. The number of alleles per polymorphic locus ranged from two to eight, and the total number of alleles was 57. From the 11 isolated loci, nine were successfully amplified in the other four Pilosocereus species. CONCLUSIONS: The results showed that these markers will be useful for studies of genetic diversity, patterns of gene flow, and population genetic structure in P. machrisii, as well as across other congeneric species.