Historical biogeography and phylogeny of Cucurbita: Insights from ancestral area reconstruction and niche evolution

Fluctuation of ecological niches and geographic range shifts along chile pepper's domestication gradient

Domestication is an ongoing well-described process. However, while many have studied the changes domestication causes in plant genetics, few have explored its impact on the portion of the geographic landscape in which the plants exist. Therefore, the goal of this study was to understand how the process of domestication changed the geographic space suitable for chile pepper (Capsicum annuum) in its center of origin (domestication). C. annuum is a major crop species globally whose center of domestication, Mexico, has been well-studied. It provides a unique opportunity to explore the degree to which ranges of different domestication classes diverged and how these ranges might be altered by climate change. To this end, we created ecological niche models for four domestication classes (wild, semiwild, landrace, modern cultivar) based on present climate and future climate scenarios for 2050, 2070, and 2090. Considering present environment, we found substantial overlap in the geographic niches of all the domestication classes. Yet, environmental and geographic aspects of the current ranges did vary among classes. Wild and commercial varieties could grow in desert conditions, while landraces could not. With projections into the future, habitat was lost asymmetrically, with wild, semiwild, and landraces at greater risk of territorial declines than modern cultivars. Further, we identified areas where future suitability overlap between landraces and wilds is expected to be lost. While range expansion is widely associated with domestication, we found little support of a constant niche expansion (either in environmental or geographical space) throughout the domestication gradient in chile peppers in Mexico. Instead, particular domestication transitions resulted in loss, followed by capturing or recapturing environmental or geographic space. The differences in environmental characterization among domestication gradient classes and their future potential range shifts increase the need for conservation efforts to preserve landraces and semiwild genotypes.

More phylogenetically diverse polycultures inconsistently suppress insect herbivore populations

Because the diet of many herbivorous insects is restricted to closely related taxa with similar chemistry, intercropping with diverse plant communities may reduce both pest populations and reliance on chemical pesticides in agroecosystems. We tested whether the effectiveness of intercropping against herbivorous insects depends on the phylogenetic relatedness of neighboring crops, using butternut squash (Cucurbita moschata) as a focal crop species in a series of different intercropping combinations. We found that increased phylogenetic divergence of neighboring plants could reduce abundance of herbivorous insects, but the effect was only detectable mid-season. In addition, we tested two hypothesized mechanisms for a negative association between phylogenetic distance of neighboring plants and reduced herbivore populations: one, we tested using Y-tube olfactometer and choice cage trials whether diverse volatile cues impede host-plant location by the dominant pest of butternut squash in our experiment, striped cucumber beetle Acalymma vittatum. Two, we recorded predator and parasitoid abundance relative to crop phylodiversity to test whether diverse crops support larger natural-enemy populations that can better control pest species. Our results, however, did not support either hypothesis. Striped cucumber beetles preferentially oriented toward non-host-plant volatiles, and predator populations more often decreased with phylodiversity than increased. Thus, the mechanisms driving associations in the field between phylogenetic divergence and herbivore populations remain unclear.

The effect of squash domestication on a belowground tritrophic interaction

The domestication of plants has commonly resulted in the loss of plant defense metabolites, with important consequences for the plants' interactions with herbivores and their natural enemies. Squash domestication started 10'000 years ago and has led to the loss of cucurbitacins, which are highly toxic triterpenes. The banded cucumber beetle (Diabrotica balteata), a generalist herbivore, is adapted to feed on plants from the Cucurbitaceae and is known to sequester cucurbitacins, supposedly for its own defense. However, the evidence for this is inconclusive. In this study we tested the impact of squash domestication on the chemical protection of D. balteata larvae against a predatory rove beetle (Dalotia coriaria). We found that cucurbitacins do not defend the larvae against this common soil dwelling predator. In fact, D. balteata larvae were less attacked when they fed on cucurbitacin-free roots of domesticated varieties compared to high-cucurbitacin roots of wild plants. This study appears to be the first to look at the consequences of plant domestication on belowground tritrophic interactions. Our results challenge the generalized assumption that sequestered cucurbitacins protect this herbivore against natural enemies, and instead reveals an opposite effect that may be due to a tradeoff between coping with cucurbitacins and avoiding predation.

Contrasting Metabolic Fingerprints and Seed Protein Profiles of Cucurbita foetidissima and C. radicans Fruits from Feral Plants Sampled in Central Mexico

Cucurbita foetidissima and C. radicans are scarcely studied wild pumpkin species that grow in arid and semi-arid areas of Mexico and the United States. This study describes the morphological, proximal composition, metabolic finger-prints and seed protein profiles of C. foetidissima and C. radicans fruits collected in the wild during a one-year period in different locations of central-western Mexico. The results obtained complement the limited information concerning the fruit composition of C. foetidissima and greatly expand information in this respect regarding C. radicans. Morphology and proximal composition of their fruits varied significantly. Different metabolic fingerprints and seed protein profiles were detected between them and also with the chemical composition of domesticated Cucurbita fruits. The neutral lipids in seed, pulp and peels were rich in wax content and in unsaturated compounds, probably carotenoids and tocopherols, in addition to tri-, di- and mono-acylglycerols. The tri- and diacylglycerol profiles of their seed oils were different from commercial seed oils and between each other. They also showed unusual fatty acid compositions. Evidence of a possible alkaloid in the pulp and peel of both species was obtained in addition to several putative cucurbitacins. An abundance of phenolic acids was found in all fruit parts, whereas flavonoids were only detected in the peels. Unlike most cucurbits, globulins were not the main protein fraction in the seeds of C. radicans, whereas the non-structural carbohydrate and raffinose oligosaccharide content in their fruit parts was lower than in other wild cucurbit species. These results emphasize the significantly different chemical composition of these two marginally studied Cucurbita species, which was more discrepant in C. radicans, despite the notion regarding C. foetidissima as an aberrant species with no affinity to any other Cucurbita species.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Genome-wide annotation reveals that the gene birth-death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

Recent and Historical Gene Flow in Cultivars, Landraces, and a Wild Taxon of Cucurbita pepo in Mexico

Gene flow among crops and their wild relatives is an active study area in evolutionary biology and horticulture, because genetic exchange between them may impact their evolutionary trajectories and increase the genetic variation of the cultivated lineages. Mexico is a center of diversity for the genus Cucurbita that includes pumpkins, squash and gourds. Gene flow between domesticated and wild species has been reported as common in Cucurbita; but gene flow among populations of C. pepo ssp. pepo from Mexico and its wild relative has not been studied. We used 2,061 SNPs, derived from tunable genotyping by sequencing (tGBS) to estimate gene flow among 14 Mexican traditional landraces of C. pepo ssp. pepo, also including individuals from five improved cultivars of C. pepo ssp. pepo and C. pepo ssp. ovifera var. ovifera, and individuals of their wild relative C. pepo ssp. fraterna. We found moderate to high levels of genetic diversity, and low to moderate genetic differentiation. In the test of introgression between lineages, we found that all possible arrangements for ancestral and derived sites between the lineages showed similar frequencies; thus, incomplete lineage sorting, but also gene flow, might be taking place in C. pepo. Overall, our results suggest that gene flow between these subspecies and cultigens, incomplete lineage sorting and the retention of ancestral characters shaped the evolutionary trajectory of C. pepo in its area of origin and diversification. In addition, we found evidence of the use of Mexican landraces as genetic material for the improvement of commercial cultivars. The landraces of Mexico are an important source of genetic diversity for C. pepo, which has been preserved both by management practices of small farmers and by the natural gene flow that exists between the different crop fields of the region.

The domestication of Cucurbita argyrosperma as revealed by the genome of its wild relative

Despite their economic importance and well-characterized domestication syndrome, the genomic impact of domestication and the identification of variants underlying the domestication traits in Cucurbita species (pumpkins and squashes) is currently lacking. Cucurbita argyrosperma, also known as cushaw pumpkin or silver-seed gourd, is a Mexican crop consumed primarily for its seeds rather than fruit flesh. This makes it a good model to study Cucurbita domestication, as seeds were an essential component of early Mesoamerican diet and likely the first targets of human-guided selection in pumpkins and squashes. We obtained population-level data using tunable Genotype by Sequencing libraries for 192 individuals of the wild and domesticated subspecies of C. argyrosperma across Mexico. We also assembled the first high-quality wild Cucurbita genome. Comparative genomic analyses revealed several structural variants and presence/absence of genes related to domestication. Our results indicate a monophyletic origin of this domesticated crop in the lowlands of Jalisco. We found evidence of gene flow between the domesticated and wild subspecies, which likely alleviated the effects of the domestication bottleneck. We uncovered candidate domestication genes that are involved in the regulation of growth hormones, plant defense mechanisms, seed development, and germination. The presence of shared selected alleles with the closely related species Cucurbita moschata suggests domestication-related introgression between both taxa.

Disentangling Domestication from Food Production Systems in the Neotropics

The Neolithic Revolution narrative associates early-mid Holocene domestications with the development of agriculture that fueled the rise of late Holocene civilizations. This narrative continues to be influential, even though it has been deconstructed by archaeologists and geneticists in its homeland. To further disentangle domestication from reliance on food production systems, such as agriculture, we revisit definitions of domestication and food production systems, review the late Pleistocene–early Holocene archaeobotanical record, and quantify the use, management and domestication of Neotropical plants to provide insights about the past. Neotropical plant domestication relies on common human behaviors (selection, accumulation and caring) within agroecological systems that focus on individual plants, rather than populations—as is typical of agriculture. The early archaeobotanical record includes numerous perennial and annual species, many of which later became domesticated. Some of this evidence identifies dispersal with probable cultivation, suggesting incipient domestication by 10,000 years ago. Since the Pleistocene, more than 6500, 1206 and 6261 native plant species have been used in Mesoamerica, the Central Andes and lowland South America, respectively. At least 1555, 428 and 742 are managed outside and inside food production systems, and at least 1148, 428 and 600 are cultivated, respectively, suggesting at least incipient domestication. Full native domesticates are more numerous in Mesoamerica (251) than the Andes (124) and the lowlands (45). This synthesis reveals that domestication is more common in the Neotropics than previously recognized and started much earlier than reliance on food production systems. Hundreds of ethnic groups had, and some still have, alternative strategies that do involve domestication, although they do not rely principally on food production systems, such as agriculture.

Climatic-niche evolution follows similar rules in plants and animals

Climatic niches are essential in determining where species can occur and how they will respond to climate change. However, it remains unclear if climatic-niche evolution is similar in plants and animals or is intrinsically different. For example, previous authors have proposed that plants have broader environmental tolerances than animals but are more sensitive to climate change. Here, we test ten predictions about climatic-niche evolution in plants and animals, using phylogenetic and climatic data for 19 plant clades and 17 vertebrate clades (2,087 species total). Surprisingly, we find that for all ten predictions, plants and animals show similar patterns. For example, in both groups, climatic niches change at similar mean rates and species have similar mean niche breadths, and niche breadths show similar relationships with latitude across groups. Our results suggest that there are general 'rules' of climatic-niche evolution that span plants and animals, despite the fundamental differences in their biology. These results may help to explain why plants and animals have similar responses to climate change and why they often have shared species richness patterns, biogeographic regions, biomes and biodiversity hotspots.

Phylogeographic and population genetic analyses of Cucurbita moschata reveal divergence of two mitochondrial lineages linked to an elevational gradient

PREMISE

Domestication usually involves local adaptation to environmental conditions. Cucurbita species are a promising model for studying these processes. Cucurbita moschata is the third major crop in the genus because of its economic value and because it displays high landrace diversity, but research about its genetic diversity, population structure, and phylogeography is limited. We aimed at understanding how geography and elevation shape the distribution of genetic diversity in C. moschata landraces in Mexico.

METHODS

We sampled fruits from 24 localities throughout Mexico. We assessed 11 nuclear microsatellite loci, one mtDNA region, and three cpDNA regions but found no variation in cpDNA. We explored genetic structure with cluster analysis, and phylogeographic relationships with haplotype network analysis.

RESULTS

Mitochondrial genetic diversity was high, and nuclear genetic differentiation among localities was intermediate compared to other domesticated Cucurbita. We found high levels of inbreeding. We recovered two mitochondrial lineages: highland (associated with the Trans-Mexican Volcanic Belt) and lowland. Nuclear microsatellites show that localities from the Yucatan Peninsula constitute a well-differentiated group.

CONCLUSIONS

Mexico is an area of high diversity for C. moschata, and these landraces represent important plant genetic resources. In Mexico this species is characterized by divergence processes linked to an elevational gradient, which could be related to adaptation and may be of value for applications in agriculture. The Isthmus of Tehuantepec may be a partial barrier to gene flow. Morphological variation, agricultural management, and cultural differences may be related to this pattern of genetic structure, but further studies are needed.

Evolutionary Dynamics of Transferred Sequences Between Organellar Genomes in Cucurbita

Twenty-nine DNA regions of plastid origin have been previously identified in the mitochondrial genome of Cucurbita pepo (pumpkin; Cucurbitaceae). Four of these regions harbor homolog sequences of rbcL, matK, rpl20-rps12 and trnL-trnF, which are widely used as molecular markers for phylogenetic and phylogeographic studies. We extracted the mitochondrial copies of these regions based on the mitochondrial genome of C. pepo and, along with published sequences for these plastome markers from 13 Cucurbita taxa, we performed phylogenetic molecular analyses to identify inter-organellar transfer events in the Cucurbita phylogeny and changes in their nucleotide substitution rates. Phylogenetic reconstruction and tree selection tests suggest that rpl20 and rbcL mitochondrial paralogs arose before Cucurbita diversification whereas the mitochondrial matK and trnL-trnF paralogs emerged most probably later, in the mesophytic Cucurbita clade. Nucleotide substitution rates increased one order of magnitude in all the mitochondrial paralogs compared to their original plastid sequences. Additionally, mitochondrial trnL-trnF sequences obtained by PCR from nine Cucurbita taxa revealed higher nucleotide diversity in the mitochondrial than in the plastid copies, likely related to the higher nucleotide substitution rates in the mitochondrial region and loss of functional constraints in its tRNA genes.

Tracing back the origin of pumpkins (Cucurbita pepo ssp. pepo L.) in Mexico

Cucurbita pepo is an economically important crop, which consists of cultivated C. pepo ssp. pepo, and two wild taxa (C. pepo ssp. fraterna and C. pepo ssp. ovifera). We aimed at understanding the domestication and the diversity of C. pepo in Mexico. We used two chloroplast regions and nine nuclear microsatellite loci to assess the levels of genetic variation and structure for C. pepo ssp. pepo's landraces sampled in 13 locations in Mexico, five improved varieties, one C. pepo ssp. fraterna population and ornamental C. pepo ssp. ovifera. We tested four hypotheses regarding the origin of C. pepo ssp. pepo's ancestor through approximate Bayesian computation: C. pepo ssp. ovifera as the ancestor; C. pepo ssp. fraterna as the ancestor; an unknown extinct lineage as the ancestor; and C. pepo ssp. pepo as hybrid from C. pepo ssp. ovifera and C. pepo ssp. fraterna ancestors. Cucurbita pepo ssp. pepo showed high genetic variation and low genetic differentiation. Cucurbita pepo ssp. fraterna and C. pepo ssp. pepo shared two chloroplast haplotypes. The three subspecies were well differentiated for microsatellite loci. Cucurbita pepo ssp. fraterna was probably C. pepo ssp. pepo's wild ancestor, but subsequent hybridization between taxa complicate defining C. pepo ssp. pepo's ancestor.

The Genome of Cucurbita argyrosperma (Silver-Seed Gourd) Reveals Faster Rates of Protein-Coding Gene and Long Noncoding RNA Turnover and Neofunctionalization within Cucurbita

Whole-genome duplications are an important source of evolutionary novelties that change the mode and tempo at which genetic elements evolve within a genome. The Cucurbita genus experienced a whole-genome duplication around 30 million years ago, although the evolutionary dynamics of the coding and noncoding genes in this genus have not yet been scrutinized. Here, we analyzed the genomes of four Cucurbita species, including a newly assembled genome of Cucurbita argyrosperma, and compared the gene contents of these species with those of five other members of the Cucurbitaceae family to assess the evolutionary dynamics of protein-coding and long intergenic noncoding RNA (lincRNA) genes after the genome duplication. We report that Cucurbita genomes have a higher protein-coding gene birth-death rate compared with the genomes of the other members of the Cucurbitaceae family. C. argyrosperma gene families associated with pollination and transmembrane transport had significantly faster evolutionary rates. lincRNA families showed high levels of gene turnover throughout the phylogeny, and 67.7% of the lincRNA families in Cucurbita showed evidence of birth from the neofunctionalization of previously existing protein-coding genes. Collectively, our results suggest that the whole-genome duplication in Cucurbita resulted in faster rates of gene family evolution through the neofunctionalization of duplicated genes.

Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits

Cucurbita pepo contains two cultivated subspecies, each of which encompasses four fruit-shape morphotypes (cultivar groups). The Pumpkin, Vegetable Marrow, Cocozelle, and Zucchini Groups are of subsp. pepo and the Acorn, Crookneck, Scallop, and Straightneck Groups are of subsp. ovifera. Recently, a de novo assembly of the C. pepo subsp. pepo Zucchini genome was published, providing insights into its evolution. To expand our knowledge of evolutionary processes within C. pepo and to identify variants associated with particular morphotypes, we performed whole-genome resequencing of seven of these eight C. pepo morphotypes. We report for the first time whole-genome resequencing of the four subsp. pepo (Pumpkin, Vegetable Marrow, Cocozelle, green Zucchini, and yellow Zucchini) morphotypes and three of the subsp. ovifera (Acorn, Crookneck, and Scallop) morphotypes. A high-depth resequencing approach was followed, using the BGISEQ-500 platform that enables the identification of rare variants, with an average of 33.5X. Approximately 94.5% of the clean reads were mapped against the reference Zucchini genome. In total, 3,823,977 high confidence single-nucleotide polymorphisms (SNPs) were identified. Within each accession, SNPs varied from 636,918 in green Zucchini to 2,656,513 in Crookneck, and were distributed homogeneously along the chromosomes. Clear differences between subspecies pepo and ovifera in genetic variation and linkage disequilibrium are highlighted. In fact, comparison between subspecies pepo and ovifera indicated 5710 genes (22.5%) with Fst > 0.80 and 1059 genes (4.1%) with Fst = 1.00 as potential candidate genes that were fixed during the independent evolution and domestication of the two subspecies. Linkage disequilibrium was greater in subsp. ovifera than in subsp. pepo, perhaps reflective of the earlier differentiation of morphotypes within subsp. ovifera. Some morphotype-specific genes have been localized. Our results offer new clues that may provide an improved understanding of the underlying genomic regions involved in the independent evolution and domestication of the two subspecies. Comparisons among SNPs unique to particular subspecies or morphotypes may provide candidate genes responsible for traits of high economic importance.