Comparative transcriptome analysis reveals lineage- and environment-specific adaptations in cacti from the Brazilian Atlantic Forest.

MAIN CONCLUSION: Natural selection influenced adaptive divergence between Cereus fernambucensis and Cereus insularis, revealing key genes governing abiotic stress responses and supporting neoteny in C. insularis. Uncovering the molecular mechanisms driving adaptive divergence in traits related to habitat adaptation remains a central challenge. In this study, we focused on the cactus clade, which includes Cereus sericifer F.Ritter, Cereus fernambucensis Lem., and Cereus insularis Hemsley. These allopatric species inhabit distinct relatively drier regions within the Brazilian Atlantic Forest, each facing unique abiotic conditions. We leveraged whole transcriptome data and abiotic variables datasets to explore lineage-specific and environment-specific adaptations in these species. Employing comparative phylogenetic methods, we identified genes under positive selection (PSG) and examined their association with non-synonymous genetic variants and abiotic attributes through a PhyloGWAS approach. Our analysis unveiled signatures of selection in all studied lineages, with C. fernambucensis northern populations and C. insularis showing the most PSGs. These PSGs predominantly govern abiotic stress regulation, encompassing heat tolerance, UV stress response, and soil salinity adaptation. Our exclusive observation of gene expression tied to early developmental stages in C. insularis supports the hypothesis of neoteny in this species. We also identified genes associated with abiotic variables in independent lineages, suggesting their role as environmental filters on genetic diversity. Overall, our findings suggest that natural selection played a pivotal role in the geographic range of these species in response to environmental and biogeographic transitions.

Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence.

The senescence of fruit is a complex physiological process, with various cell types within the pericarp, making it highly challenging to elucidate their individual roles in fruit senescence. In this study, a single-cell expression atlas of the pericarp of pitaya (Hylocereus undatus) is constructed, revealing exocarp and mesocarp cells undergoing the most significant changes during the fruit senescence process. Pseudotime analysis establishes cellular differentiation and gene expression trajectories during senescence. Early-stage oxidative stress imbalance is followed by the activation of resistance in exocarp cells, subsequently senescence-associated proteins accumulate in the mesocarp cells at late-stage senescence. The central role of the early response factor HuCMB1 is unveiled in the senescence regulatory network. This study provides a spatiotemporal perspective for a deeper understanding of the dynamic senescence process in plants.

Host-pathogen interaction between pitaya and Neoscytalidium dimidiatum reveals the mechanisms of immune response associated with defense regulators and metabolic pathways.

BACKGROUND: Understanding how plants and pathogens regulate each other's gene expression during their interactions is key to revealing the mechanisms of disease resistance and controlling the development of pathogens. Despite extensive studies on the molecular and genetic basis of plant immunity against pathogens, the influence of pitaya immunity on N. dimidiatum metabolism to restrict pathogen growth is poorly understood, and how N. dimidiatum breaks through pitaya defenses. In this study, we used the RNA-seq method to assess the expression profiles of pitaya and N. dimidiatum at 4 time periods after interactions to capture the early effects of N. dimidiatum on pitaya processes. RESULTS: The study defined the establishment of an effective method for analyzing transcriptome interactions between pitaya and N. dimidiatum and to obtain global expression profiles. We identified gene expression clusters in both the host pitaya and the pathogen N. dimidiatum. The analysis showed that numerous differentially expressed genes (DEGs) involved in the recognition and defense of pitaya against N. dimidiatum, as well as N. dimidiatum's evasion of recognition and inhibition of pitaya. The major functional groups identified by GO and KEGG enrichment were responsible for plant and pathogen recognition, phytohormone signaling (such as salicylic acid, abscisic acid). Furthermore, the gene expression of 13 candidate genes involved in phytopathogen recognition, phytohormone receptors, and the plant resistance gene (PG), as well as 7 effector genes of N. dimidiatum, including glycoside hydrolases, pectinase, and putative genes, were validated by qPCR. By focusing on gene expression changes during interactions between pitaya and N. dimidiatum, we were able to observe the infection of N. dimidiatum and its effects on the expression of various defense components and host immune receptors. CONCLUSION: Our data show that various regulators of the immune response are modified during interactions between pitaya and N. dimidiatum. Furthermore, the activation and repression of these genes are temporally coordinated. These findings provide a framework for better understanding the pathogenicity of N. dimidiatum and its role as an opportunistic pathogen. This offers the potential for a more effective defense against N. dimidiatum.

Evidences of differential methylation in the genome during development in the cactophilic Drosophila species.

DNA methylation with 5-methylcytosine (5mC) has been reported in the genome of several eukaryotes, with marked differences between vertebrates and invertebrates. DNA methylation is poorly understood as its role in evolution in insects. Drosophila gouveai (cluster Drosophila buzzatii) presents larvae that develop obligatorily in necrotic tissues of cacti in nature, with the distribution of populations in South America, and plasticity of phenotypes in insect-plant interaction. We characterize organisms at developmental stages and analyze variations at multiple methylation-sensitive loci in pupae, and adult flies using methylation sensitive amplification polymorphism. We obtained 326 loci with CCGG targets in the genome of D. gouveai. Genomic regions with molecular lengths from 100 to 700 pb were most informative about methylation states. Multiple loci show differences in methylation-sensitive sites (MSL) concerning developmental stages, such as in pupae (MSL = 40), female reproductive tissue (MSL = 76), and male reproductive tissues (MSL = 58). Our results are the first evidence of genome-wide methylation in D. gouveai organisms.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications.

BACKGROUND: Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported. RESULTS: This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3). CONCLUSIONS: This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

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.

Identification of HuSWEET Family in Pitaya (Hylocereus undatus) and Key Roles of HuSWEET12a and HuSWEET13d in Sugar Accumulation.

The sugar composition and content of fruit have a significant impact on their flavor and taste. In pitaya, or dragon fruit, sweetness is a crucial determinant of fruit taste and consumer preference. The sugars will eventually be exported transporters (SWEETs), a novel group of sugar transporters that have various physiological functions, including phloem loading, seed filling, nectar secretion, and fruit development. However, the role of SWEETs in sugar accumulation in pitaya fruit is not yet clear. Here, we identified 19 potential members (HuSWEET genes) of the SWEET family in pitaya and analyzed their conserved motifs, physiochemical characteristics, chromosomal distribution, gene structure, and phylogenetic relationship. Seven highly conserved α-helical transmembrane domains (7-TMs) were found, and the HuSWEET proteins can be divided into three clades based on the phylogenetic analysis. Interestingly, we found two HuSWEET genes, HuSWEET12a and HuSWEET13d, that showed strong preferential expressions in fruits and an upward trend during fruit maturation, suggesting they have key roles in sugar accumulation in pitaya. This can be further roughly demonstrated by the fact that transgenic tomato plants overexpressing HuSWEET12a/13d accumulated high levels of sugar in the mature fruit. Together, our result provides new insights into the regulation of sugar accumulation by SWEET family genes in pitaya fruit, which also set a crucial basis for the further functional study of the HuSWEETs.

Are 150 km of open sea enough? Gene flow and population differentiation in a bat-pollinated columnar cactus.

Genetic differentiations and phylogeographical patterns are controlled by the interplay between spatial isolation and gene flow. To assess the extent of gene flow across an oceanic barrier, we explored the effect of the separation of the peninsula of Baja California on the evolution of mainland and peninsular populations of the long-lived columnar cactus Stenocereus thurberi. We analyzed twelve populations throughout the OPC distribution range to assess genetic diversity and structure using chloroplast DNA sequences. Genetic diversity was higher (Hd = 0.81), and genetic structure was lower (GST = 0.143) in mainland populations vs peninsular populations (Hd = 0.71, GST = 0.358 respectively). Genetic diversity was negatively associated with elevation but positively with rainfall. Two mainland and one peninsular ancestral haplotypes were reconstructed. Peninsular populations were as isolated among them as with mainland populations. Peninsular haplotypes formed a group with one mainland coastal population, and populations across the gulf shared common haplotypes giving support to regular gene flow across the Gulf. Gene flow is likely mediated by bats, the main pollinators and seed dispersers. Niche modeling suggests that during the Last Glacial Maximum (c. 130 ka), OPC populations shrank to southern locations. Currently, Stenocereus thurberi populations are expanding, and the species is under population divergence despite ongoing gene flow. Ancestral populations are located on the mainland and although vicariant peninsular populations cannot be ruled out, they are likely the result of gene flow across the seemingly formidable barrier of the Gulf of California. Still, unique haplotypes occur in the peninsula and the mainland, and peninsular populations are more structured than those on the mainland.

Stem transcriptome screen for selection in wild and cultivated pitahaya (Selenicereus undatus): an epiphytic cactus with edible fruit.

Dragon fruit, pitahaya or pitaya are common names for the species in the Hylocereus group of Selenicereus that produce edible fruit. These Neotropical epiphytic cacti are considered promising underutilized crops and are currently cultivated around the world. The most important species, S. undatus, has been managed in the Maya domain for centuries and is the focus of this article. Transcriptome profiles from stems of wild and cultivated plants of this species were compared. We hypothesized that differences in transcriptomic signatures could be associated with genes related to drought stress. De novo transcriptome assembly and the analysis of differentially expressed genes (DEGs) allowed us to identify a total of 9,203 DEGs in the Hunucmá cultivar relative of wild Mozomboa plants. Of these, 4,883 represent up-regulated genes and 4,320, down-regulated genes. Additionally, 6,568 DEGs were identified from a comparison between the Umán cultivar and wild plants, revealing 3,286 up-regulated and 3,282 down-regulated genes. Approximately half of the DEGs are shared by the two cultivated plants. Differences between the two cultivars that were collected in the same region could be the result of differences in management. Metabolism was the most representative functional category in both cultivars. The up-regulated genes of both cultivars formed a network related to the hormone-mediated signaling pathway that includes cellular responses to auxin stimulus and to hormone stimulus. These cellular reactions have been documented in several cultivated plants in which drought-tolerant cultivars modify auxin transport and ethylene signaling, resulting in a better redistribution of assimilates.

Identification of RT-qPCR reference genes suitable for gene function studies in the pitaya canker disease pathogen Neoscytalidium dimidiatum.

Neoscytalidium dimidiatum is the main causal agent of pitaya canker. Most studies of virulence and pathogenicity genes have measured expression levels using real-time quantitative polymerase chain reaction (RT-qPCR). Suitable reference genes are essential for ensuring that estimates of gene expression levels by RT-qPCR are accurate. However, no reference genes can be robustly applied across all contexts and species. No studies to date have evaluated the most effective reference genes for normalizing gene expression levels estimated by RT-qPCR in N. dimidiatum. In this study, RT-qPCR data for individual candidate reference genes were analyzed using four different methods: the delta Ct method and the geNorm, NormFinder, and BestKeeper algorithms. We evaluated the utility of eight candidate reference genes (18S rRNA, Actin (1), Actin (2), Actin, GAPDH (1), GAPDH (2), UBQ, and Tubulin) for normalizing expression levels estimated by RT-qPCR in N. dimidiatum at different developmental stages, at different temperatures, and during interaction with pitaya. All candidate reference genes were suitable for gene expression analysis except for Actin (2). Tubulin and Actin (1) were the most stably expressed reference genes under different temperatures. Actin (1) and Actin were the most stably expressed reference genes in N. dimidiatum at different developmental stages. Tubulin and UBQ were the most stably expressed reference genes during interaction with pitaya. Actin and 18s rRNA were the most stably expressed across all experimental conditions. Subsequently, Tubulin and UBQ were further investigated in analyses of pectinase expression during the pitaya-N. dimidiatum interaction. Our results provide insights that will aid future RT-qPCR studies of gene expression in N. dimidiatum.

Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects.

Hypogeococcus pungens is a species complex native to southern South America that is composed of at least five putative species, each one specialized in the use of different host plants. Two of these undescribed species were registered as invasive in Central and North America: Hyp-C is a cactophagous mealybug that became an important pest that threatens endemic cactus species in Puerto Rico, and Hyp-AP feeds on Amaranthaceae and Portulacaceae hosts, but does not produce severe damage to the host plants. We quantified genomic variation and investigated the demographic history of both invasive species by means of coalescent-based simulations using high throughput sequencing data. We also evaluated the incidence of host plant infestation produced by both species and used an ecological niche modeling approach to assess potential distribution under current and future climatic scenarios. Our genetic survey evinced the footprints of strong effective population size reduction and signals of genetic differentiation among populations within each species. Incidence of plant attacks varied between species and among populations within species, with some host plant species preferred over others. Ecological niche modeling suggested that under future climatic scenarios both species would expand their distribution ranges in Puerto Rico. These results provide valuable information for the design of efficient management and control strategies of the Puerto Rican cactus pest and shed light on the evolutionary pathways of biological invasions.

Genome-wide identification of GRF gene family and their contribution to abiotic stress response in pitaya (Hylocereus polyrhizus).

Growth-regulating factors (GRFs) are plant-specific transcription factors identified in many land plants. Recently, their indispensable roles in stress response are highlighted. In present work, 11 HpGRFs were cloned in pitaya. Segmental duplication is considered essential for the expansion of HpGRFs. A phylogenetic tree suggested that GRFs could be divided into eight categories, among which G-I was a Caryophyllales-specific one. The categorization was further evidenced by differences in the gene structure, collinearity, protein domain of HpGRFs. Five miR396 hairpins giving rise to two types of matured miR396s were identified in pitaya via sRNA-Seq in combination with bioinformatic analysis. Parallel analysis of RNA ends proved that HpGRFs except HpGRF5 were degraded by miR396-directed cleavages at the regions which code the conserved WRC motifs of HpGRFs. Multiple cis-regulatory elements were discovered in the promoters of HpGRFs. Among the elements, most are involved in stress and phytohormone response as well as plant growth, indicating a crosstalk between them. Expression analysis showed the responsive patterns of the miR396-GRF module under abiotic stresses. To conclude, our work systematically identified the miR396-targeted HpGRFs in pitaya and confirmed their involvement in stress response, providing novel insights into the comprehensive understanding of the stress resistance of pitaya.

The evidence for anthocyanins in the betalain-pigmented genus Hylocereus is weak.

Here we respond to Zhou (BMC Genomics 21:734, 2020) "Combined Transcriptome and Metabolome analysis of Pitaya fruit unveiled the mechanisms underlying peel and pulp color formation" published in BMC Genomics. Given the evolutionary conserved anthocyanin biosynthesis pathway in betalain-pigmented species, we are open to the idea that species with both anthocyanins and betalains might exist. However, in absence of LC-MS/MS spectra, apparent lack of biological replicates, and no comparison to authentic standards, the findings of Zhou (BMC Genomics 21:734, 2020) are not a strong basis to propose the presence of anthocyanins in betalain-pigmented pitaya. In addition, our re-analysis of the datasets indicates the misidentification of important genes and the omission of key flavonoid and anthocyanin synthesis genes ANS and DFR. Finally, our re-analysis of the RNA-Seq dataset reveals no correlation between anthocyanin biosynthesis gene expression and pigment status.

Remarkably rapid, recent diversification of Cochemiea and Mammillaria in the Baja California, Mexico region.

PREMISE: The Cactaceae of northwestern Mexico and the southwestern United States constitute a major component of the angiosperm biodiversity of the region. The Mammilloid clade, (Cactaceae, tribe Cacteae), composed of the genera Cochemiea, Coryphantha, Cumarinia, Mammillaria, and Pelecyphora is especially species rich. We sought to understand the timing, geographical and climate influences correlated with expansion of the Mammilloid clade, through the Sonoran Desert into Baja California. METHODS: We reconstructed the historical biogeography of the Mammilloid clade, using Bayesian and maximum likelihood methods, based on a strongly supported molecular phylogeny. We also estimated divergence times, the timing of emergence of key characters, and diversification rates and rate shifts of the Mammilloid clade. RESULTS: We found that the most recent common ancestor of Cochemiea arrived in the Cape region of Baja California from the Sonoran Desert region approximately 5 million years ago, coinciding with the timing of peninsular rifting from the mainland, suggesting dispersal and vicariance as causes of species richness and endemism. The diversification rate for Cochemiea is estimated to be approximately 12 times that of the mean background diversification rate for angiosperms. Divergence time estimation shows that many of the extant taxa in Cochemiea and Baja California Mammillaria emerged from common ancestors 1 million to 200,000 years ago, having a mid-Pleistocene origin. CONCLUSIONS: Cochemiea and Mammillaria of the Baja California region are examples of recent, rapid diversification. Geological and climatic forces at multiple spatial and temporal scales are correlated with the western distributions of the Mammilloid clade.

Transcriptional responses are oriented towards different components of the rearing environment in two Drosophila sibling species.

BACKGROUND: The chance to compare patterns of differential gene expression in related ecologically distinct species can be particularly fruitful to investigate the genetics of adaptation and phenotypic plasticity. In this regard, a powerful technique such as RNA-Seq applied to ecologically amenable taxa allows to address issues that are not possible in classic model species. Here, we study gene expression profiles and larval performance of the cactophilic siblings Drosophila buzzatii and D. koepferae reared in media that approximate natural conditions and evaluate both chemical and nutritional components of the diet. These closely related species are complementary in terms of host-plant use since the primary host of one is the secondary of the other. D. koepferae is mainly a columnar cactus dweller while D. buzzatii prefers Opuntia hosts. RESULTS: Our comparative study shows that D. buzzatii and D. koepferae have different transcriptional strategies to face the challenges posed by their natural resources. The former has greater transcriptional plasticity, and its response is mainly modulated by alkaloids of its secondary host, while the latter has a more canalized genetic response, and its transcriptional plasticity is associated with the cactus species. CONCLUSIONS: Our study unveils a complex pleiotropic genetic landscape in both species, with functional links that relate detox responses and redox mechanisms with developmental and neurobiological processes. These results contribute to deepen our understanding of the role of host plant shifts and natural stress driving ecological specialization.

Pitaya Genome and Multiomics Database (PGMD): A Comprehensive and Integrative Resource of Selenicereus undatus.

Pitaya (Selenicereus) is a kind of novel fruit with a delicious taste and superior horticulture ornamental value. The potential economic impact of the pitaya lies in its diverse uses not only as agricultural produce and processed foods but also in industrial and medicinal products. It is also an excellent plant material for basic and applied biological research. A comprehensive database of pitaya would facilitate studies of pitaya and the other Cactaceae plant species. Here, we constructed pitaya genome and multiomics database, which is a collection of the most updated and high-quality pitaya genomic assemblies. The database contains various information such as genomic variation, gene expression, miRNA profiles, metabolite and proteomic data from various tissues and fruit developmental stages of different pitaya cultivars. In PGMD, we also uploaded videos on the flowering process and planting tutorials for practical usage of pitaya. Overall, these valuable data provided in the PGMD will significantly facilitate future studies on population genetics, molecular breeding and function research of pitaya.

Effectiveness of Two Universal Angiosperm Probe Sets Tested In Silico for Caryophyllids Taxa with Emphasis on Cacti Species.

In angiosperms, huge advances in massive DNA sequencing technologies have impacted phylogenetic studies. Probe sets have been developed with the purpose of recovering hundreds of orthologous loci of targeted DNA sequences (TDS) across different plant lineages. We tested in silico the effectiveness of two universal probe sets in the whole available genomes of Caryophyllids, emphasizing phylogenetic issues in cacti species. A total of 870 TDS (517 TDS from Angiosperm v.1 and 353 from Angiosperms353) were individually tested in nine cacti species and Amaranthus hypochondriacus (external group) with ≥17 Gbp of available DNA data. The effectiveness was measured by the total number of orthologous loci recovered and their length, the percentage of loci discarded by paralogy, and the proportion of informative sites (PIS) in the alignments. The results showed that, on average, Angiosperms353 was better than Angiosperm v.1 for cacti species, since the former obtained an average of 275.6 loci that represent 123,687 bp, 2.48% of paralogous loci, and 4.32% of PIS in alignments, whereas the latter recovered 148.4 loci (37,683 bp), 10.38% of paralogous loci, and 3.49% of PIS. We recommend the use of predesigned universal probe sets for Caryophyllids, since these recover a high number of orthologous loci that resolve phylogenetic relationships.

A target Capture Probe Set Useful for Deep- and Shallow-Level Phylogenetic Studies in Cactaceae.

The molecular phylogenies of Cactaceae have enabled us to better understand their systematics, biogeography, and diversification ages. However, most of the phylogenetic relationships within Cactaceae major groups remain unclear, largely due to the lack of an appropriate set of molecular markers to resolve its contentious relationships. Here, we explored the genome and transcriptome assemblies available for Cactaceae and identified putative orthologous regions shared among lineages of the subfamily Cactoideae. Then we developed a probe set, named Cactaceae591, targeting both coding and noncoding nuclear regions for representatives from the subfamilies Pereskioideae, Opuntioideae, and Cactoideae. We also sampled inter- and intraspecific variation to evaluate the potential of this panel to be used in phylogeographic studies. We retrieved on average of 547 orthologous regions per sample. Targeting noncoding nuclear regions showed to be crucial to resolving inter- and intraspecific relationships. Cactaceae591 covers 13 orthologous genes shared with the Angiosperms353 kit and two plastid regions largely used in Cactaceae studies, enabling the phylogenies generated by our panel to be integrated with angiosperm and Cactaceae phylogenies, using these sequences. We highlighted the importance of using coalescent-based species tree approaches on the Cactaceae591 dataset to infer accurate phylogenetic trees in the presence of extensive incomplete lineage sorting in this family.

Evolutionary Genetics of Cacti: Research Biases, Advances and Prospects.

Here, we present a review of the studies of evolutionary genetics (phylogenetics, population genetics, and phylogeography) using genetic data as well as genome scale assemblies in Cactaceae (Caryophyllales, Angiosperms), a major lineage of succulent plants with astonishing diversity on the American continent. To this end, we performed a literature survey (1992-2021) to obtain detailed information regarding key aspects of studies investigating cactus evolution. Specifically, we summarize the advances in the following aspects: molecular markers, species delimitation, phylogenetics, hybridization, biogeography, and genome assemblies. In brief, we observed substantial growth in the studies conducted with molecular markers in the past two decades. However, we found biases in taxonomic/geographic sampling and the use of traditional markers and statistical approaches. We discuss some methodological and social challenges for engaging the cactus community in genomic research. We also stressed the importance of integrative approaches, coalescent methods, and international collaboration to advance the understanding of cactus evolution.

Metagenomic analysis of the intestinal microbiome in goats on cactus and Salicornia-based diets.

Background: The Peruvian coast is characterized by its arid and saline soils, the cactus being an alternative for arid soils and Salicornia for saline soils. Therefore, it is necessary to develop nutrition based on the intestinal microbiota in goats. Aim: To identify the intestinal microbiota in goats through a metagenomic analysis. Methods: In this study, goats and kids were randomly selected and fed cacti and Salicornia as potential forage species compared to native grass to study the changes in the microbiota using massive sequencing using the 16S rRNA gene as a marker. Results: The sequencing results showed the taxonomic levels of Bacteroidetes and Firmicutes at the phylum level as the most abundant in creole goats' microbiome, varying from 18% to 36% and 47% to 66%, respectively. At the genus level, variants of the genus Ruminococcaceae stand out, related to cellulose degradation, as the most dominant in all samples, followed by Christensenellaceae, Rikenellaceae, and Prevotellaceae. Also, the genus Akkermansia appeared in greater abundance in kids fed with cactus, being necessary for being related to the intestinal mucosa's health and avoiding the adhesion of pathogens to the intestinal epithelium. Conclusion: These microbiota changes based on diets with high fiber content are necessary to understand the adaptation of this species to favorable dietary changes.