# Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin

Priscila Ambrósio Moreira, Cédric Mariac, Leila Zekraoui, Marie Couderc, Doriane Picanço Rodrigues, Charles Roland Clément, Yves Vigouroux

<figure><img src="image_1.png" /></figure>

## To cite this version:

Priscila Ambrósio Moreira, Cédric Mariac, Leila Zekraoui, Marie Couderc, Doriane Picanço Rodrigues, et al.. Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin. Evolutionary Applications, 2017, 10 (6), pp.577-589. 10.1111/eva.12474. hal-03340254

HAL Id: hal-03340254

https://hal.umontpellier.fr/hal-03340254v1

Submitted on 16 Sep 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin

Article in Evolutionary Applications · March 2017

DOI: 10.1111/eva.12474

CITATIONS

9

READS

247

7 authors, including:

<figure><img src="image_3.png" /></figure>

Priscila Ambrósio Moreira

Instituto Nacional de Pesquisas da Amazônia

18 PUBLICATIONS

207 CITATIONS

SEE PROFILE

<figure><img src="image_4.png" /></figure>

Mariac Cedric

Institute of Research for Development

129 PUBLICATIONS

1,965 CITATIONS

SEE PROFILE

<figure><img src="image_5.png" /></figure>

Leila Zekraoui

Institute of Research for Development

55 PUBLICATIONS

443 CITATIONS

SEE PROFILE

<figure><img src="image_6.png" /></figure>

Doriane Picanço-Rodrigues

Federal University of Amazonas

39 PUBLICATIONS

560 CITATIONS

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

<figure><img src="image_7.png" /></figure>

Comparative Study of Heavy Metal Contamination at Common Biomedical Waste Treatment and Disposal Sites (Incineration and Deep Burial) in Mumbai, Maharashtra, India
View project

<figure><img src="image_8.png" /></figure>

Studying the adaptive potential of Ugandan Coffea canephora (Robusta Coffee) diversity to drought stress
View project
# Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin

Priscila Ambrósio Moreira<fnref n="1" /> ID

Doriane Picanço Rodrigues<fnref n="1" />,<fnref n="3" />

Cédric Mariac<fnref n="2" />

Charles R. Clement<fnref n="1" />,<fnref n="4" />

Leila Zekraoui<fnref n="2" />

Yves Vigouroux<fnref n="2" />

Marie Couderc<fnref n="2" />

### Correspondence

Priscila Ambrósio Moreira, Post-Graduate Program in Botany, Instituto Nacional de Pesquisas da Amazônia (INPA), Manaus, Amazonas, Brazil.
Email: pri.ambrosio@hotmail.com

### Funding information

Conselho Nacional de Desenvolvimento Científico e Tecnológico, Grant/Award Number: CNPq-473422/2012-3; Fundação de Apoio à Pesquisa do Estado do Amazonas, Grant/Award Number: FAPEAM 062.03.137/2012; Agence Nationale de la Recherche, Grant/Award Number: ANR-13-BVS7-0017; Agropolis Fondation; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Grant/Award Number: CAPES-99999.010075/2014-03

### Abstract

Local people's perceptions of cultivated and wild agrobiodiversity, as well as their management of hybridization are still understudied in Amazonia. Here we analyze domesticated treegourd (Crescentia cujete), whose versatile fruits have technological, symbolic, and medicinal uses. A wild relative (C. amazonica) of the cultivated species grows spontaneously in Amazonian flooded forests. We demonstrated, using whole chloroplast sequences and nuclear microsatellites, that the two species are strongly differentiated. Nonetheless, they hybridize readily throughout Amazonia and the proportions of admixture correlate with fruit size variation of cultivated trees. New morphotypes arise from hybridization, which are recognized by people and named as local varieties. Small hybrid fruits are used to make the important symbolic rattle (maracá), suggesting that management of hybrid trees is an ancient human practice in Amazonia. Effective conservation of Amazonian agrobiodiversity needs to incorporate this interaction between wild and cultivated populations that is managed by smallholder families. Beyond treegourd, our study clearly shows that hybridization plays an important role in tree crop phenotypic diversification and that the integration of molecular analyses and farmers' perceptions of diversity help disentangle crop domestication history.

### KEYWORDS

agrobiodiversity, Crescentia amazonica, Crescentia cujete, introgression, plant domestication, wetlands

## 1 | INTRODUCTION

Amazonia is an important center of plant domestication (Clement, 1999; Meyer, Duval, & Jensen, 2012). Its great biological and cultural diversity (Balée, 2013) make it an especially interesting area to study the role of human societies in plant domestication and diversification (Balée, 2013; Clement, de Cristo-Araújo, D'Eeckenbrugge, Pereira, & Picanço-Rodrigues, 2010). The distinction between wild and cultivated is one of the basic questions of plant domestication (Lévi-Strauss, 1950; Pickersgill, 2013; Terrell et al., 2003). The distinction,

however, often goes unnoticed, given the lack of understanding of how local people perceive biological diversity in traditional societies (Caillon & Degeorges, 2007). People's perceptions of cultivated and wild diversity, as well as their management practices that deal with plant hybridization, are still understudied in Amazonia, especially for tree species (Moreira, Lins, Dequigiovanni, Veasey, & Clement, 2015; Rollo et al., 2016; Smith & Fausto, 2016). Hybridization between related cultivated and wild plants may be favored or discouraged by local farmers (Jarvis & Hodgkin, 1999). It can promote domestication and diversification (Gompert & Buerkle, 2016; Miller & Gross, 2011),

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2017 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd

<footnote n="1">¹Post-Graduate Program in Botany, Instituto Nacional de Pesquisas da Amazônia (INPA), Manaus, Amazonas, Brazil</footnote>

<footnote n="2">²Institut de Recherche pour le Développement, Université de Montpellier (IRD), UMR DIADE, Montpellier, France</footnote>

<footnote n="3">³Laboratório de Evolução Aplicada, Universidade Federal do Amazonas (UFAM), Manaus, Amazonas, Brazil</footnote>

<footnote n="4">⁴Coordenação de Tecnologia e Inovação, INPA, Manaus, Amazonas, Brazil</footnote>
because hybrids often present interesting traits that can be selected and maintained (Goldschmidt, 2013; Miller & Gross, 2011; Zohary & Spiegel-Roy, 1975), but may not be adaptive in natural environments (Ellstrand, 2003). This is especially true of introgressive hybrids, as backcrossing to one parent maintains its useful characteristics with minor influence of the other parent (Ellstrand, 2003; Harrison & Larson, 2014). It follows that hybridization and introgression between cultivated and wild plants, as well as the human practices that maintain diversity, are important for effective agrobiodiversity conservation. The distinction between wild and cultivated and its linkages with natural ecosystems are essential for a broader understanding of agriculture (Aumeeruddy-Thomas, Hmimsa, Ater, & Khadari, 2014). More efforts are necessary for its recognition and implementation by public agricultural and forestry policies (Michon, Nasi, & Balent, 2013).

Crescentia spp. (Bignoniaceae) are excellent candidates to study hybridization and domestication associated with floodplains in Amazonia. *Crescentia cujete* Linnaeus (1753), known as treegourd or calabash tree, is an important tree crop for Amazonian smallholders (Lima & Saragoussi, 2000; Wittmann & Wittmann, 2011). Its versatile fruits, called *cuia* in Portuguese, are traditionally used as storage vessels, drinking cups, scoops to bail water from canoes, traps for fishing, diving masks, bags, body ornaments, ritualistic musical instruments, and, more recently, as “ecological” cups; they also have medicinal applications (Acostupa, Bardales, & Teco, 2013; Bennett, 1992; Bustamante, Hidalgo, & Frausin, 2011; Heiser, 1993; Morton, 1968; Patiño, 1967; Steward, 1948). *C. cujete* presents an ample variation in fruit shapes and sizes (Aguirre-Dugua, Eguiarte, González-Rodríguez, & Casas, 2012; Arango-Ulloa, Bohorquez, Duque, & Maass, 2009; Gentry, 1980) that support the wide range of uses. A wild relative (*C. amazonica* Ducke 1937) occurs in flooded forests in the Orinoco and Amazon Basins, as well as smaller rivers of the Guianas (Díaz, 2009; Gentry, 1980; Godoy, Petts, & Salo, 1999; Wittmann et al., 2006). *Crescentia amazonica* fruits are smaller with thinner rinds that float in the water and are dispersed by fish (Waldhoff, Ulrich, & Furch, 1996).

The relationship between the two species is largely speculative. *Crescentia amazonica* was hypothesized to be the wild progenitor from which treegourd was domesticated (Ducke, 1946). Alternatively, treegourd was domesticated in Mesoamerica and later distributed to Amazonia (Gentry, 1980). In this case, it was hypothesized that *C. amazonica* was feralized *C. cujete* (Gentry, 1980). Their chloroplast genetic diversity does not support the possibility of *C. amazonica* being derived from *C. cujete* (Moreira et al., 2016). Nonetheless, a close relationship between *C. cujete* and *C. amazonica* is recognized by local human populations. In Guyana, people “called spirits” when *C. amazonica* was found, because they recognized it is a sort of “shadow” of the domesticated *C. cujete* (van Andel, 2000). As the two species co-occur in Amazonia and *Crescentia* species are hypothesized to be largely interfertile (Gentry, 1980), gene flow might be abundant, but has not yet been shown at the molecular level.

In this study, we asked whether (i) hybridization plays a significant role in shaping genetic and morphologic diversity in *Crescentia* species, and whether (ii) hybrid and introgressed individuals are managed

by Amazonian smallholders. To address these questions, we combined a genetic study based on chloroplast (single nucleotide polymorphisms—SNP) and nuclear (simple sequence repeats—SSR) markers with local farmer interviews. Using these datasets, we analyzed (i) the genetic differences between *C. amazonica* and *C. cujete* in the Amazon Basin, and gene flow between them; and (ii) the relationship between genetic and morphological diversity, and how people use and perceive this diversity.

## 2 | METHODS

### 2.1 | Field sampling and interviews

We visited rural and peri-urban villages in 36 municipalities distributed along the major rivers of the Brazilian Amazon Basin (Figure 1). The broad geographical sampling followed two criteria for village selection: dependence on river resources and treegourd use in daily life. Data were collected after an informed consent invitation that was read collectively in each village and signed by a local representative. This research followed the International Society for Ethnobiology's code of ethics (International Society of Ethnobiology, 2006) and was approved by the Committee for Ethics in Research with Human Beings of the National Research Institute for Amazonia (CEP INPA, proc. no. 408.611, 2013). We collected leaves for genetic analyses of each treegourd (N = 469) found in domestic areas in the villages (Table S1). We considered domestic areas to include homegardens, swiddens, ports near the river, football fields, trails, and old homestead sites. Clones propagated from the same individual were avoided in order to better assess available diversity. Fruits were measured and photographed and their shape was categorized according to Arango-Ulloa et al. (2009). We performed semi-structured interviews about use, management, and history of most trees collected, and in nine municipalities we practiced participant observation of daily activities and urban farmers' markets. To map wild treegourd distribution, we surveyed herbarium records of *C. amazonica*, and photographs of its fruit were presented to farmers in all villages visited to stimulate their memory of its presence in local flooded forests. We collected leaves of *C. amazonica* (N = 32) growing spontaneously in seven areas of flooded forests in six municipalities along the Solimões-Amazonas and lower Madeira rivers (Figure 1). Collection was authorized by the Brazilian System for Authorization and Information in Biodiversity, Chico Mendes Institute for Biodiversity Conservation, proc. no. 25052-1, 2012, and transportation by the Brazilian Institute for the Environment and Renewable Natural Resources, proc. no. 14BR015576/DF, 2014.

### 2.2 | Genetic analysis

DNA was extracted from dried leaves (N = 234) using the CTAB 5% protocol (Doyle & Doyle, 1990) with minor modifications. For nuclear SSR analysis, we genotyped all 234 samples, among which 184 were also analyzed for chloroplast SNPs (Table S1). Collected leaves without enough DNA or with low quality were excluded. To develop the
<figure><img src="image_9.png" /><figcaption>**FIGURE 1** Geographical distribution of *Crescentia cujete* and *Crescentia amazonica* treegourds collected for this study along major rivers of Brazilian Amazonia (N = 234). Proportions of admixture identified by Structure at K = 2 are indicated. Proportions above 0.9 were considered pure and below 0.9 were classified as admixed. *Crescentia amazonica* records in northern South America are from the Global Biodiversity Information Facility</figcaption></figure>

nuclear SSRs, a barcoded library of *C. cujete* genomic DNA was sequenced (Moreira et al., 2016) using an Illumina MiSeqv.3 (San Diego, California, USA). QDD software 3.1.2 (Meglécz et al., 2009) was used to identify nuclear SSR motifs and design primers from 113,865 merged reads. The parameters used to select the nSSR primers were as follows: one primer pair for each read to avoid repeated regions of the genome; avoidance of mononucleotide microsatellites; preference for perfect microsatellites with ≥8 repeats; and avoidance of primers that are very close (≤20 bp) to the target SSR. A total of 1,436 SSR were identified, of which 1,068 were perfect SSR motifs with 819 di-, 191 tri-, 47 tetra-, 10 penta-, and 1 hexa-repeat motifs, and 368 compound motifs. The primers were designed for the perfect SSRs (Table S2). We performed preliminary tests of amplification of 15 SSRs using *C. cujete* samples (N = 3). Five of the primer pairs failed to amplify, even using different temperatures and DNA concentrations, and were discarded. The remaining 10 SSR primers were labeled with fluorescence (FAM, NED, HEX; Applied Biosystems, Foster City, California, USA) and genotyped in an ABI 3130xL Genetic Analyzer (Applied Biosystems) using GS-500 LIZ as the size standard (Applied Biosystems). Although all 10 SSR were polymorphic, at least for cultivated samples (N = 221), two

loci were excluded (SSR2 and SSR9) because they failed to amplify in 60% of the samples. We kept the remaining eight SSRs (Table S3) for hybridization analysis, as they also amplified successfully for *C. amazonica*. Locus amplifications were made in simplex and multiplexed for fragment analysis, using the PCR kit (Qiagen, n.206143, Hilden, Germany) with the following program: 95°C for 15 min; 38 cycles, each of 94°C for 30 s, T~a~ °C for 1.30 min, and 72°C for 1 min; and a final step of 60°C for 30 min. Fragment size and allele identification were determined using GeneMapper (Applied Biosystems). We also obtained SNPs observed in the whole chloroplast using previously described approaches (Moreira et al., 2016; Scarcelli et al., 2016).

## 2.3 | Diversity analyses

The eight nuclear SSRs were used to genotype 234 treegourds, among which 221 were from domestic areas and 13 from flooded forests (Table 1). Relationships among individuals were assessed with Structure 2.3, using the admixture model (Pritchard, Stephens, & Donnelly, 2000). We varied the number of genetic clusters (K) from K = 1 to 20, with 100.000 burn-in, 100.000 iterations, and five different runs for each K
<table><caption>**TABLE 1** Summary of nuclear (SSR) admixture proportions (N = 234) and chloroplast (SNPs) haplotypes (N = 184) of the *Crescentia cujete* and *Crescentia amazonica* collections analyzed in this study and the habitats they were collected in. A. Admixture proportions for *C. cujete* and *C. amazonica* in columns and chloroplast haplotypes in lines. Not confirmed means the chloroplast was not analyzed in plants that were genotyped with nSSR. B. Habitats in which pure and admixed treegourds were collected</caption><thead><tr><th></th><th></th><th></th><th>Pure cujete</th><th>Cujete admixed</th><th>Hybrids</th><th>Amazonica admixed</th><th>Pure amazonica</th></tr><tr><th></th><th>nSSR (N = 234)</th><th>N</th><th>175</th><th>25</th><th>11</th><th>5</th><th>18</th></tr></thead><tbody><tr><td rowspan="3">A)</td><td>Haplotypes<br>(N = 184)</td><td>C. cujete</td><td>170</td><td>142</td><td>16</td><td>8</td><td>1</td><td>3</td></tr><tr><td></td><td>C. amazonica</td><td>14</td><td>0</td><td>0</td><td>1</td><td>2</td><td>11<sup>a</sup></td></tr><tr><td></td><td>Not confirmed</td><td>50</td><td>33</td><td>9</td><td>2</td><td>2</td><td>4</td></tr><tr><td rowspan="2">B)</td><td>Domestic areas (N = 221)</td><td></td><td>175</td><td>25</td><td>11</td><td>5</td><td>5</td></tr><tr><td>Flooded forest (N = 13)</td><td></td><td>0</td><td>0</td><td>0</td><td>0</td><td>13<sup>b</sup></td></tr></tbody></table>

value. The ad hoc ΔK (Evanno, Regnaut, & Goudet, 2005) was used to identify the most likely number of clusters in the matrix. We considered an individual to belong to a given cluster if its proportion of admixture was <0.10, that is, with more than 0.90 of the individual's SSR profile attributable to the given cluster. In this study, we used the term "admixed" for ancestry between 0.90 and 0.60, and the term "hybrid" for ancestry from 0.60 to 0.40 (Table S1). We also calculated a hybridization index using Introgress 1.2.3 (Gompert & Buerkle, 2010) and compared it with the Structure admixture proportions. The whole chloroplast sequences of 174 domestic treegourds and 10 from flooded forests were analyzed (Table S1). We built a haplotype network based on 250 SNPs using the median joining algorithm (Bandelt, Forster, & Röhl, 1999). The network was visualized using the samples with <4% of missing data, according to software requirements in POPART 1.7 (Leigh & Bryant, 2015). The nuclear and chloroplast comparison defined paternal and maternal introgression, respectively, and determined the final botanical identification (170 *C. cujete* and 14 *C. amazonica*). To assess the impact of the uneven number of samples of the two species, we performed a complementary Structure analysis using the same sample size for both species (14 *C. amazonica*, 14 *C. cujete*; Fig. S1), with the *C. cujete* samples chosen at random. Finally, we assessed the potential impact of null alleles on admixture inferences (Fig. S2), by coding any missing data as a homozygote recessive allele (Falush, Stephens, & Pritchard, 2007).

Nuclear genetic diversity of *C. cujete* and *C. amazonica* was explored with principal components analysis (PCA) executed with stats R package (R Core Team, 2015). Nuclear genetic diversity and species differentiation were estimated using hierfstat (Goudet, 2005). Pairwise $F_{ST}$ were calculated and statistically assessed using 1,000 bootstraps (Nei, 1987). The significance of $F_{IS}$ was measured as deviation from Hardy-Weinberg equilibrium using pegas R (Paradis, 2010). Chloroplast diversity was estimated using DNAsp 5.10.1 (Librado & Rozas, 2009) and Arlequin 3.5 (Excoffier & Lischer, 2010). We examined the relationship between *C. amazonica* admixture proportions and *C. cujete* fruit diameters (N = 61) with simple regression in R package (R Core Team, 2015).

# 3 | RESULTS

## 3.1 | Genetic structure revealed with nuclear SSRs

The Structure analyses identified two clusters as the most likely structure (K = 2, Figure 2a), based on the ad hoc ΔK approach (Evanno et al., 2005; Fig. S1). The two clusters correspond to domesticated *C. cujete* and wild *C. amazonica*, the two accepted botanical species already described in the Amazon Basin. We also found a significant amount of admixture between them (Figure 2a, Table 1), especially along the Solimões-Amazonas River (Figures 1 and 2b). The frequency of individuals with admixture proportions suggesting hybridization and introgression was similar in both groups. Among the 234 samples, 200 were assigned primarily to *C. cujete*, of which 24 (12%) were admixed to some degree, and 23 were assigned to *C. amazonica*, of which five (21%) were admixed (Table 1). Eleven plants presented admixture proportions between 40% and 60% and were classified as hybrids. These proportions were similar when the *C. cujete* sample size was reduced at random to be equal to the *C. amazonica* sample size ($r^2$ = .99, $p$ < 10<sup>-15</sup>, Fig. S1). At the second possible grouping (K = 3, Fig. S3), similar admixture proportions were also found ($r^2$ = .98, $p$ < 10<sup>-15</sup>), while the original *C. cujete* group was subdivided without relation to geography or fruit morphology (Fig. S3). The hybridization index calculated with Introgress was highly correlated with the admixture proportions calculated with Structure (Fig. S4, $r^2$ = .83, $p$ < 10<sup>-15</sup>). Finally, coding all missing data as homozygous recessive alleles did not have an impact on admixture inferences (Fig. S2, $r^2$ = .96, $p$ < 10<sup>-15</sup>).

## 3.2 | Comparison between chloroplast SNPs and nuclear SSR diversity

The chloroplast analyses also clearly identified the two botanical species (Figure 3), with 63 mutational differences between the *C. cujete* and *C. amazonica* chloroplast sequences. In the haplotype network,

<footnote><sup>a</sup>Crescentia amazonica with "pure amazonica" assignment in Structure at K = 2 is predominant in flooded forests (N = 10), but one was found cultivated (N = 1). Admixed *C. amazonica* were only found in domestic areas.</footnote>

<footnote><sup>b</sup>All of these are likely to be *C. amazonica*, but three were not analyzed for their chloroplast haplotypes.</footnote>
<figure><img src="image_10.png" /><figcaption>**FIGURE 2** Structure analysis of 234 treegourd samples collected in Brazilian Amazonia. The y-axis shows the proportion of assignment to the groups at K = 2 (red—Crescentia *cujete* and blue—Crescentia *amazonica*). (a) Samples are ordered by their proportion of admixture: admixed if >0.1, hybrids if 0.4 to 0.6, pure if >0.9 of assignment to the group. (b) Samples were ordered by their geographical location along the main rivers: the Negro, Solimões, and Amazonas rivers are ordered west to east; the Branco River is ordered north to south; the Madeira River is ordered south to north. (c) Samples are ordered by seven fruit shapes (see Figure 4) and fruit size, with size increasing from left to right</figcaption></figure>

samples with *C. cujete* nuclear assignments showed exclusively *C. cujete* chloroplast haplotypes (Figure 3, Table 1). For individuals with *C. amazonica* nuclear assignments (N = 23), 77% had *C. amazonica* haplotypes (N = 13) and 17% had *C. cujete* haplotypes (N = 4). Among the hybrids, most had *C. cujete* haplotypes (88%, N = 8), and one had a *C. amazonica* haplotype. The flooded forests harbored exclusively pure *C. amazonica* samples (N = 10), while domestic areas harbored

pure and admixed samples of both species (Table 1, Figure 5c). The genetic diversities of the pure samples of both species were lower than their admixed samples (Table S4). The rarified allele count ($A_r$) varied from 1.1 to 5.4, with the lowest value in *C. amazonica*, slightly higher within its admixed samples and the highest value in admixed *C. cujete*. The mean expected heterozygosity ($H_s$) of pure *C. cujete* was 0.31 and that of pure *C. amazonica* was 0.09, while their admixed
<figure><img src="image_11.png" /><figcaption>**FIGURE 3** Chloroplast haplotype median joining network of *Crescentia cujete* and *Crescentia amazonica* from Brazilian Amazonia based on 250 chloroplast SNPs. Nuclear ancestry of each sample (N = 184) was evaluated using Structure (Table 1). The size of the circle reflects the number of individuals presenting the same haplotype. Numbers of mutations are indicated as hatch marks and numbers between haplotypes. Colors represent percentage of nuclear admixture, where *C. amazonica* >0.9 (dark blue); admixed *C. amazonica* <0.9 and >0.6 (light blue); hybrids <0.6 and >0.4 (violet); admixed *C. cujete* <0.9 and >0.6 (orange); and *C. cujete* >0.9 (red). A clear chloroplast difference is observed between *C. amazonica* and *C. cujete*, although 17% (N = 4) of the *C. amazonica* ancestry samples (N = 23, Table 1) have *C. cujete* haplotypes</figcaption></figure>

samples showed mean *H*<sub>s</sub> values of 0.58 and 0.33, respectively (Table S4). *C. amazonica* showed extremely low diversity, with fixed alleles at six loci. These low diversity values for *C. amazonica* are probably due to marker development from *C. cujete* with consequent poor transfe- ral due to the significant divergence between species. Although such bias might lead to imprecise estimation of *C. amazonica* diversity, it does not have an impact on the identification of hybrids, as admix- ture is based on allele frequency differences and not diversity per se (Pritchard et al., 2000).

Differentiation between *C. cujete* and *C. amazonica* was very high (*F*<sub>ST</sub> = 0.74, IC<sub>95%</sub> = 0.59-0.80 excluding admixed samples; *F*<sub>ST</sub> = 0.56, IC<sub>95%</sub> = 0.41-0.66 with all samples). *F*<sub>IS</sub> (Table S4) was not significant for *C. cujete*, but significant for *C. amazonica* (*F*<sub>IS</sub> = 0.44, *p* < .05). In the chloroplasts, we found 92 SNPs in 14 individuals of *C. amazonica* with 14 haplotypes, and 178 SNPs in 170 individuals of *C. cujete* with 71 haplotypes. Nucleotide diversity (π) was 3 × 10<sup>-3</sup> and 1.2 × 10<sup>-3</sup>, for *C. amazonica* and *C. cujete*, respectively. Chloroplast differentiation between species was high (*F*<sub>ST</sub> = 0.89, *p* < .05).
<figure><img src="image_12.png" /><figcaption>**FIGURE 4** Diversity of fruit shapes and sizes in Amazonian treegourds. Seven fruit shapes of *Crescentia cujete* and one shape of *Crescentia amazonica* were found in the Brazilian Amazon Basin. Shape classification follows Arango-Ulloa et al. (2009). Smaller fruits of types 1 and 2 have local names that are shown in italics in the corresponding photograph. Box plots of the variation in diameter (cm) of each fruit shape, with domestic (●) and wild individuals (▲). The colors represent chloroplast haplotypes (red—C. cujete; blue—C. amazonica; gray—not confirmed). In fruit shape 2, fruits smaller than 10 cm can be C. amazonica or admixed C. cujete, in which case they correspond to the *maraca* fruit type</figcaption></figure>

## 3.3 | Admixture proportions correlated with fruit size

There is ample morphological diversity of treegourd fruits in the Brazilian Amazon Basin (Figure 4). We recorded seven types of *C. cujete* fruit shapes in domestic areas (N = 167) and one fruit type of *C. amazonica* in flooded forests (N = 10). In domestic areas, two of them account for 86.4% of the plants: 63% flattened (type 1) and 23.4% oblong (type 2). The other five types were rare: 0.6% was cuneate (type 3), 3.6% elongated (type 4), 1.8% globular (type 5), 4.8% rounded-drop (type 6), and 1.2% oblong-drop (type 7). The two *C. amazonica* samples found in domestic areas showed type 2 fruits. While the fruits of *C. amazonica* in the flooded forest have small diameters ranging from 2.5 to 5.4 cm (median: 4.85 cm), in domestic areas diameters were slightly larger, ranging from 4 to 8.4 cm (median: 6.2 cm) (Figure 4). Among *C. cujete*, there is great size variation, especially within types 1, 2, and 5, with variation from 5.2 to 29 cm. Smaller fruits of *C. cujete* were correlated with *C. amazonica* admixture ($r^2$ = .34, *p* = 1.26 × 10⁻⁶, Figure 5), especially the oblong (type 2) fruit ($r^2$ = .62, *p* = 4.81 × 10⁻⁶). This relation remains even if we exclude the most extreme admixed samples (wild ancestry higher than 0.8, $r^2$ = .23, *p* = 1.3 × 10⁻⁴). The smaller *C. cujete* fruits,

between 3 and 10 cm in diameter, are those with the highest admixture proportions (Figure 5). Some *C. amazonica* ancestry could be observed in all shapes found in domestic areas, except in the rare type 7 (Figure 2c).

## 3.4 | Use, management and perception of fruit diversity

We recorded 11 current domestic uses of *C. cujete* fruits that are related to fruit morphology and were categorized in levels of specialization. Five uses are related with fruit size, without restriction of shape. The small gourds that are mainly used for food consumption (*xibé*, a meal of water and manioc flour), *acai* (juice from *Euterpe oleracea* or *E. precatoria*), water, or some particular medicine were reported in 31% of the municipalities. In the upper Negro River, these are called *cuiupi*. Large gourds are preferred for use as baskets to store seeds, seedlings, and cleaning products or to serve food (fish, *tapioca*), and were reported in 66% of the municipalities. Large gourds are also used as a unit of measurement of manioc during the preparation of flour and were reported in 47% of the municipalities, for the daily bath in the river (39%), to bail water from inside the canoe (39%), and pulp of broken fallen fruits as fodder for livestock (16%).
<figure><img src="image_13.png" /><figcaption>**FIGURE 5** Admixture proportions of treegourds based on eight nSSR and its relationship with fruit size and propagation method, and the dispersal of diversity in the principal component analysis. (a) Fruit diameter of genetically confirmed *Crescentia cujete* samples (N = 61) as a function of the proportion of admixture with *Crescentia amazonica*, where diameter = 17.9–2.07*admixture* (r² = .34, *p* = 1.26 × 10⁻⁶). (b) Admixture proportions of *C. cujete* (N = 80) samples propagated spontaneously by seeds and cultivated by stems and by seeds. Black dots indicate the average size and white circles are outliers that highlight the active maintenance of admixed individuals by people. (c) Principal components analysis of the genetic relationships among *Crescentia cujete* (N = 170; red (●)) and *C. amazonica* (N = 14; blue (▲)) individuals. The proportion of the variance explained by each PC is shown in parentheses along each axis. The gradient of admixture is visible along PC1 and the admixed individuals correspond to the smaller fruit size varieties *paranã* (flattened type), *maracá*, and *cuiupi* (oblong types) found exclusively in domestic areas</figcaption></figure>

Four uses require greater specialization because they need a com-
bination of size and fruit shape. Larger (20–24 cm) oblong fruits (type
2) are used to make a kind of bag with fiber handles (called coió, com-
boró, mocó), recorded in 30% of the municipalities, and used to trans-
port things to the swiddens or to fish. Flattened fruits (type 1) with
relatively smaller sizes (13–15 cm) are used in different contexts. They
are involved in *tacacá* commerce, a typical Amazonian soup served in
treegourd bowls (called *paranã*), whose manufacture was recorded in
11% of the municipalities. They are also used to make bowls to pre-
pare blessings (14 cm) during religious events, recorded in 8% of the
municipalities, and more rarely as parts of clothing sewn for regional
festivals (2.8% of the municipalities).

Finally, there is an opportunistic use of fruit diversity, recorded in
16% of the municipalities, which incorporates different fruit shapes
and sizes into the repertoire of manufactured objects sold as hand-
icrafts. In this repertoire, even the wild species, *C. amazonica*, with

its small fruits, is included, as recorded in four municipalities along
the Solimões-Amazonas River. Its fruits are designed as small cups,
painted and sold to cosmetic and food companies as "ecological cups,"
and are worth U$ 6 per 100 cups for handicrafters (price in 2014).
Other wild fruit type handicraft designs include the reproduction of
ancient musical instruments called *maracá*, and also as toys for chil-
dren (*carrapeta*). The use of these artifacts was recorded in domestic
contexts in 8% of the municipalities, all along the Solimões-Amazonas
River, where admixture is relatively frequent and with higher levels
(Figure 2b).

While large fruits have a general name (*cuia* and *coité*) with vari-
ations that include shape information (e.g., long gourd), the smaller
fruits are distinguished with different local names that are related to
flooded environments: *cuia do igapó*, *cuiupi*, *paranã*, and *maracá*. “*Cuia
do igapó*” is the wild tree that grows spontaneously in the *igapó*, a local
word for a forest environment that is periodically flooded by black
or clear water (Irion, Junk, & Mello, 1997), and has smaller size fruit ("it is the same as cuia, but is small, is from the forest, is not planted"). The other smaller gourds mentioned depend on human intervention to survive and are found exclusively in domestic contexts. *Cuiupi* was cited along the Madeira River as similar to the wild type ("there is cuiupi in the lakes, near the river, but I prefer the big one"), while along the upper Negro River it is the general term for small gourds. *Paranã* is a local word to describe secondary river channels, often between the edge of the forest and the white-water floodplain system rich in nutrients and sediments (Irion et al., 1997), the habitat of *C. amazonica*. *Maracá* develops increased fruit size in domestic contexts ("maracá is the cuia do igapó, but it is bigger") and its fruit production is restricted to flooding season ("the fruit fails, does not give all the time"). Based on genetic analyses and people's designation of varieties, *cuia do igapó* refers to *C. amazonica*, while *cuiupi* and *paranã* are *C. cujete* with moderate admixture proportions. *Maracá* is a special case, as it is applied not only to admixed *C. cujete*, but also to *C. amazonica*.

Cultivation by seeds, but also by stems, maintains highly admixed samples (Figure 5b), as many desirable gourds have smaller sizes, which is partly due to admixture effects (Figure 2c). Most treegourds are propagated by stem cuttings (61% of the sampled trees), and less frequently by seeds (37%). Stem cuttings have a purpose: They ensure a fast way to produce fruits ("by seeds takes more time") and the maintenance of fruit morphology, avoiding those fruits that break easily ("by stems it is better for the cuiupi not to become soft and break"). This is also the preferred method among traders who cultivate large clonal areas in upland environments for *tacacá* bowls, the flattened smaller gourd (type 1) with moderate levels of admixture called *paranã*. Seed propagation can be a spontaneous event when flood water brings seeds (28%) or spontaneously propagated through discarded fruit pulp near the house (71%). Seedlings, despite the lack of guarantee in fruit morphology, can survive floods better than cuttings in floodplain landscapes. People stated that: "it is not correct to plant seeds, but stems do not take here," "when it grows by seeds, from the pulp, it grows smaller." In the floodplain, stem cuttings require greater effort ("to make big cuia takes much work"), and success is not guaranteed ("it is not all the stems that work, here I've tried hard"). Thus, seed propagation is a way to deal with the high annual flooding events in the floodplain. Another motivation is that seedlings are more suitable to produce smaller fruit sizes that are also useful, especially for handicraft purposes, as observed in the social movement of craftswomen along the middle Amazonas River.

# 4 | DISCUSSION

## 4.1 | Identification of geographically widespread admixture

The tenuous nature of reproductive barriers among *Crescentia* species (Gentry, 1980) was confirmed by the large amount of admixture observed between *C. cujete* and *C. amazonica*. Our inference of admixture proportions was very robust even with unequal sample sizes and was correlated with the hybridization index of Gompert and Buerkle

(2010), even with our small set of nuclear makers. The main reason is certainly the very high differentiation between the two species ($F_{ST} = 0.74$ with nSSR; $F_{ST} = 0.89$ with cSNP). Although we could not rule out the existence of null nSSR alleles in *C. amazonica*, they do not influence the high differentiation observed between species. Different histories of gene flow might result in similar patterns of admixture proportions (Barton & Hewitt, 1985; Gompert & Buerkle, 2016). In this case, one question is whether these admixture proportions reveal hybridization after secondary contact or a long-term divergence. The high number of substitutions between *C. cujete* and *C. amazonica* chloroplast sequences suggests ancient divergence between the two species. As *C. cujete* is a species with cultivated populations, the origin of this domestication is unlikely to be older than other domestications in the Americas, which started by 11,000 years (Piperno, 2011). The chloroplast sequence divergence that we found suggests that the two species diverged earlier. Therefore, *C. cujete* was likely introduced by humans into South America (Gentry, 1980), and the admixture observed is secondary contact.

The hybrids are concentrated along the east-west axis of the Amazon Basin, the Solimões and Amazonas rivers (Figure 1). Although proximity to floodplains is an important parameter for the occurrence of hybrids, admixed individuals are also found beyond *C. amazonica*'s known distribution. One example is the occurrence of admixture along the Negro River, where no *C. amazonica* has been collected to date (Figures 1 and 2b). This pattern might result from social networks and propagule exchange of admixed plants by humans between riv- ers, such as between the Orinoco and Negro basins (Hornborg, 2005; Lathrap, 2010).

## 4.2 | Flooded forests are a source for cultivated treegourd phenotypic diversity

Fruit size variation of domesticated *C. cujete* in Amazonia is partly shaped by admixture between wild and cultivated plants. Note here that we do not have a common garden experiment to evaluate the fruit phenotype. Establishing such a common garden will be difficult because treegourd is of minor economic importance outside local communities, and consequently, there is no conservation or breeding in a research institution to allow working in already available common gardens. However, most of the cultivated plants measured shared a cultivated environment in homegardens. If size was simply associated with variability of the environment, we would not detect a significant association with introgression. The absence of a common environ- ment certainly adds more variability in size, but consequently also ren- ders significant correlations with admixture more difficult to detect.

As larger fruit size is an expected feature of tree domestication syn- dromes (Meyer et al., 2012; Miller & Gross, 2011), hybridization and introgression create variation in fruit size (Cornille, Giraud, Smulders, Roldán-Ruiz, & Gladieux, 2014) that can be managed (Aumeeruddy-Thomas et al., 2014; Cornille et al., 2012; García-Marin, Hernández-Xolocotzi, & Castillo, 1986; Hughes et al., 2007; Zerega, Ragone, & Motley, 2004). The great diversity of fruit sizes in Amazonian homegardens was also observed in the floodplains of the Orinoco
River and the Caribbean regions of Colombia (Arango-Ulloa et al., 2009). Similarly, in the Yucatán Peninsula of Mexico, large propagated fruits and spontaneous smaller fruits were reported in homegardens (Aguirre-Dugua et al., 2012). This suggests that perceptions of hybridization are used to manage fruit size and shape across the Neotropics.

Our results showed that pollen gene flow occurs in both direc-
tions between these *Crescentia* species. Bat pollination observed in
both species (Fleming, Geiselman, & Kress, 2009) certainly favors this
pollen flow between villages and flooded forests. However, the hy-
brid plants were restricted to human-managed areas. Hybrid and in-
trogressed *C. cujete/C. amazonica* are certainly selected against in the
flooded forest and favored in human areas, as adaptation to natural
environments is likely to be reduced by hybridization with domesti-
cated populations (Ellstrand, 2003).

4.3 | Traditional communities manage hybridization

Use of larger treegourds is widely distributed throughout the Amazon
Basin, while the use of smaller ones is reported in less than half of
the municipalities visited. Small fruits are much appreciated for the
manufacture of handicrafts. However, the record of domestic artifacts
made of small fruits suggests ancient use and not only contemporary
commercial handicraft demand. The historical use of smaller treegourd
fruits was mentioned by many ethnographers, especially in rituals and
for medicinal purposes. *Maracá* is a symbol in spiritual practices of
different Native Amazonian cultures, as well as an ancient rattle, and
*paranã*, the flattened round type, was traditionally used for bowls in
ceremonial rituals (Lévi-Strauss, 2004; Ribeiro, 1995; Steward, 1948).
Both of these objects made from smaller treegourds are considered
by oral histories as the first things to be in the world in different cul-
tures of Mesoamerica, the Antilles, and South America (Heiser, 1993),
such as the Taíno from the Dominican Republic (Martin, 1999), the
Guarani from southern Brazil (Montardo, 2002), and the Tukano from
Amazonian Brazil and Colombia (Hugh-Jones, 2009). These smaller
fruits come from admixed individuals (Figures 4 and 5), which are
maintained and dispersed by people. In the Brazilian Amazon Basin,
they are recognized as local varieties (*cuiupi*, *paranã*, *maracá*, and *cuia
do igapó*), whose names and associated traditional ecological knowl-
edge are explicitly related to *C. amazonica* habitat, or wild-cultivated
hybridization and its morphological consequences. The association
with flooded environments is also present in the Tupi origin of the
name *cuiupi* (from *kuy*´y) that refers to gourds of the water (Ferreira,
2004). The name *maracá* (*mbara'ka*) refers to the small rattles played in
order to talk to lakes and heal sick people (Andía, 2015). *Maracá* was
also used by Ducke (1946) when he described *C. amazonica* collected
in flooded forests of the Solimões River in 1937. All of this confirms
that the small fruited varieties used to manufacture important objects
in the Amazon Basin, such as *maracás*, are the result of human selec-
tion of hybrid and admixed trees, and highlight that treegourd diver-
sity is partly dependent on hybridization between homegarden and
flooded forest genepools.

Although there is an historical use of these admixed treegourds,
there is no evidence of ancient cultivation of *C. amazonica*. Our

documentation of cultivation of *C. amazonica* along the Solimões-
Amazonas River, especially in Santarém handicrafter villages, ap-
pears to be a recent practice. Treegourd handicrafts have been
famous since before the Colonial period (Medina, 1934; Patiño,
1967; Rodrigues Ferreira, 1933) and gained prominence recently as
a Brazilian Cultural Heritage (IPHAN, 2015). Handicrafts are moti-
vated not only by social and economic demands, as highlighted by
Santos (1982) and Carvalho (2011), but also ecological pressures, as
severe flooding in these areas influenced people to use seedlings as
an alternative way to produce treegourds for sale. As a result, we
observed a high frequency of admixture in these handicrafter villages
located in the middle Amazonas River (Figure 2b). The most common
cultivation practice of *C. cujete*, however, is vegetative propaga-
tion, not only in Amazonia, but also in Mexico (Aguirre-Dugua et al.,
2012). This is the traditional way to maintain useful fruit phenotypes,
a practice that allows management of hybridization also (Miller &
Gross, 2011). The admixed treegourds have been dispersed by hu-
mans along Amazonian rivers and potentially over larger geographical
areas, which might create a complex pattern of geographical admix-
ture, as observed for several other Neotropical fruit trees and their
wild populations [*Spondias* (Miller, 2008); *Inga* (Dawson et al., 2008);
*Chrysophyllum* (Petersen, Parker, & Potter, 2014)], as well as in the
Old World genus *Prunus* (Delplancke et al., 2012). Therefore, human
activity not only maintains, but promotes congener interaction, as
expected with other crops (Anderson, 2005; Riesenberg & Wendel,
1993). Hybrids are perceived and propagated in a dynamic way, so
that hybridization is managed according to people's needs. This is
in agreement with local farmers' practices and experimentation ob-
served worldwide (Brush, 2000; García-Marin et al., 1986; Hughes
et al., 2007; Jarvis & Hodgkin, 1999).

5 | CONCLUSIONS AND PERSPECTIVES

We provided evidence that variation in fruit size of the two *Crescentia* species found in the Brazilian Amazon Basin is related to their admixture proportions. New morphotypes that arise from hybridization are clearly recognized by people and named as local varieties (*maracá*, *cuiupi*, *paranã*), whose symbolism is emblematic for Amazonian cultures. Beyond treegourd, our study clearly shows that hybridization plays an important role in crop phenotypic diversification. We also showed that the integration of molecular analyses and farmers' perceptions of diversity can help disentangle crop domestication history. The specific traditional uses suggest that admixture management is an ancient human practice, also used in current traditional communities. We found that treegourd phenotype diversity depends partially on gene flow between homegardens and flooded forests. These results highlight the linkages between agriculture and forest ecosystems necessary for effective conservation of Amazonian agrobiodiversity. This is especially important as traditional ecological knowledge and floodplain conservation are neglected by development models for Amazonia (Castello et al., 2013; Posey & Balick, 2006).
ACKNOWLEDGEMENTS

This research was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-473422/2012-3), the Fundação de Apoio à Pesquisa do Estado do Amazonas (FAPEAM 062.03.137/2012), the Agence Nationale de la Recherche (ANR-13-BVS7-0017), and the ARCAD project funded by the Agropolis Fondation. P.A.M. thanks the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior for a scholarship (CAPES-99999.010075/2014-03). We thank the Instituto de Desenvolvimento Agrário do Amazonas for logistical support and farmer families for their support, kindness, and consent for this research.

DATA SHARING

1. Geographical coordinates are available in Table S1.

2. Vouchers of *C. amazonica* (Lat -4.32/Lon -59.71 and Lat -2.11/Lon -54.72) were deposited in the National Research Institute of Amazonia Herbarium (numbers 255.829 and 266.725).

3. The sequence alignments and microsatellite genotypes are available in Dryad doi: https://doi.org/10.5061/dryad.t84p3.

AUTHOR CONTRIBUTIONS

PAM, CRC, and YV planned the study. PAM carried out the field collections and interviews. PAM, LZ, MC, and DPR performed the molecular work. PAM, YV, and CM carried out the genetic analysis. PAM, CRC, and YV wrote the manuscript.

REFERENCES

Acostupa, R. J. H., Bardales, J. J. A., & Teco, R. M. V. (2013). Uso de las plantas medicinales en la comunidad El Chino del área de conservación regional comunal Tamshiyacu-Tahuayo, Loreto, Perú. *Conocimiento Amazónico*, 4(2), 77-86.

Aguirre-Dugua, X., Eguiarte, L. E., González-Rodríguez, A., & Casas, A. (2012). Round and large: morphological and genetic consequences of artificial selection on the gourd tree *Crescentia cujete* by the Maya of the Yucatan Peninsula, Mexico. *Annals of Botany*, 109(7), 1297-1306. doi:10.1093/aob/mcs068

Andel, T. R. van (2000). *Non-timber forest products of the north-west district of Guyana Part I*. Georgetown, Guyana: Tropenbos-Guyana Series.

Anderson, E. (2005). *Plants, man and life*. Mineola, NY: Dover.

Andía, J. J. R. (2015). El arte de hablar con los cerros: instrumentos musicales, entidades no humanas, cuerpos y géneros en los Andes peruanos septentrionales. In B. Brabec de Mori, M. Lewy, & M. A. García (Eds.), *Sudamérica y sus mundos audibles. Cosmologías y prácticas sonoras de los pueblos indígenas* (pp. 255-269). Berlin: Ibero-Amerikanisches Institut.

Arango-Ulloa, J., Bohorquez, A., Duque, M. C., & Maass, B. L. (2009). Diversity of the calabash tree (*Crescentia cujete* L.) in Colombia. *Agroforestry Systems*, 76(3), 543-553. doi:10.1007/s10457-009-9207-0

Aumeeruddy-Thomas, Y., Hmimsa, Y., Ater, M., & Khadari, B. (2014). Beyond the divide between wild and domesticated: Spatiality, domesticity and practices pertaining to fig (*Ficus carica* L.) and olive (*Olea europaea* L.) agroecosystems among Jbala communities in northern Morocco. In A. Chevalier, E. Marinova, & L. Peña-Chocarro (Eds.), *Plants and people: Choices and diversity through time* (pp. 191-210). Havertown, PA: Oxbow Books.

Balée, W. (2013). *Cultural forests of the Amazon: A historical ecology of people and their landscapes*. Tuscaloosa, AL: The University of Alabama Press.

Bandelt, H. J., Forster, P., & Röhl, A. (1999). Median-joining networks for inferring intraspecific phylogenies. *Molecular Biology and Evolution*, 16(1), 37-48. doi:10.1093/oxfordjournals.molbev.a026036

Barton, N. H., & Hewitt, G. M. (1985). Analysis of hybrid zones. *Annual Review of Ecology and Systematics*, 16, 113-148. Retrieved from http://www.jstor.org/stable/2097045.

Bennett, B. C. (1992). Plants and people of the amazonian rainforests. *BioScience*, 42(8), 599-607. Retrieved from http://www.jstor.org/stable/131192.

Brush, S. B. (2000). *Genes in the field: On-farm conservation of crop diversity*. Boca Raton, FL: Lewis.

Bustamante, G. G. F., Hidalgo, A. F., & Frausin, V. (2011). Cultura material na alimentação Saterê Mawé, Maués-AM, Brasil. In T. de J. P. Fraxe, A. C. Witkoski, & H. dos S. Pereira Pereira (Eds.), *Amazônia: Cultura material e imaterial* (pp. 225-246). São Paulo: Annablume.

Caillon, S., & Degeorges, P. (2007). Biodiversity: Negotiating the border between nature and culture. *Biodiversity & Conservation*, 26, 2919-2931. doi:10.1007/s10531-007-9149-7

Carvalho, L. G. de (2011). Artesanato e mudança social: Sobre projetos e comunidades em Santarém. In L. G. de Carvalho, L. A. Gennari, E. F. da Silva, & P. P. Peralta (Eds.), *O artesanato de cuias em perspectiva - Santarém* (pp. 19-46). Rio de Janeiro, RJ: IPHAN.

Castello, L., Mcgrath, D. G., Hess, L. L., Coe, M. T., Lefebvre, P. A., Petry, P., ...Arantes, C. C. (2013). The vulnerability of Amazon freshwater ecosystems. *Conservation Letters*, 6(4), 217-229. doi:10.1111/conl.12008

Clement, C. R. (1999). 1492 and the loss of amazonian crop genetic resources. I. The relation between domestication and human population decline. *Economic Botany*, 53(2), 188-202.

Clement, C. R., de Cristo-Araújo, M., D'Eckenbrugge, G. C., Pereira, A. A., & Picanço-Rodrigues, D. (2010). Origin and domestication of native amazonian crops. *Diversity*, 2, 72-106. doi:10.3390/d2010072

Cornille, A., Giraud, T., Smulders, M. J. M., Roldán-Ruiz, I., & Gladieux, P. (2014). The domestication and evolutionary ecology of apples. *Trends in Genetics*, 30(2), 57-65. doi:10.1016/j.tig.2013.10.002

Cornille, A., Gladieux, P., Smulders, M. J. M., Rolda, I., Cam, B. Le, Nersesyan, A., ... Giraud, T. (2012). New insight into the history of domesticated apple: Secondary contribution of the european wild apple to the genome of cultivated varieties. *PLoS Genetics*, 8(5), e1002703. doi:10.1371/journal.pgen.1002703

Dawson, I. K., Hollingsworth, P. M., Doyle, J. J., Kresovich, S., Weber, J. C., Sotelo Montes, C., ... Pennington, R. T. (2008). Origins and genetic conservation of tropical trees in agroforestry systems: A case study from the peruvian Amazon. *Conservation Genetics*, 9(2), 361-372. doi:10.1007/s10592-007-9348-5

Delplancke, M., Alvarez, N., Espindola, A., Joly, H., Benoit, L., Brouck, E., & Arrigo, N. (2012). Gene flow among wild and domesticated almond species: Insights from chloroplast and nuclear markers. *Evolutionary Applications*, 5(4), 317-329. doi:10.1111/j.1752-4571.2011.00223.x

Díaz, W. (2009). Composición florística de las comunidades vegetales aledañas al tercer puente sobre el río Orinoco, Venezuela. *Boletín Del Centro de Investigaciones Biológicas*, 43(3), 337-354.

Doyle, J., & Doyle, J. L. (1990). Isolation of plant DNA from fresh tissue. *Focus*, 12, 13-15.

Ducke, A. (1946). Plantas de cultura precolombiana na Amazônia brasileira. Notas sobre as espécies ou formas espontâneas que supostamente lhes teriam dado origem. *Boletim Técnico Do Instituto Agronômico Do Norte*, 8, 1-25.

Ellstrand, N. C. (2003). *Dangerous liaisons? When cultivated plants mate with their wild relatives*. Baltimore, MD: The Johns Hopkins University Press.

Evanno, G., Regnaut, S., & Goudet, J. (2005). Detecting the number of clusters of individuals using the software structure: A simulation study. *Molecular Ecology*, 14(8), 2611-2620. doi:10.1111/j.1365-294X.2005.02553.x
Excoffier, L., & Lischer, H. E. L. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. *Molecular Ecology Resources*, 10(3), 564–567.
doi:10.1111/j.1755-0998.2010.02847.x

Falush, D., Stephens, M., & Pritchard, J. K. (2007). Inference of population structure using multilocus genotype data: Dominant markers and null alleles. *Molecular Ecology Notes*, 7(4), 574–578.
doi:10.1111/j.1471-8286.2007.01758.x

Ferreira, M. C. (2004). Dicionário morfológico Tupi-Guarani, 2nd ed.. São Paulo: Edicon.

Fleming, T. H., Geiselman, C., & Kress, W. J. (2009). The evolution of bat pollination: A phylogenetic perspective. *Annals of Botany*, 104, 1017–1043.
doi:10.1093/aob/mcp197

García-Marin, P. C., Hernández-Xolocotzi, E., & Castillo, A. C. (1986). Variacion morfologica, manejo agricola y grados de domesticacion de *Opuntia* spp. en el bajio guanajuatense. *Agrociencia*, 65, 7-45.

Gentry, A. H. (1980). *Bignoniaceae - Part 1 (Crescentieae and Tourrettieae)*. *Flora Neotropica*. Bronx, NY: New York Botanical Garden Press.

Godoy, J. R., Petts, G., & Salo, J. (1999). Riparian flooded forests of the Orinoco and Amazon basins: A comparative review. *Biodiversity & Conservation*, 8(4), 551–586.
doi:10.1023/A:1008846531941

Goldschmidt, E. E. (2013). The evolution of fruit tree productivity: A review. *Economic Bota*, 67(1), 51–62.

Gompert, Z., & Buerkle, C. A. (2010). INTROGRESS: A software package for mapping components of isolation in hybrids. *Molecular Ecology Resources*, 10, 378–384.
doi:10.1111/j.1755-0998.2009.02733.x

Gompert, Z., & Buerkle, C. A. (2016). What, if anything, are hybrids: Enduring truths and challenges associated with population structure and gene flow. *Evolutionary Applications*, 9(7), 909–923.
doi:10.1111/eva.12380

Goudet, J. (2005). HIERFSTAT, a package for R to compute and test hierarchical F -statistics. *Molecular Ecology Notes*, 5, 184–186.
doi:10.1111/j.1471-8278

Harrison, R. G., & Larson, E. L. (2014). Hybridization, introgression, and the nature of species boundaries. *Journal of Heredity*, 105(S1), 795–809.
doi:10.1093/jhered/esu033

Heiser, C. B. (1993). *The gourd book*. Norman, OK: University of Oklahoma Press.

Hornborg, A. (2005). Ethnogenesis, regional integration, and ecology in pre-historic Amazonia: Toward a system perspective. *Current Anthropology*, 46(4), 589–620.
doi:10.1086/431530

Hughes, C. E., Govindarajulu, R., Robertson, A., Filer, D. L., Harris, S. A., & Bailey, C. D. (2007). Serendipitous backyard hybridization and the origin of crops. *Proceedings of the National Academy of Sciences of the United States of America*, 104(36), 14389–14394.
doi:10.1073/pnas.0702193104

Hugh-Jones, S. (2009). The fabricated body: Objects and ancestors in northwest Amazonia. In F. Santos-Granero (Ed.), *The occult life of things: Native amazonian theories of materiality and personhood* (pp. 33–59).
Tucson, AZ: The University of Arizona Press.

IPHAN. (2015). Modos de fazer cuias do baixo Amazonas. Livro do Registro de Saberes: Instituto do Patrimônio Histórico e Artístico Nacional. Brasília, Brazil. Retrieved from http://portal.iphan.gov.br/pagina/detalhes/1055/ Irion, G., Junk, W. J., & Mello, J. A. S. N. (1997). The large central Amazonian river floodplains near Manaus: Geological, climatological, hydrological and geomorphological aspects. In W. J. Junk (Ed.), *The central Amazon floodplain: Ecology of a pulsing system* (pp. 23–46). Berlin: Springer.

Jarvis, D. I., & Hodgkin, T. (1999). Wild relatives and crop cultivars: Detecting natural introgression and farmer selection of new genetic combinations in agroecosystems. *Molecular Ecology*, 8(s1), S159–S173.
doi:10.1046/j.1365-294X.1999.00799.x

Lathrap, D. W. (2010). La cultura del bosque tropical. In S. R. Panduro (Ed.), *El alto Amazonas* (pp. 79–97). Lima-Iquitos: Instituto Cultural Rvna & Chataro.

Leigh, J. W., & Bryant, D. (2015). popart: Full feature software for haplotype network construction. *Methods in Ecology and Evolution*, 6(9), 1110–1116.

Lévi-Strauss, C. (1950). The use of wild plants in tropical South America. In J. H. Steward (Ed.), *Handbook of south American Indians: Physical anthropology, linguistics and cultural geography of south American Indians* (pp. 465–486). Washington, Columbia: Smithsonian Institution Press.

Lévi-Strauss, C. (2004). Mitológicas II: Do mel às cinzas. São Paulo: Cosac Naify.

Librado, P., & Rozas, J. (2009). DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. *Bioinformatics*, 25(11), 1451–1452.
doi:10.1093/bioinformatics/btp187

Lima, R., & Saragoussi, M. (2000). Floodplain home gardens on the central Amazon in Brazil. In W. J. Junk, J. J. Ohly, M. T. F. Piedade, & M. G. M. Soares (Eds.), *The central Amazon floodplain: Actual use and options for a sustainable management* (pp. 243–268). Leiden: Backhuys Publishers.

Martin, R. (1999). Ceremonial offerings and religious practices among Taíno indians: An archaeological investigation of gourd use in Taíno culture. *IU South Bend Undergraduate Research Journal*, 2, 71–76.

Medina, J. T. (1934). The discovery of the Amazon according to the account of the Friar Gaspar de Carvajal. *American Geographical Society*, 17, 4–467.

Meglécz, E., Costedoat, C., Dubut, V., Gilles, A., Malausa, T., Pech, N., & Martin, J. F. (2009). QDD: A user-friendly program to select microsatellite markers and design primers from large sequencing projects. *Bioinformatics*, 26(3), 403–404.
doi:10.1093/bioinformatics/btp670

Meyer, R. S., Duval, A. E., & Jensen, H. R. (2012). Patterns and processes in crop domestication: An historical review and quantitative analysis of 203 global food crops. *New Phytologist*, 196, 29–48.
doi:10.1111/j.1469-8137.2012.04253.x

Michon, G., Nasi, R., & Balent, G. (2013). Public policies and management of rural forests: Lasting alliance or fool's dialogue? *Ecology and Society*, 18(1), 30–42.
doi:10.5751/ES-05706-180130

Miller, A. J. (2008). Characterization of a domesticated tree lineage (*Spondias purpurea*, Anacardiaceae) based on nuclear and chloroplast sequence data. *Journal of the Torrey Botanical Society*, 135(4), 463–474.
doi:10.3159/08-RA-020R.1

Miller, A. J., & Gross, B. L. (2011). From forest to field: Perennial fruit crop domestication. *American Journal of Botany*, 98(9), 1389–1414.
doi:10.3732/ajb.1000522

Montardo, D. L. O. (2002). *Através do mbaraka: Música e xamanismo guanani*. São Paulo: Department of Social Anthropology, Universidade de São Paulo.

Moreira, P. A., Lins, J., Dequigiovanni, G., Veasey, E. A., & Clement, C. R. (2015). The domestication of annatto (*Bixa orellana*) from *Bixa urucurana* in Amazonia. *Economic Botany*, 69(2), 127–135.
doi:10.1007/s12231-015-9304-0

Moreira, P. A., Mariac, C., Scarcelli, N., Couderc, M., Rodrigues, D. P., Clement, C. R., & Vigouroux, Y. (2016). Chloroplast sequence of treegourd (*Crescentia cujete*) to study phylogeography and domestication. *Applications in Plant Sciences*, 4(10), 1600048.
doi:10.3732/apps.1600048

Morton, J. F. (1968). The calabash (*Crescentia cujete*) in folk medicine. *Economic Botany*, 22(3), 273–280.
doi:10.1007/BF02861961

Nei, M. (1987). *Molecular evolutionary genetics*. New York, NY: Columbia University Press.

Paradis, E. (2010). pegas: An R package for population genetics with an integrated-modular approach. *Bioinformatics*, 26(3), 419–420.
doi:10.1093/bioinformatics/btp696

Patiño, V. M. (1967). Plantas cultivadas y animales domésticos en América equinoccial III: Fibras, medicinas, misceláneas. Cali, Vale do Cauca: Imprenta Departamental.

Petersen, J. J., Parker, I. M., & Potter, D. (2014). Domestication of the neotropical tree *Chrysophyllum cainito* from a geographically limited yet genetically diverse gene pool in Panama. *Ecology and Evolution*, 4(5), 539–553.
doi:10.1002/ece3.948
Pickersgill, B. (2013). Some current topics in plant domestication: An overview with particular reference to Amazonia. *Tipiti: Journal of the Society for the Anthropology of Lowland South America, 11*(2), 16-29.

Piperno, D. R. (2011). The origins of plant cultivation and domestication in the New World tropics: Patterns, process, and new developments. *Current Anthropology*, 52(S4), S453-S470. doi:10.1086/659998

Posey, A. D., & Balick, J. M. (2006). *Human impacts on Amazonia: The role of traditional ecological knowledge in conservation and development*. Chichester, West Sussex: Columbia University Press.

Pritchard, J. K., Stephens, M., & Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics*, 155(2), 945-959.

R Foundation for Statistical Computing (2015). *R: A language and environment for statistical computing*, reference index version 3.1.3. Vienna, Austria: R Foundation for Statistical Computing.

Ribeiro, B. G. (1995). *Os indios das águas pretas: Modo de produção e equipa-mento produtivo* (B. G. Ribeiro, Ed.). São Paulo: Editora da Universidade de São Paulo & Companhia das Letras.

Riesenberg, L. H., & Wendel, J. F. (1993). Introgression and its consequences in plants. In R. G. Harrison (Ed.), *Hybrid zones and the evolutionary process* (pp. 70-109). New York, NY: Oxford University Press.

Rodrigues Ferreira, A. (1933). Memoria sobre as cuyas. *Revista Nacional de Educação*, 1(16), 58-63.

Rollo, A., Lojka, B., Honys, D., Mandák, B., Wong, J. A. C., Santos, C., & Ribeiro, M. M. (2016). Genetic diversity and hybridization in the two species *Inga ingoides* and *Inga edulis*: Potential applications for agroforestry in the Peruvian Amazon. *Annals of Forest Science*, 73(2), 425-435. doi:10.1007/s13595-015-0535-0

Santos, A. M. de S. (1982). Aritapera: uma comunidade de pequenos produtores na várzea amazônica (Santarém-PA). *Boletim Do Museu Emílio Goeldi*, 83, 1-52.

Scarcelli, N., Mariac, C., Couvreur, T. L. P., Faye, A., Richard, D., Sabot, F., & Vigouroux, Y. (2016). Intra-individual polymorphism in chloroplasts from NGS data: Where does it come from and how to handle it? *Molecular Ecology Resources*, 16(2), 434-445. doi:10.1111/1755-0998.12462

Smith, M., & Fausto, C. (2016). Sociality and diversity of pequi (Caryocar brasiliense Caryocaraceae) among the Kuikuro of the Upper Xingu river (Brazil). *Boletim Do Museu Paraense Emilio Goeldi*, 11(1), 87-113. doi:10.1590/1981.81222016000100006

Steward, J. H. (1948). *Handbook of south american indians: The tropical forest tribes*. (J. H., Steward, Ed.) Washington, Columbia: Bureau of American Ethnology, Smithsonian Institution Press.

Terrell, J. E., Hart, J. P., Barut, S., Cellinese, N., Curet, A., Denham, T., & Staller, J. E. (2003). Domesticated landscapes: The subsistence ecology of plant and animal domestication. *Journal of Archaeological Method and Theory*, 10(4), 323-368.

Waldhoff, D., Ulrich, S.-P., & Furch, B. (1996). Value of fruits and seeds from the floodplain forests of central Amazonia as food resource for fish. *Ecotropica*, 2, 143-156.

Wittmann, F., Schöngart, J., Montero, J. C., Motzer, T., Junk, W. J., Piedade, M. T. F., & Worbes, M. (2006). Tree species composition and diversity gradients in white-water forests across the Amazon basin. *Journal of Biogeography*, 33(8), 1334-1347. doi:10.1111/j.1365-2699.2006.01495.x

Wittmann, F., & Wittmann, A. de O. (2011). Use of amazonian floodplain trees. In W. J. Junk, F. Wittmann, P. Parolin, M. T. F. Piedade, & J. Schöngart (Eds.), *Amazonian floodplain forests: Ecophysiology, biodiversity and sustainable management, ecological studies* (pp. 389-418). London: Springer. doi: 10.1007/978-90-481-8725-6_19

Zerega, N. J. C., Ragone, D., & Motley, T. (2004). Complex origins of breadfruit (*Artocarpus altilis*, Moraceae): Implications for human migrations in Oceania. *American Journal of Botany*, 91(5), 760-766. doi:10.3732/ajb.91.5.760

Zohary, D., & Spiegel-Roy, P. (1975). Beginnings of fruit growing in the old world. *Science, New Series*, 187(4174), 319-327.

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the supporting information tab for this article.

**How to cite this article:** Ambrósio Moreira P, Mariac C, Zekraoui L, et al. Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin. *Evol Appl.* 2017;00:1-13. https://doi.org/10.1111/eva.12474