Genetic analyses and dispersal patterns unveil the Amazonian origin of guava domestication.

Guava (Psidium guajava L.) is a semi-domesticated fruit tree of moderate importance in the Neotropics, utilized for millennia due to its nutritional and medicinal benefits, but its origin of domestication remains unknown. In this study, we examine genetic diversity and population structure in 215 plants from 11 countries in Mesoamerica, the Andes, and Amazonia using 25 nuclear microsatellite loci to propose an origin of domestication. Genetic analyses reveal one gene pool in Mesoamerica (Mexico) and four in South America (Brazilian Amazonia, Peruvian Amazonia and Andes, and Colombia), indicating greater differentiation among localities, possibly due to isolation between guava populations, particularly in the Amazonian and Andean regions. Moreover, Mesoamerican populations show high genetic diversity, with moderate genetic structure due to gene flow from northern South American populations. Dispersal scenarios suggest that Brazilian Amazonia is the probable origin of guava domestication, spreading from there to the Peruvian Andes, northern South America, Central America, and Mexico. These findings present the first evidence of guava domestication in the Americas, contributing to a deeper understanding of its evolutionary history.

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.

Population structure in diverse pepper (Capsicum spp.) accessions.

BACKGROUND: Peppers, bell and chile, are a culturally and economically important worldwide. Domesticated Capsicum spp. are distributed globally and represent a complex of valuable genetic resources. OBJECTIVES: Explore population structure and diversity in a collection of 467 peppers representing eight species, spanning the spectrum from highly domesticated to wild using 22,916 SNP markers distributed across the twelve chromosomes of pepper. RESULTS: These species contained varied levels of genetic diversity, which also varied across chromosomes; the species also differ in the size of genetic bottlenecks they have experienced. We found that levels of diversity negatively correlate to levels of domestication, with the more diverse being the least domesticated.

Environment of origin and domestication affect morphological, physiological, and agronomic response to water deficit in chile pepper (Capsicum sp.).

Global climate change is having a significant effect on agriculture by causing greater precipitation variability and an increased risk of drought. To mitigate these effects, it is important to identify specific traits, adaptations, and germplasm that improve tolerance to soil water deficit. Local varieties, known as landraces, have undergone generations of farmer-mediated selection and can serve as sources of variation, specifically for tolerance to abiotic stress. Landraces can possess local adaptations, where accessions adapted to a particular environment will outperform others grown under the same conditions. We explore adaptations to water deficit in chile pepper landraces from across an environmental gradient in Mexico, a center of crop domestication and diversity, as well in improved varieties bred for the US. In the present study, we evaluated 25 US and Mexico accessions in a greenhouse experiment under well-watered and water deficit conditions and measured morphological, physiological, and agronomic traits. Accession and irrigation regime influenced plant biomass and height, while branching, CO2 assimilation, and fruit weight were all influenced by an interaction between accession and irrigation. A priori group contrasts revealed possible adaptations to water deficit for branching, CO2 assimilation, and plant height associated with geographic origin, domestication level, and pepper species. Additionally, within the Mexican landraces, the number of primary branches had a strong relationship with precipitation from the environment of origin. This work provides insight into chile pepper response to water deficit and adaptation to drought and identifies possibly tolerant germplasm.

Genetic diversity, gene flow, and differentiation among wild, semiwild, and landrace chile pepper (Capsicum annuum) populations in Oaxaca, Mexico.

PREMISE: Capsicum annuum (Solanaceae) was originally domesticated in Mexico, where wild (C. annuum var. glabriusculum) and cultivated (C. annuum var. annuum) chile pepper populations (>60 landraces) are common, and wild-resembling individuals (hereafter semiwild) grow spontaneously in anthropogenic environments. Here we analyze the role of elevation and domestication gradients in shaping the genetic diversity in C. annuum from the state of Oaxaca, Mexico. METHODS: We collected samples of 341 individuals from 28 populations, corresponding to wild, semiwild (C. annuum var. glabriusculum) and cultivated C. annuum, and closely related species Capsicum frutescens and C. chinense. From the genetic variation of 10 simple sequence repeat (SSR) loci, we assessed the population genetic structure, inbreeding, and gene flow through variance distribution analyses, genetic clustering, and connectivity estimations. RESULTS: Genetic diversity (HE ) did not differ across domestication levels. However, inbreeding coefficients were higher in semiwild and cultivated chiles than in wild populations. We found evidence for gene flow between wild populations and cultivated landraces along the coast. Genetic structure analysis revealed strong differentiation between most highland and lowland landraces. CONCLUSIONS: Gene flow between wild and domesticated populations may be mediated by backyards and smallholder farms, while mating systems may facilitate gene flow between landraces and semiwild populations. Domestication and elevation may overlap in their influence on genetic differentiation. Lowland Gui'ña dani clustered with highland landraces perhaps due to the social history of the Zapotec peoples. In situ conservation may play an important role in preserving semiwild populations and private alleles found in landraces.

Germination response of diverse wild and landrace chile peppers (Capsicum spp.) under drought stress simulated with polyethylene glycol.

Responses to drought within a single species may vary based on plant developmental stage, drought severity, and the avoidance or tolerance mechanisms employed. Early drought stress can restrict emergence and seedling growth. Thus, in areas where water availability is limited, rapid germination leading to early plant establishment may be beneficial. Alternatively, germination without sufficient water to support the seedling may lead to early senescence, so reduced germination under low moisture conditions may be adaptive at the level of the population. We studied the germination response to osmotic stress of diverse chile pepper germplasm collected in southern Mexico from varied ecozones, cultivation systems, and of named landraces. Drought stress was simulated using polyethylene glycol solutions. Overall, survival time analysis revealed delayed germination at the 20% concentration of PEG across all ecozones. The effect was most pronounced in the genotypes from hotter, drier ecozones. Additionally, accessions from wetter and cooler ecozones had the fastest rate of germination. Moreover, accessions of the landraces Costeño Rojo and Tusta germinated more slowly and incompletely if sourced from a drier ecozone than a wetter one, indicating that slower, reduced germination under drought stress may be an adaptive avoidance mechanism. Significant differences were also observed between named landraces, with more domesticated types from intensive cultivation systems nearly always germinating faster than small-fruited backyard- or wild-types, perhaps due to the fact that the smaller-fruited accessions may have undergone less selection. Thus, we conclude that there is evidence of local adaptation to both ecozone of origin and source cultivation system in germination characteristics of diverse chile peppers.

Genome-wide genotyping of a novel Mexican Chile Pepper collection illuminates the history of landrace differentiation after Capsicum annuum L. domestication.

Studies of genetic diversity among phenotypically distinct crop landraces improve our understanding of fruit evolution and genome structure under domestication. Chile peppers (Capsicum spp. L.) are economically valuable and culturally important species, and extensive phenotypic variation among landraces exists in southern Mexico, a center of C. annuum diversity. We collected 103 chile pepper seed accessions from 22 named landraces across 27 locations in southern Mexico. We genotyped these accessions with genotyping by sequencing (GBS), yielding 32,623 filtered single-nucleotide polymorphisms. Afterward, we genotyped 32 additional C. annuum accessions from a global collection for comparison to the Mexican collection. Within the Mexican collection, genetic assignment analyses showed clear genetic differentiation between landraces and clarified the unique nature of the Tusta landrace. Further clustering analyses indicated that the largest fresh-use Chile de Agua and dry-use Costeño landraces were part of separate clades, indicating that these two landraces likely represent distinct populations. The global accessions showed considerable admixture and limited clustering, which may be due to the collapse of use-type divisions outside of Central America. The separation of the Mexican landraces in part by fruit morphology related to use highlights the relevance of this use-type morphological diversity for plant breeders and the utility of fruit development variation for evolutionary biologists.

Appropriate homoplasy metrics in linked SSRs to predict an underestimation of demographic expansion times.

BACKGROUND: Homoplasy affects demographic inference estimates. This effect has been recognized and corrective methods have been developed. However, no studies so far have defined what homoplasy metrics best describe the effects on demographic inference, or have attempted to estimate such metrics in real data. Here we study how homoplasy in chloroplast microsatellites (cpSSR) affects inference of population expansion time. cpSSRs are popular markers for inferring historical demography in plants due to their high mutation rate and limited recombination. RESULTS: In cpSSRs, homoplasy is usually quantified as the probability that two markers or haplotypes that are identical by state are not identical by descent (Homoplasy index, P). Here we propose a new measure of multi-locus homoplasy in linked SSR called Distance Homoplasy (DH), which measures the proportion of pairwise differences not observed due to homoplasy, and we compare it to P and its per cpSSR locus average, which we call Mean Size Homoplasy (MSH). We use simulations and analytical derivations to show that, out of the three homoplasy metrics analyzed, MSH and DH are more correlated to changes in the population expansion time and to the underestimation of that demographic parameter using cpSSR. We perform simulations to show that Approximate Bayesian Computation (ABC) can be used to obtain reasonable estimates of MSH and DH. Finally, we use ABC to estimate the expansion time, MSH and DH from a chloroplast SSR dataset in Pinus caribaea. To our knowledge, this is the first time that homoplasy has been estimated in population genetic data. CONCLUSIONS: We show that MSH and DH should be used to quantify how homoplasy affects estimates of population expansion time. We also demonstrate how ABC provides a methodology to estimate homoplasy in population genetic data.

Phylogenetic relationships and species delimitation in pinus section trifoliae inferrred from plastid DNA.

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.

Phylogeography of Pinus subsection Australes in the Caribbean Basin.

BACKGROUND AND AIMS: Four species of Pinus subsection Australes occur in the Caribbean Basin: P. caribaea, P. cubensis, P. maestrensis and P. occidentalis. This study analyses the phylogeography of these species to assess possible colonization events from Central America to the islands and subsequent population expansions during glacial periods driven by both drier climate and larger emerged land areas. METHODS: Allele size data were obtained for plastid microsatellites for 314 individuals from 24 populations, covering the distribution range of subsection Australes in the Caribbean Basin. KEY RESULTS: In total, 113 plastid haplotypes were identified. The highest genetic diversity was found in populations of P. caribaea. Overall, Caribbean Basin populations fit the isolation by distance model. Significant phylogeographical structure was found (R(ST) = 0·671 > permuted R(ST) = 0·101; P < 0·0001). The haplotype network and a Bayesian analysis of population structure (BAPS) indicated different Central American origins for P. caribaea var. bahamensis and P. caribaea var. caribaea plastids, with Central America populations in northern and south-eastern groups. Sudden expansion times for BAPS clusters were close to three glacial maxima. CONCLUSIONS: Central America contains ancestral plastid haplotypes. Population expansion has played a major role in the distribution of genetic diversity in P. caribaea var. hondurensis. Two colonization events gave rise to the P. caribaea var. bahamensis and P. caribaea var. caribaea lineages. Plastid variation in the eastern species (P. cubensis, P. maestrensis and P. occidentalis) evolved independently from that in P. caribaea var. caribaea. Incomplete lineage sorting between P. cubensis and P. maestrensis is apparent. Inferred expansion times for P. caribaea var. bahamensis and for the eastern lineages correspond to glacial maxima, whereas those for P. caribaea var. hondurensis correspond to the beginning of the temperature decrease that led to Marine Isotope Stage 8.